Error while trying ANNOVAR
0
0
Entering edit mode
2.1 years ago

I am trying to annotate a set of mice variants using annovar. I have used the command below:

perl table_annovar.pl C131T137.avinput mousedb/ -buildver mm10 -out myanno -remove -protocol refGene  -operation g

But I am getting the following error
Error:

FASTA sequence file mousedb/mm10_refGeneMrna.fa does not exist. Please use 'annotate_variation.pl --downdb refGene mousedb' download the database.
Error running system command: <annotate_variation.pl -geneanno -buildver mm10 -dbtype refGene -outfile myanno.refGene -exonsort C131T137.avinput mousedb/>

I downloaded the mousedb as suggested in the command above but still there is the above error. Where am I going wrong?

annovar • 935 views
ADD COMMENT
0
Entering edit mode

Output of ls inside that mousedb folder?

ADD REPLY
0
Entering edit mode

I changed the name of the ref genome file to mm10_refGeneMrna.fa and the program ran. I hope I did it correctly. But the output of the annotated file tells whether the variant is exonic/intergenic etc. How do I get the sift/polyphen scores for each variant?

ADD REPLY

Login before adding your answer.

Traffic: 2374 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