Not able to filter 1000G (hg19_ALL.sites.2015_08.txt) in Annovar
0
1
Entering edit mode
9.2 years ago

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 annotate_variation.pl line 2277, <DB> line ###.

This error occurs in every line of the database.

On line 2277 from annotate_variation.pl seems 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 am I doing wrong?

I downloaded all databases using:

perl annotate_variation.pl -downdb -webfrom annovar --buildver hg19 1000g2015aug humandb

My commands are:

convert2annovar.pl -format vcf4 <vcf> -includeinfo -outfile <vcf_avinput> -allsample
table_annovar.pl <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>

Thanks

annovar • 4.6k views
ADD COMMENT
0
Entering edit mode

Hi, I have the same problem, how did you fix it?

ADD REPLY
0
Entering edit mode

Hi biotefarosero,

please do not ask question in old threads, same here A: Annovar Annotation With 1000G. Please open a new question, provide the necessary background information.

ADD REPLY

Login before adding your answer.

Traffic: 1607 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6