loftee icon indicating copy to clipboard operation
loftee copied to clipboard

DBD::SQLite::db prepare failed: no such table: phylocsf_summary

Open NLSVTN opened this issue 5 years ago • 0 comments

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.

NLSVTN avatar Feb 20 '20 21:02 NLSVTN