I have whole genome sequencing paired end data that I am trying to map to the human genome primary sequence hg38.
I mapped using bwa-mem after index the genome with bwa index command.
However, the mapping rate is very low here is the outcome of samtools flagstat downstream I am interested in variant calling using varscan2 or other variant caller.
681299381 + 0 in total (QC-passed reads + QC-failed reads)
0 + 0 secondary
6190495 + 0 supplementary
0 + 0 duplicates
217048997 + 0 mapped (31.86% : N/A)
675108886 + 0 paired in sequencing
337554443 + 0 read1
337554443 + 0 read2
161718456 + 0 properly paired (23.95% : N/A)
178076376 + 0 with itself and mate mapped
Do you have any idea how to improve the mapping rate, I mean is this rate of mapping expected in general or am I doing something wrong? I usually do exome or rnaseq where mapping rate is in the 90s.
0 answers
No answers yet.
Log in to answer this question.
-
Variation in mapping percentage with genome
written by onkar •
-
Separating unmapped reads from bam files is not working using samtools -f 4
written by sorrymouse
-
Bowtie2 and BWA-MEM giving very different results in metagenomic data
written by Antonio Camargo
-
BWA gives wrong number of input sequences
written by agata88
-
Comparing alignment stats BWAmem vs. Soap2 aligner
written by pinn
-
QC-passed reads in BWA aln and BWA mem
written by kai •
-
Cannot align reads to plasmid
written by David
-
Different alignments rates with bwa mem (0%) and bowtie2 (82%)
written by James Ashmore
-
when bwa mem run with different thread, output.sam is different.
written by Peavil •
-
High Mapping but No Reads Properly Paired
written by pld
You can take a sample of reads that are not mapping and check them using blast at NCBI. This should help you eliminate potential contamination issues.