Showing docs for version 3.7-0
MVLikelihoodRatio
Likelihood of being a Mendelian Violation
Category
Annotation Modules
VCF Field
INFO (variant-level)
Type
RodRequiringAnnotation
Header definition line
INFO=<ID=MVLR,Number=1,Type=Float,Description="Mendelian violation likelihood ratio: L[MV] - L[No MV]">
Overview
This annotation uses the likelihoods of the genotype calls to assess whether a site is transmitted from parents to offspring according to Mendelian rules. The output is the likelihood of the site being a Mendelian violation, which can be tentatively interpreted either as an indication of error (in the genotype calls) or as a possible mutation. The higher the output value, the more likely there is to be a Mendelian violation. Note that only positive values indicating likely MVs will be annotated; if the value for a given site is negative (indicating that there is no violation) the annotation is not written to the file.
Statistical notes
This annotation considers all possible combinations of all possible genotypes (homozygous-reference, heterozygous, and homozygous-variant) for each member of a trio, which amounts to 27 possible combinations. Using the Phred-scaled genotype likelihoods (PL values) from each individual, the likelihood of each combination is calculated, and the result contributes to the likelihood of the corresponding case (mendelian violation or non-violation) depending on which set it belongs to. See the method document on statistical tests for a more detailed explanation of this statistical test.
Caveats
- The calculation assumes that the organism is diploid.
- This annotation requires a valid pedigree file.
- When multiple trios are present, the annotation is simply the maximum of the likelihood ratios, rather than the strict 1-Prod(1-p_i) calculation, as this can scale poorly for uncertain sites and many trios.
- This annotation can only be used from the VariantAnnotator. If you attempt to use it from the UnifiedGenotyper, the run will fail with an error message to that effect. If you attempt to use it from the HaplotypeCaller, the run will complete successfully but the annotation will not be added to any variants.
Related annotations
- PossibleDeNovo annotates the existence of a de novo mutation in at least one of a set of families/trios.