I have a fasta file that corresponds to my reference genome and a vcf file that corresponds to the SNP call of my data. I would like to get the sequence of each one of the S