Question: define PFAM dir in get_homologues perl script
gravatar for dago
5.6 years ago by
dago2.6k wrote:

I think this is really a thrivial question but I am stuck on it.

I would like to use get_homolohues to perform a comaprative genomic analysis. It requires some dependencies, such as MCL, Blast, hmmer. I am using Biolinux 8 and most of them are installed already.

I follow the installation instructions. When checking the binaries needed I got the following:

 perl -v

Checking required binaries and data sources, all set in :
        EXE_BLASTP : OK (path:/usr/local/lib/get_homologues-x86_64-20140930/bin/ncbi-blast-2.2.27+/bin/blastp)
        EXE_BLASTN : OK (path:/usr/local/lib/get_homologues-x86_64-20140930/bin/ncbi-blast-2.2.27+/bin/blastn)
      EXE_FORMATDB : OK (path:/usr/local/lib/get_homologues-x86_64-20140930/bin/ncbi-blast-2.2.27+/bin/makeblastdb)
           EXE_MCL : OK (path:/usr/local/lib/get_homologues-x86_64-20140930//bin/mcl-02-063/shmcl/mcl)
      EXE_MAKEHASH : OK (path:/usr/local/lib/get_homologues-x86_64-20140930//bin/COGsoft/COGmakehash/COGmakehash )
     EXE_READBLAST : OK (path:/usr/local/lib/get_homologues-x86_64-20140930//bin/COGsoft/COGreadblast/COGreadblast )
        EXE_COGLSE : OK (path:/usr/local/lib/get_homologues-x86_64-20140930//bin/COGsoft/COGlse/COGlse )
        EXE_COGTRI : OK (path:/usr/local/lib/get_homologues-x86_64-20140930//bin/COGsoft/COGtriangles/COGtriangles )
       EXE_HMMPFAM : Pfam-A.hmm file needs to be installed (

        EXE_INPARA : OK (path:/usr/local/lib/get_homologues-x86_64-20140930/
         EXE_ORTHO : OK (path:/usr/local/lib/get_homologues-x86_64-20140930/
         EXE_HOMOL : OK (path:/usr/local/lib/get_homologues-x86_64-20140930/
    EXE_SPLITBLAST : OK (path:/usr/local/lib/get_homologues-x86_64-20140930/
  EXE_SPLITHMMPFAM : OK (path:/usr/local/lib/get_homologues-x86_64-20140930/

So I just needed to install PFAM. I dowloaded as follow: 

wget -N

gunzip Pfam-A.hmm.gz


I then changed the directory in the appropriate perl scipt to my actual dir where I have the Pfam-A.hmm, which is


However, when I check again the binaries needed I still get the same message as above, telling me that Pfam-A is not installed.

Anyone can help me with this?

uname -a                                                                                                                                           [ 3:20PM]
Linux bl8vbox 3.13.0-39-generic #66-Ubuntu SMP Tue Oct 28 13:30:27 UTC 2014 x86_64 x86_64 x86_64 GNU/Linux


I downloaded again the Pfam-A.hmm.gz and this time it worked.

I run the with the option -D for Pfam search and I do not get any error. However, when I open the file resulted from the analysis XXX.fasta.pfam I do not have any data, but this message:

Error: Unrecognized format, trying to open hmm file /usr/local/lib/get_homologues-x86_64-20140930/db/Pfam-A.hmm for reading.

Any idea?


hmmer pfam biolinux perl • 2.8k views
ADD COMMENTlink modified 5.6 years ago by b.contreras.moreira10 • written 5.6 years ago by dago2.6k

Did you follow all the instructions given here? Did you install HMMER3? Did you run the commands given below:

Generate binary files for Pfam-A.hmm using 

/your/path/hmmer-3.0/src/hmmpress Pfam-A.hmm

Edit file lib/ and modify the relevant line to:

if( ! defined($ENV{"PFAMDB"}) ){ $ENV{"PFAMDB"} = "db/Pfam-A.hmm"; }
ADD REPLYlink modified 5.6 years ago • written 5.6 years ago by komal.rathi3.6k

Thanks a lot.

Hmmer Is already installed in Biolinux. I followed the instruction as in the manual. I

ADD REPLYlink written 5.6 years ago by dago2.6k
gravatar for b.contreras.moreira
5.6 years ago by
b.contreras.moreira10 wrote:

Hi Dago,

get_homologues is bundled with all the required sources/binaries, so after running install you should be ready to go if all you get are OKs. The script will also take care of pressing the PfamA.hmm file located at the $ENV{"PFAMDB"} location, as mentioned above. I do recommend proceeding this way, as thus the software is using compatible hmmpress and then hmmscan binaries. If you choose to press the hmms yourself then you should check that the hmmpress binary in biolinux is compatible with the bundled hmmscan binary,

hope this helps,


ADD COMMENTlink written 5.6 years ago by b.contreras.moreira10
Please log in to add an answer.


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