I have called variants in WGS data with GATK. Now I would like to filter out alla heterozygous SNPs in the exome and which genes they are present in.
The output should be structured:
chrom | position | GT (nucleotides) | gene
1 | 100 | A/T | geneA
1 | 200 | G/T | geneA
2 | 100 | A/C | geneB
Thanks for any help!