why my WGS bam file (after alignment) shows more reads than the starting fastq file
I have one paired ended (150x2) WGS fastq file.
fastqc QC total read = ~ 459 M
fastqc QC after trimming total reads = ~ 458 M
after the alignment with bwa-mem and realigning around the indel, I used QUALIMAP for QC and the report shows ~723 M reads!
Why do I get more reads after the alignment?
1 answer
Log in to answer this question.
More posts like this
-
Different flagstat number after using MergeBamAlignment
written by ThomasLam •
-
BWA MEM Alignment mapped small part of the reads!!
written by Bunit •
-
Is it ok to trim assembled paired-end fastq files with cutadapt 3.4 ?
written by Bengisu •
-
Strange alignments results
written by antmantras •
-
miRNA fastqc sequence length distribution with UMI
written by maria2019
-
What is the right illumina universal adapter sequence for trimming paired-end reads?
written by maria2019
-
why does qualimap2 shows different results for number of mismatches for bowtie2 and BWAMEM?
written by maria2019
-
QC-passed reads in BWA aln and BWA mem
written by kai •
-
High Mapping but No Reads Properly Paired
written by pld
-
Whole Genome Sequence (WGS) throughput (FastQC, Picard)
written by toshnam
I am having the same issue. The total number of reads is even larger in the bam file generated with filted fastq file than the raw fastq file. Do you know why it happened like this? Thanks