Hi there, I am mapping single-end reads to a draft reference genome using BWA. The workflow I was given says that a 90% map is good, but does not have anymore information than that. The example flagstat output is given as something like this (i've truncated it for easier reading):
1182945 + 0 in total (QC-passed reads + QC-failed reads)
0 + 0 secondary
2498 + 0 supplementary
0 + 0 duplicates
1093076 + 0 mapped (92.40% : N/A)
However, I get a huge range in percent mapped (eg: up to almost 100% and way down to 0.6%). Is there information about expected % mapping? I'm using draft genomes that are known to have some mistakes, so can I lower this acceptability threshold down to say 80%?
•
written
by
BeetleCheers •
0 answers
No answers yet.
Log in to answer this question.
-
Samtools flagstat
written by michelafrancesconi8 •
-
"Secondary", "Supplementary", "Duplicates" and "Paired in sequencing" in samtools flagstat
written by pasha64t •
-
Separating unmapped reads from bam files is not working using samtools -f 4
written by sorrymouse
-
Interpreting mapping contaminants
written by cecilio11 •
-
Why do I have more reads in my .bam file than in my .fastq input file?
written by francois •
-
QC-passed reads in BWA aln and BWA mem
written by kai •
-
Samtools: How can I extract properly-paired QC-passed reads instead of extracting only properly-pai…
written by bioinfo8
-
when bwa mem run with different thread, output.sam is different.
written by Peavil •
-
High Mapping but No Reads Properly Paired
written by pld
-
What does secondary (and some other fields) mean in flagstat output?
written by Niek De Klein
What 'goodness" are you trying to measure here? A low mapping % might mean a poor quality run, or it might mean your reference genome doesn't match your sample.