Category

Genome Variant Analysis


Usage

java -jar GenomeAnalysisTK.jar -T GenotypeGVCFs -R reference.fasta --variant sample1.g.vcf --variant sample2.g.vcf -o output.vcf


Manual

Argument name(s)Default valueSummary
Required Inputs
--variant
 -V
NAOne or more input gVCF files
Optional Inputs
--dbsnp
 -D
nonedbSNP file
Optional Outputs
--out
 -o
stdoutFile to which variants should be written
Optional Parameters
--group
 -G
[StandardAnnotation]One or more classes/groups of annotations to apply to variant calls
--heterozygosity
 -hets
0.001Heterozygosity value used to compute prior likelihoods for any locus
--heterozygosity_stdev
 -heterozygosityStandardDeviation
0.01Standard deviation of eterozygosity for SNP and indel calling.
--indel_heterozygosity
 -indelHeterozygosity
1.25E-4Heterozygosity for indel calling
--sample_ploidy
 -ploidy
2Ploidy per sample. For pooled data, set to (Number of samples in each pool * Sample Ploidy).
--standard_min_confidence_threshold_for_calling
 -stand_call_conf
10.0The minimum phred-scaled confidence threshold at which variants should be called
Optional Flags
--annotateNDA
 -nda
falseAnnotate number of alleles observed
--includeNonVariantSites
 -allSites
falseInclude loci found to be non-variant after genotyping
--useNewAFCalculator
 -newQual
falseUse new AF model instead of the so-called exact model
Advanced Parameters
--annotation
 -A
[]One or more specific annotations to recompute. The single value 'none' removes the default annotations
--input_prior
 -inputPrior
[]Input prior for calls
--max_alternate_alleles
 -maxAltAlleles
6Maximum number of alternate alleles to genotype
--max_genotype_count
 -maxGT
1024Maximum number of genotypes to consider at any site
--max_num_PL_values
 -maxNumPLValues
100Maximum number of PL values to output


Share your experience or ask a question