******************************************************************************** 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/351/351.seqs.fa ALPHABET= ACGT Sequence name Weight Length Sequence name Weight Length ------------- ------ ------ ------------- ------ ------ 8760 1.0000 500 9266 1.0000 500 36981 1.0000 500 54757 1.0000 500 37779 1.0000 500 9947 1.0000 500 43478 1.0000 500 9912 1.0000 500 49596 1.0000 500 50133 1.0000 500 41624 1.0000 500 44566 1.0000 500 35158 1.0000 500 45852 1.0000 500 35669 1.0000 500 48314 1.0000 500 44963 1.0000 500 36843 1.0000 500 45746 1.0000 500 49110 1.0000 500 48130 1.0000 500 37659 1.0000 500 50127 1.0000 500 7073 1.0000 500 44962 1.0000 500 55100 1.0000 500 44195 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/351/351.seqs.fa -oc motifs/351 -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= 27 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= 13500 N= 27 strands: + sample: seed= 0 seqfrac= 1 Letter frequencies in dataset: A 0.280 C 0.241 G 0.215 T 0.264 Background letter frequencies (from dataset with add-one prior applied): A 0.280 C 0.241 G 0.215 T 0.264 ******************************************************************************** ******************************************************************************** MOTIF 1 MEME width = 21 sites = 16 llr = 204 E-value = 8.1e-001 ******************************************************************************** -------------------------------------------------------------------------------- Motif 1 Description -------------------------------------------------------------------------------- Simplified A 1::13316111:11:146:11 pos.-specific C 1:11843:16981:34342:: probability G 4113::7:1::1:3121128: matrix T 4986:3:483:186633:629 bits 2.2 2.0 1.8 1.6 * * * Relative 1.3 * * * Entropy 1.1 ** * * *** ** (18.4 bits) 0.9 ** * ******* * ** 0.7 ** * ********* **** 0.4 *************** **** 0.2 **************** **** 0.0 --------------------- Multilevel TTTTCCGATCCCTTTCAATGT consensus G GAACT T GCTCC sequence T T -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 sites sorted by position p-value -------------------------------------------------------------------------------- Sequence name Start P-value Site ------------- ----- --------- --------------------- 44195 59 9.51e-09 AGTTTTCACC TTTTCCCATCCCCTTTACTGT TAGCGTATTT 9266 208 1.22e-08 ATCGATCATT GTTGCAGTTCCCTATCCACGT AAAGCCTCCC 43478 474 3.99e-08 CAGTACGCTA TTTTCAGAACCCTTCCAACGT TTGATC 44962 419 4.98e-08 TGCAGCTTTC TTGTCCGATCCCTTTCGCTTT GAAGAGCTAG 55100 267 6.88e-08 CCACGTGGAA GGTGATGTTCCCTTTCTATGT TTTTGAATCG 44963 130 1.41e-07 TGCCTACGAT GTCGCAGTTTCCTTTTAATTT CGACGTGGCG 41624 213 3.97e-07 TGAGAAAGAA CTTTACATTCCCTTCGTATGT CAGGCAAGGA 9912 177 4.34e-07 AAAACAAGTT GTTCACCATTCCTGTCCAGGT TCGAAACGAA 45746 478 6.69e-07 AAAATTCGTC GTTGACGTTCACTGCCTGTGT CC 35158 406 7.92e-07 TTTTTTCTCC ATTTCAGATCCCTTCAACTAT TCTCACAATT 9947 91 1.51e-06 CAAGCCCTCG TTCGCTGTTTCCATTTCAGGT TCCCAAACCT 35669 112 1.76e-06 GAAGTAGGTC CTTTCCGATTCTTGGGTATGT TCTGGCGAAA 36981 366 1.76e-06 AAGAATTTCT TTTTCAGACCCGCTTTCACGT ATTCTTTTTT 48130 161 2.04e-06 CAGAACAGTG GTTTCCCAGTACTTCGGCTGT CGACGGTGCA 7073 407 2.55e-06 TGAAACCTCA TTTACTGTTACCTATAACTGT GCCGTCCTGT 48314 461 5.88e-06 TAGGGTCTGG TTTTCTCATCCGTGTCACGTA TAATTCGCTT -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 block diagrams -------------------------------------------------------------------------------- SEQUENCE NAME POSITION P-VALUE MOTIF DIAGRAM ------------- ---------------- ------------- 44195 9.5e-09 58_[+1]_421 9266 1.2e-08 207_[+1]_272 43478 4e-08 473_[+1]_6 44962 5e-08 418_[+1]_61 55100 6.9e-08 266_[+1]_213 44963 1.4e-07 129_[+1]_350 41624 4e-07 212_[+1]_267 9912 4.3e-07 176_[+1]_303 45746 6.7e-07 477_[+1]_2 35158 7.9e-07 405_[+1]_74 9947 1.5e-06 90_[+1]_389 35669 1.8e-06 111_[+1]_368 36981 1.8e-06 365_[+1]_114 48130 2e-06 160_[+1]_319 7073 2.6e-06 406_[+1]_73 48314 5.9e-06 460_[+1]_19 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 in BLOCKS format -------------------------------------------------------------------------------- BL MOTIF 1 width=21 seqs=16 44195 ( 59) TTTTCCCATCCCCTTTACTGT 1 9266 ( 208) GTTGCAGTTCCCTATCCACGT 1 43478 ( 474) TTTTCAGAACCCTTCCAACGT 1 44962 ( 419) TTGTCCGATCCCTTTCGCTTT 1 55100 ( 267) GGTGATGTTCCCTTTCTATGT 1 44963 ( 130) GTCGCAGTTTCCTTTTAATTT 1 41624 ( 213) CTTTACATTCCCTTCGTATGT 1 9912 ( 177) GTTCACCATTCCTGTCCAGGT 1 45746 ( 478) GTTGACGTTCACTGCCTGTGT 1 35158 ( 406) ATTTCAGATCCCTTCAACTAT 1 9947 ( 91) TTCGCTGTTTCCATTTCAGGT 1 35669 ( 112) CTTTCCGATTCTTGGGTATGT 1 36981 ( 366) TTTTCAGACCCGCTTTCACGT 1 48130 ( 161) GTTTCCCAGTACTTCGGCTGT 1 7073 ( 407) TTTACTGTTACCTATAACTGT 1 48314 ( 461) TTTTCTCATCCGTGTCACGTA 1 // -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 position-specific scoring matrix -------------------------------------------------------------------------------- log-odds matrix: alength= 4 w= 21 n= 12960 bayes= 9.66 E= 8.1e-001 -216 -94 80 73 -1064 -1064 -178 183 -1064 -94 -178 162 -216 -194 54 109 -16 164 -1064 -1064 16 86 -1064 -8 -216 5 168 -1064 100 -1064 -1064 73 -216 -194 -178 162 -216 138 -1064 24 -116 186 -1064 -1064 -1064 175 -78 -208 -216 -94 -1064 162 -116 -1064 22 124 -1064 38 -178 124 -116 86 -20 -8 42 5 -78 -8 100 64 -178 -1064 -1064 -36 -20 124 -216 -1064 180 -49 -216 -1064 -1064 183 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 position-specific probability matrix -------------------------------------------------------------------------------- letter-probability matrix: alength= 4 w= 21 nsites= 16 E= 8.1e-001 0.062500 0.125000 0.375000 0.437500 0.000000 0.000000 0.062500 0.937500 0.000000 0.125000 0.062500 0.812500 0.062500 0.062500 0.312500 0.562500 0.250000 0.750000 0.000000 0.000000 0.312500 0.437500 0.000000 0.250000 0.062500 0.250000 0.687500 0.000000 0.562500 0.000000 0.000000 0.437500 0.062500 0.062500 0.062500 0.812500 0.062500 0.625000 0.000000 0.312500 0.125000 0.875000 0.000000 0.000000 0.000000 0.812500 0.125000 0.062500 0.062500 0.125000 0.000000 0.812500 0.125000 0.000000 0.250000 0.625000 0.000000 0.312500 0.062500 0.625000 0.125000 0.437500 0.187500 0.250000 0.375000 0.250000 0.125000 0.250000 0.562500 0.375000 0.062500 0.000000 0.000000 0.187500 0.187500 0.625000 0.062500 0.000000 0.750000 0.187500 0.062500 0.000000 0.000000 0.937500 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 regular expression -------------------------------------------------------------------------------- [TG]TT[TG][CA][CAT][GC][AT]T[CT]CCT[TG][TC][CT][ACT][AC]TGT -------------------------------------------------------------------------------- Time 6.02 secs. ******************************************************************************** ******************************************************************************** MOTIF 2 MEME width = 12 sites = 12 llr = 133 E-value = 2.5e+001 ******************************************************************************** -------------------------------------------------------------------------------- Motif 2 Description -------------------------------------------------------------------------------- Simplified A 1::91::1163: pos.-specific C 9:a181::9:52 probability G :3::::a::322 matrix T :8::19:9:217 bits 2.2 * 2.0 * * 1.8 * * 1.6 * * **** Relative 1.3 * ** **** Entropy 1.1 ********* (16.0 bits) 0.9 ********* 0.7 ********** * 0.4 ********** * 0.2 ************ 0.0 ------------ Multilevel CTCACTGTCACT consensus G GA sequence -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 2 sites sorted by position p-value -------------------------------------------------------------------------------- Sequence name Start P-value Site ------------- ----- --------- ------------ 43478 412 6.66e-08 TCGGCCCGAG CTCACTGTCACT CCAACCACCT 54757 80 2.49e-07 CAAATCCTGT CTCACTGTCAAT AAGGACAAAA 44566 396 3.09e-07 AAAGTAGTGG CTCACTGTCAGT GTCTTCGTAT 36981 214 4.86e-07 GCAGGAGGCA CTCACTGTCACC ATCGTCAAGA 48130 370 5.64e-06 CTGTGGCGGA CTCATTGTCTCT CGTTTTTGTT 36843 360 5.64e-06 GAGAAAACAG ATCACTGTCAGT GCAAGCGAAC 35669 273 6.72e-06 GACCCTAATC CTCAATGTCACC TGTGCTCATT 8760 445 7.43e-06 ACGTCGGATC CGCACTGTCTTT GCGAATCCGT 49110 479 7.86e-06 CATCTTCTTG CTCACTGTAGAT CAAAAAGACC 44963 192 8.26e-06 TTCAACCAAC CTCACCGTCGCG CAATCAGTGT 7073 109 1.21e-05 CGTAAGATTT CGCCCTGTCGAT GACGCCGTCC 50127 119 1.21e-05 CTTGCGGATG CGCACTGACACG ATTGACCGCG -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 2 block diagrams -------------------------------------------------------------------------------- SEQUENCE NAME POSITION P-VALUE MOTIF DIAGRAM ------------- ---------------- ------------- 43478 6.7e-08 411_[+2]_77 54757 2.5e-07 79_[+2]_409 44566 3.1e-07 395_[+2]_93 36981 4.9e-07 213_[+2]_275 48130 5.6e-06 369_[+2]_119 36843 5.6e-06 359_[+2]_129 35669 6.7e-06 272_[+2]_216 8760 7.4e-06 444_[+2]_44 49110 7.9e-06 478_[+2]_10 44963 8.3e-06 191_[+2]_297 7073 1.2e-05 108_[+2]_380 50127 1.2e-05 118_[+2]_370 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 2 in BLOCKS format -------------------------------------------------------------------------------- BL MOTIF 2 width=12 seqs=12 43478 ( 412) CTCACTGTCACT 1 54757 ( 80) CTCACTGTCAAT 1 44566 ( 396) CTCACTGTCAGT 1 36981 ( 214) CTCACTGTCACC 1 48130 ( 370) CTCATTGTCTCT 1 36843 ( 360) ATCACTGTCAGT 1 35669 ( 273) CTCAATGTCACC 1 8760 ( 445) CGCACTGTCTTT 1 49110 ( 479) CTCACTGTAGAT 1 44963 ( 192) CTCACCGTCGCG 1 7073 ( 109) CGCCCTGTCGAT 1 50127 ( 119) CGCACTGACACG 1 // -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 2 position-specific scoring matrix -------------------------------------------------------------------------------- log-odds matrix: alength= 4 w= 12 n= 13203 bayes= 10.5501 E= 2.5e+001 -175 193 -1023 -1023 -1023 -1023 22 151 -1023 205 -1023 -1023 171 -153 -1023 -1023 -175 179 -1023 -166 -1023 -153 -1023 180 -1023 -1023 222 -1023 -175 -1023 -1023 180 -175 193 -1023 -1023 106 -1023 22 -66 -16 105 -37 -166 -1023 -53 -37 134 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 2 position-specific probability matrix -------------------------------------------------------------------------------- letter-probability matrix: alength= 4 w= 12 nsites= 12 E= 2.5e+001 0.083333 0.916667 0.000000 0.000000 0.000000 0.000000 0.250000 0.750000 0.000000 1.000000 0.000000 0.000000 0.916667 0.083333 0.000000 0.000000 0.083333 0.833333 0.000000 0.083333 0.000000 0.083333 0.000000 0.916667 0.000000 0.000000 1.000000 0.000000 0.083333 0.000000 0.000000 0.916667 0.083333 0.916667 0.000000 0.000000 0.583333 0.000000 0.250000 0.166667 0.250000 0.500000 0.166667 0.083333 0.000000 0.166667 0.166667 0.666667 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 2 regular expression -------------------------------------------------------------------------------- C[TG]CACTGTC[AG][CA]T -------------------------------------------------------------------------------- Time 11.87 secs. ******************************************************************************** ******************************************************************************** MOTIF 3 MEME width = 21 sites = 10 llr = 150 E-value = 2.2e+002 ******************************************************************************** -------------------------------------------------------------------------------- Motif 3 Description -------------------------------------------------------------------------------- Simplified A 29923:8851::235:334:: pos.-specific C 1:17:3214::a:31::72:3 probability G 71:167:1:2::83:8::355 matrix T ::::1:::17a::1427:152 bits 2.2 2.0 ** 1.8 ** 1.6 ** Relative 1.3 ** * *** * Entropy 1.1 ** ** *** *** * (21.7 bits) 0.9 ******** **** *** * 0.7 ************* *** ** 0.4 ************* **** ** 0.2 ********************* 0.0 --------------------- Multilevel GAACGGAAATTCGAAGTCAGG consensus A AACC CG ACTTAAGTC sequence G C T -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 3 sites sorted by position p-value -------------------------------------------------------------------------------- Sequence name Start P-value Site ------------- ----- --------- --------------------- 45746 65 2.12e-10 TTGGATCGGA GAACGGAATTTCGATGTCAGG ATCTCCCATT 43478 386 5.69e-09 CGCTCACCGC GAACGGAACTTCACATTCGGC CCGAGCTCAC 44195 345 2.99e-08 TGGCTGCCAG GAAGTGCAATTCGCAGTCGTG GGAATCATAT 54757 332 4.30e-08 CAGATTGAGA GAACAGAAAGTCATTGTCCGG AATGGCCATA 37779 267 6.10e-08 GACAAAAGCT AAACACAACTTCGGAGTCTTT GACGTCCAAA 45852 256 7.84e-08 CGCCGTCACC GACCGCAAATTCGCAGAACTC GTAGGCGCGA 36981 327 1.17e-07 CGACAAAAAA GAAAGGAACGTCGATTAAATG AAATCTCGAA 50127 29 1.36e-07 TATCGGGGCG AGACGGAAATTCGGCGAAAGG TTAGTGCCAG 9947 136 1.97e-07 TATTGACTTG GAACGCCGAATCGGAGTCGTT CTTCGCAACT 44566 194 2.42e-07 TGCAAGTCTC CAAAAGACCTTCGATGTCAGC CTCGCACGTT -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 3 block diagrams -------------------------------------------------------------------------------- SEQUENCE NAME POSITION P-VALUE MOTIF DIAGRAM ------------- ---------------- ------------- 45746 2.1e-10 64_[+3]_415 43478 5.7e-09 385_[+3]_94 44195 3e-08 344_[+3]_135 54757 4.3e-08 331_[+3]_148 37779 6.1e-08 266_[+3]_213 45852 7.8e-08 255_[+3]_224 36981 1.2e-07 326_[+3]_153 50127 1.4e-07 28_[+3]_451 9947 2e-07 135_[+3]_344 44566 2.4e-07 193_[+3]_286 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 3 in BLOCKS format -------------------------------------------------------------------------------- BL MOTIF 3 width=21 seqs=10 45746 ( 65) GAACGGAATTTCGATGTCAGG 1 43478 ( 386) GAACGGAACTTCACATTCGGC 1 44195 ( 345) GAAGTGCAATTCGCAGTCGTG 1 54757 ( 332) GAACAGAAAGTCATTGTCCGG 1 37779 ( 267) AAACACAACTTCGGAGTCTTT 1 45852 ( 256) GACCGCAAATTCGCAGAACTC 1 36981 ( 327) GAAAGGAACGTCGATTAAATG 1 50127 ( 29) AGACGGAAATTCGGCGAAAGG 1 9947 ( 136) GAACGCCGAATCGGAGTCGTT 1 44566 ( 194) CAAAAGACCTTCGATGTCAGC 1 // -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 3 position-specific scoring matrix -------------------------------------------------------------------------------- log-odds matrix: alength= 4 w= 21 n= 12960 bayes= 10.5905 E= 2.2e+002 -49 -127 170 -997 168 -997 -110 -997 168 -127 -997 -997 -49 154 -110 -997 10 -997 148 -140 -997 32 170 -997 151 -27 -997 -997 151 -127 -110 -997 83 73 -997 -140 -148 -997 -11 141 -997 -997 -997 192 -997 205 -997 -997 -49 -997 189 -997 10 32 48 -140 83 -127 -997 60 -997 -997 189 -40 10 -997 -997 141 10 154 -997 -997 51 -27 48 -140 -997 -997 122 92 -997 32 122 -40 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 3 position-specific probability matrix -------------------------------------------------------------------------------- letter-probability matrix: alength= 4 w= 21 nsites= 10 E= 2.2e+002 0.200000 0.100000 0.700000 0.000000 0.900000 0.000000 0.100000 0.000000 0.900000 0.100000 0.000000 0.000000 0.200000 0.700000 0.100000 0.000000 0.300000 0.000000 0.600000 0.100000 0.000000 0.300000 0.700000 0.000000 0.800000 0.200000 0.000000 0.000000 0.800000 0.100000 0.100000 0.000000 0.500000 0.400000 0.000000 0.100000 0.100000 0.000000 0.200000 0.700000 0.000000 0.000000 0.000000 1.000000 0.000000 1.000000 0.000000 0.000000 0.200000 0.000000 0.800000 0.000000 0.300000 0.300000 0.300000 0.100000 0.500000 0.100000 0.000000 0.400000 0.000000 0.000000 0.800000 0.200000 0.300000 0.000000 0.000000 0.700000 0.300000 0.700000 0.000000 0.000000 0.400000 0.200000 0.300000 0.100000 0.000000 0.000000 0.500000 0.500000 0.000000 0.300000 0.500000 0.200000 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 3 regular expression -------------------------------------------------------------------------------- [GA]AA[CA][GA][GC][AC]A[AC][TG]TC[GA][ACG][AT][GT][TA][CA][AGC][GT][GCT] -------------------------------------------------------------------------------- Time 17.86 secs. ******************************************************************************** ******************************************************************************** SUMMARY OF MOTIFS ******************************************************************************** -------------------------------------------------------------------------------- Combined block diagrams: non-overlapping sites with p-value < 0.0001 -------------------------------------------------------------------------------- SEQUENCE NAME COMBINED P-VALUE MOTIF DIAGRAM ------------- ---------------- ------------- 8760 6.59e-03 2_[+2(9.50e-05)]_430_[+2(7.43e-06)]_\ 44 9266 4.08e-05 207_[+1(1.22e-08)]_272 36981 3.87e-09 213_[+2(4.86e-07)]_101_\ [+3(1.17e-07)]_18_[+1(1.76e-06)]_114 54757 2.58e-07 79_[+2(2.49e-07)]_240_\ [+3(4.30e-08)]_148 37779 3.15e-04 266_[+3(6.10e-08)]_213 9947 2.73e-06 90_[+1(1.51e-06)]_24_[+3(1.97e-07)]_\ 344 43478 1.04e-12 385_[+3(5.69e-09)]_5_[+2(6.66e-08)]_\ 50_[+1(3.99e-08)]_6 9912 7.68e-04 176_[+1(4.34e-07)]_303 49596 5.17e-01 500 50133 8.40e-01 500 41624 6.56e-04 212_[+1(3.97e-07)]_267 44566 6.42e-07 193_[+3(2.42e-07)]_181_\ [+2(3.09e-07)]_93 35158 3.41e-03 405_[+1(7.92e-07)]_74 45852 7.18e-04 255_[+3(7.84e-08)]_224 35669 1.73e-04 111_[+1(1.76e-06)]_140_\ [+2(6.72e-06)]_216 48314 2.33e-02 460_[+1(5.88e-06)]_19 44963 2.50e-06 129_[+1(1.41e-07)]_41_\ [+2(8.26e-06)]_297 36843 1.06e-02 359_[+2(5.64e-06)]_129 45746 7.21e-09 64_[+3(2.12e-10)]_392_\ [+1(6.69e-07)]_2 49110 5.18e-02 478_[+2(7.86e-06)]_10 48130 1.47e-04 160_[+1(2.04e-06)]_188_\ [+2(5.64e-06)]_119 37659 6.93e-01 500 50127 4.03e-05 28_[+3(1.36e-07)]_69_[+2(1.21e-05)]_\ 370 7073 2.82e-04 108_[+2(1.21e-05)]_286_\ [+1(2.55e-06)]_73 44962 6.11e-04 418_[+1(4.98e-08)]_61 55100 3.66e-04 266_[+1(6.88e-08)]_213 44195 1.53e-08 58_[+1(9.51e-09)]_265_\ [+3(2.99e-08)]_135 -------------------------------------------------------------------------------- ******************************************************************************** ******************************************************************************** Stopped because nmotifs = 3 reached. ******************************************************************************** CPU: seaotter.hsd1.wa.comcast.net ********************************************************************************