Question: Erorr running HHblits
1
gravatar for max_19
12 months ago by
max_19130
max_19130 wrote:

Hi all,

I'm trying to run HHblits but running into an error. I have several multiple sequence alignments that I have trimmed and concatenated and then converted the resulting .phy file to a .fasta (so HHblits can accept it).

Downloaded the database using:

wget http://www.gwdg.de/~compbiol/uniclust/2018_08/uniclust30_2018_08_hhsuite.tar.gz
tar xzvf uniclust30_2018_08_hhsuite.tar.gz

When I run HHblits though, I keep getting this error message:

$ hhblits -i sample.fasta -o hhblits_results.hhr -d ./databases/uniclust30_2018_08

- 13:51:18.331 ERROR: In /dev/shm/ebuser/avx2/HHsuite/3.0-beta.3/iccifort-2018.3/hh-suite-3.0-beta.3/src/hhdatabase.cpp:35: FFindexDatabase:
- 13:51:18.331 ERROR:   could not open file './databases/uniclust30_2018_08_cs219.ffdata'

my ./databases/uniclust30_2018_08 contains the following files:

uniclust30_2018_08_a3m_db uniclust30_2018_08.cs219 uniclust30_2018_08_hhm_db uniclust30_2018_08_md5sum uniclust30_2018_08_a3m_db.index uniclust30_2018_08_cs219.ffdata uniclust30_2018_08_hhm_db.index uniclust30_2018_08_a3m.ffdata uniclust30_2018_08_cs219.ffindex uniclust30_2018_08_hhm.ffdata uniclust30_2018_08_a3m.ffindex uniclust30_2018_08.cs219.sizes uniclust30_2018_08_hhm.ffindex

Thank you

protein homology hmm • 744 views
ADD COMMENTlink modified 10 months ago by sofimukhtar10 • written 12 months ago by max_19130

I have a similar question as the author of this post at - https://github.com/soedinglab/hh-suite/issues/167

Which link did you download your hhblits database from, and how are they different for the hhblits runs / results? uniclust30_2018_08_hhsuite.tar.gz = 25G, @ http://gwdu111.gwdg.de/~compbiol/uniclust/2018_08/uniclust30_2018_08_hhsuite.tar.gz OR uniclust30_2018_08.tar.gz = 6.6G, @ http://gwdu111.gwdg.de/~compbiol/uniclust/2018_08/uniclust30_2018_08.tar.gz

Since you are running hhblits successfully, I thought you might know :) Could you help please? Thanks!

ADD REPLYlink modified 10 months ago • written 10 months ago by Anand Rao290
3
gravatar for Joe
12 months ago by
Joe16k
United Kingdom
Joe16k wrote:

In the older versions at least, you had to provide the base name of the file itself, not just the directory.

With hhsearch its sufficient to use the /path/to/databasename/databasename.ffdata. With hhblits I believe it needs more than a single file, so you provide the basename of the files.

In your case it may look something like ./databases/uniclust30_2018_08/uniclust30_2018_08.

This is also assuming you’re in the correct directory, since you’re using relative filepaths.

ADD COMMENTlink written 12 months ago by Joe16k
1

That fixed it! I just needed to provide the basename of the files. Thanks very much!

ADD REPLYlink written 12 months ago by max_19130
1

No problem. Learned from bitter experience :P

ADD REPLYlink written 12 months ago by Joe16k
0
gravatar for sofimukhtar1
10 months ago by
sofimukhtar10 wrote:

try to give the full path of the query sequence as well as database

ADD COMMENTlink written 10 months ago by sofimukhtar10

I have a similar question as the author of this post at - https://github.com/soedinglab/hh-suite/issues/167

Which link did you download your hhblits database from, and how are they different for the hhblits runs / results? uniclust30_2018_08_hhsuite.tar.gz = 25G, @ http://gwdu111.gwdg.de/~compbiol/uniclust/2018_08/uniclust30_2018_08_hhsuite.tar.gz OR uniclust30_2018_08.tar.gz = 6.6G, @ http://gwdu111.gwdg.de/~compbiol/uniclust/2018_08/uniclust30_2018_08.tar.gz

Since you are running hhblits successfully, I thought you might know :) Could you help please? Thanks!

ADD REPLYlink written 10 months ago by Anand Rao290
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1623 users visited in the last hour