![]() |
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
how to get reads with "mate mapped to a different chr" | gene_x | Bioinformatics | 12 | 07-18-2018 04:43 AM |
SNP call by samtools, but ALT is "X" QUAL is "0" | illuminaGA | Bioinformatics | 1 | 12-16-2015 11:25 PM |
Samtools "is recognized as '*'" "truncated file" error | axiom7 | Bioinformatics | 3 | 11-26-2014 02:53 AM |
samtools flagstat reports an odd number of reads "with itself and mate mapped" | syfo | Bioinformatics | 0 | 04-24-2013 03:43 AM |
Tophat error "read #1 does not have a matching mate in the same order" | wangli | RNA Sequencing | 16 | 11-22-2012 04:21 AM |
![]() |
|
Thread Tools |
![]() |
#1 |
Junior Member
Location: Toronto Join Date: Feb 2021
Posts: 1
|
![]()
Here is my flagstat output after running bwa mem on my paired end data.
25730940 + 0 in total (QC-passed reads + QC-failed reads) 111474 + 0 secondary 0 + 0 supplementary 0 + 0 duplicates 25300974 + 0 mapped (98.33% : N/A) 25619466 + 0 paired in sequencing 12809733 + 0 read1 12809733 + 0 read2 24045298 + 0 properly paired (93.86% : N/A) 25182938 + 0 with itself and mate mapped 6562 + 0 singletons (0.03% : N/A) 1071166 + 0 with mate mapped to a different chr 297377 + 0 with mate mapped to a different chr (mapQ>=5) So from my understanding, 93.86% of my reads are properly paired which means the read and its mate mapped with decent distance to the same chromosome. So does "with itself and mate mapped" simply mean that both reads did map but maybe not all of them are within good distance from each other? Does it also mean both did map to the same chromosome? It is unclear to me what that category means but it does make sense to me that the number of reads would be higher. |
![]() |
![]() |
![]() |
Tags |
samtools |
Thread Tools | |
|
|