How to check if a read in a Bam file is paired or not?
0
0
Entering edit mode
6.0 years ago

From the flag bit 8, we can get obtain that the next segment in the template is unmapped so consequently, the read is not paired but this rule is not always true like in the case of flag value 129 so what's a reliable way to check if a read contains a pair or not?

sam bam reads paired reads single paired reads • 2.2k views
ADD COMMENT
1
Entering edit mode

https://broadinstitute.github.io/picard/explain-flags.html Enjoy :P

"Mate unmapped" means the read is paired, just the mate doesn't map. What exactly do you want to know? Flag 1: read is paired (this is always set in paired-end data), flag 2: both map properly (i.e. FR orientation, same reference sequence; this is what you want in most(?) cases, but is useless if you want to link contigs for instance).

ADD REPLY

Login before adding your answer.

Traffic: 2813 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6