meme.lipocalin.tcm

EM算法的改进

--------------------------------------------------------------------------------
	Motif 1 position-specific scoring matrix
--------------------------------------------------------------------------------
log-odds matrix: alength= 20 w= 28 n= 768 bayes= 7.65181 E= 1.5e-004 
  -179   -268   -226   -205      7   -247    380   -213   -160   -202   -113   -104   -202     10    220   -133   -151   -225   -138    188 
  -216   -385    335    -74   -371   -240   -199   -326   -279   -356   -267    -75    167   -146   -194   -181   -222   -327   -393   -369 
   -98    154   -277   -331   -377    372   -289   -326   -338   -381   -261   -219   -262   -255   -206   -136   -211   -296   -366   -401 
   120   -212   -167   -110     -1   -150     -3   -110    -66   -138    -48   -113   -127      6    261    -52    -65   -115   -142    193 
  -242    486   -509   -468   -350   -398   -408   -149   -499   -283   -212   -438   -390   -371   -342   -318   -259    -26   -466   -456 
   -85   -296   -258   -236   -310   -212   -220   -256   -247   -261   -204   -234    440   -128   -144     88   -127   -243   -385   -379 
  -124   -370    238    171   -329   -166   -121   -272   -136   -285   -178    -50   -167      4    -78    117   -109   -258   -360   -304 
  -104   -200   -381   -337    167   -306   -235     90   -348    -55    -14   -316   -248   -226   -204   -223   -110    297   -243   -172 
  -293   -397   -419   -332   -477   -365   -258   -331    325   -382   -273   -289   -326   -156    160   -305   -272   -392   -394   -440 
   -93   -286   -279   -244   -272   -231   -226   -150   -254   -208   -156   -263    435   -135   -148   -123   -137     64   -362   -353 
   -77   -211   -273   -202   -214   -253   -181     59     33   -130    -66   -223   -196    -96    -38   -167    -84    301   -304   -294 
  -192   -391    297    121   -368   -198   -153   -321   -213   -345   -243    113   -236    -74   -148   -134   -178   -317   -388   -346 
  -246   -359    145   -180   -354   -182    -86   -301   -251   -371   -268    356   -246   -124   -178    -92   -162   -348   -357   -344 
  -311   -286   -454   -446    429   -407   -331   -186   -472   -134   -134   -414   -322   -362   -348   -287   -346   -251   -188    -83 
  -299   -415    366   -112   -405   -316   -248   -363   -403   -411   -320    -90   -365   -257   -282   -267   -316   -376   -410   -408 
  -163   -274   -276   -214    -91   -253   -156    -96     62     74    -49   -207   -216    -71    -14   -165   -140   -141    467   -132 
   178   -161   -279   -269   -296   -125   -243   -245   -276   -283   -168   -173   -170   -163   -158    312     59   -196   -343   -340 
   124   -254    -99    -67   -235   -111    -42   -173    -28   -191    -86    169   -112     44    283    -24    -43   -161   -259   -219 
  -279   -277   -412   -397    399   -373   -138   -188   -413   -145   -129   -344   -299   -283   -271   -254   -299   -242   -103    168 
   290   -154   -205   -165   -212   -108    218   -169   -176   -201   -101   -151   -163    -63    -74    -22    -84   -109   -277   -211 
  -175   -369   -278   -341   -423    394   -306   -394   -350   -448   -312   -229   -304   -286   -217   -196   -294   -377   -379   -430 
   148   -221   -133    -68   -218   -132    -79   -137    124   -171    -64    -80   -100     31     44     19    188   -116   -259   -222 
  -381   -381   -432   -434   -127   -378   -344   -357   -437   -238   -232   -390   -373   -296   -268   -367   -354   -351    568   -199 
  -221   -269   -316   -305    101   -302    300   -202   -318   -184   -123   -209   -251   -121   -143   -194   -218   -229    367    261 
    22   -435    -99    328   -422   -272   -243   -326   -226   -357   -253   -226   -244    -62   -173   -200   -214   -311   -453   -414 
  -175   -244   -408   -375   -169   -342   -312    324   -369    -12    249   -332   -301   -254   -248   -250   -146    183   -306   -291 
   341   -132   -343   -301   -283   -113   -294   -208   -334   -248   -145   -289   -258   -221   -205    -24   -122   -116   -332   -368 
  -326   -410   -459   -401   -505   -404   -333   -352    339   -433   -317   -335   -358   -263     42   -362   -320   -436   -421   -487 
--------------------------------------------------------------------------------

--------------------------------------------------------------------------------
	Motif 1 position-specific probability matrix
--------------------------------------------------------------------------------
letter-probability matrix: alength= 20 w= 28 nsites= 3 E= 1.5e-004 
 0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.333333  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.333333  0.000000  0.000000  0.000000  0.000000  0.333333 
 0.000000  0.000000  0.666667  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.333333  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000 
 0.000000  0.333333  0.000000  0.000000  0.000000  0.666667  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000 
 0.333333  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.333333  0.000000  0.000000  0.000000  0.000000  0.333333 
 0.000000  0.666667  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.333333  0.000000  0.000000 
 0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.666667  0.000000  0.000000  0.333333  0.000000  0.000000  0.000000  0.000000 
 0.000000  0.000000  0.333333  0.333333  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.333333  0.000000  0.000000  0.000000  0.000000 
 0.000000  0.000000  0.000000  0.000000  0.333333  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.666667  0.000000  0.000000 
 0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.666667  0.000000  0.000000  0.000000  0.000000  0.000000  0.333333  0.000000  0.000000  0.000000  0.000000  0.000000 
 0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.666667  0.000000  0.000000  0.000000  0.000000  0.333333  0.000000  0.000000 
 0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.333333  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.666667  0.000000  0.000000 
 0.000000  0.000000  0.333333  0.333333  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.333333  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000 
 0.000000  0.000000  0.333333  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.666667  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000 
 0.000000  0.000000  0.000000  0.000000  1.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000 
 0.000000  0.000000  1.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000 
 0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.333333  0.333333  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.333333  0.000000 
 0.333333  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.666667  0.000000  0.000000  0.000000  0.000000 
 0.333333  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.333333  0.000000  0.000000  0.333333  0.000000  0.000000  0.000000  0.000000  0.000000 
 0.000000  0.000000  0.000000  0.000000  0.666667  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.333333 
 0.666667  0.000000  0.000000  0.000000  0.000000  0.000000  0.333333  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000 
 0.000000  0.000000  0.000000  0.000000  0.000000  1.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000 
 0.333333  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.333333  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.333333  0.000000  0.000000  0.000000 
 0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  1.000000  0.000000 
 0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.333333  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.333333  0.333333 
 0.333333  0.000000  0.000000  0.666667  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000 
 0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.333333  0.000000  0.000000  0.333333  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.333333  0.000000  0.000000 
 1.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000 
 0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  1.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000 
--------------------------------------------------------------------------------

--------------------------------------------------------------------------------
	Motif 1 regular expression
--------------------------------------------------------------------------------
[HRY][DP][GC][ARY][CV][PS][DES][VF][KR][PV][VK][DEN][ND]FD[KLW][SA][ANR][FY][AH]G[AKT]W[HWY][EA][IMV]AK
--------------------------------------------------------------------------------




Time  2.27 secs.

********************************************************************************


********************************************************************************
MOTIF  2	width =   42   sites =   2   llr = 218   E-value = 9.9e-004
********************************************************************************
--------------------------------------------------------------------------------
	Motif 2 Description
--------------------------------------------------------------------------------
Simplified        A  :::::5::::::5::::::::::::5::::5:::::::::::
pos.-specific     C  :::::::::::::::::a::::::::::::::::::::::::
probability       D  :::::::a::::::::::5:5:a:::::5:::::::::::::
matrix            E  :::::::::::::::::::::5::::::::::::::::::5:
                  F  :5:::::::::::::::::::::::::::5::::::::::::
                  G  ::::::::::::::5::::::::::5:::::::::::::::a
                  H  ::::::::::::::::::::5:::::a::5::::::::::::
                  I  ::::::::::::5a::::::::::::::5:::5:::::::::
                  K  :::::::::a::::::::5::::aa::::::::::5:a::::
                  L  ::::a::::::::::::::::::::::::::::a:::::a::
                  M  ::::::::::::::::::::::::::::::::::::::::::
                  N  ::5:::::5:a:::5:5:::::::::::::::::::::::::
                  P  :5:::::::::::::::::::5::::::::::::::::::::
                  Q  :::::::::::::::::::::::::::5::::::::::::::
                  R  :::::::::::::::::::::::::::::::::::5::::::
                  S  :::::5:::::::::::::::::::::5::::::a:a:::::
                  T  ::::::a:::::::::::::::::::::::::::::::::5:
                  V  a::a::::::::::::::::::::::::::5:5:::::a:::
                  W  ::5::::::::::::::::::::::::::::a::::::::::
                  Y  ::::::::5::a:::a5::a::::::::::::::::::::::

         bits    5.9                                *          
                 5.3                  *        *    *          
                 4.7       *      *   *        *    *          
                 4.1 * **  **  ** * * * *  *   *  * *  * * *  *
Information      3.5 ***** ****** * *** ****** ** * * ******* *
content          2.9 ******************************************
(157.6 bits)     2.3 ******************************************
                 1.8 ******************************************
                 1.2 ******************************************
                 0.6 ******************************************
                 0.0 ------------------------------------------

Multilevel           VFNVLATDNKNYAIGYNCDYDEDKKAHQDFAWILSKSKVLEG
consensus             PW  S  Y   I N Y K HP   G SIHV V  R    T 
sequence                                                       
                                                               
                                                               
--------------------------------------------------------------------------------

--------------------------------------------------------------------------------
	Motif 2 sites sorted by position p-value
--------------------------------------------------------------------------------
Sequence name             Start   P-value                               Site                 
-------------             ----- ---------            ------------------------------------------
ICYA_MANSE                  102  8.32e-48 FKFGQRVVNL VPWVLATDYKNYAINYNCDYHPDKKAHSIHAWILSKSKVLEG NTKEVVDNVL
BBP_PIEBR                    98  3.83e-47 LTYGGVTKEN VFNVLSTDNKNYIIGYYCKYDEDKKGHQDFVWVLSRSKVLTG EAKTAVENYL
--------------------------------------------------------------------------------

--------------------------------------------------------------------------------
	Motif 2 block diagrams
--------------------------------------------------------------------------------
SEQUENCE NAME            POSITION P-VALUE  MOTIF DIAGRAM
-------------            ----------------  -------------
ICYA_MANSE                        8.3e-48  101_[2]_46
BBP_PIEBR                         3.8e-47  97_[2]_34
--------------------------------------------------------------------------------

--------------------------------------------------------------------------------