Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
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
Old 02-03-2021, 05:37 AM   #1
Junior Member
Location: Toronto

Join Date: Feb 2021
Posts: 1
Default Samtools Flagstat, what is "with itself and mate mapped" mean?

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.
kdb.chau is offline   Reply With Quote


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 10:09 AM.

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