Variant annotations

GATK_TeamGATK_Team
edited January 9 in Dictionary

Variant annotations can be produced by HaplotypeCaller, Mutect2, VariantAnnotator and GenotypeGVCFs. The available annotations are listed under Annotations in the Tool Documentation.

Note that some annotations values calculated by the different tools may be different for the same original data. There are a few things that generally account for differences in annotation values, linked to the normal behaviors of the tools. In all cases, you can end up looking at different sets or numbers of reads, which causes some of the annotation values to be different. It's usually not a cause for alarm. Remember that many of these annotations should be interpreted relatively, not absolutely.

Local realignment

HaplotypeCaller and Mutect2 apply a read realignment step that can modify local coverage counts and other annotation values that take coverage into account. In contrast, VariantAnnotator will calculate annotations either based on the pileup if you give it the original BAM file, or it will calculate summary metrics based on existing VCF fields. GenotypeGVCFs only calculates annotations based on existing VCF fields (that's pretty much its raison d'être).

Read filtering

Some tools apply different read filters by default.

Downsampling

Some tools apply downsampling in order to ensure good performance, but they may do so to different depths of coverage by default compared to others.


See related forum discussions here and here.

Post edited by Sheila on
Sign In or Register to comment.