SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
SAM flag idioms kgulukota Bioinformatics 14 08-23-2016 01:32 AM
bitwise flag in sam format and others aurelielaugraud Bioinformatics 48 11-25-2012 12:01 AM
SAM format flag eoh001 Bioinformatics 1 05-14-2012 12:48 PM
Looking process to convert gff3 format into ace format or sam format andylai Bioinformatics 1 05-17-2011 03:09 AM
Flag=4 in SAM Rachelly Bioinformatics 2 12-22-2010 03:54 AM

Reply
 
Thread Tools
Old 11-08-2012, 01:08 PM   #1
mimozz
Junior Member
 
Location: Los angeles

Join Date: Nov 2012
Posts: 2
Default Confusion about one flag in sam format

Hello,

I met one problem in the Tophat mapping. Does anyone know what the flag difference between 97 and 99 is?

The annotation of 99 is: read paired, read mapped in proper pair, mate reverse strand, first in pair

The annotation of 97 is: read paired, mate reverse strand, first in pair

Seems that the only difference is “read mapped in proper pair”. Is it due to the estimated insertion length? Can 97 be counted as a concordant pair?

Many thanks.
mimozz is offline   Reply With Quote
Old 11-08-2012, 02:02 PM   #2
swbarnes2
Senior Member
 
Location: San Diego

Join Date: May 2008
Posts: 912
Default

"Proper pair" likely means that both reads are on the same chromosome, point in towards each other, and maybe that their distance is some appropriate number (I don't know how bowtie determines that, maybe the settings you give it in the command line)

So any of those things could be off, and give you a 97 flag.

So if the reads point away from each other, or are on different chromosomes, no, you should not count them as concordant.
swbarnes2 is offline   Reply With Quote
Old 11-09-2012, 03:18 AM   #3
maubp
Peter (Biopython etc)
 
Location: Dundee, Scotland, UK

Join Date: Jul 2009
Posts: 1,543
Default

The meaning of "read mapped in proper pair" is down to the aligner/assembler, and would consider the expected read pair orientation (typically --> <--, also <-- --> for mate pair, and --> --> or <-- <-- for native Roche 454 reads if not flipped to look like traditional paired reads) if mapped to the same reference, and the expected fragment/template length (i.e. separation of the reads, field TLEN aka ISIZE).

Last edited by maubp; 11-09-2012 at 03:19 AM. Reason: typo
maubp is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 03:13 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO