ohmmbuild

 

Wiki

The master copies of EMBOSS documentation are available at http://emboss.open-bio.org/wiki/Appdocs on the EMBOSS Wiki.

Please help by correcting and extending the Wiki pages.

Function

Build HMM

Description

EMBASSY HMMER is a port of the original hmmer v2.2.1 applications written by Sean Eddy.

Algorithm

Please read the Userguide.pdf distributed with the original HMMER and included in the EMBASSY HMMER distribution under the DOCS directory.

Usage

Here is a sample session with ohmmbuild


% ohmmbuild globins50.msf globin.hmm -auto 

hmmbuild - build a hidden Markov model from an alignment
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).
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Training alignment:                ../../data/hmm/globins50.msf
Number of sequences:               50
Search algorithm configuration:    Local  (hmmsw)
S/W aggregate entry probability:   0.50
S/W aggregate exit probability:    0.50
Model construction strategy:       Fast/ad hoc (gapmax 0.50)
Null model used:                   (default)
Prior used:                        (default)
Prior strategy:                    Dirichlet
Sequence weighting method:         G/S/C tree weights
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Determining effective sequence number    ... done. [13]
Weighting sequences heuristically        ... done.
Constructing model architecture          ... done.
Saving count vector file                 ... done. [globins50.out3]
Converting counts to probabilities       ... done.
Setting model name, etc.                 ... 
Constructed a profile HMM (length 147)
Average score:      280.62 bits
Minimum score:       96.03 bits
Maximum score:      340.28 bits
Std. deviation:      58.76 bits

Finalizing model configuration           ... done.
Saving model to file                     ... done. [globin.hmm]
Saving annotated alignment               ... done. [globins50.out2]

Go to the input files for this example
Go to the output files for this example

Command line arguments

Build HMM
Version: EMBOSS:6.4.0.0

   Standard (Mandatory) qualifiers:
  [-sequences]         seqset     (Aligned) protein sequence set filename and
                                  optional format, or reference (input USA)
  [-outfile]           outfile    [*.ohmmbuild] HMMER hidden markov model
                                  output file

   Additional (Optional) qualifiers: (none)
   Advanced (Unprompted) qualifiers:
   -priorfile          infile     HMMER dirichlet prior file (optional)
   -nullfile           infile     HMMER null model file (optional)
   -pamfile            infile     HMMER matrix file (optional)
   -strategy           menu       [L] Select preference (Values: M (multiple);
                                  L (local); G (global))
   -append             boolean    [N] Append to file
   -amino              boolean    [N] Treat as protein
   -nucleic            boolean    [N] Treat as dna
   -archpri            float      [0.85] Architecture prior (Any numeric
                                  value)
   -binary             boolean    [N] Write HMM as binary
   -cstrategy          menu       [F] Select strategy (Values: F (Fast); H
                                  (Hand))
   -fast               boolean    [N] Work in fast mode
   -gapmax             float      [0.5] Fast mode control (Any numeric value)
   -hand               boolean    [N] Specify model by hand
   -idlevel            float      [0.62] Cutoff ID threhold (Any numeric
                                  value)
   -[no]efficiency     boolean    [Y] Be efficient
   -pamweight          float      [20.0] Weighting for PAM (Any numeric value)
   -swentry            float      [0.5] Probability control for local entries
                                  (Any numeric value)
   -swexit             float      [0.5] Probability control for exits (Any
                                  numeric value)
   -more               boolean    [N] Verbosity
   -weighting          menu       [G] Select weighting (Values: B (Blosum); G
                                  (Gerstein/Sonnhammer/Chothia); K
                                  (Krogh/Mitchison); N (None); V (Voronoi))
   -resavefile         outfile    [*.ohmmbuild] Selex alignment output file
                                  (optional)
   -cfile              outfile    [*.ohmmbuild] HMMER emission and transition
                                  count output file (optional)

   Associated qualifiers:

   "-sequences" associated qualifiers
   -sbegin1            integer    Start of each sequence to be used
   -send1              integer    End of each sequence to be used
   -sreverse1          boolean    Reverse (if DNA)
   -sask1              boolean    Ask for begin/end/reverse
   -snucleotide1       boolean    Sequence is nucleotide
   -sprotein1          boolean    Sequence is protein
   -slower1            boolean    Make lower case
   -supper1            boolean    Make upper case
   -sformat1           string     Input sequence format
   -sdbname1           string     Database name
   -sid1               string     Entryname
   -ufo1               string     UFO features
   -fformat1           string     Features format
   -fopenfile1         string     Features file name

   "-outfile" associated qualifiers
   -odirectory2        string     Output directory

   "-resavefile" associated qualifiers
   -odirectory         string     Output directory

   "-cfile" associated qualifiers
   -odirectory         string     Output directory

   General qualifiers:
   -auto               boolean    Turn off prompts
   -stdout             boolean    Write first file to standard output
   -filter             boolean    Read first file from standard input, write
                                  first file to standard output
   -options            boolean    Prompt for standard and additional values
   -debug              boolean    Write debug output to program.dbg
   -verbose            boolean    Report some/full command line options
   -help               boolean    Report command line options and exit. More
                                  information on associated and general
                                  qualifiers can be found with -help -verbose
   -warning            boolean    Report warnings
   -error              boolean    Report errors
   -fatal              boolean    Report fatal errors
   -die                boolean    Report dying program messages
   -version            boolean    Report version number and exit

Qualifier Type Description Allowed values Default
Standard (Mandatory) qualifiers
[-sequences]
(Parameter 1)
seqset (Aligned) protein sequence set filename and optional format, or reference (input USA) Readable set of sequences Required
[-outfile]
(Parameter 2)
outfile HMMER hidden markov model output file Output file <*>.ohmmbuild
Additional (Optional) qualifiers
(none)
Advanced (Unprompted) qualifiers
-priorfile infile HMMER dirichlet prior file (optional) Input file Required
-nullfile infile HMMER null model file (optional) Input file Required
-pamfile infile HMMER matrix file (optional) Input file Required
-strategy list Select preference
M (multiple)
L (local)
G (global)
L
-append boolean Append to file Boolean value Yes/No No
-amino boolean Treat as protein Boolean value Yes/No No
-nucleic boolean Treat as dna Boolean value Yes/No No
-archpri float Architecture prior Any numeric value 0.85
-binary boolean Write HMM as binary Boolean value Yes/No No
-cstrategy list Select strategy
F (Fast)
H (Hand)
F
-fast boolean Work in fast mode Boolean value Yes/No No
-gapmax float Fast mode control Any numeric value 0.5
-hand boolean Specify model by hand Boolean value Yes/No No
-idlevel float Cutoff ID threhold Any numeric value 0.62
-[no]efficiency boolean Be efficient Boolean value Yes/No Yes
-pamweight float Weighting for PAM Any numeric value 20.0
-swentry float Probability control for local entries Any numeric value 0.5
-swexit float Probability control for exits Any numeric value 0.5
-more boolean Verbosity Boolean value Yes/No No
-weighting list Select weighting
B (Blosum)
G (Gerstein/Sonnhammer/Chothia)
K (Krogh/Mitchison)
N (None)
V (Voronoi)
G
-resavefile outfile Selex alignment output file (optional) Output file <*>.ohmmbuild
-cfile outfile HMMER emission and transition count output file (optional) Output file <*>.ohmmbuild
Associated qualifiers
"-sequences" associated seqset qualifiers
-sbegin1
-sbegin_sequences
integer Start of each sequence to be used Any integer value 0
-send1
-send_sequences
integer End of each sequence to be used Any integer value 0
-sreverse1
-sreverse_sequences
boolean Reverse (if DNA) Boolean value Yes/No N
-sask1
-sask_sequences
boolean Ask for begin/end/reverse Boolean value Yes/No N
-snucleotide1
-snucleotide_sequences
boolean Sequence is nucleotide Boolean value Yes/No N
-sprotein1
-sprotein_sequences
boolean Sequence is protein Boolean value Yes/No N
-slower1
-slower_sequences
boolean Make lower case Boolean value Yes/No N
-supper1
-supper_sequences
boolean Make upper case Boolean value Yes/No N
-sformat1
-sformat_sequences
string Input sequence format Any string  
-sdbname1
-sdbname_sequences
string Database name Any string  
-sid1
-sid_sequences
string Entryname Any string  
-ufo1
-ufo_sequences
string UFO features Any string  
-fformat1
-fformat_sequences
string Features format Any string  
-fopenfile1
-fopenfile_sequences
string Features file name Any string  
"-outfile" associated outfile qualifiers
-odirectory2
-odirectory_outfile
string Output directory Any string  
"-resavefile" associated outfile qualifiers
-odirectory string Output directory Any string  
"-cfile" associated outfile qualifiers
-odirectory string Output directory Any string  
General qualifiers
-auto boolean Turn off prompts Boolean value Yes/No N
-stdout boolean Write first file to standard output Boolean value Yes/No N
-filter boolean Read first file from standard input, write first file to standard output Boolean value Yes/No N
-options boolean Prompt for standard and additional values Boolean value Yes/No N
-debug boolean Write debug output to program.dbg Boolean value Yes/No N
-verbose boolean Report some/full command line options Boolean value Yes/No Y
-help boolean Report command line options and exit. More information on associated and general qualifiers can be found with -help -verbose Boolean value Yes/No N
-warning boolean Report warnings Boolean value Yes/No Y
-error boolean Report errors Boolean value Yes/No Y
-fatal boolean Report fatal errors Boolean value Yes/No Y
-die boolean Report dying program messages Boolean value Yes/No Y
-version boolean Report version number and exit Boolean value Yes/No N

Input file format

ohmmbuild reads any normal sequence USAs.

Input files for usage example

File: globins50.msf



   MSF:  171  Type: P    Check:  4694   .. 

 Name: GLB2_MORMR      oo  Len:  171  Check:  6522  Weight: 2.7687
 Name: GLBZ_CHITH      oo  Len:  171  Check:  6733  Weight: 2.9329
 Name: HBA2_BOSMU      oo  Len:  171  Check:  5006  Weight: 0.6394
 Name: HBA2_GALCR      oo  Len:  171  Check:  6652  Weight: 0.5183
 Name: HBA4_SALIR      oo  Len:  171  Check:  5434  Weight: 1.9511
 Name: HBAD_CHLME      oo  Len:  171  Check:  6621  Weight: 1.2927
 Name: HBAD_PASMO      oo  Len:  171  Check:  8118  Weight: 1.2927
 Name: HBAZ_HORSE      oo  Len:  171  Check:  8382  Weight: 1.6223
 Name: HBA_AILME       oo  Len:  171  Check:  5402  Weight: 0.4145
 Name: HBA_ANSSE       oo  Len:  171  Check:  3688  Weight: 0.8315
 Name: HBA_COLLI       oo  Len:  171  Check:  4420  Weight: 0.8557
 Name: HBA_ERIEU       oo  Len:  171  Check:  5528  Weight: 0.8390
 Name: HBA_FRAPO       oo  Len:  171  Check:  4136  Weight: 0.5014
 Name: HBA_MACFA       oo  Len:  171  Check:  5986  Weight: 0.2233
 Name: HBA_MACSI       oo  Len:  171  Check:  6064  Weight: 0.2233
 Name: HBA_MESAU       oo  Len:  171  Check:  5499  Weight: 0.6722
 Name: HBA_PAGLA       oo  Len:  171  Check:  6189  Weight: 0.5388
 Name: HBA_PHACO       oo  Len:  171  Check:  5129  Weight: 0.5014
 Name: HBA_PONPY       oo  Len:  171  Check:  5894  Weight: 0.3907
 Name: HBA_PROLO       oo  Len:  171  Check:  5810  Weight: 0.4145
 Name: HBA_TRIOC       oo  Len:  171  Check:  6427  Weight: 0.6883
 Name: HBB1_VAREX      oo  Len:  171  Check:  7239  Weight: 1.1252
 Name: HBB2_TRICR      oo  Len:  171  Check:  7790  Weight: 1.9629
 Name: HBB2_XENTR      oo  Len:  171  Check:  9537  Weight: 1.4685
 Name: HBBL_RANCA      oo  Len:  171  Check:  7490  Weight: 1.4685
 Name: HBB_CALAR       oo  Len:  171  Check:  6568  Weight: 0.4226
 Name: HBB_COLLI       oo  Len:  171  Check:  5043  Weight: 0.7672
 Name: HBB_EQUHE       oo  Len:  171  Check:  6101  Weight: 0.6734
 Name: HBB_LARRI       oo  Len:  171  Check:  5673  Weight: 0.7672
 Name: HBB_MANSP       oo  Len:  171  Check:  7148  Weight: 0.4226
 Name: HBB_ORNAN       oo  Len:  171  Check:  6249  Weight: 0.6615
 Name: HBB_RABIT       oo  Len:  171  Check:  7043  Weight: 0.5259
 Name: HBB_SPECI       oo  Len:  171  Check:  3680  Weight: 0.5422
 Name: HBB_SPETO       oo  Len:  171  Check:  4246  Weight: 0.5422
 Name: HBB_SUNMU       oo  Len:  171  Check:  5601  Weight: 0.6734
 Name: HBB_TACAC       oo  Len:  171  Check:  7133  Weight: 0.6615
 Name: HBB_TRIIN       oo  Len:  171  Check:  4009  Weight: 0.8445
 Name: HBB_TUPGL       oo  Len:  171  Check:  7197  Weight: 0.7375
 Name: HBB_URSMA       oo  Len:  171  Check:  7200  Weight: 0.4695
 Name: HBE_PONPY       oo  Len:  171  Check:  5902  Weight: 1.0101
 Name: HBF1_URECA      oo  Len:  171  Check:  6462  Weight: 2.9329
 Name: LGB1_PEA        oo  Len:  171  Check:  4791  Weight: 2.0005
 Name: LGB1_VICFA      oo  Len:  171  Check:  7426  Weight: 2.0005
 Name: MYG_ESCGI       oo  Len:  171  Check:  9170  Weight: 0.7432
 Name: MYG_HORSE       oo  Len:  171  Check:  1290  Weight: 0.7432
 Name: MYG_LYCPI       oo  Len:  171  Check:  1107  Weight: 0.8773


  [Part of this file has been deleted for brevity]

GLB2_MORMR      KLMSMICILL SSAY...... .
GLBZ_CHITH      ATFDAFFGAV FAK....... M
HBA2_BOSMU      KFLANVSTVL TSKYR..... .
HBA2_GALCR      KFMASVSTVL TSKYR..... .
HBA4_SALIR      KFLQQLALAL AEKYR..... .
HBAD_CHLME      KFLSAVAAVL AEKYR..... .
HBAD_PASMO      KFMSAVASVL AEKYR..... .
HBAZ_HORSE      KFLSIVSSVL TEKYR..... .
HBA_AILME       KFFSAVSTVL TSKYR..... .
HBA_ANSSE       KFLCAVATVL TAKYR..... .
HBA_COLLI       KFVLAVGTVL TAKYR..... .
HBA_ERIEU       KFLATVATVL TSKYR..... .
HBA_FRAPO       KFLCAVGNVL TAKYR..... .
HBA_MACFA       KFLASVSTVL TSKYR..... .
HBA_MACSI       KFLASVSTVL TSKYR..... .
HBA_MESAU       KFFASVSTVL TSKYR..... .
HBA_PAGLA       KFFSAVSTVL TSKYR..... .
HBA_PHACO       KFLCAVGTVL TAKYR..... .
HBA_PONPY       KFLASVSTVL TSKYR..... .
HBA_PROLO       KFFTSVSTVL TSKYR..... .
HBA_TRIOC       KFLCAVGNVL SAKYR..... .
HBB1_VAREX      KLVNVVSHSL ARRYH..... .
HBB2_TRICR      KLMKGVSHAL GHGY...... .
HBB2_XENTR      KLNATLVAAL SHGYF..... .
HBBL_RANCA      KFIAVLVDAL SHSYH..... .
HBB_CALAR       KVVAGVANAL AHKYH..... .
HBB_COLLI       KLVRVVAHAL ARKYH..... .
HBB_EQUHE       KVVAGVANAL AHKYH..... .
HBB_LARRI       KLVRVVAHAL ARKYH..... .
HBB_MANSP       KVVAGVANAL AHKYH..... .
HBB_ORNAN       KLVSGVAHAL GHKYH..... .
HBB_RABIT       KVVAGVANAL AHKYH..... .
HBB_SPECI       KVVAGVANAL AHKYH..... .
HBB_SPETO       KVVAGVANAL SHKYH..... .
HBB_SUNMU       KVVAGVANAL AHKYH..... .
HBB_TACAC       KLVSGVSHAL AHKYH..... .
HBB_TRIIN       KVVAGVANAL AHKYH..... .
HBB_TUPGL       KVVAGVANAL AHKYH..... .
HBB_URSMA       KVVAGVANAL AHKYH..... .
HBE_PONPY       KLVSAVAIAL AHKYH..... .
HBF1_URECA      DAAGVLVAAM .......... K
LGB1_PEA        VAYDGLATAI KKAMKT.... A
LGB1_VICFA      VAYDGLATAI K....A.... A
MYG_ESCGI       KALELFRKDI AAKYKELGFQ G
MYG_HORSE       KALELFRNDI AAKYKELGFQ G
MYG_LYCPI       KALELFRNDI AAKYKELGFQ G
MYG_MOUSE       KALELFRNDI AAKYKELGFQ G
MYG_MUSAN       GAFKIICSDI EKEYKAANFQ G
MYG_PROGU       KALELFRNDI AAKYKELGFQ G
MYG_SAISC       KALELFRNDM AAKYKELGFQ G

Output file format

ohmmbuild outputs a graph to the specified graphics device. outputs a report format file. The default format is ...

Output files for usage example

File: globin.hmm

HMMER2.0
NAME  GLB2_MORMR
ACC   
DESC  
LENG  147
ALPH  Amino
RF    yes
CS    yes
MAP   yes
COM   ohmmbuild ../../data/hmm/globins50.msf globin.hmm -auto
NSEQ  50
DATE  Fri Jul 15 12:00:00 2011
CKSUM 8939
XT      -8455     -4      0      *  -8455     -4      0      * 
NULT      -4  -8455
NULE     595  -1558     85    338   -294    453  -1158    197    249    902  -1085   -142    -21   -313     45    531    201    384  -1998   -644 
HMM        A      C      D      E      F      G      H      I      K      L      M      N      P      Q      R      S      T      V      W      Y    
         m->m   m->i   m->d   i->m   i->i   d->m   d->d   b->m   m->e
         -104      *  -3846
     1  -1478  -1780   -402   -384  -1798    746  -1522    195  -1557  -1660   1115  -1800    -97  -1422  -1853     13  -1419   2556  -2187  -1765     4
     P   -143   -495    230     41   -385    395    124   -625    208   -465   -725    274    391     47     92    364    113   -364   -299   -254 
     P  -1004  -1008  -9940  -2163   -365   -701  -1378  -1104  -8190 
     2  -3899  -3461  -5717   -754   -329  -5403  -3143  -1795  -4957   2999  -1134  -4739  -5080  -4188  -4625  -4609  -3788  -2395   3057  -1547    16
     L   -149   -500    233     43   -381    399    106   -626    210   -466   -720    275    394     45     96    359    117   -369   -294   -249 
     L     -9  -8898  -9940   -894  -1115   -701  -1378  -8294  -8185 
     3  -1995  -2562  -4756  -4903    948  -2899  -4002  -4323  -4709  -4557  -3803  -3536  -3665  -4271  -4506   2845   2539  -3465  -4401  -3785    17
     S   -149   -500    233     43   -381    399    106   -626    210   -466   -720    275    394     45     96    359    117   -369   -294   -249 
     S     -9  -8898  -9940   -894  -1115   -701  -1378  -8294  -8180 
     4   1974  -2894   1467   -816  -3309    781  -1179  -3058    963  -3017  -2101  -1116    189   -730  -1296   -776    170  -2597  -3194  -2516    18
     D   -149   -500    233     43   -381    399    106   -626    210   -466   -720    275    394     45     96    359    117   -369   -294   -249 
     D     -9  -8898  -9940   -894  -1115   -701  -1378  -8294  -8175 
     5   1383  -2883    585   1876  -3204    405  -1043  -2954    287  -2899  -1973    -61  -2477    990  -1132   -921   -518   -602  -3067  -2384    19
     A   -149   -500    233     43   -381    399    106   -626    210   -466   -720    275    394     45     96    359    117   -369   -294   -249 
     A     -9  -8964 -10006   -894  -1115   -701  -1378  -8294  -8170 
     6  -3444  -5495   2619   2744  -5608  -3223  -2646  -5588  -3024  -5423  -4830    874  -3800   1490  -3998  -3057  -3535  -5032  -5625  -4510    20
     E   -149   -500    233     43   -381    399    106   -626    210   -466   -720    275    394     45     96    359    117   -369   -294   -249 
     E     -9  -8964 -10006   -894  -1115   -701  -1378  -8294  -8165 
     7  -2536  -2957  -3306  -2413  -3119  -3548  -1859    874   3060  -1039  -2114  -2413   -446  -1563    710  -2671  -2388  -2506   3242  -2796    21
     K   -149   -500    233     43   -381    399    106   -626    210   -466   -720    275    394     45     96    359    117   -369   -294   -249 
     K     -9  -8964 -10006   -894  -1115   -701  -1378  -8294  -8160 
     8   1087  -2878  -1251    152  -3199     57  -1039  -2949   1245  -2894  -1968   1023   -478   1653  -1127    491    694  -2500  -3062  -2379    22
     N   -149   -500    233     43   -381    399    106   -626    210   -466   -720    275    394     45     96    359    117   -369   -294   -249 
     N     -9  -8964 -10006   -894  -1115   -701  -1378  -8294  -8155 
     9    998   -413  -1742   -198  -2292   -825  -1300  -1037    714   1062    420   2005  -2677   -395  -1494  -1541   -425   -661  -2537  -2029    23
     K   -149   -500    233     43   -381    399    106   -626    210   -466   -720    275    394     45     96    359    117   -369   -294   -249 
     K     -9  -8964 -10006   -894  -1115   -701  -1378  -8294  -8150 
    10  -3221  -2748  -5841  -5461  -3140  -5557  -5075   2551  -5338  -1071  -1986  -5198  -5343  -5155  -5383  -4856  -3201   2966  -4560   1098    24
     I   -149   -500    233     43   -381    399    106   -626    210   -466   -720    275    394     45     96    359    117   -369   -294   -249 
     I     -9  -8964 -10006   -894  -1115   -701  -1378  -8294  -8145 
    11      0     29  -1335   -786  -2969  -2418  -1085  -1608   1448    -65  -1823    700  -2509   1617    182    529   1205   -899  -2941  -2302    25


  [Part of this file has been deleted for brevity]

     F   -149   -500    233     43   -381    399    106   -626    210   -466   -720    275    394     45     96    359    117   -369   -294   -249 
     F    -13  -8968 -10010   -894  -1115   -701  -1378  -8294  -7340 
   132  -1696  -3203   2107   1686  -3513    508  -1292  -3276    -99  -3215  -2304   -545  -2698   2093  -1462    197   -974  -2821  -3384  -2676   150
     E   -149   -500    233     43   -381    399    106   -626    210   -466   -720    275    394     45     96    359    117   -369   -294   -249 
     E    -13  -8968 -10010   -894  -1115   -701  -1378  -8294  -7331 
   133   -206  -3207    301  -1311  -3800   -540  -1754  -3505   3357  -3525  -2667  -1586  -3037  -1345  -1850  -1908  -2007    106  -3732  -3091   151
     K   -149   -500    233     43   -381    399    106   -626    210   -466   -720    275    394     45     96    359    117   -369   -294   -249 
     K    -13  -8968 -10010   -894  -1115   -701  -1378  -8294  -7322 
   134   1654  -1721  -4245  -3631   2645  -3462  -2390   -910  -3244   1144   -791  -3131  -3508  -2856  -3053  -2564     87    891  -2223  -1892   152
     L   -149   -500    233     43   -381    399    106   -626    210   -466   -720    275    394     45     96    359    117   -369   -294   -249 
     L    -13  -8968 -10010   -894  -1115   -701  -1378  -8294  -7313 
   135   -137  -1432  -3948  -3313   1581  -3155  -2026   -114  -2909   1428   2225   -648  -3203  -2531  -2710  -2240  -1547   1624  -1889   1109   153
     M   -149   -500    233     43   -381    399    106   -626    210   -466   -720    275    394     45     96    359    117   -369   -294   -249 
     M    -13  -8968 -10010   -894  -1115   -701  -1378  -8294  -7304 
   136   1325   1101   1112    721  -3197   -363  -1036  -2947    702  -1945  -1965   -289  -2471    375   -257   1183  -1037  -2498  -3059  -2377   154
     S   -149   -500    233     43   -381    399    106   -626    210   -466   -720    275    394     45     96    359    117   -369   -294   -249 
     S    -13  -8968 -10010   -894  -1115   -701  -1378  -8294  -7295 
   137   1277  -1469  -3536  -2925  -1430   1509  -1928    601  -2600    216   1186  -1524  -3132   -175  -2535   -508    -85   1083  -1920  -1568   155
     M   -149   -500    233     43   -381    399    106   -626    210   -466   -720    275    394     45     96    359    117   -369   -294   -249 
     M    -14  -8968 -10010   -894  -1115   -701  -1378  -8294  -7285 
   138  -3610  -3121  -6153  -5690   1786  -5837  -5020   1076  -5518   1386  -1290  -5502  -5406  -4901  -5339  -5127  -3551   2808  -4105  -3986   156
     I   -149   -500    233     43   -381    399    106   -626    210   -466   -720    275    394     45     96    359    117   -369   -294   -249 
     I    -14  -8968 -10010   -894  -1115   -701  -1378  -8294  -7276 
   139   2030   2259  -3344  -2779    473   -514  -2054  -1307  -2515  -1647   -982  -2529  -3154  -2266    681   1110  -1625    881  -2225  -1864   157
     C   -149   -500    233     43   -381    399    106   -626    210   -466   -720    275    394     45     96    359    117   -369   -294   -249 
     C    -14  -8968 -10011   -894  -1115   -701  -1378  -8294  -7267 
   140    502  -2761   -339   -762  -3031   -384   1988    -18   -326  -1254  -1863   1897  -2499   -629  -1172    471   1548  -2356  -2975  -2324   158
     I   -149   -500    233     43   -381    399    106   -626    210   -466   -720    275    394     45     96    359    117   -369   -294   -249 
     I    -14  -8968 -10011   -894  -1115   -701  -1378  -8294  -7257 
   141   2557  -2132   1224  -2787  -2497  -3049  -2499  -1705  -2776   -586  -1671  -2638  -3386  -2536  -2949  -1127  -1982   1773  -2964  -2581   159
     L   -149   -500    233     43   -381    399    106   -626    210   -466   -720    275    394     45     96    359    117   -369   -294   -249 
     L    -14  -8969 -10011   -894  -1115   -701  -1378  -8294  -7248 
   142  -4284  -3705  -6735  -6190  -2138  -6510  -5341   1840  -6011   2724   1723  -6217  -5667  -4910  -5578  -5867  -4158    115  -3948  -4121   160
     L   -149   -500    233     43   -381    399    106   -626    210   -466   -720    275    394     45     96    359    117   -369   -294   -249 
     L    -14  -8884  -9926   -894  -1115   -701  -1378  -8294  -7238 
   143   1959  -2667  -1304    160    207   -440  -1064  -2658    652  -2686  -1806  -1057  -2476   -630  -1169    822   1391  -2278  -2923  -2284   161
     S   -149   -500    233     43   -381    399    106   -626    210   -466   -720    275    394     45     96    359    117   -369   -294   -249 
     S    -14  -8825  -9867   -894  -1115   -701  -1378  -8294  -7229 
   144   1365  -2863  -1251   1036  -3195  -2373   3163  -2932    818  -2872  -1957  -1010  -2467   -558    298    930  -1351  -2492  -3031  -2370   162
     S   -149   -500    233     43   -381    399    106   -626    210   -466   -720    275    394     45     96    359    117   -369   -294   -249 
     S    -14  -8825  -9868   -894  -1115   -701  -1378  -8294  -7219 
   145    360  -3058  -1454    144  -3448    -36  -1205  -3163   3161  -3086  -2193  -1228  -2682   -756   -274   -567  -1600  -2731  -3220  -2602   163
     A   -149   -500    233     43   -381    399    106   -626    210   -466   -720    275    394     45     96    359    117   -369   -294   -249 
     A    -15  -8733  -9776   -894  -1115   -701  -1378  -8294  -7209 
   146  -5001  -4245  -5636  -5890   -650  -5388  -2113  -3919  -5374  -3258    935  -4373  -5399  -4454  -4899  -4809  -4973  -4162  -1378   4785   164
     Y   -149   -500    233     43   -381    399    106   -626    210   -466   -720    275    394     45     96    359    117   -369   -294   -249 
     Y    -15  -8575  -9617   -894  -1115   -701  -1378  -8294  -7200 
   147  -3220  -3923  -3873  -2549   -201  -3790   3575  -4055   2069  -3728  -3058  -2457  -3723  -1277   2773  -3101  -2911  -3829  -3495  -3253   165
     -      *      *      *      *      *      *      *      *      *      *      *      *      *      *      *      *      *      *      *      * 
     -      *      *      *      *      *      *      *  -8294      0 
//

File: globins50.out2

#=AU HMMER 2.1.1
#=SQ GLB2_MORMR 0.7013 GLB2_MORMR - 0..0::0 -
#=SQ GLBZ_CHITH 0.7742 GLBZ_CHITH - 0..0::0 -
#=SQ HBA2_BOSMU 0.1490 HBA2_BOSMU - 0..0::0 -
#=SQ HBA2_GALCR 0.1303 HBA2_GALCR - 0..0::0 -
#=SQ HBA4_SALIR 0.4958 HBA4_SALIR - 0..0::0 -
#=SQ HBAD_CHLME 0.3491 HBAD_CHLME - 0..0::0 -
#=SQ HBAD_PASMO 0.3491 HBAD_PASMO - 0..0::0 -
#=SQ HBAZ_HORSE 0.4250 HBAZ_HORSE - 0..0::0 -
#=SQ HBA_AILME  0.1053 HBA_AILME - 0..0::0 -
#=SQ HBA_ANSSE  0.2173 HBA_ANSSE - 0..0::0 -
#=SQ HBA_COLLI  0.2297 HBA_COLLI - 0..0::0 -
#=SQ HBA_ERIEU  0.2146 HBA_ERIEU - 0..0::0 -
#=SQ HBA_FRAPO  0.1349 HBA_FRAPO - 0..0::0 -
#=SQ HBA_MACFA  0.0579 HBA_MACFA - 0..0::0 -
#=SQ HBA_MACSI  0.0579 HBA_MACSI - 0..0::0 -
#=SQ HBA_MESAU  0.1699 HBA_MESAU - 0..0::0 -
#=SQ HBA_PAGLA  0.1331 HBA_PAGLA - 0..0::0 -
#=SQ HBA_PHACO  0.1349 HBA_PHACO - 0..0::0 -
#=SQ HBA_PONPY  0.1014 HBA_PONPY - 0..0::0 -
#=SQ HBA_PROLO  0.1053 HBA_PROLO - 0..0::0 -
#=SQ HBA_TRIOC  0.1835 HBA_TRIOC - 0..0::0 -
#=SQ HBB1_VAREX 0.2879 HBB1_VAREX - 0..0::0 -
#=SQ HBB2_TRICR 0.5211 HBB2_TRICR - 0..0::0 -
#=SQ HBB2_XENTR 0.4019 HBB2_XENTR - 0..0::0 -
#=SQ HBBL_RANCA 0.4019 HBBL_RANCA - 0..0::0 -
#=SQ HBB_CALAR  0.1119 HBB_CALAR - 0..0::0 -
#=SQ HBB_COLLI  0.1918 HBB_COLLI - 0..0::0 -
#=SQ HBB_EQUHE  0.1787 HBB_EQUHE - 0..0::0 -
#=SQ HBB_LARRI  0.1918 HBB_LARRI - 0..0::0 -
#=SQ HBB_MANSP  0.1119 HBB_MANSP - 0..0::0 -
#=SQ HBB_ORNAN  0.1693 HBB_ORNAN - 0..0::0 -
#=SQ HBB_RABIT  0.1284 HBB_RABIT - 0..0::0 -
#=SQ HBB_SPECI  0.1323 HBB_SPECI - 0..0::0 -
#=SQ HBB_SPETO  0.1323 HBB_SPETO - 0..0::0 -
#=SQ HBB_SUNMU  0.1725 HBB_SUNMU - 0..0::0 -
#=SQ HBB_TACAC  0.1693 HBB_TACAC - 0..0::0 -
#=SQ HBB_TRIIN  0.1918 HBB_TRIIN - 0..0::0 -
#=SQ HBB_TUPGL  0.1918 HBB_TUPGL - 0..0::0 -
#=SQ HBB_URSMA  0.1221 HBB_URSMA - 0..0::0 -
#=SQ HBE_PONPY  0.2395 HBE_PONPY - 0..0::0 -
#=SQ HBF1_URECA 0.7874 HBF1_URECA - 0..0::0 -
#=SQ LGB1_PEA   0.5206 LGB1_PEA - 0..0::0 -
#=SQ LGB1_VICFA 0.5206 LGB1_VICFA - 0..0::0 -
#=SQ MYG_ESCGI  0.1987 MYG_ESCGI - 0..0::0 -
#=SQ MYG_HORSE  0.1987 MYG_HORSE - 0..0::0 -
#=SQ MYG_LYCPI  0.2397 MYG_LYCPI - 0..0::0 -
#=SQ MYG_MOUSE  0.2397 MYG_MOUSE - 0..0::0 -
#=SQ MYG_MUSAN  0.6181 MYG_MUSAN - 0..0::0 -
#=SQ MYG_PROGU  0.2044 MYG_PROGU - 0..0::0 -


  [Part of this file has been deleted for brevity]

GLB2_MORMR  KLMSMICILLSSAY-......
GLBZ_CHITH  ATFDAFFGAVFAK--m.....
HBA2_BOSMU  KFLANVSTVLTSKYR......
HBA2_GALCR  KFMASVSTVLTSKYR......
HBA4_SALIR  KFLQQLALALAEKYR......
HBAD_CHLME  KFLSAVAAVLAEKYR......
HBAD_PASMO  KFMSAVASVLAEKYR......
HBAZ_HORSE  KFLSIVSSVLTEKYR......
HBA_AILME   KFFSAVSTVLTSKYR......
HBA_ANSSE   KFLCAVATVLTAKYR......
HBA_COLLI   KFVLAVGTVLTAKYR......
HBA_ERIEU   KFLATVATVLTSKYR......
HBA_FRAPO   KFLCAVGNVLTAKYR......
HBA_MACFA   KFLASVSTVLTSKYR......
HBA_MACSI   KFLASVSTVLTSKYR......
HBA_MESAU   KFFASVSTVLTSKYR......
HBA_PAGLA   KFFSAVSTVLTSKYR......
HBA_PHACO   KFLCAVGTVLTAKYR......
HBA_PONPY   KFLASVSTVLTSKYR......
HBA_PROLO   KFFTSVSTVLTSKYR......
HBA_TRIOC   KFLCAVGNVLSAKYR......
HBB1_VAREX  KLVNVVSHSLARRYH......
HBB2_TRICR  KLMKGVSHALGHGY-......
HBB2_XENTR  KLNATLVAALSHGYF......
HBBL_RANCA  KFIAVLVDALSHSYH......
HBB_CALAR   KVVAGVANALAHKYH......
HBB_COLLI   KLVRVVAHALARKYH......
HBB_EQUHE   KVVAGVANALAHKYH......
HBB_LARRI   KLVRVVAHALARKYH......
HBB_MANSP   KVVAGVANALAHKYH......
HBB_ORNAN   KLVSGVAHALGHKYH......
HBB_RABIT   KVVAGVANALAHKYH......
HBB_SPECI   KVVAGVANALAHKYH......
HBB_SPETO   KVVAGVANALSHKYH......
HBB_SUNMU   KVVAGVANALAHKYH......
HBB_TACAC   KLVSGVSHALAHKYH......
HBB_TRIIN   KVVAGVANALAHKYH......
HBB_TUPGL   KVVAGVANALAHKYH......
HBB_URSMA   KVVAGVANALAHKYH......
HBE_PONPY   KLVSAVAIALAHKYH......
HBF1_URECA  DAAGVLVAAM-----k.....
LGB1_PEA    VAYDGLATAIKKAMKta....
LGB1_VICFA  VAYDGLATAIK----aa....
MYG_ESCGI   KALELFRKDIAAKYKelgfqg
MYG_HORSE   KALELFRNDIAAKYKelgfqg
MYG_LYCPI   KALELFRNDIAAKYKelgfqg
MYG_MOUSE   KALELFRNDIAAKYKelgfqg
MYG_MUSAN   GAFKIICSDIEKEYKaanfqg
MYG_PROGU   KALELFRNDIAAKYKelgfqg
MYG_SAISC   KALELFRNDMAAKYKelgfqg

File: globins50.out3

M 0.00 0.00 0.52 0.52 0.00 1.87 0.00 0.77 0.00 0.00 0.70 0.00 0.70 0.00 0.00 0.92 0.00 6.37 0.00 0.00 
M 0.00 0.00 0.00 0.52 0.24 0.00 0.00 0.00 0.00 10.15 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.48 0.00 
M 0.00 0.00 0.00 0.00 1.04 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 6.91 4.43 0.00 0.00 0.00 
M 4.86 0.00 2.25 0.00 0.00 1.89 0.00 0.00 1.47 0.00 0.00 0.00 0.86 0.00 0.00 0.26 0.79 0.00 0.00 0.00 
M 3.31 0.00 1.12 3.79 0.00 1.52 0.00 0.00 0.67 0.00 0.00 0.41 0.00 1.04 0.00 0.17 0.34 0.62 0.00 0.00 
M 0.00 0.00 4.46 5.94 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.04 0.00 1.56 0.00 0.00 0.00 0.00 0.00 0.00 
M 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.56 7.55 0.52 0.00 0.00 0.52 0.00 0.95 0.00 0.00 0.00 1.90 0.00 
M 2.59 0.00 0.00 0.62 0.00 1.14 0.00 0.00 2.01 0.00 0.00 1.37 0.52 1.96 0.00 1.41 1.38 0.00 0.00 0.00 
M 2.35 0.17 0.00 0.52 0.00 0.52 0.00 0.24 1.32 3.30 0.43 3.14 0.00 0.20 0.00 0.00 0.40 0.40 0.00 0.00 
M 0.00 0.00 0.00 0.00 0.00 0.00 0.00 4.51 0.00 0.35 0.00 0.00 0.00 0.00 0.00 0.00 0.00 7.10 0.00 1.04 
M 0.92 0.29 0.00 0.00 0.00 0.00 0.00 0.19 2.45 1.46 0.00 1.02 0.00 1.91 0.70 1.48 2.15 0.43 0.00 0.00 
M 3.05 0.00 0.85 0.17 0.00 1.14 0.00 0.52 0.00 0.00 0.00 1.62 0.00 0.70 0.00 3.34 1.09 0.52 0.00 0.00 
M 1.63 0.00 0.00 0.00 0.21 0.00 0.79 2.90 0.00 3.66 0.00 0.00 0.00 0.00 0.00 0.00 1.12 2.70 0.00 0.00 
M 0.00 0.00 0.00 0.00 1.94 0.00 0.00 0.00 0.00 0.52 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 9.76 0.77 
M 0.62 0.00 1.22 0.48 0.00 7.27 0.00 0.00 0.00 0.79 0.00 0.00 0.00 0.40 0.00 0.86 1.36 0.00 0.00 0.00 
M 0.62 0.00 0.00 0.00 0.00 0.00 0.17 1.74 9.68 0.00 0.00 0.79 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
M 0.00 0.00 0.00 0.00 0.00 0.00 0.00 3.93 0.00 0.56 0.24 0.00 0.00 0.00 0.00 0.00 0.00 8.27 0.00 0.00 
M 0.48 0.00 0.00 1.90 0.00 2.02 0.00 0.00 1.56 1.54 0.00 0.00 0.00 0.00 0.00 0.56 0.18 0.00 0.00 0.70 
M 1.07 0.00 0.00 0.52 0.00 4.09 0.00 0.00 0.70 0.00 0.43 0.00 0.50 0.52 0.00 0.89 0.24 0.00 0.00 0.00 
M 0.35 0.79 3.19 0.00 0.00 0.00 2.38 0.00 1.54 0.00 0.00 3.33 0.00 0.65 0.00 0.00 0.00 0.00 0.00 0.00 
M 3.64 0.00 0.00 0.70 0.00 0.21 0.00 0.97 0.00 2.28 0.00 0.00 0.00 0.35 0.00 0.50 0.00 2.87 0.00 0.70 
M 1.67 0.00 1.96 3.61 0.00 1.59 0.00 0.00 0.19 0.00 0.00 0.34 1.25 0.79 0.00 0.20 0.62 0.00 0.00 0.00 
M 1.83 0.00 1.16 4.40 0.00 1.08 0.00 0.00 0.89 0.29 0.00 0.00 0.00 0.80 0.00 0.73 0.00 0.00 0.00 0.00 
M 1.03 0.38 1.58 0.00 0.35 0.00 1.29 1.88 0.00 1.09 0.00 0.70 0.00 0.00 0.00 0.00 0.00 1.57 0.00 2.09 
M 1.31 0.00 0.00 0.00 0.00 10.40 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.77 0.52 0.00 0.00 0.00 
M 3.68 0.00 0.79 0.50 0.00 3.08 0.80 0.00 0.00 0.00 0.00 0.00 0.00 1.90 0.00 0.35 0.43 1.48 0.00 0.00 
M 0.00 0.00 1.65 6.92 0.00 0.77 0.00 0.00 1.04 0.00 0.00 0.62 0.00 1.21 0.00 0.79 0.00 0.00 0.00 0.00 
M 6.91 0.52 0.00 0.00 0.00 0.52 0.00 3.08 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.19 0.69 1.09 0.00 0.00 
M 0.00 0.00 0.00 0.00 0.79 0.00 0.00 0.00 0.00 11.69 0.52 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
M 1.19 0.00 0.21 1.69 0.79 2.05 0.00 1.29 0.00 0.62 0.00 0.00 0.00 0.99 0.00 0.90 0.40 0.70 0.35 0.77 
M 0.77 0.00 0.00 0.00 0.00 0.53 0.00 0.00 1.49 0.00 0.00 0.00 0.00 0.00 8.77 0.40 0.00 0.00 0.00 0.00 
M 0.00 0.00 0.00 0.00 0.70 0.00 0.00 0.00 0.00 6.61 2.75 0.00 0.00 0.00 0.00 0.00 0.34 0.77 0.00 0.79 
M 0.00 0.00 0.00 0.00 7.03 0.00 0.00 0.52 0.00 4.41 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
M 0.21 0.00 0.00 0.62 0.17 0.00 0.35 1.79 1.86 0.85 0.00 0.00 0.00 0.43 0.00 0.43 1.69 3.57 0.00 0.00 
M 1.78 0.00 0.00 0.00 0.00 1.51 0.00 0.29 0.00 0.00 0.00 0.24 0.00 0.62 0.00 1.61 1.49 4.43 0.00 0.00 
M 0.00 0.00 0.77 0.00 1.06 0.00 1.50 0.00 0.00 0.00 0.00 0.52 0.00 0.00 0.00 0.00 0.70 0.00 0.00 7.40 
M 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 11.96 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
M 0.70 0.00 0.00 1.90 0.00 0.79 0.00 0.00 0.00 0.00 0.00 0.00 0.18 1.72 0.00 0.62 2.00 0.00 4.05 0.00 
M 0.70 0.00 0.79 0.00 0.00 0.00 0.00 0.77 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.14 8.56 0.00 0.00 0.00 
M 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 3.74 2.07 0.00 0.00 0.00 5.62 0.52 0.00 0.00 0.00 0.00 0.00 
M 2.76 0.00 0.48 0.81 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.62 0.00 0.00 4.05 0.00 3.25 0.00 0.00 0.00 
M 0.77 0.00 0.00 0.00 3.92 0.00 0.62 0.00 1.29 0.00 0.00 0.00 0.00 0.29 0.00 0.00 0.00 0.00 0.00 5.07 
M 0.00 0.00 0.00 0.00 11.96 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
M 0.19 0.00 2.57 0.77 0.00 0.00 1.31 0.00 0.00 0.00 0.00 0.00 5.34 0.00 0.00 1.78 0.00 0.00 0.00 0.00 
M 0.17 0.00 0.52 0.00 0.00 0.00 4.22 0.00 3.39 0.00 0.00 0.00 0.00 0.77 0.00 2.48 0.40 0.00 0.00 0.00 
M 0.00 0.00 0.00 0.00 12.50 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.50 0.00 
M 0.50 0.00 3.25 0.00 0.00 4.05 0.00 0.00 2.60 0.00 0.00 0.00 0.00 0.00 0.00 1.83 0.00 0.77 0.00 0.00 
M 0.00 0.00 2.74 0.00 1.04 1.48 1.05 0.00 0.00 0.00 0.00 2.17 0.00 0.00 0.00 1.28 0.00 0.00 0.00 0.00 
M 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.77 10.15 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.45 0.00 0.00 
M 0.50 0.00 0.00 0.00 0.00 0.00 0.96 0.00 2.94 0.00 0.00 0.21 0.00 0.22 0.00 5.91 0.70 0.00 0.00 0.00 


  [Part of this file has been deleted for brevity]

T 13.00 0.00 0.00 0.00 0.00 0.00 0.00 
T 13.00 0.00 0.00 0.00 0.00 0.00 0.00 
T 13.00 0.00 0.00 0.00 0.00 0.00 0.00 
T 13.00 0.00 0.00 0.00 0.00 0.00 0.00 
T 13.00 0.00 0.00 0.00 0.00 0.00 0.00 
T 11.96 0.00 1.04 0.00 0.00 0.00 0.00 
T 11.96 0.00 0.00 0.00 0.00 1.04 0.00 
T 13.00 0.00 0.00 0.00 0.00 0.00 0.00 
T 13.00 0.00 0.00 0.00 0.00 0.00 0.00 
T 13.00 0.00 0.00 0.00 0.00 0.00 0.00 
T 13.00 0.00 0.00 0.00 0.00 0.00 0.00 
T 13.00 0.00 0.00 0.00 0.00 0.00 0.00 
T 13.00 0.00 0.00 0.00 0.00 0.00 0.00 
T 13.00 0.00 0.00 0.00 0.00 0.00 0.00 
T 13.00 0.00 0.00 0.00 0.00 0.00 0.00 
T 12.30 0.00 0.70 0.00 0.00 0.00 0.00 
T 12.30 0.00 0.00 0.00 0.00 0.00 0.70 
T 12.30 0.00 0.00 0.00 0.00 0.00 0.70 
T 12.30 0.00 0.00 0.00 0.00 0.00 0.70 
T 12.30 0.00 0.00 0.00 0.00 0.00 0.70 
T 12.30 0.00 0.00 0.00 0.00 0.00 0.70 
T 12.30 0.00 0.00 0.00 0.00 0.00 0.70 
T 12.30 0.00 0.00 0.00 0.00 0.00 0.70 
T 10.74 0.00 1.56 0.00 0.00 0.00 0.70 
T 10.74 0.00 0.00 0.00 0.00 0.70 1.56 
T 11.44 0.00 0.00 0.00 0.00 1.56 0.00 
T 13.00 0.00 0.00 0.00 0.00 0.00 0.00 
T 13.00 0.00 0.00 0.00 0.00 0.00 0.00 
T 13.00 0.00 0.00 0.00 0.00 0.00 0.00 
T 13.00 0.00 0.00 0.00 0.00 0.00 0.00 
T 13.00 0.00 0.00 0.00 0.00 0.00 0.00 
T 13.00 0.00 0.00 0.00 0.00 0.00 0.00 
T 13.00 0.00 0.00 0.00 0.00 0.00 0.00 
T 13.00 0.00 0.00 0.00 0.00 0.00 0.00 
T 13.00 0.00 0.00 0.00 0.00 0.00 0.00 
T 13.00 0.00 0.00 0.00 0.00 0.00 0.00 
T 13.00 0.00 0.00 0.00 0.00 0.00 0.00 
T 13.00 0.00 0.00 0.00 0.00 0.00 0.00 
T 13.00 0.00 0.00 0.00 0.00 0.00 0.00 
T 13.00 0.00 0.00 0.00 0.00 0.00 0.00 
T 13.00 0.00 0.00 0.00 0.00 0.00 0.00 
T 13.00 0.00 0.00 0.00 0.00 0.00 0.00 
T 13.00 0.00 0.00 0.00 0.00 0.00 0.00 
T 13.00 0.00 0.00 0.00 0.00 0.00 0.00 
T 13.00 0.00 0.00 0.00 0.00 0.00 0.00 
T 13.00 0.00 0.00 0.00 0.00 0.00 0.00 
T 12.21 0.00 0.00 0.00 0.00 0.00 0.00 
T 11.69 0.00 0.00 0.00 0.00 0.00 0.00 
T 11.69 0.00 0.00 0.00 0.00 0.00 0.00 
T 10.92 0.00 0.00 0.00 0.00 0.00 0.00 
T 9.70 0.00 0.00 0.00 0.00 0.00 0.00 

Data files

None.

Notes

None.

References

None.

Warnings

None.

Diagnostic Error Messages

None.

Exit status

It always exits with status 0.

Known bugs

None.

See also

Program name Description
ehmmalign Align sequences to an HMM profile
ehmmbuild Build a profile HMM from an alignment
ehmmcalibrate Calibrate HMM search statistics
ehmmconvert Convert between profile HMM file formats
ehmmemit Generate sequences from a profile HMM
ehmmfetch Retrieve an HMM from an HMM database
ehmmindex Create a binary SSI index for an HMM database
ehmmpfam Search one or more sequences against an HMM database
ehmmsearch Search a sequence database with a profile HMM
libgen Generate discriminating elements from alignments
ohmmalign Align sequences with an HMM
ohmmcalibrate Calibrate a hidden Markov model
ohmmconvert Convert between HMM formats
ohmmemit Extract HMM sequences
ohmmfetch Extract HMM from a database
ohmmindex Index an HMM database
ohmmpfam Align single sequence with an HMM
ohmmsearch Search sequence database with an HMM

Author(s)

This program is an EMBOSS conversion of a program written by Sean Eddy as part of his HMMER package.

Please report all bugs to the EMBOSS bug team (emboss-bug © emboss.open-bio.org) not to the original author.

History

Target users

This program is intended to be used by everyone and everything, from naive users to embedded scripts.