meme.crp0.tcm
来自「EM算法的改进」· TCM 代码 · 共 391 行 · 第 1/2 页
TCM
391 行
13 -1081 -182 101 ---------------------------------------------------------------------------------------------------------------------------------------------------------------- Motif 1 position-specific probability matrix--------------------------------------------------------------------------------letter-probability matrix: alength= 4 w= 19 nsites= 18 E= 2.8e-004 0.000000 0.222222 0.000000 0.777778 0.000000 0.055556 0.500000 0.444444 0.000000 0.222222 0.000000 0.777778 0.222222 0.000000 0.777778 0.000000 0.666667 0.222222 0.000000 0.111111 0.222222 0.166667 0.277778 0.333333 0.166667 0.555556 0.166667 0.111111 0.222222 0.222222 0.388889 0.166667 0.666667 0.000000 0.166667 0.166667 0.055556 0.277778 0.333333 0.333333 0.222222 0.111111 0.388889 0.277778 0.000000 0.111111 0.000000 0.888889 0.055556 0.944444 0.000000 0.000000 0.833333 0.000000 0.111111 0.055556 0.000000 0.944444 0.000000 0.055556 0.888889 0.000000 0.055556 0.055556 0.333333 0.444444 0.000000 0.222222 0.333333 0.111111 0.000000 0.555556 0.333333 0.000000 0.055556 0.611111 ---------------------------------------------------------------------------------------------------------------------------------------------------------------- Motif 1 regular expression--------------------------------------------------------------------------------[TC][GT][TC][GA][AC][TGA]C[GAC]A[GTC][GTA]TCACA[CAT][TA][TA]--------------------------------------------------------------------------------Time 3.17 secs.****************************************************************************************************************************************************************MOTIF 2 width = 11 sites = 5 llr = 63 E-value = 4.9e+003********************************************************************************-------------------------------------------------------------------------------- Motif 2 Description--------------------------------------------------------------------------------Simplified A :::6aa:::::pos.-specific C 2::::::2:::probability G 8:a4::::a:2matrix T :a::::a8:a8 bits 2.4 * * 2.1 * * 1.9 * * 1.6 *** *** ** Information 1.4 *** *** ** content 1.2 *** *******(18.2 bits) 0.9 *********** 0.7 *********** 0.5 *********** 0.2 *********** 0.0 -----------Multilevel GTGAAATTGTTconsensus C G C Gsequence ---------------------------------------------------------------------------------------------------------------------------------------------------------------- Motif 2 sites sorted by position p-value--------------------------------------------------------------------------------Sequence name Strand Start P-value Site ------------- ------ ----- --------- -----------uxu1 + 10 9.04e-07 CCCATGAGA GTGAAATTGTT GTGATGTGGTlac - 95 9.04e-07 . GTGAAATTGTT ATCCGCTCAClac + 77 2.07e-06 CGTATGTTGT GTGGAATTGTG AGCGGATAACgale - 31 2.07e-06 GAATAAATTA GTGGAATCGTT TACACAAGAAtnaa - 94 2.97e-06 T CTGAAATTGTT TAAATGTGAA---------------------------------------------------------------------------------------------------------------------------------------------------------------- Motif 2 block diagrams--------------------------------------------------------------------------------SEQUENCE NAME POSITION P-VALUE MOTIF DIAGRAM------------- ---------------- -------------uxu1 9e-07 9_[+2]_85lac 2.1e-06 76_[+2]_7_[-2]gale 2.1e-06 30_[-2]_64tnaa 3e-06 93_[-2]_1---------------------------------------------------------------------------------------------------------------------------------------------------------------- Motif 2 in BLOCKS format--------------------------------------------------------------------------------BL MOTIF 2 width=11 seqs=5uxu1 ( 10) GTGAAATTGTT 1 lac ( 95) GTGAAATTGTT 1 lac ( 77) GTGGAATTGTG 1 gale ( 31) GTGGAATCGTT 1 tnaa ( 94) CTGAAATTGTT 1 //---------------------------------------------------------------------------------------------------------------------------------------------------------------- Motif 2 position-specific scoring matrix--------------------------------------------------------------------------------log-odds matrix: alength= 4 w= 11 n= 1710 bayes= 8.66584 E= 4.9e+003 -897 3 203 -897 -897 -897 -897 172 -897 -897 235 -897 98 -897 103 -897 172 -897 -897 -897 172 -897 -897 -897 -897 -897 -897 172 -897 3 -897 139 -897 -897 235 -897 -897 -897 -897 172 -897 -897 3 139 ---------------------------------------------------------------------------------------------------------------------------------------------------------------- Motif 2 position-specific probability matrix--------------------------------------------------------------------------------letter-probability matrix: alength= 4 w= 11 nsites= 5 E= 4.9e+003 0.000000 0.200000 0.800000 0.000000 0.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.600000 0.000000 0.400000 0.000000 1.000000 0.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 1.000000 0.000000 0.200000 0.000000 0.800000 0.000000 0.000000 1.000000 0.000000 0.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.200000 0.800000 ---------------------------------------------------------------------------------------------------------------------------------------------------------------- Motif 2 regular expression--------------------------------------------------------------------------------[GC]TG[AG]AAT[TC]GT[TG]--------------------------------------------------------------------------------Time 5.27 secs.****************************************************************************************************************************************************************SUMMARY OF MOTIFS********************************************************************************-------------------------------------------------------------------------------- Combined block diagrams: non-overlapping sites with p-value < 0.0001--------------------------------------------------------------------------------SEQUENCE NAME COMBINED P-VALUE MOTIF DIAGRAM------------- ---------------- -------------ce1cg 1.17e-04 16_[-1(5.52e-06)]_28_[+1(5.44e-07)]_23ara 1.94e-05 57_[+1(2.02e-08)]_29bglr1 1.81e-03 78_[+1(3.88e-06)]_8crp 2.53e-03 65_[+1(2.43e-06)]_21cya 1.15e-03 49_[-1(9.07e-06)]_37deop2 8.27e-04 9_[+1(2.00e-06)]_31_[-1(1.19e-06)]_27gale 8.39e-07 30_[-2(2.07e-06)]_3_[+1(6.87e-07)]_42ilv 1.31e-02 41_[+1(1.92e-05)]_45lac 4.37e-06 11_[+1(9.07e-06)]_46_[+2(2.07e-06)]_7_[-2(9.04e-07)]male 1.62e-03 16_[+1(1.47e-06)]_70malk 6.73e-03 60_[-1(2.06e-05)]_26malt 1.76e-05 40_[-1(6.87e-07)]_46ompa 9.97e-04 50_[+1(8.61e-07)]_36tnaa 1.54e-05 73_[+1(1.06e-05)]_1_[-2(2.97e-06)]_1uxu1 3.88e-04 9_[+2(9.04e-07)]_85pbr322 5.63e-03 52_[-1(6.54e-06)]_34trn9cat 3.22e-01 105tdc 1.58e-03 78_[+1(6.54e-06)]_8--------------------------------------------------------------------------------****************************************************************************************************************************************************************Stopped because nmotifs = 2 reached.********************************************************************************CPU: chromo.imb.uq.edu.au********************************************************************************
⌨️ 快捷键说明
复制代码Ctrl + C
搜索代码Ctrl + F
全屏模式F11
增大字号Ctrl + =
减小字号Ctrl + -
显示快捷键?