Hello,

I have a VCF file generated from GATK and the reference genome used for the analysis. I'm interested in detecting amino acid changes resulting from the variants listed in the VCF file without going for annotation like VEP or other tools. I understand that this involves interpreting the variants within coding regions and translating the resulting nucleotide changes into amino acids.

Specifically, I would like to:

1. Parse the VCF file to extract variant information.
2. Use the reference genome to determine the sequence context around each variant.
3. Identify the codon changes and translate them to detect amino acid changes.