[TOC] [Prev] [Next]

Accurate Statistics for Smith-Waterman


If the statistical estimates calculated from the database search are accurate, then (1) there should be good agreement between the actual distribution of similarity scores and the predicted distribution (Fig. 4, below) and (2) the highest scoring unrelated sequence should have an expectation value between 0.5 - 2.0. In the table below, the highest scoring unrelated sequence has an expectation value of 4.4.



In the table below, >>...<< indicate high-scoring, unrelated sequences; ** indicate related sequences with expectations values < 11 that were detected with SSEARCH but not FASTA, ktup=2.

The best scores are:                                  s-w Z-score E(43470)
GTB1_MOUSE GLUTATHIONE S-TRANSFERASE GT8.7 (EC 2.5   1490 2006.4   0
GTM1_HUMAN GLUTATHIONE S-TRANSFERASE HB SUBUNIT 4    1235 1661.9   0
GT2_CHICK GLUTATHIONE S-TRANSFERASE 2 (EC 2.5.1.18    954 1282.1   0
GT26_FASHE GLUTATHIONE S-TRANSFERASE 26 KD (EC 2.5    703  942.9   0
GTP_PIG GLUTATHIONE S-TRANSFERASE P (EC 2.5.1.18)     371  494.8 4.1e-21
GTP_MOUSE GLUTATHIONE S-TRANSFERASE P (EC 2.5.1.18    361  481.2 2.3e-20
GTP_HUMAN GLUTATHIONE S-TRANSFERASE P (EC 2.5.1.18    356  474.4 5.6e-20
GTP_CAEEL GLUTATHIONE S-TRANSFERASE P (EC 2.5.1.18    325  432.6 1.2e-17
GTA2_MOUSE GLUTATHIONE S-TRANSFERASE YA CHAIN (EC     229  302.2 2.2e-10
SC11_OMMSL S-CRYSTALLIN SL11 (MAJOR LENS POLYPEPTI    217  289.4 1.1e-09
GTA1_MOUSE GLUTATHIONE S-TRANSFERASE GT41A (EC 2.5    218  287.4 1.5e-09
GTC_MOUSE GLUTATHIONE S-TRANSFERASE YC (EC 2.5.1.1    215  283.4 2.4e-09
SC3_OCTDO S-CRYSTALLIN 3 (OL3).                       210  278.2 4.8e-09
GTH2_HUMAN GLUTATHIONE S-TRANSFERASE A2-2 (EC 2.5.    199  261.7 3.9e-08
GT28_SCHMA GLUTATHIONE S-TRANSFERASE 28 KD (EC 2.5    196  258.1 6.2e-08
GT5A_MOUSE GLUTATHIONE S-TRANSFERASE GST 5.7 (EC 2    183  240.1 6.3e-07
GT28_SCHJA GLUTATHIONE S-TRANSFERASE 28 KD (EC 2.5    169  221.9 6.5e-06
GT2_DROME GLUTATHIONE S-TRANSFERASE 2 (EC 2.5.1.18    164  213.4 1.9e-05
SC1_OCTVU S-CRYSTALLIN 1.                             159  209.0 3.4e-05
GTAC_CHICK GLUTATHIONE S-TRANSFERASE, CL-3 SUBUNIT    144  187.1 0.0006
SC18_OMMSL S-CRYSTALLIN SL18.                         131  163.0 0.01
GT1_MUSDO GLUTATHIONE S-TRANSFERASE 1 (EC 2.5.1.18    122  158.3 0.02
GTTR_RAT GLUTATHIONE S-TRANSFERASE YRS-YRS (EC 2.5    123  158.1 0.02
**GT1_MAIZE GLUTATHIONE S-TRANSFERASE I (EC 2.5.1.18  120  155.3 0.03 **
ARP_TOBAC AUXIN-REGULATED PROTEIN (STR246C PROTEIN    117  151.0 0.06
**GT32_MAIZE GLUTATHIONE S-TRANSFERASE III (EC 2.5.1  115  148.2 0.08 **
GT1_DROME GLUTATHIONE S-TRANSFERASE 1-1 (EC 2.5.1.    100  128.5  1.0
GT_PROMI GLUTATHIONE S-TRANSFERASE GST-6.0 (EC 2.5     97  124.7  1.7
**DCMA_METSP DICHLOROMETHANE DEHALOGENASE (EC 4.5.1.   98  122.7  2.6 **
GTY2_ISSOR GLUTATHIONE S-TRANSFERASE Y-2 (EC 2.5.1     94  121.3  2.6
ARP2_TOBAC AUXIN-INDUCED PROTEIN PGNT35/PCNT111.       93  118.4  3.8
GTT1_RAT GLUTHATHIONE S-TRANSFERASE 5 (EC 2.5.1.18     93  117.8  4.1
>>MOD5_YEAST TRNA ISOPENTENYLTRANSFERASE (EC 2.5.1.8  100  117.2  4.4 <<
GT2_WHEAT GLUTATHIONE S-TRANSFERASE 2 (EC 2.5.1.18     92  114.5  6.2
>>MYSP_MOUSE MYOSIN HEAVY CHAIN, PERINATAL SKELETAL    81  113.5  7.0 <<
LIGE_PSEPA BETA-ETHERASE (BETA-ARYL ETHER CLEAVING     91  113.5  7.0
>>YFHE_ECOLI HYPOTHETICAL 20.1 KD PROTEIN IN HSCA 5'   86  113.5  7.1 <<
**EF1G_HUMAN ELONGATION FACTOR 1-GAMMA (EF-1-GAMMA)    94  113.3  7.2 **
GT_ECOLI GLUTATHIONE S-TRANSFERASE (EC 2.5.1.18).      88  112.7  7.9
>>ABF2_YEAST ARS-BINDING FACTOR 2 PRECURSOR.           87  112.2  8.3 <<
>>KKQ1_YEAST PROBABLE SERINE/THREONINE-PROTEIN KINASE  92  110.7 10.1 <<
**EF1G_RABIT ELONGATION FACTOR 1-GAMMA (EF-1-GAMMA).   92  110.6 10.2 **
ARP3_TOBAC AUXIN-INDUCED PROTEIN PCNT103.              87  110.3 10.6

[TOC] [Prev] [Next]


wrp@virginia.edu