如何在conda中运行ensembl-vep



我已经这样安装了:

conda install ensembl-vep=105.0-0

然后安装了这样的人工缓存:

vep_install -a cf -s homo_sapiens -y GRCh38 -c /mnt/gpfs/live/rd01__/ritd-ag-project-rd018o-mdflo13/refs/vep —CONVERT

但我无法用任何命令运行它,例如

vep --dir_cache "/mnt/gpfs/live/rd01__/ritd-ag-project-rd018o-mdflo13/conda/envs/bioinfo/share/ensembl-vep-105.0-0" 
-i "/mnt/gpfs/live/rd01__/ritd-ag-project-rd018o-mdflo13/data/test/manual/results/variants/cohort.norm_recalibrated.vcf" 
-o "/mnt/gpfs/live/rd01__/ritd-ag-project-rd018o-mdflo13/data/test/manual/results/variants/vep_output.txt”

这会给出一条关于下载缓存的错误消息:

IMPORTANT INFORMATION:
The VEP can read gene data from either a local cache or local/remote databases.

或者这个:

vep --cache 
-i "/mnt/gpfs/live/rd01__/ritd-ag-project-rd018o-mdflo13/data/test/manual/results/variants/cohort.norm_recalibrated.vcf" 
-o "/mnt/gpfs/live/rd01__/ritd-ag-project-rd018o-mdflo13/data/test/manual/results/variants/vep_output.txt”

这给出了错误:

MSG: ERROR: Cache directory /mnt/gpfs/home/skgtmdf/.vep/homo_sapiens not found

我想没有人能为我指明正确的方向吗?

我找到了答案。您需要同时添加--cache和--dir_cache参数:
vep --cache --dir_cache "/mnt/gpfs/live/rd01__/ritd-ag-project-rd018o-mdflo13/conda/envs/bioinfo/share/ensembl-vep-105.0-0" 
-i "/mnt/gpfs/live/rd01__/ritd-ag-project-rd018o-mdflo13/data/test/manual/results/variants/cohort.norm_recalibrated.vcf" 
-o "/mnt/gpfs/live/rd01__/ritd-ag-project-rd018o-mdflo13/data/test/manual/results/variants/vep_output.txt"

最新更新