Allele-specific rank sum test for mapping qualities of REF versus ALT reads (AS_MQRankSum)
Category Variant Annotations
OverviewAllele specific Rank Sum Test for mapping qualities of REF versus ALT reads
This variant-level annotation compares the mapping qualities of the reads supporting the reference allele with those supporting each alternate allele. To be clear, it does so separately for each alternate allele.
The ideal result is a value close to zero, which indicates there is little to no difference. A negative value indicates that the reads supporting the alternate allele have lower mapping quality scores than those supporting the reference allele. Conversely, a positive value indicates that the reads supporting the alternate allele have higher mapping quality scores than those supporting the reference allele.
Finding a statistically significant difference in quality either way suggests that the sequencing and/or mapping process may have been biased or affected by an artifact. In practice, we only filter out low negative values when evaluating variant quality because the idea is to filter out variants for which the quality of the data supporting the alternate allele is comparatively low. The reverse case, where it is the quality of data supporting the reference allele that is lower (resulting in positive ranksum scores), is not really informative for filtering variants.
The value output for this annotation is the u-based z-approximation from the Mann-Whitney-Wilcoxon Rank Sum Test for mapping qualities (MAPQ of reads supporting REF vs. MAPQ of reads supporting ALT). See the method document on statistical tests for a more detailed explanation of the ranksum test.
- The mapping quality rank sum test can not be calculated for sites without a mixture of reads showing both the reference and alternate alleles.
- Uninformative reads are not used in these annotations.
- MappingQualityRankSumTest outputs a version of this annotation that includes all alternate alleles in a single calculation.
- RMSMappingQuality gives an estimation of the overal read mapping quality supporting a variant call.
GATK version 22.214.171.124-1-gf548ccd-SNAPSHOT built at Mon, 2 Aug 2021 14:22:45 -0700.
Please sign in to leave a comment.