[ensembl-dev] Problems using VEP cache with github version
Hans Vasquez-Gross
havasquezgross at ucdavis.edu
Tue Mar 28 22:34:04 BST 2017
Hi ENSEMBL,
I have successfully been using ensembl tool release 75 for a while. Now,
I'm planning to start using the github version because of the speed
increase as well as new features.
However, I am having trouble running VEP. I get the following error.
ernie-sears:ensembl-vep_git_beta havasquezgross$ ./vep.pl --species
triticum_aestivum -i Josh_mut_v4.mapspart2.HetMinCov3HetMinPer15HomMinCov2.vcf
-o Josh_mut_v4.mapspart2.HetMinCov3HetMinPer15HomMinCov2.vep.vcf --fork 4
--db_version 34 --assembly TGACv1 --force_overwrite
WARNING: Chromosome TGACv1_scaffold_001462_1AL not found in annotation
sources or synonyms on line 1
WARNING: Chromosome TGACv1_scaffold_001462_1AL not found in annotation
sources or synonyms on line 2
WARNING: Chromosome TGACv1_scaffold_001462_1AL not found in annotation
sources or synonyms on line 3
WARNING: Chromosome TGACv1_scaffold_001462_1AL not found in annotation
sources or synonyms on line 4
WARNING: Chromosome TGACv1_scaffold_001462_1AL not found in annotation
sources or synonyms on line 5
WARNING: Chromosome TGACv1_scaffold_001462_1AL not found in annotation
sources or synonyms on line 6
However, looking at the VEP cache, this folder is actually present. Do you
have any suggestions why this error could be showing up?
ernie-sears:TGACv1_scaffold_001462_1AL havasquezgross$ pwd
/Users/havasquezgross/.vep/triticum_aestivum/34_TGACv1/
TGACv1_scaffold_001462_1AL
ernie-sears:TGACv1_scaffold_001462_1AL havasquezgross$ ls -1
1-1000000.gz
1-1000000_var.gz
I also have a symlink in the .vep/triticum_aestivum folder which points the
directory 34_TGACv1 to 34 as well, so I can run with the old release 75
version which I have been using for years. Running with release 75 doesn't
print this same warning and produces results.
time ./variant_effect_predictor.pl -species triticum_aestivum -i
v41.mapspart2.HetMinCov3HetMinPer15HomMinCov2.vcf -o v41.mapspart2.
HetMinCov3HetMinPer15HomMinCov2.vep.vcf --fork 4 --offline --db_version 34
--force_overwrite --compress gzcat --vcf
[...clip...]
2017-03-27 15:22:32 - Writing output
2017-03-27 15:22:32 - Processed 239683 total variants (585 vars/sec, 781
vars/sec total)
2017-03-27 15:22:34 - Wrote stats summary to v41.mapspart2.
HetMinCov3HetMinPer15HomMinCov2.vep.vcf_summary.html
2017-03-27 15:22:34 - Finished!
real 5m9.791s
user 17m47.120s
sys 5m17.075s
Thank you,
-Hans
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.ensembl.org/pipermail/dev_ensembl.org/attachments/20170328/d9a25440/attachment.html>
More information about the Dev
mailing list