It looks like you're new here. If you want to get involved, click one of these buttons!
Is there a way to select ONLY homozygote calls from a vcf file using SelectVariants?
I understand that the GT for homozygotes is 1/1, whereas the genotype for heterozygotes is 0/1. But when I ude the following command, I get an empty vcf file (i.e. only with a header):
$ java -Xmx4g -jar GenomeAnalysisTK.jar -SelectVariants -R reference.fasta --variant input.vcf -select 'GT == "1/1";' -o output.vcf
I'd appreciate your help on this matter.