qualimap2 mean mapping quality

I’ve done a contrast experiment to see the difference between the bam with BQSR and the bam without BQSR.
I use qualimap to evaluate both bams.
This is the confusing part. Using hap.py and the giab na12878 truth vcf, shows the bam with BQSR is better than the one without BQSR, but the mean mapping quality with bqsr is only 5.86.
The mean mapping quality without bqsr is 27.47.
I calculate use the bam:
samtools view aln.bam |
awk '{sum+=$5} END { print "Mean MAPQ =",sum/NR}': the result is around 56.
I wonder what is the algorithm theory of the qualimap used to cauculate the mean mapping quality.



Source link