Question: Understanding Annovar Output
0
gravatar for newbinf
4 weeks 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 • 111 views
ADD COMMENTlink modified 4 weeks ago by h.mon18k • written 4 weeks ago by newbinf0
2
gravatar for h.mon
4 weeks ago by
h.mon18k
Brazil
h.mon18k 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 4 weeks ago by h.mon18k
0
gravatar for newbinf
4 weeks 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 4 weeks 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: 1462 users visited in the last hour