Question: Understanding Annovar Output
0
gravatar for newbinf
3 months ago by
newbinf0
newbinf0 wrote:

I have a vcf output from Annovar (and wAnnovar), and need some help understanding all of the headers. I know i am given the reference and the alternative allele. However, I want to know how many reads in my input file were reference and alternative.

I'm looking at the vcf file, and see some trailing information. What do these mean? Do these value tell me how many reads per allele?

ADP=5;WT=0;HET=1;HOM=0;NC=0

GT:GQ:SDP:DP:RD:AD:FREQ:PVAL:RBQ:ABQ:RDF:RDR:ADF:ADR 0/1:6:5:5:3:2:40%:2.2222E-1:75:48:3:0:2:0

annovar vcf • 161 views
ADD COMMENTlink modified 3 months ago by h.mon20k • written 3 months ago by newbinf0
2
gravatar for h.mon
3 months ago by
h.mon20k
Brazil
h.mon20k wrote:

There is an ANNOVAR intro here:

Introduction to VCF file and some of its complications

And you can get a description of each attribute from the VCF headers, see this Pierre Lindenbaum answer.

The VCF specification is quite lengthy, and you have to read it if you want to understand VCF files in depth. However, as VCF files are a bit complicated there is plenty of discussions around, for example here at BioStars:

Regarding vcf output

VCF format: what does AD 1,0 denote?

Vcf files fields: AD vs DP4

ADD COMMENTlink written 3 months ago by h.mon20k
0
gravatar for newbinf
3 months ago by
newbinf0
newbinf0 wrote:

I realize that DP:RD:AD is total reads: reference allele reads: alternative allele reads. Please correct me if I am wrong.

How do I isolate these values. Or are there columns in the VCF file output that contain these values already?

ADD COMMENTlink written 3 months ago by newbinf0
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: 1584 users visited in the last hour