t/data/hmmpfam.out

   1 hmmpfam - search a single seq against HMM database
   2 HMMER 2.1.1 (Dec 1998)
   3 Copyright (C) 1992-1998 Washington University School of Medicine
   4 HMMER is freely distributed under the GNU General Public License (GPL).
   5 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
   6 HMM file:                 pfam
   7 Sequence file:            /home/birney/src/wise2/example/road.pep
   8 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
   9 Query:  roa1_drome
  10
  11 Scores for sequence family classification (score includes all domains):
  12 Model    Description                                    Score    E-value  N
  13 -------- -----------                                    -----    ------- ---
  14 SEED                                                    146.1    6.3e-40   2
  15
  16 Parsed for domains:
  17 Model    Domain  seq-f seq-t    hmm-f hmm-t      score  E-value
  18 -------- ------- ----- -----    ----- -----      -----  -------
  19 SEED       1/2      33   103 ..     1    77 []    71.2  2.2e-17
  20 SEED       2/2     124   193 ..     1    77 []    75.5  1.1e-18
  21
  22 Alignments of top-scoring domains:
  23 SEED: domain 1 of 2, from 33 to 103: score 71.2, E = 2.2e-17
  24                    *->lfVgNLppdvteedLkdlFskfGpivsikivrDiiekpketgkskGf
  25                       lf+g+L + +t+e Lk++F+k G iv++ +++D     + t++s+Gf
  26   roa1_drome    33    LFIGGLDYRTTDENLKAHFEKWGNIVDVVVMKD-----PRTKRSRGF 74
  27
  28                    aFVeFeseedAekAlealnG.kelggrklrv<-*
  29                    +F+++  ++  + A +    +++++gr+++
  30   roa1_drome    75 GFITYSHSSMIDEAQK--SRpHKIDGRVVEP    103
  31
  32 SEED: domain 2 of 2, from 124 to 193: score 75.5, E = 1.1e-18
  33                    *->lfVgNLppdvteedLkdlFskfGpivsikivrDiiekpketgkskGf
  34                       lfVg L  d +e+ ++d+F++fG iv+i+iv+D     ketgk +Gf
  35   roa1_drome   124    LFVGALKDDHDEQSIRDYFQHFGNIVDINIVID-----KETGKKRGF 165
  36
  37                    aFVeFeseedAekAlealnGkelggrklrv<-*
  38                    aFVeF++++ ++k +     ++l+g+ + v
  39   roa1_drome   166 AFVEFDDYDPVDKVVL--KQHQLNGKMVDV    193
  40
  41 //