This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Not able to filter 1000G (hg19_ALL.sites.2015_08.txt) in Annovar

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

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

0 answers

No answers yet.

Log in to answer this question.