sambamba-depth: All files must be coordinate-sorted error?
1
0
Entering edit mode
11 months ago
MobiusT ▴ 10

I am trying to use sambamba to an indexed whole genome dataset.

basename_X $dataset
f="$(basename_X -- $dataset)"
sambamba depth base -L genomic.bed $dataset> ./reads/"$f.txt"

However, the above script gives the following error.

sambamba-depth: All files must be coordinate-sorted

What is the reason and how can I solve it?

read linux sambamba • 594 views
ADD COMMENT
1
Entering edit mode
11 months ago

What is the reason

All files must be coordinate-sorted

how can I solve it?

samtools sort

ADD COMMENT
0
Entering edit mode

samtools sort $filename works in that case, right?

ADD REPLY
1
Entering edit mode

samtools sort -o sorted.bam input.bam

ADD REPLY

Login before adding your answer.

Traffic: 1412 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6