- Notifications
You must be signed in to change notification settings - Fork219
Open
Description
Hi,
I want to convert the vcf file to the maf file, but an error occurs when I run the following code:perl vcf2maf.pl --input-vcf 10084_10084.vcf --output-maf 10084_10084.maf --vep-path /home/data/t230420/miniconda3/envs/vep/bin --vep-data /home/data/t230420 --ref-fasta /home/data/t230420/hg18.fa --vep-forks 1 --tumor-id 10084_10084 --filter-vcf 0
.
The error is:
WARNING: Couldn't retrieve bps around 23:152850324 from reference FASTA: /home/data/t230420/hg18.faWARNING: Couldn't retrieve bps around 23:152880620 from reference FASTA: /home/data/t230420/hg18.faWARNING: Couldn't retrieve bps around 23:152892914 from reference FASTA: /home/data/t230420/hg18.faWARNING: Couldn't retrieve bps around 23:152901148 from reference FASTA: /home/data/t230420/hg18.faWARNING: Couldn't retrieve bps around 23:153313264 from reference FASTA: /home/data/t230420/hg18.faWARNING: Couldn't retrieve bps around 25:214201 from reference FASTA: /home/data/t230420/hg18.faWARNING: Couldn't retrieve bps around 25:1946325 from reference FASTA: /home/data/t230420/hg18.faWARNING: Couldn't retrieve bps around 25:1948777 from reference FASTA: /home/data/t230420/hg18.faWARNING: Couldn't retrieve bps around 25:154685032 from reference FASTA: /home/data/t230420/hg18.faWARNING: Couldn't retrieve bps around 25:154745988 from reference FASTA: /home/data/t230420/hg18.faWARNING: Couldn't retrieve bps around 26:5443 from reference FASTA: /home/data/t230420/hg18.faSTATUS: Running VEP and writing to: ./10084_10084.vep.vcfPossible precedence issue with control flow operator at /home/data/t230420/miniconda3/envs/vep/lib/site_perl/5.26.2/Bio/DB/IndexedBase.pm line 845.-------------------- EXCEPTION --------------------MSG: ERROR: No cache found for homo_sapiens, version 93STACK Bio::EnsEMBL::VEP::CacheDir::dir /home/data/t230420/miniconda3/envs/vep/share/ensembl-vep-93.2-0/modules/Bio/EnsEMBL/VEP/CacheDir.pm:328STACK Bio::EnsEMBL::VEP::CacheDir::init /home/data/t230420/miniconda3/envs/vep/share/ensembl-vep-93.2-0/modules/Bio/EnsEMBL/VEP/CacheDir.pm:227STACK Bio::EnsEMBL::VEP::CacheDir::new /home/data/t230420/miniconda3/envs/vep/share/ensembl-vep-93.2-0/modules/Bio/EnsEMBL/VEP/CacheDir.pm:111STACK Bio::EnsEMBL::VEP::AnnotationSourceAdaptor::get_all_from_cache /home/data/t230420/miniconda3/envs/vep/share/ensembl-vep-93.2-0/modules/Bio/EnsEMBL/VEP/AnnotationSourceAdaptor.pm:115STACK Bio::EnsEMBL::VEP::AnnotationSourceAdaptor::get_all /home/data/t230420/miniconda3/envs/vep/share/ensembl-vep-93.2-0/modules/Bio/EnsEMBL/VEP/AnnotationSourceAdaptor.pm:91STACK Bio::EnsEMBL::VEP::BaseRunner::get_all_AnnotationSources /home/data/t230420/miniconda3/envs/vep/share/ensembl-vep-93.2-0/modules/Bio/EnsEMBL/VEP/BaseRunner.pm:175STACK Bio::EnsEMBL::VEP::Runner::init /home/data/t230420/miniconda3/envs/vep/share/ensembl-vep-93.2-0/modules/Bio/EnsEMBL/VEP/Runner.pm:123STACK Bio::EnsEMBL::VEP::Runner::run /home/data/t230420/miniconda3/envs/vep/share/ensembl-vep-93.2-0/modules/Bio/EnsEMBL/VEP/Runner.pm:194STACK toplevel /home/data/t230420/miniconda3/envs/vep/bin/vep:224Date (localtime) = Thu Sep 12 09:31:55 2024Ensembl API version = 93---------------------------------------------------ERROR: Failed to run the VEP annotator! Command: /home/data/t230420/miniconda3/envs/vep/bin/perl /home/data/t230420/miniconda3/envs/vep/bin/vep --species homo_sapiens --assembly GRCh37 --offline --no_progress --no_stats --buffer_size 5000 --sift b --ccds --uniprot --hgvs --symbol --numbers --domains --gene_phenotype --canonical --protein --biotype --uniprot --tsl --pubmed --variant_class --shift_hgvs 1 --check_existing --total_length --allele_number --no_escape --xref_refseq --failed 1 --vcf --flag_pick_allele --pick_order canonical,tsl,biotype,rank,ccds,length --dir /home/data/t230420 --fasta /home/data/t230420/hg18.fa --format vcf --input_file 10084_10084.vcf --output_file ./10084_10084.vep.vcf --polyphen b --af --af_1kg --af_esp --af_gnomad --regulatory
It's been determined that the human reference genome I'm going to use is hg18.
The vep version used in the current environment is 93.2.
What can I do to solve this problem?
Thank you.
Metadata
Metadata
Assignees
Labels
No labels