******************************************************************************** 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/403/403.seqs.fa ALPHABET= ACGT Sequence name Weight Length Sequence name Weight Length ------------- ------ ------ ------------- ------ ------ 1946 1.0000 500 20999 1.0000 500 2126 1.0000 500 22692 1.0000 500 2278 1.0000 500 23546 1.0000 500 25825 1.0000 500 261172 1.0000 500 261705 1.0000 500 3006 1.0000 500 3157 1.0000 500 31979 1.0000 500 32971 1.0000 500 3987 1.0000 500 4115 1.0000 500 6681 1.0000 500 6848 1.0000 500 7921 1.0000 500 9465 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/403/403.seqs.fa -oc motifs/403 -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= 19 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= 9500 N= 19 strands: + sample: seed= 0 seqfrac= 1 Letter frequencies in dataset: A 0.269 C 0.230 G 0.238 T 0.262 Background letter frequencies (from dataset with add-one prior applied): A 0.269 C 0.230 G 0.238 T 0.262 ******************************************************************************** ******************************************************************************** MOTIF 1 MEME width = 12 sites = 16 llr = 158 E-value = 8.3e-001 ******************************************************************************** -------------------------------------------------------------------------------- Motif 1 Description -------------------------------------------------------------------------------- Simplified A ::84a:4:9216 pos.-specific C 11:::1:1:::: probability G 8826:9391693 matrix T 11::::3::2:2 bits 2.1 1.9 * 1.7 * * 1.5 ** ** * Relative 1.3 * * ** ** * Entropy 1.1 ****** ** * (14.2 bits) 0.8 ****** ** * 0.6 ****** **** 0.4 ************ 0.2 ************ 0.0 ------------ Multilevel GGAGAGAGAGGA consensus A T G sequence G -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 sites sorted by position p-value -------------------------------------------------------------------------------- Sequence name Start P-value Site ------------- ----- --------- ------------ 31979 465 6.09e-08 TCAGAGCAGA GGAGAGAGAGGA CTCATTGCAC 32971 333 5.98e-07 GGGAATGCGA GGAGAGTGAGGG GTGTATTGGA 2126 290 8.35e-07 GTCGGCCGTC GGAGAGAGAAGA GTTTTTCTTA 3987 125 4.19e-06 TCACGTAGCA GGAGACGGAGGA GGCAGCAGCT 2278 60 7.16e-06 TCCGTGTGTG GTAGAGAGAGGT TAGAGGGTGT 1946 120 7.16e-06 TTCCGAGTAG GGGGAGTGATGA TTGGAGAGTT 7921 458 8.05e-06 GTCTGGAGGT GGAAAGAGGGGA GGAAGACCGA 3157 75 9.64e-06 GGGAGAATAC GGAGAGGCAGGA ACTTTGTAAG 6681 239 1.27e-05 ACAACTTTTA GCAAAGAGAAGA AGGTTGCACG 22692 341 1.37e-05 GGTGGAGTGT TGAGAGTGATGA GTGGTTGATG 261705 42 1.66e-05 TGAGATGGGT GCAGAGAGATGG AGAGGAGCGT 25825 40 1.79e-05 CACTACTGAT GGAAACGGAGGT GTCAAAAACC 23546 198 3.81e-05 ATCTCCCCGC GTGAAGGGAGGG GTGGGGCCGA 261172 264 4.07e-05 AACCCAGAAG GGGAAGTGAGAG GCTGGAAATG 3006 402 4.62e-05 ACTTCGCCAA CGAAAGAGAGAA AACACAGAAG 4115 19 4.96e-05 CACTTAAAGC TGAAAGTGAAGT GTACGAGAAG -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 block diagrams -------------------------------------------------------------------------------- SEQUENCE NAME POSITION P-VALUE MOTIF DIAGRAM ------------- ---------------- ------------- 31979 6.1e-08 464_[+1]_24 32971 6e-07 332_[+1]_156 2126 8.4e-07 289_[+1]_199 3987 4.2e-06 124_[+1]_364 2278 7.2e-06 59_[+1]_429 1946 7.2e-06 119_[+1]_369 7921 8.1e-06 457_[+1]_31 3157 9.6e-06 74_[+1]_414 6681 1.3e-05 238_[+1]_250 22692 1.4e-05 340_[+1]_148 261705 1.7e-05 41_[+1]_447 25825 1.8e-05 39_[+1]_449 23546 3.8e-05 197_[+1]_291 261172 4.1e-05 263_[+1]_225 3006 4.6e-05 401_[+1]_87 4115 5e-05 18_[+1]_470 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 in BLOCKS format -------------------------------------------------------------------------------- BL MOTIF 1 width=12 seqs=16 31979 ( 465) GGAGAGAGAGGA 1 32971 ( 333) GGAGAGTGAGGG 1 2126 ( 290) GGAGAGAGAAGA 1 3987 ( 125) GGAGACGGAGGA 1 2278 ( 60) GTAGAGAGAGGT 1 1946 ( 120) GGGGAGTGATGA 1 7921 ( 458) GGAAAGAGGGGA 1 3157 ( 75) GGAGAGGCAGGA 1 6681 ( 239) GCAAAGAGAAGA 1 22692 ( 341) TGAGAGTGATGA 1 261705 ( 42) GCAGAGAGATGG 1 25825 ( 40) GGAAACGGAGGT 1 23546 ( 198) GTGAAGGGAGGG 1 261172 ( 264) GGGAAGTGAGAG 1 3006 ( 402) CGAAAGAGAGAA 1 4115 ( 19) TGAAAGTGAAGT 1 // -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 position-specific scoring matrix -------------------------------------------------------------------------------- log-odds matrix: alength= 4 w= 12 n= 9291 bayes= 9.17913 E= 8.3e-001 -1064 -188 177 -107 -1064 -88 165 -107 159 -1064 -35 -1064 70 -1064 124 -1064 189 -1064 -1064 -1064 -1064 -88 188 -1064 70 -1064 7 25 -1064 -188 198 -1064 180 -1064 -193 -1064 -52 -1064 139 -48 -111 -1064 188 -1064 106 -1064 7 -48 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 position-specific probability matrix -------------------------------------------------------------------------------- letter-probability matrix: alength= 4 w= 12 nsites= 16 E= 8.3e-001 0.000000 0.062500 0.812500 0.125000 0.000000 0.125000 0.750000 0.125000 0.812500 0.000000 0.187500 0.000000 0.437500 0.000000 0.562500 0.000000 1.000000 0.000000 0.000000 0.000000 0.000000 0.125000 0.875000 0.000000 0.437500 0.000000 0.250000 0.312500 0.000000 0.062500 0.937500 0.000000 0.937500 0.000000 0.062500 0.000000 0.187500 0.000000 0.625000 0.187500 0.125000 0.000000 0.875000 0.000000 0.562500 0.000000 0.250000 0.187500 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 regular expression -------------------------------------------------------------------------------- GGA[GA]AG[ATG]GAGG[AG] -------------------------------------------------------------------------------- Time 3.21 secs. ******************************************************************************** ******************************************************************************** MOTIF 2 MEME width = 12 sites = 19 llr = 174 E-value = 1.5e-002 ******************************************************************************** -------------------------------------------------------------------------------- Motif 2 Description -------------------------------------------------------------------------------- Simplified A :84286638218 pos.-specific C a228:4372391 probability G :::::::1:3:: matrix T ::4:2:1::3:1 bits 2.1 * 1.9 * * 1.7 * * 1.5 * * Relative 1.3 ** ** * * Entropy 1.1 ** *** * ** (13.2 bits) 0.8 ** *** ** ** 0.6 ** ****** ** 0.4 ********* ** 0.2 ********* ** 0.0 ------------ Multilevel CATCAAACACCA consensus CAA CCACG sequence C T -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 2 sites sorted by position p-value -------------------------------------------------------------------------------- Sequence name Start P-value Site ------------- ----- --------- ------------ 25825 412 1.96e-07 GCATAGCACC CAACAAACACCA GCCCCACCGG 261705 204 1.35e-06 ATGCGACCGT CATCAAACAACA CTGATGCAGC 9465 90 2.91e-06 CATTTGGAGG CCTCAAACAGCA GAGCAGGGAG 4115 485 5.04e-06 CAGTCCACGC CCCCAAACACCA AACC 3006 381 5.04e-06 GTCTCCAAAG CACCACACAACA CTTCGCCAAC 261172 11 6.05e-06 ACTCTCAGCC CATCAACAAGCA ATGCATGCCT 6681 285 6.91e-06 TTACACTAGA CAACACAAAGCA GTAGCAGCGA 7921 390 8.23e-06 GGCGACCATG CAACACACCTCA TCTCTGCATC 32971 42 9.07e-06 GTTCTCACAC CACCACACCCCA AGGATCCAAC 2126 477 1.27e-05 CTCCTTCAGC CCTCAAACCCCA GCATCATCCA 23546 417 2.16e-05 TTTCACAGAA CATCACCCAGCC AGCACGCCGG 6848 470 4.25e-05 ACTGATGCCG CAAATAACATCA TTGGTATTGA 22692 443 4.25e-05 CAAGTCAGAG CAAAAACAATCA ATGGATCCAG 2278 248 5.39e-05 CATACCAAAA CAACAACCAGAA TCAATCCACT 20999 82 9.18e-05 AGTATGATCC CATCAATGATCA AGATAGCATG 3987 481 1.05e-04 CTCACTGCTT CACCACCCCTCC CGCCCGAC 3157 287 1.37e-04 GATGTGTTCA CCTAAATCAACA GGAGGATATC 31979 180 1.53e-04 CTATTATTTC CAAATCCAACCA GCGTATAAAA 1946 432 1.78e-04 ACTCTTTCTT CATCTAAAACCT AATTTACCCA -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 2 block diagrams -------------------------------------------------------------------------------- SEQUENCE NAME POSITION P-VALUE MOTIF DIAGRAM ------------- ---------------- ------------- 25825 2e-07 411_[+2]_77 261705 1.3e-06 203_[+2]_285 9465 2.9e-06 89_[+2]_399 4115 5e-06 484_[+2]_4 3006 5e-06 380_[+2]_108 261172 6.1e-06 10_[+2]_478 6681 6.9e-06 284_[+2]_204 7921 8.2e-06 389_[+2]_99 32971 9.1e-06 41_[+2]_447 2126 1.3e-05 476_[+2]_12 23546 2.2e-05 416_[+2]_72 6848 4.2e-05 469_[+2]_19 22692 4.2e-05 442_[+2]_46 2278 5.4e-05 247_[+2]_241 20999 9.2e-05 81_[+2]_407 3987 0.00011 480_[+2]_8 3157 0.00014 286_[+2]_202 31979 0.00015 179_[+2]_309 1946 0.00018 431_[+2]_57 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 2 in BLOCKS format -------------------------------------------------------------------------------- BL MOTIF 2 width=12 seqs=19 25825 ( 412) CAACAAACACCA 1 261705 ( 204) CATCAAACAACA 1 9465 ( 90) CCTCAAACAGCA 1 4115 ( 485) CCCCAAACACCA 1 3006 ( 381) CACCACACAACA 1 261172 ( 11) CATCAACAAGCA 1 6681 ( 285) CAACACAAAGCA 1 7921 ( 390) CAACACACCTCA 1 32971 ( 42) CACCACACCCCA 1 2126 ( 477) CCTCAAACCCCA 1 23546 ( 417) CATCACCCAGCC 1 6848 ( 470) CAAATAACATCA 1 22692 ( 443) CAAAAACAATCA 1 2278 ( 248) CAACAACCAGAA 1 20999 ( 82) CATCAATGATCA 1 3987 ( 481) CACCACCCCTCC 1 3157 ( 287) CCTAAATCAACA 1 31979 ( 180) CAAATCCAACCA 1 1946 ( 432) CATCTAAAACCT 1 // -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 2 position-specific scoring matrix -------------------------------------------------------------------------------- log-odds matrix: alength= 4 w= 12 n= 9291 bayes= 9.12593 E= 1.5e-002 -1089 212 -1089 -1089 155 -13 -1089 -1089 45 -13 -1089 68 -36 178 -1089 -1089 164 -1089 -1089 -73 123 68 -1089 -1089 110 46 -1089 -132 -3 157 -218 -1089 155 -13 -1089 -1089 -77 46 14 0 -235 204 -1089 -1089 164 -113 -1089 -231 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 2 position-specific probability matrix -------------------------------------------------------------------------------- letter-probability matrix: alength= 4 w= 12 nsites= 19 E= 1.5e-002 0.000000 1.000000 0.000000 0.000000 0.789474 0.210526 0.000000 0.000000 0.368421 0.210526 0.000000 0.421053 0.210526 0.789474 0.000000 0.000000 0.842105 0.000000 0.000000 0.157895 0.631579 0.368421 0.000000 0.000000 0.578947 0.315789 0.000000 0.105263 0.263158 0.684211 0.052632 0.000000 0.789474 0.210526 0.000000 0.000000 0.157895 0.315789 0.263158 0.263158 0.052632 0.947368 0.000000 0.000000 0.842105 0.105263 0.000000 0.052632 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 2 regular expression -------------------------------------------------------------------------------- C[AC][TAC][CA]A[AC][AC][CA][AC][CGT]CA -------------------------------------------------------------------------------- Time 6.18 secs. ******************************************************************************** ******************************************************************************** MOTIF 3 MEME width = 21 sites = 11 llr = 161 E-value = 1.6e-002 ******************************************************************************** -------------------------------------------------------------------------------- Motif 3 Description -------------------------------------------------------------------------------- Simplified A 2282:54:32:54:95:5:1a pos.-specific C 48:1a21a1:7:19:16:95: probability G 5:23:3::5213:11:21:5: matrix T :::5::5:26225::5251:: bits 2.1 * * 1.9 * * * 1.7 * * * * * 1.5 * * * ** * * Relative 1.3 ** * * ** * * Entropy 1.1 ** * * * ** * * (21.1 bits) 0.8 ** * * * ** * * * 0.6 *** **** ** ********* 0.4 *** **** ************ 0.2 ********************* 0.0 --------------------- Multilevel GCATCATCGTCATCAACACCA consensus C G GA A GA T T G sequence -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 3 sites sorted by position p-value -------------------------------------------------------------------------------- Sequence name Start P-value Site ------------- ----- --------- --------------------- 4115 171 3.69e-09 CCGTTGATGC CCATCATCATCATCATCATCA TCATCGTCAT 9465 10 1.40e-08 TACACCACC CCAACATCATCATGATCACGA ACTGTTAGCT 2278 319 1.74e-08 ATAATCTCAA AAATCAACGGCATCAACACCA CAATACTAGT 7921 42 2.65e-08 CAAGATTGGA CCATCAACTTCACCACCTCCA ACGAAGCCCC 6848 250 5.78e-08 CAGCTCAAAT ACGGCATCGTCGACATTTCCA GCTTCAAGAA 261172 194 5.78e-08 TAATCTAGTA GAATCAACGTGGACAACTCGA TGGATTCCAC 261705 221 6.32e-08 CAACACTGAT GCAGCCTCAACAACAACACAA AACACACAAC 6681 88 8.25e-08 GAGATCTAGT GCATCGTCTTTTTCAAGACGA CATTATGGAG 32971 59 2.34e-07 CCCCAAGGAT CCAACGCCGACATCGTCTCCA TTCAAAAACT 3987 16 6.38e-07 TTGTCTTGCG GCGGCGACGGCGACATGGCGA GCTGAGTAGC 31979 93 6.38e-07 GCCTTTCGGT GCACCCTCCTTTTCAATTCGA GAGATTCATC -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 3 block diagrams -------------------------------------------------------------------------------- SEQUENCE NAME POSITION P-VALUE MOTIF DIAGRAM ------------- ---------------- ------------- 4115 3.7e-09 170_[+3]_309 9465 1.4e-08 9_[+3]_470 2278 1.7e-08 318_[+3]_161 7921 2.7e-08 41_[+3]_438 6848 5.8e-08 249_[+3]_230 261172 5.8e-08 193_[+3]_286 261705 6.3e-08 220_[+3]_259 6681 8.3e-08 87_[+3]_392 32971 2.3e-07 58_[+3]_421 3987 6.4e-07 15_[+3]_464 31979 6.4e-07 92_[+3]_387 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 3 in BLOCKS format -------------------------------------------------------------------------------- BL MOTIF 3 width=21 seqs=11 4115 ( 171) CCATCATCATCATCATCATCA 1 9465 ( 10) CCAACATCATCATGATCACGA 1 2278 ( 319) AAATCAACGGCATCAACACCA 1 7921 ( 42) CCATCAACTTCACCACCTCCA 1 6848 ( 250) ACGGCATCGTCGACATTTCCA 1 261172 ( 194) GAATCAACGTGGACAACTCGA 1 261705 ( 221) GCAGCCTCAACAACAACACAA 1 6681 ( 88) GCATCGTCTTTTTCAAGACGA 1 32971 ( 59) CCAACGCCGACATCGTCTCCA 1 3987 ( 16) GCGGCGACGGCGACATGGCGA 1 31979 ( 93) GCACCCTCCTTTTCAATTCGA 1 // -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 3 position-specific scoring matrix -------------------------------------------------------------------------------- log-odds matrix: alength= 4 w= 21 n= 9120 bayes= 10.721 E= 1.6e-002 -57 66 93 -1010 -57 183 -1010 -1010 160 -1010 -39 -1010 -57 -134 19 79 -1010 212 -1010 -1010 102 -34 19 -1010 43 -134 -1010 105 -1010 212 -1010 -1010 2 -134 93 -53 -57 -1010 -39 128 -1010 166 -139 -53 102 -1010 19 -53 43 -134 -1010 105 -1010 198 -139 -1010 175 -1010 -139 -1010 75 -134 -1010 79 -1010 147 -39 -53 75 -1010 -139 79 -1010 198 -1010 -153 -157 98 93 -1010 189 -1010 -1010 -1010 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 3 position-specific probability matrix -------------------------------------------------------------------------------- letter-probability matrix: alength= 4 w= 21 nsites= 11 E= 1.6e-002 0.181818 0.363636 0.454545 0.000000 0.181818 0.818182 0.000000 0.000000 0.818182 0.000000 0.181818 0.000000 0.181818 0.090909 0.272727 0.454545 0.000000 1.000000 0.000000 0.000000 0.545455 0.181818 0.272727 0.000000 0.363636 0.090909 0.000000 0.545455 0.000000 1.000000 0.000000 0.000000 0.272727 0.090909 0.454545 0.181818 0.181818 0.000000 0.181818 0.636364 0.000000 0.727273 0.090909 0.181818 0.545455 0.000000 0.272727 0.181818 0.363636 0.090909 0.000000 0.545455 0.000000 0.909091 0.090909 0.000000 0.909091 0.000000 0.090909 0.000000 0.454545 0.090909 0.000000 0.454545 0.000000 0.636364 0.181818 0.181818 0.454545 0.000000 0.090909 0.454545 0.000000 0.909091 0.000000 0.090909 0.090909 0.454545 0.454545 0.000000 1.000000 0.000000 0.000000 0.000000 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 3 regular expression -------------------------------------------------------------------------------- [GC]CA[TG]C[AG][TA]C[GA]TC[AG][TA]CA[AT]C[AT]C[CG]A -------------------------------------------------------------------------------- Time 9.12 secs. ******************************************************************************** ******************************************************************************** SUMMARY OF MOTIFS ******************************************************************************** -------------------------------------------------------------------------------- Combined block diagrams: non-overlapping sites with p-value < 0.0001 -------------------------------------------------------------------------------- SEQUENCE NAME COMBINED P-VALUE MOTIF DIAGRAM ------------- ---------------- ------------- 1946 1.17e-02 119_[+1(7.16e-06)]_132_\ [+1(1.37e-05)]_225 20999 1.45e-01 81_[+2(9.18e-05)]_407 2126 1.61e-04 289_[+1(8.35e-07)]_175_\ [+2(1.27e-05)]_12 22692 7.74e-04 340_[+1(1.37e-05)]_90_\ [+2(4.25e-05)]_46 2278 1.84e-07 59_[+1(7.16e-06)]_176_\ [+2(5.39e-05)]_59_[+3(1.74e-08)]_161 23546 7.74e-03 197_[+1(3.81e-05)]_207_\ [+2(2.16e-05)]_72 25825 8.15e-05 39_[+1(1.79e-05)]_360_\ [+2(1.96e-07)]_77 261172 3.65e-07 10_[+2(6.05e-06)]_171_\ [+3(5.78e-08)]_49_[+1(4.07e-05)]_225 261705 4.49e-08 41_[+1(1.66e-05)]_150_\ [+2(1.35e-06)]_5_[+3(6.32e-08)]_221_[+3(3.54e-05)]_17 3006 6.75e-04 380_[+2(5.04e-06)]_9_[+1(4.62e-05)]_\ 87 3157 5.26e-03 74_[+1(9.64e-06)]_414 31979 1.61e-07 92_[+3(6.38e-07)]_351_\ [+1(6.09e-08)]_24 32971 4.07e-08 41_[+2(9.07e-06)]_5_[+3(2.34e-07)]_\ 253_[+1(5.98e-07)]_89_[+3(5.00e-05)]_46 3987 5.28e-06 15_[+3(6.38e-07)]_88_[+1(4.19e-06)]_\ 364 4115 3.01e-08 18_[+1(4.96e-05)]_140_\ [+3(3.69e-09)]_293_[+2(5.04e-06)]_4 6681 2.00e-07 87_[+3(8.25e-08)]_130_\ [+1(1.27e-05)]_34_[+2(6.91e-06)]_204 6848 4.94e-05 249_[+3(5.78e-08)]_199_\ [+2(4.25e-05)]_19 7921 5.48e-08 41_[+3(2.65e-08)]_327_\ [+2(8.23e-06)]_56_[+1(8.05e-06)]_31 9465 1.29e-06 9_[+3(1.40e-08)]_59_[+2(2.91e-06)]_\ 399 -------------------------------------------------------------------------------- ******************************************************************************** ******************************************************************************** Stopped because nmotifs = 3 reached. ******************************************************************************** CPU: seaotter.hsd1.wa.comcast.net ********************************************************************************