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)