I'm working with a colleague of mine to compare the genomes of two intra-generic vertebrate speceis. We're interested in identifying genome-wide differences and at the moment are looking at the effects of single nucleotide and multi nucleotide differences in genes. We've identified variants by making a whole genome alignment and saving the variants in a vcf file. Our problem is that when we try to annotate the vcf file with SNPEff all the single nucleotide differences within gene boundaries are labeled as FRAME_SHIFT which makes no sense.
Has anyone else run into this issue?
Has anyone else run into this issue?
Comment