Variant calling with GATK HaplotypeCaller following best practices. Covers germline SNP/indel calling, GVCF workflow for cohorts, joint genotyping, and variant quality score recalibration (VQSR).
SkillScore100Maintained
Automatic
◆Research pack
Last verified 8/2/2026
What it changes
Passed 24 of 24 test cases with the skill, and 21 of 24 without it.
Without skill88%
With skill100%
Checks it passed that the model alone missed
Sets --truth-sensitivity-filter-level 99.5.
Runs gatk HaplotypeCaller followed by gatk GenotypeGVCFs.
Sets HapMap resource options as --resource:hapmap,known=false,training=true,truth=true,prior=15.0 hapmap.vcf.gz.
3 fixed21 passed either way0 still failing0 broke
Pass-rate lift
+13 pp
Tokens
+233%
Agent turns
0%
Response time
-50%
Model gemini-3.6-flashMeasured 8/2/2026v124 of 24 cases passed with skill
Tokens, turns and time describe resource use, not a dollar price.
See it in action
One of 3 test cases the skill fixed. Same model, same prompt, with and without the skill.
1 of 3
Prompt
In the standard GATK VQSR workflow for human whole genome SNP filtering, what value should be passed to `--truth-sensitivity-filter-level` in `gatk ApplyVQSR`?
What was checked, without and with the skill
Checked
Without skill
With skill
Sets --truth-sensitivity-filter-level 99.5.
Missed
Passed
case-12 · gemini-3.6-flash · 8/2/2026 · answers are capped and redacted
What's inside
The instructions your model reads when the skill loads.
SKILL.md387 linessingle file
GATK Variant Calling
GATK HaplotypeCaller is the gold standard for germline variant calling. This skill covers the GATK Best Practices workflow.
Prerequisites
BAM files should be preprocessed:
Mark duplicates
Base quality score recalibration (BQSR) - optional but recommended