loftee
loftee copied to clipboard
DBD::SQLite::db prepare failed: no such table: phylocsf_summary
I am trying to run loftee as a part of VEP with the following command:
/vep/ensembl-tools-release-95/vep -i SF0637658_WES_CIDR.vcf --format vcf --json --everything --allele_number --no_stats --cache --offline --minimal --verbose --assembly GRCh38 --dir_cache /vep/vep_cache --fasta /vep/homo_sapiens/95_GRCh38/hg38.fa --plugin LoF,loftee_path:/vep/loftee_grch38,gerp_bigwig:/vep/loftee_data_grch38/gerp_conservation_scores.homo_sapiens.GRCh38.bw,human_ancestor_fa:/vep/loftee_data_grch38/human_ancestor.fa.gz,filter_position:0.05,min_intron_size:15,conservation_file:/vep/loftee_data_grch38/loftee.sql,run_splice_predictions:0 --dir_plugins /vep/loftee_grch38 -o vep_output
The error that I am getting seems to be all about conservation_file:
DBD::SQLite::db prepare failed: no such table: phylocsf_summary at /var/lib/spark/vep/loftee_grch38/LoF.pm line 553, <ANONIO> line 15156.
I tried supplying as conservation file phylocsf_gerp.sql.gz, phylocsf_gerp.sql, loftee.sql but none really works.
How could it be fixed?
I have also installed PhyloCSF and placed it on the PATH. Not sure whether it is needed. The version of VEP that I am using is 99. I checked out the correct grch38 loftee branch, and not just master. My OS is the following:
(py37) -bash-4.2$ lsb_release -a
LSB Version: :core-4.1-amd64:core-4.1-noarch
Distributor ID: RedHatEnterpriseServer
Description: Red Hat Enterprise Linux Server release 7.7 (Maipo)
Release: 7.7
Codename: Maipo
I believe it is the problem of loftee.sql database. From README it seems to me that the only thing that I need to do here is to use loftee.sql, but maybe there is smth undocumented regarding that.