HMMER: Failed to open target sequence database for reading
1
0
Entering edit mode
6.9 years ago
Cayro • 0

Good Afternoon

I'm trying to run a protein search through HMMER (phmmer) against PFam database in command line. However the process is interrupted and the following message appears: Failed to open target database sequence HMMER/database/Pfam-A.fasta for reading. I would like to know what's happening.

HMMER phmmer Database • 8.2k views
ADD COMMENT
1
Entering edit mode

Hi cayro

Did you prepare your database with hmmpress?

hmmpress creates four binary files. The database cannot be read because it needs to have those four auxiliary binary files generated by hmmpress.

ADD REPLY
1
Entering edit mode

hmmpress is for HMM files that you're going to search with hmmscan; not for FASTA files that'll be searched with phmmer.

ADD REPLY
0
Entering edit mode

How could I run this database? Because I tried with the uniprot database in fasta format and it worked. Would I have to change the format of my bank?

ADD REPLY
0
Entering edit mode

I tried this approach, however the following message appears:

Error: File format problem in trying to open HMM file database/Pfam/Pfam-A.fasta. Format tag is '>S6ENI5_9CLOT/170-197': unrecognized. Current H3 format is 'HMMER3/f'. Previous H2/H3 formats also supported.

ADD REPLY
1
Entering edit mode

Pfam-A.fasta is a FASTA sequence file. It is not a profile HMM database. (Have you looked at it? If not, look at it.)

ADD REPLY
0
Entering edit mode

It's a fasta sequence. I can't understand why it does not work. Should I change the format of my database?

ADD REPLY
0
Entering edit mode

Hi, Did you ever find a solution to this problem? I am getting the same exact error message.

ADD REPLY
1
Entering edit mode
6.9 years ago

It means that your file HMMER/database/Pfam-A.fasta can't be opened for reading. It either doesn't exist, or you don't have read permission for it.

Also, if you are trying to search a protein sequence against Pfam, you generally want to be using hmmscan or hmmsearch, not phmmer. phmmer is for searching sequence(s) against a sequence database (like UniProt or NCBI NR).

ADD COMMENT

Login before adding your answer.

Traffic: 4961 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6