Hello,
It will be very helpful if somebody can please suggest me if the this a poor alignment or not and I need to exclude QC-failed reads.
I did a tophat alignment from paired end reads from Illumina and executed Samtools flagstat on the aligned BAM file,it shows this stats below:
117590453 + 26766923 in total (QC-passed reads + QC-failed reads)
17600285 + 3418020 duplicates
101724522 + 8960862 mapped (86.51%:33.48%)
117590453 + 26766923 paired in sequencing
58968738 + 13380987 read1
58621715 + 13385936 read2
93403124 + 5491156 properly paired (79.43%:20.51%)
93590388 + 5502264 with itself and mate mapped
8134134 + 3458598 singletons (6.92%:12.92%)
0 + 0 with mate mapped to a different chr
0 + 0 with mate mapped to a different chr (mapQ>=5)
Please suggest.
Thanks in advance
It will be very helpful if somebody can please suggest me if the this a poor alignment or not and I need to exclude QC-failed reads.
I did a tophat alignment from paired end reads from Illumina and executed Samtools flagstat on the aligned BAM file,it shows this stats below:
117590453 + 26766923 in total (QC-passed reads + QC-failed reads)
17600285 + 3418020 duplicates
101724522 + 8960862 mapped (86.51%:33.48%)
117590453 + 26766923 paired in sequencing
58968738 + 13380987 read1
58621715 + 13385936 read2
93403124 + 5491156 properly paired (79.43%:20.51%)
93590388 + 5502264 with itself and mate mapped
8134134 + 3458598 singletons (6.92%:12.92%)
0 + 0 with mate mapped to a different chr
0 + 0 with mate mapped to a different chr (mapQ>=5)
Please suggest.
Thanks in advance