Question: Parsing through VCF annotated by Funcotator
gravatar for abbysue
11 months ago by
abbysue10 wrote:

I have a variant set annotated by GATK's Funcotator... creates an extra field in the VCF called "Funcotator" with several annotations separated by pipe characters. I want to get some statistics on % intronic vs exonic and missense/nonsense mutations. Is the easiest way to do this using grep? Or does GATK's SelectVariants work?

funcotator annotated gatk vcf • 702 views
ADD COMMENTlink modified 9 months ago by jonnsmith0 • written 11 months ago by abbysue10
gravatar for jonnsmith
9 months ago by
jonnsmith0 wrote:

There are a few ways you can do this. I'd recommend grabbing a script that I wrote to turn Functoator VCF annotations into CSVs.

I did this so that I could more easily debug things when I was writing it.

The script is here:

That should work for you (though it may not work correctly if the VCF header is very large). If you have more troubles you can post over on the GATK forums or open a question ticket on GitHub.

ADD COMMENTlink written 9 months ago by jonnsmith0
Please log in to add an answer.


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