there is an error after samtools index
1
0
Entering edit mode
4.1 years ago
543541656 ▴ 20

Hi, there is an error occuring after the following command

[E::hts_idx_push] Unsorted positions on sequence #4: 27130750 followed by 27130362

samtools index: failed to create index for "1_sorted.bam"

and my samtools sort command is

samtools sort 1_raw.bam -n -m 2G -@ 4 -o 1_sorted.bam

[bam_sort_core] merging from 8 files and 4 in-memory blocks...

I want to know where is the problem,who can give we some advice on the error??

DNA-SEQ • 1.9k views
ADD COMMENT
1
Entering edit mode

(remove the -n option); It is because the -n option is sorted by name.

 -n         Sort by read name

Index, requires the bam file sorted by Coordinates.

See the manual of samtools: http://www.htslib.org/doc/samtools-index.html

Index a coordinate-sorted BGZIP-compressed SAM, BAM or CRAM file for fast random access.

ADD REPLY
0
Entering edit mode

thank you,I have solved the problem according to your requirements

ADD REPLY
0
Entering edit mode
4.1 years ago

Is there a reason you are sitting by name? It looks like you can't index a bam so sorted.

ADD COMMENT
0
Entering edit mode

no,I just copy someone else's code,because it is the first time to use that for me.

ADD REPLY

Login before adding your answer.

Traffic: 3148 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