double @HD header for bam
0
0
Entering edit mode
4 months ago
ZheFrench ▴ 570

The @HD should be unique , no ? I don't get why it can be sorted and unsorted at the same time in a two lines header fashion.

samtools view -H sorted_test.bam |grep '@HD'
@HD VN:1.3  SO:coordinate
@HD VN:1.5  SO:unsorted GO:query
bam • 334 views
ADD COMMENT
0
Entering edit mode

How did that happen? Were two headers joined together at some step manually? Strictly speaking that breaks the SAM format spec since only one @HD line is expected though it is optional.

ADD REPLY

Login before adding your answer.

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