Question: Phastcons 46-Way Alignments To Annotate Variants That Fall Within Conserved Genomic Regions
gravatar for ivivek_ngs
7.1 years ago by
Seattle,WA, USA
ivivek_ngs5.0k wrote:

Dear All,

I am trying to annotate variants that fall within the conserved genomic regions with annovar script I have downloaded by the database in my /humandb folder of the annovar

path: /scratch/GT/vdas/tools/annovar/humandb/ -rwxrw-r-x 1 vdas DPT 115424594 Jan 30 2013 hg18_phastConsElements28way.txt -rwxrw-r-x 1 vdas DPT 205139082 Jan 30 2013 hg19_phastConsElements46way.txt

But when I am running the below script

perl /scratch/GT/vdas/tools/annovar/ --filter --buildver hg19 --dbtype mce46way /scratch/GT/vdas/pietro/exome_seq/results/T_S7999/SNP_AD_annovar/T_S7999.recal_AD.snps.vcf.annovar /scratch/GT/vdas/tools/annovar/humandb/

I get the below error:

NOTICE: the --dbtype mce46way is assumed to be in generic ANNOVAR database format NOTICE: Variants matching filtering criteria are written to /scratch/GT/vdas/pietro/exome_seq/results/T_S7999/SNP_AD_annovar/T_S7999.recal_AD.snps.vcf.annovar.hg19_mce46way_dropped, other variants are written to /scratch/GT/vdas/pietro/exome_seq/results/T_S7999/SNP_AD_annovar/T_S7999.recal_AD.snps.vcf.annovar.hg19_mce46way_filtered NOTICE: Processing next batch with 27052 unique variants in 27052 input lines Error: cannot read from input database file /scratch/GT/vdas/tools/annovar/humandb/hg19_mce46way.txt: No such file or directory vdas@iit-n2:/scratch/GT/vdas/pietro/exome_seq/results/T_S7999/SNP_AD_annovar$

Even I tried to download it again using the below command: perl /scratch/GT/vdas/tools/annovar/ -downdb -webfrom annovar -buildver hg19 mce46way /scratch/GT/vdas/tools/annovar/humandb/

I get the below error: NOTICE: Web-based checking to see whether ANNOVAR new version is available ... Done NOTICE: Downloading annotation database ... Failed NOTICE: Downloading annotation database ... Failed WARNING: Some files cannot be downloaded, including,

It seems the mce46way which is downloaded already in my cluster is not being recognized by annovar script as the name it is suggesting is different? Is there any other db file that I need to download to annotate my variants within conserved regions? I dont see any wrong command in the script. Did anyone face this issue? Please let me know where I am getting wrong. Any valuable suggestions are welcome.



annovar exome-sequencing snps • 3.3k views
ADD COMMENTlink written 7.1 years ago by ivivek_ngs5.0k

Why 46? UCSC just released 100 species conservation track.

ADD REPLYlink modified 7.1 years ago • written 7.1 years ago by PoGibas4.8k

How to retrieve it and then is the annovar compatible with this new track? Did anyone use this ? If so how do you change the command so that the annovar script recognizes the new 100 track conservation track

ADD REPLYlink written 7.1 years ago by ivivek_ngs5.0k

But I still dint receive the reply of the post I made. Did anyone face these kind of issues regarding the phastcons 46 database recognition? Please share

ADD REPLYlink written 7.1 years ago by ivivek_ngs5.0k
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: 1410 users visited in the last hour