Annovar Annotation With 1000G
2
2
Entering edit mode
8.8 years ago

Hello everyone,

I have recently downloaded Annovar, and succesfully filtered my files with a few databases, downloaded from Annovar. But when I try to do the exact same thing with the 1000 genomes project database, I get this error many times in the terminal, one for each line of **.

Argument "&" isn't numeric in numeric eq (==) at annotate_variation.pl line 1962, <DB> line &&.

& can be A, C, G, or T.

&& every line of the database

I have tried filtering with 1000g, 1000g2010 and 1000g2010jul databases, but all give these errors.

An example of the commands I use for filtering:

perl annotate_variation.pl --filter --dbtype CEU.sites.2010_07 testfile1.txt humandb/

Can anyone tell me if I'm doing something wrong, if this is known, and if I can do anything about it?

annovar 1000genomes • 6.0k views
1
Entering edit mode

I encounter the same issue as yours. But I still don't know how to solve this.

0
Entering edit mode

did you try to open

annotate_variation.pl script in line 1962 and see what is written in it?

as shown in the error in your question it expecting some thing to check against "numeric" but instead it find & give it a try and see annotate_variation.pl at 1962 or around it.

0
Entering edit mode

Hi, I have the same problem, how did you fix it? please help me, Argument "A" isn't numeric in numeric eq (==) at ./annotate_variation.pl line 2583, <db> line 5991546.

3
Entering edit mode
8.8 years ago

you are getting this error because annotate_variation.pl seems to be filling $end with the reference allele column (line 1962 contains an equivalence check for$start and $end), but you only would have entered in that "else" section if the line you are processing isn't recognized as a 1000g data line. since the software you're using is the same as the one we others use (unless you have edited it, of course), the only thing I can think of is that you're using wrong 1000g data files. make sure you download the 1000g data files as follows: perl annotate_variation.pl -downdb -buildver hg19 1000g2012feb humandb  and then you filter/annotate your data like this: perl annotate_variation.pl -filter -dbtype 1000g2012feb_all -buildver hg19 testfile1.txt humandb  of course, you may change the above "1000g2012feb" tag for any desired 1000g version available. ADD COMMENT 2 Entering edit mode in case you are using hg18 you should go for the appropriate databases. in the subsection number 4 of the "Prepare Database" section you'll see a list of databases and their corresponding reference: some will have both hg18 and hg19 options, some will be only available for one reference only. just download the ones you are interested in and the ones that suit your experiment's requirements. a detailed list of the available databases to download can be found here. moderator's note: if you need to comment or ask anything just press the appropriate button down here, avoiding creating answers with new questions. if you need to modify your initial question you can edit it, and if you need to ask something new just post a new one. ADD REPLY 0 Entering edit mode Excuse me: i encountered another question when i use annovar.The commmand was as follows:perl convert2annovar.pl -format vcf4 input.vcf -outfile output.avinput -allsample -withfreq -include -withzyg Then is printed the error message "Use of uninitialized value$zygosity in join or string at /home/linan/annovar/convert2annovar.pl line 2112, <VAR> line 143.
Can't use an undefined value as a symbol reference at /home/linan/annovar/convert2annovar.pl line 2112, <VAR> line 143."

i do not know why.thx in advance.

0
Entering edit mode

you should open a new question, instead of commenting a (rather old) previous answer. before you create it, have in mind that maybe the "-include" option you are using should be the "-includeinfo", and that the -withfreq and -withzyg options used together generate pairs of annotations, one with frequency and other with zygosity information. you may want to further study your command before asking for help.

0
Entering edit mode

I filtered it out & rectified the above problem. But i only get 0 values at 1000g column in the final CSV file. please help.