How to extract properly paired alignments from a bam file using samtools?
This is what samtools flagstat returns :
87410692 + 0 in total (QC-passed reads + QC-failed reads)
17796957 + 0 secondary
0 + 0 supplimentary
0 + 0 duplicates
87410692 + 0 mapped (100.00%:-nan%)
69613735 + 0 paired in sequencing
34834040 + 0 read1
34779695 + 0 read2
19690094 + 0 properly paired (28.28%:-nan%)
65530882 + 0 with itself and mate mapped
4082853 + 0 singletons (5.87%:-nan%)
4317734 + 0 with mate mapped to a different chr
2061682 + 0 with mate mapped to a different chr (mapQ>=5)
1 answer
Hint: samtools view -f2
• 0 views
•
link
updated
by
Ram
•
written
by
Ashutosh Pandey
•
written
by
Ashutosh Pandey
samtools view -f 0x02 also works
• 1 views
•
link
written
by
Kevin Blighe
Log in to answer this question.
More posts like this
-
samtools filter proper paired alignment
written by mynameliuxiang •
-
Samtools: How can I extract properly-paired QC-passed reads instead of extracting only properly-pai…
written by bioinfo8
-
BWA mem paired end vs single end shows unusual flagstat summary
written by rrdavis •
-
RNA seq samtools flagstat
written by pshubhamoy •
-
flagstat interpretation from a bam file
written by zizigolu
-
Tophat mapping less reads (and less properly paired) when more real -r is specified?
written by manekineko
-
Samtools Flagstat Output Format
written by thsay •
-
samtool flagstat question
written by alex
-
Thorough Documentation Of Samtools Flagstat
written by Biomonika (Noolean)
-
Samtools Flagstat Output
written by Nitram •
