Results of SIM with:

Sequence 1: Globin (146 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


23.8% identity in 126 residues overlap; Score: 130.0; Gap frequency: 5.6%

Globin        22 WKVIGSDKKGNGVALMTTLFADNQETIGYFKRLGDV---SQGMANDKLRGHSITLMYALQ
Myoglobin     14 WAKVEADVAGHGQDILIRLFKSHPETLEKFDRFKHLKTEAEMKASEDLKKHGVTVLTALG
                 *     *  * *      **    **   * *           *   *  *  *   ** 

Globin        79 NFIDQLDNPDDLVCVVEKFAVNHITR-KISAAEFGKINGPIKKVLASKNFGDKYANAWAK
Myoglobin     74 AILKKKGHHEAE---LKPLAQSHATKHKIPIKYLEFISEAIIHVLHSRHPGDFGADAQGA
                                    *  * *  **       *   *  ** *   **  * *   

Globin       138 LVAVVQ
Myoglobin    131 MNKALE
                       


17.7% identity in 79 residues overlap; Score: 69.0; Gap frequency: 2.5%

Globin        60 GMANDKLRGHSITLMYALQNFIDQLDNPDDLVCVVEKFAVNHITRKISAAEFGKINGPIK
Myoglobin     73 GAILKKKGHHEAELKPLAQSHATKHKIPIKYLEFISE-AIIHVLHSRHPGDFGA-DAQGA
                 *    *   *   *    *        *          *  *         **       

Globin       120 KVLASKNFGDKYANAWAKL
Myoglobin    131 MNKALELFRKDIAAKYKEL
                    *   *    *     *


23.2% identity in 56 residues overlap; Score: 48.0; Gap frequency: 7.1%

Globin        20 DSWKVIGSDKKGNGVALMTTLFADNQETIGYFKRLGDVSQGMANDKLRGHSITLMY
Myoglobin     52 EAEMKASEDLKKHGVTVLTALGAILKKKGHHEAELKPLAQSHATK----HKIPIKY
                         * *  **   * * *           *    *  *      * *   *


38.5% identity in 13 residues overlap; Score: 39.0; Gap frequency: 0.0%

Globin        54 LGDVSQGMANDKL
Myoglobin    123 FGADAQGAMNKAL
                  *   **  *  *


50.0% identity in 10 residues overlap; Score: 35.0; Gap frequency: 0.0%

Globin       134 AWAKLVAVVQ
Myoglobin     13 VWAKVEADVA
                  ***  * * 


13.3% identity in 83 residues overlap; Score: 35.0; Gap frequency: 1.2%

Globin        59 QGMANDKLRGHSITLMYALQNFIDQLDNPDDLVCVVEKFAVNHITRKISA-AEFGKINGP
Myoglobin     26 QDILIRLFKSHPETLEKFDRFKHLKTEAEMKASEDLKKHGVTVLTALGAILKKKGHHEAE
                 *         *  **                      *  *   *         *     

Globin       118 IKKVLASKNFGDKYANAWAKLVA
Myoglobin     86 LKPLAQSHATKHKIPIKYLEFIS
                  *    *     *          


23.8% identity in 21 residues overlap; Score: 35.0; Gap frequency: 0.0%

Globin         6 AAAQLTADVKKDLRDSWKVIG
Myoglobin    130 AMNKALELFRKDIAAKYKELG
                 *         **     *  *


16.7% identity in 30 residues overlap; Score: 34.0; Gap frequency: 0.0%

Globin        26 GSDKKGNGVALMTTLFADNQETIGYFKRLG
Myoglobin    124 GADAQGAMNKALELFRKDIAAKYKELGYQG
                 * *  *           *           *


24.1% identity in 29 residues overlap; Score: 32.0; Gap frequency: 0.0%

Globin         5 DAAAQLTADVKKDLRDSWKVIGSDKKGNG
Myoglobin     52 EAEMKASEDLKKHGVTVLTALGAILKKKG
                  *      * **         *   *  *


25.8% identity in 31 residues overlap; Score: 32.0; Gap frequency: 0.0%

Globin         5 DAAAQLTADVKKDLRDSWKVIGSDKKGNGVA
Myoglobin    122 DFGADAQGAMNKALELFRKDIAAKYKELGYQ
                 *  *       * *    * *    *  *  


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

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

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