Results of SIM with:

Sequence 1: Leghemoglobin (153 residues)
Sequence 2: Myoglobin (153 residues)

using the parameters:

Comparison matrix: BLOSUM30
Number of alignments computed: 10
Gap open penalty: 12
Gap extension penalty: 4



16.9% identity in 148 residues overlap; Score: 99.0; Gap frequency: 2.7%

Leghemog       2 ALTESQAALVKSSWEEFNANIPKHTHRFFILVLEIAPAAKDLFSFLKGTSEVPQNNPELQ
Myoglobin      1 VLSEGEWQLVLHVWAKVEADVAGHGQDILIRLFKSHPETLEKFDRFK---HLKTEAEMKA
                  * *    **   *    *    *     *      *     *   *             

Leghemog      62 AHAGKVFKLVYEAAIQLQVTGVVVTDATLKNLGSVHVSK-GVADAHFPVVKEAILKTIKE
Myoglobin     58 SEDLKKHGVTVLTALGAILKKKGHHEAELKPLAQSHATKHKIPIKYLEFISEAIIHVLHS
                     *        *            * ** *   *  *            ***      

Leghemog     121 VVGAKWSEELNSAWTIAYDELAIVIKKE
Myoglobin    118 RHPGDFGADAQGAMNKALELFRKDIAAK
                             *   *       *   


16.3% identity in 104 residues overlap; Score: 51.0; Gap frequency: 2.9%

Leghemog      28 RFFILVLEI-APAAKDLFSF-LKGTSEVPQNNPELQAHAGKVFKLVYEAAIQLQVTGVVV
Myoglobin     45 RFKHLKTEAEMKASEDLKKHGVTVLTALGAILKKKGHHEAELKPLAQSHATKHKIPIKYL
                 **  *  *    *  **                    *      *    *          

Leghemog      86 TDATLKNLGSVHVS-KGVADAHFPVVKEAILKTIKEVVGAKWSE
Myoglobin    105 EFISEAIIHVLHSRHPGDFGADAQGAMNKALELFRKDIAAKYKE
                            *    *   *         *        **  *


17.5% identity in 97 residues overlap; Score: 48.0; Gap frequency: 2.1%

Leghemog      27 HRFFILVLEIAPAAKDLFS-FLKGTSEVPQNNPELQAHAGKVFKLVYEAAIQLQVTGVVV
Myoglobin     12 HVWAKVEADVAGHGQDILIRLFKSHPETLEKFDRFKHLKTEAEMKASEDLKKHGVTVLTA
                 *         *    *      *   *                    *      **    

Leghemog      86 TDATLKNLGSVHVSKGVADAHFPVVKEAILKTIKEVV
Myoglobin     72 LGAILKKKGH-HEAELKPLAQSHATKHKIPIKYLEFI
                   * **  *  *       *     *  *     *  


19.7% identity in 76 residues overlap; Score: 44.0; Gap frequency: 0.0%

Leghemog      73 EAAIQLQVTGVVVTDATLKNLGSVHVSKGVADAHFPVVKEAILKTIKEVVGAKWSEELNS
Myoglobin      4 EGEWQLVLHVWAKVEADVAGHGQDILIRLFKSHPETLEKFDRFKHLKTEAEMKASEDLKK
                 *   **         *     *                *    *  *     * ** *  

Leghemog     133 AWTIAYDELAIVIKKE
Myoglobin     64 HGVTVLTALGAILKKK
                         *    ** 


20.0% identity in 60 residues overlap; Score: 43.0; Gap frequency: 1.7%

Leghemog      58 PELQAHAGK-VFKLVYEAAIQLQVTGVVVTDATLKNLGSVHVSKGVADAHFPVVKEAILK
Myoglobin     88 PLAQSHATKHKIPIKYLEFISEAIIHVLHSRHPGDFGADAQGAMNKALELFRKDIAAKYK
                 *  * ** *      *   *      *                   *   *     *  *


20.0% identity in 30 residues overlap; Score: 39.0; Gap frequency: 0.0%

Leghemog      35 EIAPAAKDLFSFLKGTSEVPQNNPELQAHA
Myoglobin     85 ELKPLAQSHATKHKIPIKYLEFISEAIIHV
                 *  * *       *          *   * 


10.2% identity in 88 residues overlap; Score: 36.0; Gap frequency: 0.0%

Leghemog      21 NIPKHTHRFFILVLEIAPAAKDLFSFLKGTSEVPQNNPELQAHAGKVFKLVYEAAIQLQV
Myoglobin     60 DLKKHGVTVLTALGAILKKKGHHEAELKPLAQSHATKHKIPIKYLEFISEAIIHVLHSRH
                    **          *          **                                

Leghemog      81 TGVVVTDATLKNLGSVHVSKGVADAHFP
Myoglobin    120 PGDFGADAQGAMNKALELFRKDIAAKYK
                  *    **                *   


21.6% identity in 51 residues overlap; Score: 36.0; Gap frequency: 0.0%

Leghemog      58 PELQAHAGKVFKLVYEAAIQLQVTGVVVTDATLKNLGSVHVSKGVADAHFP
Myoglobin     37 PETLEKFDRFKHLKTEAEMKASEDLKKHGVTVLTALGAILKKKGHHEAELK
                 **          *  **               *  **     **   *   


16.0% identity in 25 residues overlap; Score: 35.0; Gap frequency: 0.0%

Leghemog       2 ALTESQAALVKSSWEEFNANIPKHT
Myoglobin     90 AQSHATKHKIPIKYLEFISEAIIHV
                 *              **      * 


31.8% identity in 22 residues overlap; Score: 35.0; Gap frequency: 0.0%

Leghemog      84 VVTDATLKNLGSVHVSKGVADA
Myoglobin    106 FISEAIIHVLHSRHPGDFGADA
                     *    * * *     ***


 SSEARCH searches a sequence database
 using the Smith-Waterman algorithm
 version 3.1t02 March, 1998
Please cite:
 T. F. Smith and M. S. Waterman, (1981) J. Mol. Biol. 147:195-197; 
 W.R. Pearson (1991) Genomics 11:635-650

wwwtmp/.PRSS.4565.1.seq - >Leghemoglobin II
 vs wwwtmp/.PRSS.4565.2.seq - >Myoglobin shuffled sequence

       opt      E()
< 20     0     0:
  22     0     0:           one = represents 1 library sequences
  24     0     0:
  26     0     0:
  28     0     0:
  30     1     0:=
  32     1     1:*
  34     7     3:==*====
  36     7     6:=====*=
  38    13    10:=========*===
  40    10    14:==========   *
  42    11    17:===========     *
  44    20    18:=================*==
  46    26    19:==================*=======
  48    17    18:=================*
  50    14    16:============== *
  52    15    14:=============*=
  54     6    12:======     *
  56     6    10:======   *
  58     7     8:=======*
  60     7     7:======*
  62     8     5:====*===
  64     9     4:===*=====
  66     5     3:==*==
  68     4     3:==*=
  70     2     2:=*
  72     1     2:=*
  74     0     1:*
  76     1     1:*
  78     1     1:*
  80     0     1:*
  82     0     0:
  84     0     0:
  86     0     0:
  88     1     0:=
  90     0     0:
  92     0     0:
  94     0     0:
  96     0     0:
  98     0     0:
 100     0     0:
 102     0     0:
 104     0     0:
 106     0     0:
 108     0     0:
 110     0     0:
 112     0     0:
 114     0     0:
 116     0     0:
 118     0     0:
>120     0     0:
  30753 residues in   201 sequences
 unscaled statistics: mu= 61.3350  var=104.9073
 Kolmogorov-Smirnov  statistic: 0.0449 (N=24) at  38

 Smith-Waterman (3.1 March, 1997) function (/usr/molbio/share/fasta3/blosum30.mat matrix), gap-penalty: -16/-4 reg.-scaled
 Scan time:  0.370

 unshuffled s-w score: 99
For 200 sequences, a score >= 99 is expected 1.005 time