Question: SnpSift annotate reporting extra alleles not present in the input VCF
0
gravatar for fwuffy
3.8 years ago by
fwuffy50
Michigan, USA
fwuffy50 wrote:

Hello everyone

SnpSift annotate seems to be including a variant allele that is not in my input VCF, when I annotate to clinVar.

For example, here's a row from a test VCF:

1 201364336 rs74315379 G A

G>A is the only call made. G matches the reference assembly. 

Yet in the INFO field, we get:
CLNHGVS=NC_000001.11:g.201364336G>A,NC_000001.11:g.201364336G>T;

Why would SnpSift be also reporting the G>T variant?  There is only 1 data row in the test file and no mixed allele frequencies or anything else to confound the annotation. It should just be reporting records matching the A variant allele because other alleles are not present in my sample.

Sure hope I'm missing something obvious here, but it appears to be grabbing every possible variant for that position from the VCF I'm annotating to (just those two exist for that position in ClinVar), without regard to the actual variant allele.  Is there any way to make SnpSift do what I want here and just annotate to the variant that's in my sample?

Thanks

(Edits for clarity)

snp snpeff snpsift • 1.6k views
ADD COMMENTlink modified 3.1 years ago by Biostar ♦♦ 20 • written 3.8 years ago by fwuffy50
0
gravatar for Dan Gaston
3.6 years ago by
Dan Gaston7.1k
Canada
Dan Gaston7.1k wrote:

All it is adding is the ClinVar info and nomenclature that is available for that SNP ID/position. All it is doing is grabbing the info present in the Clinvar VCF file you are annotating with and adding it to the appropriate info field. My suspicion, which can be born out by going into the ClinVar VCF file and looking, is that it is not split by different reported alleles, hence why you get both in the annotation field. It is just a function of how things are reported in the ClinVar vcf.

ADD COMMENTlink written 3.6 years ago by Dan Gaston7.1k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1924 users visited in the last hour