Question: Variants with multiple functional annotations -- counting occurrences
0
gravatar for olavur
6 weeks ago by
olavur40
Tórshavn, Faroe Islands
olavur40 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 • 121 views
ADD COMMENTlink modified 6 weeks ago by Kevin Blighe17k • written 6 weeks ago by olavur40
1
gravatar for Kevin Blighe
6 weeks ago by
Kevin Blighe17k
University College London Cancer Institute
Kevin Blighe17k 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 6 weeks ago by Kevin Blighe17k

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 6 weeks ago by olavur40

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 6 weeks ago by Kevin Blighe17k
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: 1161 users visited in the last hour