meme.ino_up800.oops
来自「EM算法的改进」· OOPS 代码 · 共 335 行 · 第 1/2 页
OOPS
335 行
0.000000 0.571429 0.428571 0.000000 0.000000 0.857143 0.000000 0.142857 0.285714 0.142857 0.428571 0.142857 0.000000 1.000000 0.000000 0.000000 ---------------------------------------------------------------------------------------------------------------------------------------------------------------- Motif 1 regular expression--------------------------------------------------------------------------------TTCACATG[CG]C[GA]C--------------------------------------------------------------------------------Time 4.83 secs.****************************************************************************************************************************************************************MOTIF 2 width = 10 sites = 7 llr = 81 E-value = 1.1e+002********************************************************************************-------------------------------------------------------------------------------- Motif 2 Description--------------------------------------------------------------------------------Simplified A ::1:::9:6:pos.-specific C :a:::a:911probability G 3:1aa:1:19matrix T 7:7::::11: bits 2.5 * *** 2.3 * *** 2.0 * *** * 1.8 * *** * *Information 1.5 * *** * *content 1.3 * ***** *(16.7 bits) 1.0 ** ***** * 0.8 ** ***** * 0.5 ******** * 0.3 ********** 0.0 ----------Multilevel TCTGGCACAGconsensus G sequence ---------------------------------------------------------------------------------------------------------------------------------------------------------------- Motif 2 sites sorted by position p-value--------------------------------------------------------------------------------Sequence name Strand Start P-value Site ------------- ------ ----- --------- ----------OPI3 - 186 3.24e-07 GAAAACCAGA TCTGGCACAG ACCGTTGTCAACC1 + 232 3.24e-07 CCAGTCGTAT TCTGGCACAG TATAGCCTAGCHO1 - 559 3.24e-07 ATATTCAGTG TCTGGCACAG AAGTCTGCACINO1 - 283 5.29e-06 ACGGTCTACG GCGGGCGCAG TCGCATGTCTFAS1 + 44 6.25e-06 TACACGAGGT GCAGGCACGG TTCACTACTCFAS2 - 185 8.48e-06 TTCTTGCTTT TCTGGCACTC TTGACGGCTTCHO2 - 413 8.48e-06 TTTTGCCGTT TCTGGCATCG CCGTTCATTT---------------------------------------------------------------------------------------------------------------------------------------------------------------- Motif 2 block diagrams--------------------------------------------------------------------------------SEQUENCE NAME POSITION P-VALUE MOTIF DIAGRAM------------- ---------------- -------------OPI3 3.2e-07 185_[-2]_605ACC1 3.2e-07 231_[+2]_559CHO1 3.2e-07 558_[-2]_232INO1 5.3e-06 282_[-2]_508FAS1 6.3e-06 43_[+2]_747FAS2 8.5e-06 184_[-2]_606CHO2 8.5e-06 412_[-2]_378---------------------------------------------------------------------------------------------------------------------------------------------------------------- Motif 2 in BLOCKS format--------------------------------------------------------------------------------BL MOTIF 2 width=10 seqs=7OPI3 ( 186) TCTGGCACAG 1 ACC1 ( 232) TCTGGCACAG 1 CHO1 ( 559) TCTGGCACAG 1 INO1 ( 283) GCGGGCGCAG 1 FAS1 ( 44) GCAGGCACGG 1 FAS2 ( 185) TCTGGCACTC 1 CHO2 ( 413) TCTGGCATCG 1 //---------------------------------------------------------------------------------------------------------------------------------------------------------------- Motif 2 position-specific scoring matrix--------------------------------------------------------------------------------log-odds matrix: alength= 4 w= 10 n= 5537 bayes= 9.62571 E= 1.1e+002 -945 -945 70 114 -945 251 -945 -945 -118 -945 -30 114 -945 -945 251 -945 -945 -945 251 -945 -945 251 -945 -945 140 -945 -30 -945 -945 229 -945 -118 82 -30 -30 -118 -945 -30 229 -945 ---------------------------------------------------------------------------------------------------------------------------------------------------------------- Motif 2 position-specific probability matrix--------------------------------------------------------------------------------letter-probability matrix: alength= 4 w= 10 nsites= 7 E= 1.1e+002 0.000000 0.000000 0.285714 0.714286 0.000000 1.000000 0.000000 0.000000 0.142857 0.000000 0.142857 0.714286 0.000000 0.000000 1.000000 0.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.857143 0.000000 0.142857 0.000000 0.000000 0.857143 0.000000 0.142857 0.571429 0.142857 0.142857 0.142857 0.000000 0.142857 0.857143 0.000000 ---------------------------------------------------------------------------------------------------------------------------------------------------------------- Motif 2 regular expression--------------------------------------------------------------------------------[TG]CTGGCACAG--------------------------------------------------------------------------------Time 9.64 secs.****************************************************************************************************************************************************************SUMMARY OF MOTIFS********************************************************************************-------------------------------------------------------------------------------- Combined block diagrams: non-overlapping sites with p-value < 0.0001--------------------------------------------------------------------------------SEQUENCE NAME COMBINED P-VALUE MOTIF DIAGRAM------------- ---------------- -------------CHO1 5.44e-06 152_[+2(1.10e-05)]_396_[-2(3.24e-07)]_42_[+1(4.23e-07)]_17_[+1(1.23e-05)]_149CHO2 6.96e-05 353_[+1(2.52e-07)]_47_[-2(8.48e-06)]_378FAS1 4.61e-06 43_[+2(6.25e-06)]_41_[+1(1.85e-08)]_694FAS2 2.34e-04 184_[-2(8.48e-06)]_372_[+1(9.43e-07)]_222ACC1 2.09e-06 82_[+1(1.52e-07)]_137_[+2(3.24e-07)]_559INO1 3.95e-06 282_[-2(5.29e-06)]_327_[-1(1.85e-08)]_55_[+1(7.55e-06)]_102OPI3 3.70e-05 185_[-2(3.24e-07)]_144_[+1(3.32e-06)]_449--------------------------------------------------------------------------------****************************************************************************************************************************************************************Stopped because nmotifs = 2 reached.********************************************************************************CPU: chromo.imb.uq.edu.au********************************************************************************
⌨️ 快捷键说明
复制代码Ctrl + C
搜索代码Ctrl + F
全屏模式F11
增大字号Ctrl + =
减小字号Ctrl + -
显示快捷键?