domains / families / science / variant-annotation
Annotates VCF variants with functional consequences, population frequencies, and pathogenicity scores using bcftools annotate/csq, Ensembl VEP, SnpEff, and ANNOVAR. Use when deciding which annotation engine and version to pin, which transcript set to report on (RefSeq vs Ensembl vs MANE Select/Plus Clinical, and why VEP --pick is dangerous clinically), how to reconcile HGVS 3'-shifting with VCF left-alignment, which consequence plus NMD status governs PVS1 eligibility, which single calibrated predictor to use for PP3/BP4 (REVEL, AlphaMissense, CADD, SpliceAI deltas), or how to read gnomAD v2/v3/v4 grpmax filtering allele frequency instead of one global AF cutoff. Not for ACMG combining rules or final classification (see variant-calling/clinical-interpretation). Query and annotate gene variants from ClinVar and dbSNP databases. \n\.
The three skills cover the same high-level task but operate at different levels of abstraction and tooling. The first (aipoch) is a high-level wrapper for a Python script that queries ClinVar/dbSNP for clinical significance. The second (FreedomIntelligence) is a practical guide for using command-line tools (bcftools, VEP, SnpEff, ANNOVAR) to add annotations to VCF files. The third (GPTomics) is a deep-dive into the theoretical and practical pitfalls of variant annotation, focusing on reproducibility, transcript selection (MANE vs RefSeq), and normalization issues. No single copy contains the full procedure; the first lacks the CLI tooling details, the second lacks the deep theoretical context and specific clinical reporting standards (MANE), and the third lacks the simple 'how-to' commands for basic annotation. Combining them provides a complete workflow from basic annotation to clinical-grade reproducibility.