Call stats aggregator

It’s now easier than ever to recompute info field annotations like AC and AF! We’ve added an aggregator to compute the following metrics all at once:

AC: Allele count per allele, including reference
AF: Allele frequency per allele, including reference
AN: Total number of called alleles
GC: Genotype count per possible genotype

Sample usage to compute population-specific call stats:

  annotatevariants expr -c 'va.EUR_stats = gs.filter(g => sa.pop == "EUR").callStats(g => v)'

See more here:

1 Like