Results of SIM with:

Sequence 1: Globin (146 residues)
Sequence 2: Leghemog (153 residues)

using the parameters:

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



14.2% identity in 134 residues overlap; Score: 94.0; Gap frequency: 4.5%


Globin         9 QLTADVKKDLRDSWKVIGSDKKGNGVALMTTLFADNQETIGYFKRLGDVSQ-GMANDKLR
Leghemog       2 ALTESQAALVKSSWEEFNANIPKHTHRFFILVLEIAPAAKDLFSFLKGTSEVPQNNPELQ
                  **         **                            *  *   *     *  * 

Globin        68 GHSI-TLMYALQNFID-QLDNPDDLVCVVEKFAVNHITRKISAAEFGKINGPIKKVL---
Leghemog      62 AHAGKVFKLVYEAAIQLQVTGVVVTDATLKNLGSVHVSKGVADAHFPVVKEAILKTIKEV
                  *            *  *                 *       * *      * *     

Globin       123 ASKNFGDKYANAWA
Leghemog     122 VGAKWSEELNSAWT
                            ** 


13.2% identity in 106 residues overlap; Score: 60.0; Gap frequency: 1.9%

Globin        41 FADNQETIGYFKRLGDVSQGMANDKLRGHSITLMYALQNFIDQLDNPDDLVCVVEKFAVN
Leghemog      46 FLKGTSEVPQNNPELQAHAGKVFKLVYEAAIQLQVTGVVVTDATLKNLGSVHVSKGVADA
                 *                  *          * *        *        * *    *  

Globin       101 HITRKISAAEFGKINGPIKKVLASKNFGDKYANAWAKLVAVVQAAL
Leghemog     106 HF--PVVKEAILKTIKEVVGAKWSEELNSAWTIAYDELAIVIKKEM
                 *           *          *         *   *  *     


36.4% identity in 33 residues overlap; Score: 52.0; Gap frequency: 0.0%

Globin         3 VYDAAAQLTADVKKDLRDSWKVIGSDKKGNGVA
Leghemog      71 VYEAAIQLQVTGVVVTDATLKNLGSVHVSKGVA
                 ** ** **            *  **     ***


32.1% identity in 28 residues overlap; Score: 41.0; Gap frequency: 0.0%

Globin       115 NGPIKKVLASKNFGDKYANAWAKLVAVV
Leghemog      56 NNPELQAHAGKVFKLVYEAAIQLQVTGV
                 * *     * * *   *  *    *  *


19.7% identity in 71 residues overlap; Score: 38.0; Gap frequency: 2.8%

Globin        12 ADVKKDLRDSWKVIGSDKKGNGV--ALMTTLFADNQETIGYFKRLGDVSQGMANDKLRGH
Leghemog      37 APAAKDLFSFLKGTSEVPQNNPELQAHAGKVFKLVYEAAIQLQVTGVVVTDATLKNLGSV
                 *   ***    *        *    *     *    *        * *        *   

Globin        70 SITLMYALQNF
Leghemog      97 HVSKGVADAHF
                       *   *


19.0% identity in 21 residues overlap; Score: 38.0; Gap frequency: 0.0%

Globin       117 PIKKVLASKNFGDKYANAWAK
Leghemog     108 PVVKEAILKTIKEVVGAKWSE
                 *  *    *         *  


16.5% identity in 79 residues overlap; Score: 37.0; Gap frequency: 0.0%

Globin        40 LFADNQETIGYFKRLGDVSQGMANDKLRGHSITLMYALQNFIDQLDNPDDLVCVVEKFAV
Leghemog      18 FNANIPKHTHRFFILVLEIAPAAKDLFSFLKGTSEVPQNNPELQAHAGKVFKLVYEAAIQ
                   *        *  *       * *       *      *   *         * *    

Globin       100 NHITRKISAAEFGKINGPI
Leghemog      78 LQVTGVVVTDATLKNLGSV
                    *         *  *  


21.4% identity in 28 residues overlap; Score: 35.0; Gap frequency: 0.0%

Globin        88 DDLVCVVEKFAVNHITRKISAAEFGKIN
Leghemog     104 DAHFPVVKEAILKTIKEVVGAKWSEELN
                 *    **       *     *      *


23.1% identity in 13 residues overlap; Score: 30.0; Gap frequency: 0.0%

Globin        12 ADVKKDLRDSWKV
Leghemog     124 AKWSEELNSAWTI
                 *     *   *  


16.7% identity in 12 residues overlap; Score: 30.0; Gap frequency: 0.0%

Globin       131 YANAWAKLVAVV
Leghemog      11 VKSSWEEFNANI
                     *    *  



 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.28602.1.seq - >Globin I (Dimeric hemoglobin)
 vs wwwtmp/.PRSS.28602.2.seq - >Leghemoglobin II 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     0     0:
  32     1     1:*
  34     5     3:==*==
  36     4     6:==== *
  38    13    10:=========*===
  40    10    14:==========   *
  42    25    17:================*========
  44    16    18:================ *
  46    17    19:================= *
  48    21    18:=================*===
  50    15    16:===============*
  52    17    14:=============*===
  54    10    12:========== *
  56     7    10:=======  *
  58    10     8:=======*==
  60     5     7:===== *
  62     7     5:====*==
  64     4     4:===*
  66     1     3:= *
  68     5     3:==*==
  70     1     2:=*
  72     1     2:=*
  74     1     1:*
  76     1     1:*
  78     0     1:*
  80     0     1:*
  82     0     0:
  84     0     0:
  86     0     0:
  88     1     0:=
  90     2     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= 59.2350  var=148.6631
 Kolmogorov-Smirnov  statistic: 0.0379 (N=24) at  42

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

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