Use VariantsToTable to extract alternate allele count

I'm using the following VariantsToTable command options to extract fields from a VCF file:

/usr/lib/jvm/jre-1.8.0-openjdk/bin/java -Xmx8g -jar /home/1GenomeRef/GATK/GATK_3.5/GenomeAnalysisTK.jar \
-T VariantsToTable \
-R ref.fa \
-V file1.vcf \
-F POS -F ID -F REF -F ALT -F QUAL -F FILTER -F AC -F AN -GF GT  \
--showFiltered \
--out outputfile \

This extracts the correct information, but my original VCF file reports each sample genotype (GT field) as an alternate allele count (0/0, 0/1 or 1/1) and the new output file reports the genotype as the base (C/T, for example.) So the GT for sample1 in my original file might be "0/1" but in the new file it's recoded as "C/T."

I prefer to retain the original genotype format but do not see an option that allows me to request this. Is there an option I can use for this? Or another tool I can apply that will quickly recode the new VCF?

Thanks so much. (And I am following the Best Practices Guidelines - although we are using GATK version 3.5, this is a choice we made for the purpose of ensuring the highest possible consistency with older data called using version 3.5)

Tagged:

Best Answer

Answers

Sign In or Register to comment.