Question: Question About Samtools Flag Number 117
1
gravatar for Alby
7.3 years ago by
Alby90
Alby90 wrote:

One of the rows in my alignment data has the samtool flag of 117.

I used http://picard.sourceforge.net/explain-flags.html

to figure out what it means. The following is the summary for flag 117:

  1. read paired
  2. read unmapped
  3. read reverse strand
  4. mate reverse strand
  5. first in pair

My question is, how could a read be "unmapped" and be "reverse strand(aligned to reference)" at the same time? My guess is samtools define "unmapped" as anything that is not "properly mapped"(meaning both reads in a pair mapped correctly)?

Thanks

alignment samtools • 2.6k views
ADD COMMENTlink written 7.3 years ago by Alby90

just to clarify, what do you mean by samtool flag? The flag is definitely ambiguous. In any case, as long as the flag "read unmapped" is set, I wouldn't consider it for any downstream analysis.

ADD REPLYlink written 7.3 years ago by Arun2.3k
3
gravatar for Vikas Bansal
7.3 years ago by
Vikas Bansal2.3k
Berlin, Germany
Vikas Bansal2.3k wrote:

You can find an answer by Istvan here.

117: read paired, read unmapped, read on reverse strand, mate on reverse strand, first in pair

So the output is correct, what it probably means is that the read and its mate are both on the reverse strand therefore the mapping of this read (and its mate) are questionable, thus it got the unmapped designation.

ADD COMMENTlink written 7.3 years ago by Vikas Bansal2.3k

I am not sure if this is correct. The SAM specification states: "If 0x4 is set, no assumptions can be made about [...] bits 0x2, 0x10 and 0x100 ..." where 0x10 is reverse complement 0x20 mate is reverse complement. So to me this means if 0x4 is set the others can be random. Never the less the explanation Istvan came up with might be right, but I would not rely on it.

ADD REPLYlink written 7.2 years ago by lelle800
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 806 users visited in the last hour