Question

Sam FLAGS : how could a read be described as "read unmapped" but "read in reverse strand"?

0

Entering edit mode

6.5 years ago

johnnytam100 ▴ 110

I am so confused why an unmapped read could have a direction. Any explanation for that?

Thank you!

sam bwa • 2.4k views

ADD COMMENT • link updated 6.5 years ago by glihm ▴ 660 • written 6.5 years ago by johnnytam100 ▴ 110

score 2 · Answer 1 · 2017-10-20

2

Entering edit mode

6.5 years ago

glihm ▴ 660

If the read is unmapped, the "strandness bit" has an other significance. From SAMTOOLS documentation:

Bit 0x10 indicates whether SEQ has been reverse complemented and QUAL reversed. When bit 0x4 is unset, this corresponds to the strand to which the segment has been mapped. When 0x4 is set, this indicates whether the unmapped read is stored in its original orientation as it came off the sequencing machine.

0x4 corresponds to the mapped/unmapped bit (If set, the read is considered as unmapped. If unset, the read is considered as mapped)

0x10 corresponds to the strandness bit

So, the translation of the documentation could be "If the reads is not mapped, the strandness bit indicates if the read (unmapped) is stored in itrs original orientation from the sequecing machine".

Which means, you can have an unmapped read with an orientation without any problem. You just have to give to it the good interpretation. ;)

ADD COMMENT • link 6.5 years ago by glihm ▴ 660

0

Entering edit mode

I see. Then what does "stored in its original orientation as it came off the sequencing machine" mean?

And how should I compare a read being described as "read unmapped" and a read being described as "read unmapped" and "read in reverse strand"? In short, one of them just give me one more description as "read in reverse strand", and practically they are the same because they are unmapped anyway?

ADD REPLY • link 6.5 years ago by johnnytam100 ▴ 110

1

Entering edit mode

Actually, I don't think you've to care about unmapped reads in general. Except that if the unmapped reads are your reads of interest (filtering data for instance).

practically they are the same because they are unmapped anyway?

Exactly. ;)

The only thing you have to keep in mind is that, for people interested in oriented sequencing, they are interested in knowing why some reads are unmapped and they want to keep the info of "which was the original orientation of the read during sequencing process". But for most of uses, unmapped reads are reads that are not meeting the parameters you have chosen for your alignments.