******************************************************************************** MEME - Motif discovery tool ******************************************************************************** MEME version 4.10.0 (Release date: Wed May 21 10:35:36 2014 +1000) For further information on how to interpret these results or to get a copy of the MEME software please access http://meme.nbcr.net. This file may be used as input to the MAST algorithm for searching sequence databases for matches to groups of motifs. MAST is available for interactive use and downloading at http://meme.nbcr.net. ******************************************************************************** ******************************************************************************** REFERENCE ******************************************************************************** If you use this program in your research, please cite: Timothy L. Bailey and Charles Elkan, "Fitting a mixture model by expectation maximization to discover motifs in biopolymers", Proceedings of the Second International Conference on Intelligent Systems for Molecular Biology, pp. 28-36, AAAI Press, Menlo Park, California, 1994. ******************************************************************************** ******************************************************************************** TRAINING SET ******************************************************************************** DATAFILE= motifs/233/233.seqs.fa ALPHABET= ACGT Sequence name Weight Length Sequence name Weight Length ------------- ------ ------ ------------- ------ ------ 10835 1.0000 500 10858 1.0000 500 11217 1.0000 500 11562 1.0000 500 16584 1.0000 500 22811 1.0000 500 25110 1.0000 500 261038 1.0000 500 261093 1.0000 500 263053 1.0000 500 263118 1.0000 500 2886 1.0000 500 2992 1.0000 500 5250 1.0000 500 6170 1.0000 500 6645 1.0000 500 669 1.0000 500 6779 1.0000 500 7554 1.0000 500 7910 1.0000 500 8240 1.0000 500 ******************************************************************************** ******************************************************************************** COMMAND LINE SUMMARY ******************************************************************************** This information can also be useful in the event you wish to report a problem with the MEME software. command: meme motifs/233/233.seqs.fa -oc motifs/233 -dna -minw 12 -maxw 21 -nmotifs 3 -maxsize 500000 model: mod= zoops nmotifs= 3 evt= inf object function= E-value of product of p-values width: minw= 12 maxw= 21 minic= 0.00 width: wg= 11 ws= 1 endgaps= yes nsites: minsites= 2 maxsites= 21 wnsites= 0.8 theta: prob= 1 spmap= uni spfuzz= 0.5 global: substring= yes branching= no wbranch= no em: prior= dirichlet b= 0.01 maxiter= 50 distance= 1e-05 data: n= 10500 N= 21 strands: + sample: seed= 0 seqfrac= 1 Letter frequencies in dataset: A 0.266 C 0.229 G 0.232 T 0.273 Background letter frequencies (from dataset with add-one prior applied): A 0.266 C 0.229 G 0.232 T 0.273 ******************************************************************************** ******************************************************************************** MOTIF 1 MEME width = 12 sites = 16 llr = 164 E-value = 1.0e-002 ******************************************************************************** -------------------------------------------------------------------------------- Motif 1 Description -------------------------------------------------------------------------------- Simplified A :392341a8:38 pos.-specific C a4:77:9:1772 probability G :31::3::13:: matrix T :::1:3:::::: bits 2.1 * 1.9 * * 1.7 * ** 1.5 * * ** Relative 1.3 * * ** * * Entropy 1.1 * * * ****** (14.8 bits) 0.9 * *** ****** 0.6 * *** ****** 0.4 ************ 0.2 ************ 0.0 ------------ Multilevel CCACCACAACCA consensus A AG GA sequence G T -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 sites sorted by position p-value -------------------------------------------------------------------------------- Sequence name Start P-value Site ------------- ----- --------- ------------ 25110 480 2.11e-07 AGAGACGCAC CGACCGCAACCA CTAGACAAC 261093 428 8.10e-07 TCACATCATC CCACCACAACAA CACACCAAAT 669 260 4.50e-06 CAGATGACGA CGACAACAACAA TGACGAAGAA 7554 465 4.79e-06 ATTCCTGTCG CAGCCACAACCA ACCTACCACA 16584 227 5.53e-06 TCATATTCGA CAACCTCAGCCA ATGACTATTA 11562 458 5.53e-06 AGTTCCAGTG CAACAGCAACAA CAGTTTCCCT 2886 476 6.53e-06 TACACAAAAA CGAACACAAGCA CAAAGTATTG 10835 107 6.53e-06 GGTGCTCTCT CCAACACAACAA CCAATCTCAG 6779 479 7.73e-06 AAACAAACAA CAACCTCAAGCC AACTTCGTCA 6170 355 7.73e-06 TGGATCATCT CCATCACAAGCA AACCGATAAT 2992 479 7.73e-06 AAACAAATAA CAACCTCAAGCC AACTTTGTCA 263118 487 9.63e-06 CTGCTCCAAG CCACCGAAACCA CG 263053 80 1.36e-05 TCATAATGTG CGGCATCAACCA CCTCTGCGGT 8240 167 2.27e-05 CGCCGTCCTT CCAACGCAACAC ACCTGGGAGC 11217 441 3.32e-05 GCCACCAAGG CCACATCACGCA AATCCCATTA 6645 479 4.41e-05 CCAATCCTCC CGATAGCAGCCA TCCCGGCACA -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 block diagrams -------------------------------------------------------------------------------- SEQUENCE NAME POSITION P-VALUE MOTIF DIAGRAM ------------- ---------------- ------------- 25110 2.1e-07 479_[+1]_9 261093 8.1e-07 427_[+1]_61 669 4.5e-06 259_[+1]_229 7554 4.8e-06 464_[+1]_24 16584 5.5e-06 226_[+1]_262 11562 5.5e-06 457_[+1]_31 2886 6.5e-06 475_[+1]_13 10835 6.5e-06 106_[+1]_382 6779 7.7e-06 478_[+1]_10 6170 7.7e-06 354_[+1]_134 2992 7.7e-06 478_[+1]_10 263118 9.6e-06 486_[+1]_2 263053 1.4e-05 79_[+1]_409 8240 2.3e-05 166_[+1]_322 11217 3.3e-05 440_[+1]_48 6645 4.4e-05 478_[+1]_10 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 in BLOCKS format -------------------------------------------------------------------------------- BL MOTIF 1 width=12 seqs=16 25110 ( 480) CGACCGCAACCA 1 261093 ( 428) CCACCACAACAA 1 669 ( 260) CGACAACAACAA 1 7554 ( 465) CAGCCACAACCA 1 16584 ( 227) CAACCTCAGCCA 1 11562 ( 458) CAACAGCAACAA 1 2886 ( 476) CGAACACAAGCA 1 10835 ( 107) CCAACACAACAA 1 6779 ( 479) CAACCTCAAGCC 1 6170 ( 355) CCATCACAAGCA 1 2992 ( 479) CAACCTCAAGCC 1 263118 ( 487) CCACCGAAACCA 1 263053 ( 80) CGGCATCAACCA 1 8240 ( 167) CCAACGCAACAC 1 11217 ( 441) CCACATCACGCA 1 6645 ( 479) CGATAGCAGCCA 1 // -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 position-specific scoring matrix -------------------------------------------------------------------------------- log-odds matrix: alength= 4 w= 12 n= 10269 bayes= 11.0623 E= 1.0e-002 -1064 212 -1064 -1064 23 71 43 -1064 172 -1064 -89 -1064 -50 158 -1064 -112 23 158 -1064 -1064 50 -1064 43 20 -208 203 -1064 -1064 191 -1064 -1064 -1064 161 -187 -89 -1064 -1064 158 43 -1064 23 158 -1064 -1064 161 -29 -1064 -1064 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 position-specific probability matrix -------------------------------------------------------------------------------- letter-probability matrix: alength= 4 w= 12 nsites= 16 E= 1.0e-002 0.000000 1.000000 0.000000 0.000000 0.312500 0.375000 0.312500 0.000000 0.875000 0.000000 0.125000 0.000000 0.187500 0.687500 0.000000 0.125000 0.312500 0.687500 0.000000 0.000000 0.375000 0.000000 0.312500 0.312500 0.062500 0.937500 0.000000 0.000000 1.000000 0.000000 0.000000 0.000000 0.812500 0.062500 0.125000 0.000000 0.000000 0.687500 0.312500 0.000000 0.312500 0.687500 0.000000 0.000000 0.812500 0.187500 0.000000 0.000000 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 regular expression -------------------------------------------------------------------------------- C[CAG]AC[CA][AGT]CAA[CG][CA]A -------------------------------------------------------------------------------- Time 4.04 secs. ******************************************************************************** ******************************************************************************** MOTIF 2 MEME width = 14 sites = 19 llr = 185 E-value = 2.6e-002 ******************************************************************************** -------------------------------------------------------------------------------- Motif 2 Description -------------------------------------------------------------------------------- Simplified A 5::21:4:2:2:5: pos.-specific C 3:a1:9::232411 probability G 11:3:16:643643 matrix T 19:59::a:33:17 bits 2.1 * 1.9 * * 1.7 * * * 1.5 ** ** * Relative 1.3 ** ** * Entropy 1.1 ** **** * (14.1 bits) 0.9 ** ***** * * 0.6 ** ***** * * 0.4 *** ****** *** 0.2 ********** *** 0.0 -------------- Multilevel ATCTTCGTGGTGAT consensus C G A ATGCGG sequence CA C -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 2 sites sorted by position p-value -------------------------------------------------------------------------------- Sequence name Start P-value Site ------------- ----- --------- -------------- 7554 298 8.23e-09 GAGCGAGGAT ATCTTCGTGGGGAT TAGTTGCGTT 7910 193 1.48e-06 CTGGTGGAGC CTCGTCATGGCGAT GGTGGGGTGG 263118 288 1.48e-06 GGTACATCAA CTCTTCGTGTACGT CTTTGTGCTG 263053 425 3.35e-06 AGTGTATCAT ATCATCATGCTGAT ATACGTGTGC 10835 386 3.79e-06 TCATTGAACC ATCGTCGTCGGCGT ACCTCTAAAT 22811 9 4.79e-06 TGCCCTAA ATCTTCGTATTGAG GTGATCGTGG 2992 132 1.01e-05 ACACGCCTAC ATCATCGTATGCAT CCCCACCAGA 25110 90 1.01e-05 TCTATCCAAG GTCTTCATGCCCAT CACTAAATTC 11562 385 1.01e-05 TCTCGCTGTG TTCTTCGTGGTGAG CACTGCCTTT 6645 251 1.10e-05 TCTCAGCTTA ATCCTCGTATTGGT TGTTCATAGG 16584 70 1.45e-05 TCAATACTCC ATCGTCGTACTCGG CTCCCTCGTC 8240 280 1.71e-05 AAGGTGGTGA ATCATCATCGGCGT TGATTTCTTG 6779 79 1.87e-05 GAGGCTGGAA ATCGACGTGTCCGT TCACTGTAAC 2886 253 1.87e-05 GGAATGCTAA CGCTTCGTGCCGAT GCAACCGTTA 6170 162 2.78e-05 ATGAAGGGCT CTCTTCGTCGTCCT TGATACCACC 261093 355 2.78e-05 CCACCCTCTT CTCCACGTGGAGAT GACGTCCCGT 5250 62 3.22e-05 TTGGCTGAGA CTCTTGATGGAGAG TTTGATCCTG 261038 169 9.72e-05 TGGTGATTGC GTCGTCATGCAGGC TGTTCTCACG 669 214 1.24e-04 GAAGCGGGAG ATCTTGATGTGGTG GACGATGATA -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 2 block diagrams -------------------------------------------------------------------------------- SEQUENCE NAME POSITION P-VALUE MOTIF DIAGRAM ------------- ---------------- ------------- 7554 8.2e-09 297_[+2]_189 7910 1.5e-06 192_[+2]_294 263118 1.5e-06 287_[+2]_199 263053 3.3e-06 424_[+2]_62 10835 3.8e-06 385_[+2]_101 22811 4.8e-06 8_[+2]_478 2992 1e-05 131_[+2]_355 25110 1e-05 89_[+2]_397 11562 1e-05 384_[+2]_102 6645 1.1e-05 250_[+2]_236 16584 1.4e-05 69_[+2]_417 8240 1.7e-05 279_[+2]_207 6779 1.9e-05 78_[+2]_408 2886 1.9e-05 252_[+2]_234 6170 2.8e-05 161_[+2]_325 261093 2.8e-05 354_[+2]_132 5250 3.2e-05 61_[+2]_425 261038 9.7e-05 168_[+2]_318 669 0.00012 213_[+2]_273 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 2 in BLOCKS format -------------------------------------------------------------------------------- BL MOTIF 2 width=14 seqs=19 7554 ( 298) ATCTTCGTGGGGAT 1 7910 ( 193) CTCGTCATGGCGAT 1 263118 ( 288) CTCTTCGTGTACGT 1 263053 ( 425) ATCATCATGCTGAT 1 10835 ( 386) ATCGTCGTCGGCGT 1 22811 ( 9) ATCTTCGTATTGAG 1 2992 ( 132) ATCATCGTATGCAT 1 25110 ( 90) GTCTTCATGCCCAT 1 11562 ( 385) TTCTTCGTGGTGAG 1 6645 ( 251) ATCCTCGTATTGGT 1 16584 ( 70) ATCGTCGTACTCGG 1 8240 ( 280) ATCATCATCGGCGT 1 6779 ( 79) ATCGACGTGTCCGT 1 2886 ( 253) CGCTTCGTGCCGAT 1 6170 ( 162) CTCTTCGTCGTCCT 1 261093 ( 355) CTCCACGTGGAGAT 1 5250 ( 62) CTCTTGATGGAGAG 1 261038 ( 169) GTCGTCATGCAGGC 1 669 ( 214) ATCTTGATGTGGTG 1 // -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 2 position-specific scoring matrix -------------------------------------------------------------------------------- log-odds matrix: alength= 4 w= 14 n= 10227 bayes= 9.26464 E= 2.6e-002 99 46 -114 -237 -1089 -1089 -214 180 -1089 212 -1089 -1089 -75 -112 18 80 -133 -1089 -1089 171 -1089 196 -114 -1089 47 -1089 144 -1089 -1089 -1089 -1089 187 -34 -54 144 -1089 -1089 20 86 21 -34 -12 18 21 -1089 88 132 -1089 99 -212 67 -237 -1089 -212 18 133 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 2 position-specific probability matrix -------------------------------------------------------------------------------- letter-probability matrix: alength= 4 w= 14 nsites= 19 E= 2.6e-002 0.526316 0.315789 0.105263 0.052632 0.000000 0.000000 0.052632 0.947368 0.000000 1.000000 0.000000 0.000000 0.157895 0.105263 0.263158 0.473684 0.105263 0.000000 0.000000 0.894737 0.000000 0.894737 0.105263 0.000000 0.368421 0.000000 0.631579 0.000000 0.000000 0.000000 0.000000 1.000000 0.210526 0.157895 0.631579 0.000000 0.000000 0.263158 0.421053 0.315789 0.210526 0.210526 0.263158 0.315789 0.000000 0.421053 0.578947 0.000000 0.526316 0.052632 0.368421 0.052632 0.000000 0.052632 0.263158 0.684211 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 2 regular expression -------------------------------------------------------------------------------- [AC]TC[TG]TC[GA]T[GA][GTC][TGAC][GC][AG][TG] -------------------------------------------------------------------------------- Time 8.07 secs. ******************************************************************************** ******************************************************************************** MOTIF 3 MEME width = 20 sites = 9 llr = 137 E-value = 2.5e+001 ******************************************************************************** -------------------------------------------------------------------------------- Motif 3 Description -------------------------------------------------------------------------------- Simplified A ::2:114a131862:9a1:9 pos.-specific C :4:::81::13::861:121 probability G 71:89:3:86124:4::28: matrix T 3482:11:1:4::::::6:: bits 2.1 1.9 * * 1.7 * * 1.5 * * ** * Relative 1.3 ** * * * ** ** Entropy 1.1 * **** ** ****** ** (22.0 bits) 0.9 * **** ** ****** ** 0.6 ****** *** ****** ** 0.4 ****** *** ****** ** 0.2 ******************** 0.0 -------------------- Multilevel GCTGGCAAGGTAACCAATGA consensus TTAT G ACGGAG GC sequence -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 3 sites sorted by position p-value -------------------------------------------------------------------------------- Sequence name Start P-value Site ------------- ----- --------- -------------------- 263118 165 8.39e-11 CTTTTGAGTC TCTGGCAAGGCAACGAATGA GGTGGTGGCG 11217 165 6.97e-10 TGGTATTAGT GTTGGCCAGACAGCCAATGA TGCTCTTGAT 10858 465 1.80e-08 GCTGAACTTG TTTGGCAAGCTAGCCCATGA ACTCGATTGT 22811 399 5.67e-08 CGGGATAACG GTAGACAATGTAACCAATGA CATCACTTTT 263053 219 7.84e-08 CAATAATGTT TCTGGCGAGGAGAAGAAGGA AGTGTCAGTT 7910 230 9.10e-08 TGTTGGATTA GTTGGAGAAGTAGCGAACGA AAGATTTACA 669 330 9.10e-08 TGGCAGAGGC GGTGGCAAGACAGCCAAGCC AAAGGCAGCC 6170 232 2.61e-07 TATCCAAGAT GCATGCTAGATGACCAAAGA GGATGCATAT 16584 43 3.56e-07 TAAGAGAGCA GCTTGTGAGGGAAAGAATCA ATACTCCATC -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 3 block diagrams -------------------------------------------------------------------------------- SEQUENCE NAME POSITION P-VALUE MOTIF DIAGRAM ------------- ---------------- ------------- 263118 8.4e-11 164_[+3]_316 11217 7e-10 164_[+3]_316 10858 1.8e-08 464_[+3]_16 22811 5.7e-08 398_[+3]_82 263053 7.8e-08 218_[+3]_262 7910 9.1e-08 229_[+3]_251 669 9.1e-08 329_[+3]_151 6170 2.6e-07 231_[+3]_249 16584 3.6e-07 42_[+3]_438 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 3 in BLOCKS format -------------------------------------------------------------------------------- BL MOTIF 3 width=20 seqs=9 263118 ( 165) TCTGGCAAGGCAACGAATGA 1 11217 ( 165) GTTGGCCAGACAGCCAATGA 1 10858 ( 465) TTTGGCAAGCTAGCCCATGA 1 22811 ( 399) GTAGACAATGTAACCAATGA 1 263053 ( 219) TCTGGCGAGGAGAAGAAGGA 1 7910 ( 230) GTTGGAGAAGTAGCGAACGA 1 669 ( 330) GGTGGCAAGACAGCCAAGCC 1 6170 ( 232) GCATGCTAGATGACCAAAGA 1 16584 ( 43) GCTTGTGAGGGAAAGAATCA 1 // -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 3 position-specific scoring matrix -------------------------------------------------------------------------------- log-odds matrix: alength= 4 w= 20 n= 10101 bayes= 10.2654 E= 2.5e+001 -982 -982 152 29 -982 95 -106 70 -26 -982 -982 151 -982 -982 174 -30 -126 -982 193 -982 -126 176 -982 -129 74 -104 52 -129 191 -982 -982 -982 -126 -982 174 -129 33 -104 126 -982 -126 54 -106 70 155 -982 -6 -982 106 -982 94 -982 -26 176 -982 -982 -982 128 94 -982 174 -104 -982 -982 191 -982 -982 -982 -126 -104 -6 103 -982 -4 174 -982 174 -104 -982 -982 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 3 position-specific probability matrix -------------------------------------------------------------------------------- letter-probability matrix: alength= 4 w= 20 nsites= 9 E= 2.5e+001 0.000000 0.000000 0.666667 0.333333 0.000000 0.444444 0.111111 0.444444 0.222222 0.000000 0.000000 0.777778 0.000000 0.000000 0.777778 0.222222 0.111111 0.000000 0.888889 0.000000 0.111111 0.777778 0.000000 0.111111 0.444444 0.111111 0.333333 0.111111 1.000000 0.000000 0.000000 0.000000 0.111111 0.000000 0.777778 0.111111 0.333333 0.111111 0.555556 0.000000 0.111111 0.333333 0.111111 0.444444 0.777778 0.000000 0.222222 0.000000 0.555556 0.000000 0.444444 0.000000 0.222222 0.777778 0.000000 0.000000 0.000000 0.555556 0.444444 0.000000 0.888889 0.111111 0.000000 0.000000 1.000000 0.000000 0.000000 0.000000 0.111111 0.111111 0.222222 0.555556 0.000000 0.222222 0.777778 0.000000 0.888889 0.111111 0.000000 0.000000 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 3 regular expression -------------------------------------------------------------------------------- [GT][CT][TA][GT]GC[AG]AG[GA][TC][AG][AG][CA][CG]AA[TG][GC]A -------------------------------------------------------------------------------- Time 11.74 secs. ******************************************************************************** ******************************************************************************** SUMMARY OF MOTIFS ******************************************************************************** -------------------------------------------------------------------------------- Combined block diagrams: non-overlapping sites with p-value < 0.0001 -------------------------------------------------------------------------------- SEQUENCE NAME COMBINED P-VALUE MOTIF DIAGRAM ------------- ---------------- ------------- 10835 2.38e-04 106_[+1(6.53e-06)]_267_\ [+2(3.79e-06)]_101 10858 4.28e-04 464_[+3(1.80e-08)]_16 11217 7.44e-07 164_[+3(6.97e-10)]_256_\ [+1(3.32e-05)]_48 11562 6.37e-04 384_[+2(1.01e-05)]_59_\ [+1(5.53e-06)]_31 16584 6.86e-07 42_[+3(3.56e-07)]_7_[+2(1.45e-05)]_\ 143_[+1(5.53e-06)]_262 22811 5.48e-06 8_[+2(4.79e-06)]_322_[+2(2.78e-05)]_\ 40_[+3(5.67e-08)]_82 25110 5.05e-06 89_[+2(1.01e-05)]_376_\ [+1(2.11e-07)]_9 261038 2.10e-01 168_[+2(9.72e-05)]_318 261093 3.37e-04 354_[+2(2.78e-05)]_59_\ [+1(8.10e-07)]_61 263053 1.05e-07 79_[+1(1.36e-05)]_127_\ [+3(7.84e-08)]_186_[+2(3.35e-06)]_62 263118 6.40e-11 45_[+3(9.55e-05)]_99_[+3(8.39e-11)]_\ 103_[+2(1.48e-06)]_185_[+1(9.63e-06)]_2 2886 1.50e-03 252_[+2(1.87e-05)]_209_\ [+1(6.53e-06)]_13 2992 9.30e-04 131_[+2(1.01e-05)]_333_\ [+1(7.73e-06)]_10 5250 2.17e-02 61_[+2(3.22e-05)]_178_\ [+3(9.08e-05)]_227 6170 1.26e-06 161_[+2(2.78e-05)]_56_\ [+3(2.61e-07)]_103_[+1(7.73e-06)]_134 6645 4.99e-03 250_[+2(1.10e-05)]_214_\ [+1(4.41e-05)]_10 669 1.13e-06 259_[+1(4.50e-06)]_58_\ [+3(9.10e-08)]_151 6779 5.91e-04 78_[+2(1.87e-05)]_386_\ [+1(7.73e-06)]_10 7554 8.15e-08 222_[+3(7.00e-05)]_55_\ [+2(8.23e-09)]_153_[+1(4.79e-06)]_24 7910 2.16e-06 192_[+2(1.48e-06)]_23_\ [+3(9.10e-08)]_251 8240 3.41e-03 166_[+1(2.27e-05)]_101_\ [+2(1.71e-05)]_207 -------------------------------------------------------------------------------- ******************************************************************************** ******************************************************************************** Stopped because nmotifs = 3 reached. ******************************************************************************** CPU: seaotter.hsd1.wa.comcast.net ********************************************************************************