Question: sam flag from paired end reads
0
gravatar for Varun Gupta
3.1 years ago by
Varun Gupta1.1k
United States
Varun Gupta1.1k wrote:

HI Everyone, This might be a simple question but I am not able to understand it. For paired end reads I see 83,163 and 99,147 as the sam flags.

Looking at the description from here about sam flags When I put it as 83 I get first in pair and when i put it as 147 I get second in pair although in both the cases the read maps to reverse strand.

What do we mean by first in pair and second in pair?

Thanks

bam • 2.1k views
ADD COMMENTlink modified 3.1 years ago by John12k • written 3.1 years ago by Varun Gupta1.1k
1
gravatar for Pierre Lindenbaum
3.1 years ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum120k wrote:

in paired end experiment, the sequencer produces two files. Something like sample.R1.fastq.gz and sample.R2.fastq.gz which are the reads on 5' and 3' of the sequenced fragment. They both contain the same number of reads.

First in pair means that the read comes from sample.R1.fastq.gz Second in pair means that the read comes from sample.R2.fastq.gz

First and Second in pair flags have no meaning about the mapping of the read .

ADD COMMENTlink written 3.1 years ago by Pierre Lindenbaum120k

HI Pierre,

I understand this but then I think If the read is 147 it's second in pair although its mapped on the reverse strand while 83 would mean read is on reverse strand but it is first in pair. So that's why 83,163 go together and 99,147 go together for a paired end read

ADD REPLYlink written 3.1 years ago by Varun Gupta1.1k
1

First/second in pair refer to the order in which the reads sequenced; this is independent from alignment: even unmapped reads are first or second in pair.

ADD REPLYlink written 3.1 years ago by Charles Plessy2.7k
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: 1035 users visited in the last hour