Skip to content

TABIX module not working #1874

Open
Open
@kainamill

Description

@kainamill

Describe the issue

Hello, I am currently using ensebl-vep/113 on an HPC system to annotate somatic variants; however, when attempting to use the --custom flag to annotate with the latest version of clinvar I am coming across this error:

MSG: ERROR: Cannot use format vcf without Bio::DB::HTS::Tabix module installed

STACK Bio::EnsEMBL::VEP::AnnotationSource::File::new /gpfs/data/icelake-apps/software/gcc-12.1.0/ensembl-vep/113/modules/Bio/EnsEMBL/VEP/AnnotationSource/File.pm:178
STACK Bio::EnsEMBL::VEP::AnnotationSourceAdaptor::get_all_custom /gpfs/data/icelake-apps/software/gcc-12.1.0/ensembl-vep/113/modules/Bio/EnsEMBL/VEP/AnnotationSourceAdaptor.pm:328
STACK Bio::EnsEMBL::VEP::AnnotationSourceAdaptor::get_all /gpfs/data/icelake-apps/software/gcc-12.1.0/ensembl-vep/113/modules/Bio/EnsEMBL/VEP/AnnotationSourceAdaptor.pm:94
STACK Bio::EnsEMBL::VEP::BaseRunner::get_all_AnnotationSources /gpfs/data/icelake-apps/software/gcc-12.1.0/ensembl-vep/113/modules/Bio/EnsEMBL/VEP/BaseRunner.pm:170
STACK Bio::EnsEMBL::VEP::Runner::init /gpfs/data/icelake-apps/software/gcc-12.1.0/ensembl-vep/113/modules/Bio/EnsEMBL/VEP/Runner.pm:128
STACK Bio::EnsEMBL::VEP::Runner::run /gpfs/data/icelake-apps/software/gcc-12.1.0/ensembl-vep/113/modules/Bio/EnsEMBL/VEP/Runner.pm:200
STACK toplevel /apps/software/gcc-12.1.0/ensembl-vep/113/vep:46
Date (localtime)    = Thu May 15 09:06:34 2025
Ensembl API version = 113
---------------------------------------------------

System

  • VEP version:113

  • VEP Cache version: 113

  • Perl version: 5.36.0

  • tabix installed under: /apps/software/gcc-12.1.0/ensembl-vep/113/htslib/tabix

Full VEP command line

vep -i $sample \
        --force_overwrite \
        --fork 7 \
        --format vcf \
        --hgvs \
        --fasta $fasta \
        --assembly GRCh37 \
        --species homo_sapiens \
        --symbol \
        --cache \
        --cache_version 113 \
        --dir_cache $cache \
        --offline \
        --show_ref_allele \
        --transcript_version \
        --gene_version \
        --af_gnomade \
        --custom file=$clinvar,short_name=ClinVar,format=vcf,type=exact,coords=1,fields=CLNSIG%CLNREVSTAT%CLNDN \
        --no_intergenic \
        --vcf \
        -o ${sample%.vcf.gz}"_anno_CLIN25.vcf.gz"

Data files (if applicable)

They include:

  • clinvar file: clinvar_20250504.vcf.gz
    (clinvar file has .tbi index in same directory)

Metadata

Metadata

Assignees

Labels

No labels
No labels

Type

No type

Projects

No projects

Milestone

No milestone

Relationships

None yet

Development

No branches or pull requests

Issue actions