Question: Variants with multiple functional annotations -- counting occurrences
0
gravatar for olavur
4 months ago by
olavur60
Tórshavn, Faroe Islands
olavur60 wrote:

I annotated my VCF with SnpEff, whether they are missense variants, intron variants, and so on. Now I want to know how many variants of each different type I have. I don't know exactly how to do this because each variant may have (and often has) several functions (usually corresponding to several overlapping genes, as far as I understand).

Considering this, how would one usually go about counting the number of variants of each functional type?

annotation vcf • 206 views
ADD COMMENTlink modified 4 months ago by Kevin Blighe24k • written 4 months ago by olavur60
1
gravatar for Kevin Blighe
4 months ago by
Kevin Blighe24k
USA / Europe / Brazil
Kevin Blighe24k wrote:

I don't know exactly how to do this because each variant may have (and often has) several functions (usually corresponding to several overlapping genes, as far as I understand).

Any given variant can indeed overlap multiple genes that may, for example, share the same promoter; however, each variant can also have different functionality within the same gene's transcript isoforms. A variant can be exonic in one isoform, for example, but intronic or splicing in another.

SnpEff does output a summary file: http://snpeff.sourceforge.net/SnpEff_manual.html#outputSummary

Variant Effect Predictor also outputs a HTML stats file, which may be of use to you:

d

ADD COMMENTlink written 4 months ago by Kevin Blighe24k

I see, thank you. I wonder, though, if SnpEff counts all annotations of all variants or just the first one in the list (which is sorted by impact in decreasing order).

ADD REPLYlink written 4 months ago by olavur60

It appears to be an accumulation of all noted annotations. The author attempts to explain it on the page to which I linked you (above), but it's not something that is easy to explain!

The easiest way to test this is to just annotate a handful of variants in your VCF, and then manually check the annotated stats / results.

ADD REPLYlink written 4 months ago by Kevin Blighe24k
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: 721 users visited in the last hour