Variant Calling Nodes
Variant-calling nodes identify SNPs, indels, and structural variants from aligned reads. They expect a coordinate-sorted, de-duplicated BAM and a reference FASTA.
GATK HaplotypeCaller
Industry-standard germline short-variant caller. Performs local de novo assembly of haplotypes in active regions.
Inputs: bam: BAM, reference: FASTA
Outputs: vcf: VCF
| Parameter | Type | Default | Description |
|---|---|---|---|
intervals | BED | — | Restrict calling to regions |
gvcf | Bool | false | Emit a GVCF for joint genotyping |
ploidy | Int | 2 | Sample ploidy |
HaplotypeCaller requires a valid
@RGread group on the input BAM. Set it on your aligner node.
GATK Joint Genotyping
Combines per-sample GVCFs and jointly genotypes a cohort.
Inputs: gvcfs: List<VCF>, reference: FASTA
Outputs: vcf: VCF
bcftools call
Lightweight, fast variant caller (mpileup + call). Good for bacterial genomes and quick analyses.
Inputs: bam: BAM, reference: FASTA
Outputs: vcf: VCF
| Parameter | Type | Default | Description |
|---|---|---|---|
caller | String | -m | Multiallelic caller model |
ploidy | Int | 2 | Ploidy |
Filtering & annotation
| Node | Purpose |
|---|---|
| GATK VariantFiltration | Hard-filter variants by quality/depth. |
| VQSR | Variant quality score recalibration (cohorts). |
| bcftools filter/norm | Normalize and filter a VCF. |
| SnpEff / VEP | Functional annotation of variants. |
A complete calling branch
[Mark Duplicates BAM] → [GATK HaplotypeCaller] → [VariantFiltration] → [SnpEff] → [VCF Output]