The Perl Toolchain Summit needs more sponsors. If your company depends on Perl, please support this very important event.
hmmpfam - search a single seq against HMM database
HMMER 2.1.1 (Dec 1998)
Copyright (C) 1992-1998 Washington University School of Medicine
HMMER is freely distributed under the GNU General Public License (GPL).
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
HMM file:                 pfam
Sequence file:            /home/birney/src/wise2/example/road.pep
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query:  roa1_drome  

Scores for sequence family classification (score includes all domains):
Model    Description                                    Score    E-value  N 
-------- -----------                                    -----    ------- ---
SEED                                                    146.1    6.3e-40   2

Parsed for domains:
Model    Domain  seq-f seq-t    hmm-f hmm-t      score  E-value
-------- ------- ----- -----    ----- -----      -----  -------
SEED       1/2      33   103 ..     1    77 []    71.2  2.2e-17
SEED       2/2     124   193 ..     1    77 []    75.5  1.1e-18

Alignments of top-scoring domains:
SEED: domain 1 of 2, from 33 to 103: score 71.2, E = 2.2e-17
                   *->lfVgNLppdvteedLkdlFskfGpivsikivrDiiekpketgkskGf
                      lf+g+L + +t+e Lk++F+k G iv++ +++D     + t++s+Gf
  roa1_drome    33    LFIGGLDYRTTDENLKAHFEKWGNIVDVVVMKD-----PRTKRSRGF 74   

                   aFVeFeseedAekAlealnG.kelggrklrv<-*
                   +F+++  ++  + A +    +++++gr+++    
  roa1_drome    75 GFITYSHSSMIDEAQK--SRpHKIDGRVVEP    103  

SEED: domain 2 of 2, from 124 to 193: score 75.5, E = 1.1e-18
                   *->lfVgNLppdvteedLkdlFskfGpivsikivrDiiekpketgkskGf
                      lfVg L  d +e+ ++d+F++fG iv+i+iv+D     ketgk +Gf
  roa1_drome   124    LFVGALKDDHDEQSIRDYFQHFGNIVDINIVID-----KETGKKRGF 165  

                   aFVeFeseedAekAlealnGkelggrklrv<-*
                   aFVeF++++ ++k +     ++l+g+ + v   
  roa1_drome   166 AFVEFDDYDPVDKVVL--KQHQLNGKMVDV    193  

//