Calculation method of RMSMappingQuality
Entering edit mode
4 months ago
octpus616 ▴ 60


I am trying calculate RMS (root-mean-square) MQ (MappingQuality) from a bam file for ever sites (even no variation), I visited GATK's blog:

I note they said:


Root Mean Square of the mapping quality of reads across all samples. This annotation provides an estimation of the overall mapping quality of reads supporting a variant call, averaged over all samples in a cohort.

The raw data format for this annotation consists of a list of two entries: the sum of the squared mapping qualities and the number of reads across variant (not homRef) genotypes

So if I want to calculate a RMSMQ for no variation site (homRef), can I use following formula?

RMS_MQ = sqrt(sum(MQ_i^2)/N)      .......... (1)

where, MQ_i is the MQ of reads i coverage the site, N is the total number of reads coverage the site.

If this is correct, I note this section in GATK:

Statistical notes The root mean square is equivalent to the mean of the mapping qualities plus the standard deviation of the mapping qualities.

Its means:

RMS_MQ = mean(MQ_i) + sd(MQ_i)     ........... (2)

Is it correct? but I noted the output of formula (2) is not equal to formula (1)

NGS bam samtools vcf • 529 views
Entering edit mode
11 weeks ago


Have you solved this probelm?

Entering edit mode

Maybe, I think the formal (1) is correct for most situations (for example the earlier version of GATK). However, there are some changes in the way GATK is calculated at present, they seem to only calculate the RMS mapquality that supports the existence of variations in the lastest version. I am not sure, but you can furture check


Login before adding your answer.

Traffic: 3092 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6