Documenting my 3 hours spent trying to understand what parameter in SAM file tells read orientation in the hope it will turn out to be helpful to someone. (but if my inference below is wrong!, please post a reply!)
I have been looking up ways to find which strand a read is mapping to and in the SAM file, there seems to be two confusing parameters:
After much researching to confirm my belief of 0 or 16 corresponding to read orientation on strand, I found that the SAM specification PDF says:
I then ignored this for some time ('coz I didnt know what this meant!!) before realizing that bit value of 0X10 corresponds to a decimal value of 16 (Thanks Google!!)
So if the second column is 0, it means read is mapping to one strand, and a value of 16 means it is mapping to the other strand. Still havent figured out clearly what XS:A stands for...whether it relates to a read mapping across a splice junction or even a completely exonic read is considered...not quite sure...
these posts may be of help:
AND
I have been looking up ways to find which strand a read is mapping to and in the SAM file, there seems to be two confusing parameters:
Code:
XS:A:+ and XS:A- and 0 or 16 in the second column
Column 2 is bitwise FLAG and a bit value 0x10 means SEQ being reverse complemented
So if the second column is 0, it means read is mapping to one strand, and a value of 16 means it is mapping to the other strand. Still havent figured out clearly what XS:A stands for...whether it relates to a read mapping across a splice junction or even a completely exonic read is considered...not quite sure...
these posts may be of help:
AND
Comment