Question: Not able to filter 1000G (hg19_ALL.sites.2015_08.txt) in Annovar
gravatar for dsmarcoantonio
3.6 years ago by
dsmarcoantonio20 wrote:

I am getting this error when using the 1000G to annotate may VCF:

Scanning filter database hg19_ALL.sites.2015_08.txt...Argument & isn't numeric in numeric eq (==) at line 2277, <DB> line ###.

This error occurs in every line of the database.

On line 2277 from seens that $start==$end because the hg19_ALL.sites.2015_08.txt only has start positions. But this database comes directly from annovar site.

Anyone know how to fix it??? What are I doing wrong???


I downloaded all databases using:

perl -downdb -webfrom annovar --buildver hg19 1000g2015aug humandb

My commands are: -format vcf4 <vcf> -includeinfo -outfile <vcf_avinput> -allsample <vcf_avinput> humandb -protocol refGene,phastConsElements46way,genomicSuperDups,exac03,esp6500siv2_all,esp6500siv2_ea,esp6500siv2_aa,ALL.sites.2014_10,snp137,snp138,avsift,clinvar_20150330,ljb26_all -operation g,r,r,f,f,f,f,f,f,f,f,f,f -buildver hg19 -otherinfo -remove -nastring . -outfile <annovar.vcf>


forum • 2.1k views
ADD COMMENTlink written 3.6 years ago by dsmarcoantonio20
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: 1515 users visited in the last hour