Question: Annovar Annotation With 1000G
gravatar for daphnarisseeuw
6.4 years ago by
daphnarisseeuw30 wrote:

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 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 --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 • 4.5k views
ADD COMMENTlink modified 6.4 years ago • written 6.4 years ago by daphnarisseeuw30

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

ADD REPLYlink written 7 months ago by George.Robert10

did you try to open 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 at 1962 or around it.

ADD REPLYlink modified 6.4 years ago by Istvan Albert ♦♦ 80k • written 6.4 years ago by Medhat8.3k
gravatar for Jorge Amigo
6.4 years ago by
Jorge Amigo11k
Santiago de Compostela, Spain
Jorge Amigo11k wrote:

you are getting this error because 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 -downdb -buildver hg19 1000g2012feb humandb

and then you filter/annotate your data like this:

perl -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 COMMENTlink written 6.4 years ago by Jorge Amigo11k

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 REPLYlink modified 6.4 years ago • written 6.4 years ago by Jorge Amigo11k

Excuse me:

i encountered another question when i use annovar.The commmand was as follows:perl -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/ line 2112, <VAR> line 143.
Can't use an undefined value as a symbol reference at /home/linan/annovar/ line 2112, <VAR> line 143."

i do not know why.thx in advance.

Hoping for your reply!

ADD REPLYlink written 4.7 years ago by 89759864480

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.

ADD REPLYlink written 4.7 years ago by Jorge Amigo11k

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

ADD REPLYlink modified 15 days ago • written 20 days ago by vaish01kv0
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: 1611 users visited in the last hour