why my WGS bam file (after alignment) shows more reads than the starting fastq file
1
1
Entering edit mode
23 months ago
maria2019 ▴ 180

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?

WGS fastqc qualimap QC total read • 1.1k views
ADD COMMENT
0
Entering edit mode

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

ADD REPLY

Login before adding your answer.

Traffic: 2168 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6