******************************************************************************** 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 ------------- ------ ------ ------------- ------ ------ 10034 1.0000 500 1005 1.0000 500 10831 1.0000 500 11227 1.0000 500 11795 1.0000 500 12131 1.0000 500 2001 1.0000 500 23224 1.0000 500 24409 1.0000 500 261704 1.0000 500 263913 1.0000 500 2703 1.0000 500 4190 1.0000 500 4814 1.0000 500 5801 1.0000 500 6620 1.0000 500 6764 1.0000 500 9475 1.0000 500 bd887 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= 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.275 C 0.244 G 0.217 T 0.264 Background letter frequencies (from dataset with add-one prior applied): A 0.275 C 0.244 G 0.217 T 0.264 ******************************************************************************** ******************************************************************************** MOTIF 1 MEME width = 21 sites = 10 llr = 152 E-value = 3.3e-002 ******************************************************************************** -------------------------------------------------------------------------------- Motif 1 Description -------------------------------------------------------------------------------- Simplified A :9:::5:3123:72115::2: pos.-specific C ::21:::::1::113:1::1: probability G a1839134747a253542:29 matrix T :::6147323:::234:8a51 bits 2.2 * * 2.0 * * * 1.8 * * * * * 1.5 * * * * * * Relative 1.3 *** * ** ** * Entropy 1.1 *** * * ** ** * (21.9 bits) 0.9 *** * * * *** ** * 0.7 ******* * *** **** * 0.4 ********* *** **** * 0.2 ********************* 0.0 --------------------- Multilevel GAGTGATGGGGGAGCGATTTG consensus CG TGATTA GAGTGG A sequence T A TT G -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 sites sorted by position p-value -------------------------------------------------------------------------------- Sequence name Start P-value Site ------------- ----- --------- --------------------- bd887 74 1.68e-10 TAGAGATGAA GAGGGATAGAGGAGCGATTTG TAACGATTAT 2703 233 1.35e-08 TCGATTAAAG GAGGGATGAGGGGGCGATTAG TTGGCCATGG 2001 280 1.67e-08 GAGAGGTAGC GAGCGAGAGAGGAAGGATTTG GGGCGCAGCT 11227 172 3.29e-08 TTGATTGGAC GAGTTTGAGTAGAGTGATTTG GTTTAAAGAT 10034 52 3.61e-08 TGTGGCTTGT GACTGATGGTAGATGGCTTTG GAGAGGAGCC 12131 98 8.50e-08 TGCCAGCGAT GAGTGTTTGCGGCTTTGTTGG CACCTGCAGA 23224 1 1.08e-07 . GAGTGATGGTAGAAGAGGTAG AAGATAGTTT 1005 118 1.35e-07 GGCTTGAGCG GAGGGTGGTGGGGGCTGGTCG TCGGATTTGC 24409 205 1.56e-07 GCTCGCCATT GGCTGTTTGGGGAGTTGTTTT GCATGATTAG 6764 292 2.07e-07 GCCTCCGTAG GAGTGGTTTGGGACATATTGG TGGTTTGTTT -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 block diagrams -------------------------------------------------------------------------------- SEQUENCE NAME POSITION P-VALUE MOTIF DIAGRAM ------------- ---------------- ------------- bd887 1.7e-10 73_[+1]_406 2703 1.4e-08 232_[+1]_247 2001 1.7e-08 279_[+1]_200 11227 3.3e-08 171_[+1]_308 10034 3.6e-08 51_[+1]_428 12131 8.5e-08 97_[+1]_382 23224 1.1e-07 [+1]_479 1005 1.3e-07 117_[+1]_362 24409 1.6e-07 204_[+1]_275 6764 2.1e-07 291_[+1]_188 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 in BLOCKS format -------------------------------------------------------------------------------- BL MOTIF 1 width=21 seqs=10 bd887 ( 74) GAGGGATAGAGGAGCGATTTG 1 2703 ( 233) GAGGGATGAGGGGGCGATTAG 1 2001 ( 280) GAGCGAGAGAGGAAGGATTTG 1 11227 ( 172) GAGTTTGAGTAGAGTGATTTG 1 10034 ( 52) GACTGATGGTAGATGGCTTTG 1 12131 ( 98) GAGTGTTTGCGGCTTTGTTGG 1 23224 ( 1) GAGTGATGGTAGAAGAGGTAG 1 1005 ( 118) GAGGGTGGTGGGGGCTGGTCG 1 24409 ( 205) GGCTGTTTGGGGAGTTGTTTT 1 6764 ( 292) GAGTGGTTTGGGACATATTGG 1 // -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 position-specific scoring matrix -------------------------------------------------------------------------------- log-odds matrix: alength= 4 w= 21 n= 9120 bayes= 10.0831 E= 3.3e-002 -997 -997 220 -997 171 -997 -112 -997 -997 -29 188 -997 -997 -129 47 118 -997 -997 205 -140 86 -997 -112 60 -997 -997 47 140 13 -997 88 18 -145 -997 169 -40 -46 -129 88 18 13 -997 169 -997 -997 -997 220 -997 135 -129 -12 -997 -46 -129 120 -40 -145 30 47 18 -145 -997 120 60 86 -129 88 -997 -997 -997 -12 160 -997 -997 -997 192 -46 -129 -12 92 -997 -997 205 -140 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 position-specific probability matrix -------------------------------------------------------------------------------- letter-probability matrix: alength= 4 w= 21 nsites= 10 E= 3.3e-002 0.000000 0.000000 1.000000 0.000000 0.900000 0.000000 0.100000 0.000000 0.000000 0.200000 0.800000 0.000000 0.000000 0.100000 0.300000 0.600000 0.000000 0.000000 0.900000 0.100000 0.500000 0.000000 0.100000 0.400000 0.000000 0.000000 0.300000 0.700000 0.300000 0.000000 0.400000 0.300000 0.100000 0.000000 0.700000 0.200000 0.200000 0.100000 0.400000 0.300000 0.300000 0.000000 0.700000 0.000000 0.000000 0.000000 1.000000 0.000000 0.700000 0.100000 0.200000 0.000000 0.200000 0.100000 0.500000 0.200000 0.100000 0.300000 0.300000 0.300000 0.100000 0.000000 0.500000 0.400000 0.500000 0.100000 0.400000 0.000000 0.000000 0.000000 0.200000 0.800000 0.000000 0.000000 0.000000 1.000000 0.200000 0.100000 0.200000 0.500000 0.000000 0.000000 0.900000 0.100000 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 regular expression -------------------------------------------------------------------------------- GA[GC][TG]G[AT][TG][GAT][GT][GTA][GA]G[AG][GAT][CGT][GT][AG][TG]T[TAG]G -------------------------------------------------------------------------------- Time 3.42 secs. ******************************************************************************** ******************************************************************************** MOTIF 2 MEME width = 21 sites = 19 llr = 212 E-value = 6.0e-002 ******************************************************************************** -------------------------------------------------------------------------------- Motif 2 Description -------------------------------------------------------------------------------- Simplified A :21a346:4:44427417146 pos.-specific C 338:7436382166:41:941 probability G 121::211:1:3:2:22::22 matrix T 74:::::43142::3163::1 bits 2.2 2.0 1.8 * * 1.5 * * Relative 1.3 * * * Entropy 1.1 *** * * * (16.1 bits) 0.9 * *** * * * * ** 0.7 * *** ** * *** ** 0.4 * ********* *** ***** 0.2 * ******************* 0.0 --------------------- Multilevel TTCACAACACAACCAATACAA consensus CC ACCTC TGAATCGT CG sequence G T T G G -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 2 sites sorted by position p-value -------------------------------------------------------------------------------- Sequence name Start P-value Site ------------- ----- --------- --------------------- 4814 400 7.75e-11 AAAGAATCGA TTCACCACCCTGCCAATACCA CAGTATACCT bd887 341 1.23e-07 AGGACTATTT TTCACAATACAACAACAACAA TGGTGAATTG 1005 444 4.54e-07 ACAGCAACAC TCCACCATACAAACTCTTCAT CGTGTTCGAC 261704 120 5.12e-07 TAGGGCTTAT CTAACACTCCTGCCACTACAA GCTGCCCTAT 2001 110 6.45e-07 TCCCAAGCCG TCCAACATACTTCGTGTACAA GGATCGTTGT 9475 301 7.23e-07 CTCTATCTCA TCCACCACACCGCATACACAA TGCAACTCTT 6764 191 1.12e-06 TTTCTTCTCA CTCACACCTCAGCCACATCGA ATCAACTTGG 2703 478 1.52e-06 GACGTCGGCC GCCACAACACAAACTCGTCCA CC 11795 478 1.52e-06 AACGCATCCC TTCAAGACATCAACAATACCA TA 263913 260 2.72e-06 AGCAAAACAT TGCAAACCCTATCCAGTACCA CACGTACACT 24409 73 3.26e-06 CAGGAACGAC CACACCACAGTGACAATACCG CTTCCCGTCT 4190 415 3.56e-06 GCCTCCGGCT TCGACAACCCTAACTTTACAG CCACACCCAC 10831 158 6.42e-06 ATCAACACCT TTCAAGCTTCAAAGAGTACGG TCGCATCCTC 5801 342 8.77e-06 CAGAATCCAT CACACAACTCTTCGTCTTCAT CACATCACGA 10034 81 1.02e-05 TGGAGAGGAG CCGACGGTTCAACCACGTCAA CGTGTTGTTT 6620 380 1.10e-05 AACGTCACAA TTCAACACACTACCAAGTACC AAATGACCCC 11227 1 1.26e-05 . TGAACCAGTCATCAAATACCA TGAGATTACG 23224 319 1.66e-05 TTCACAATTT TGCACAGTCCTCCCAAGACGC GACTTTGGCC 12131 259 4.32e-05 GAAAGCCTTT TACAAGCCCCCGAAAGCACGG TGATAACCGA -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 2 block diagrams -------------------------------------------------------------------------------- SEQUENCE NAME POSITION P-VALUE MOTIF DIAGRAM ------------- ---------------- ------------- 4814 7.7e-11 399_[+2]_80 bd887 1.2e-07 340_[+2]_139 1005 4.5e-07 443_[+2]_36 261704 5.1e-07 119_[+2]_360 2001 6.5e-07 109_[+2]_370 9475 7.2e-07 300_[+2]_179 6764 1.1e-06 190_[+2]_289 2703 1.5e-06 477_[+2]_2 11795 1.5e-06 477_[+2]_2 263913 2.7e-06 259_[+2]_220 24409 3.3e-06 72_[+2]_407 4190 3.6e-06 414_[+2]_65 10831 6.4e-06 157_[+2]_322 5801 8.8e-06 341_[+2]_138 10034 1e-05 80_[+2]_399 6620 1.1e-05 379_[+2]_100 11227 1.3e-05 [+2]_479 23224 1.7e-05 318_[+2]_161 12131 4.3e-05 258_[+2]_221 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 2 in BLOCKS format -------------------------------------------------------------------------------- BL MOTIF 2 width=21 seqs=19 4814 ( 400) TTCACCACCCTGCCAATACCA 1 bd887 ( 341) TTCACAATACAACAACAACAA 1 1005 ( 444) TCCACCATACAAACTCTTCAT 1 261704 ( 120) CTAACACTCCTGCCACTACAA 1 2001 ( 110) TCCAACATACTTCGTGTACAA 1 9475 ( 301) TCCACCACACCGCATACACAA 1 6764 ( 191) CTCACACCTCAGCCACATCGA 1 2703 ( 478) GCCACAACACAAACTCGTCCA 1 11795 ( 478) TTCAAGACATCAACAATACCA 1 263913 ( 260) TGCAAACCCTATCCAGTACCA 1 24409 ( 73) CACACCACAGTGACAATACCG 1 4190 ( 415) TCGACAACCCTAACTTTACAG 1 10831 ( 158) TTCAAGCTTCAAAGAGTACGG 1 5801 ( 342) CACACAACTCTTCGTCTTCAT 1 10034 ( 81) CCGACGGTTCAACCACGTCAA 1 6620 ( 380) TTCAACACACTACCAAGTACC 1 11227 ( 1) TGAACCAGTCATCAAATACCA 1 23224 ( 319) TGCACAGTCCTCCCAAGACGC 1 12131 ( 259) TACAAGCCCCCGAAAGCACGG 1 // -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 2 position-specific scoring matrix -------------------------------------------------------------------------------- log-odds matrix: alength= 4 w= 21 n= 9120 bayes= 9.09908 E= 6.0e-002 -1089 11 -204 137 -80 37 -46 48 -138 169 -104 -1089 186 -1089 -1089 -1089 20 149 -1089 -1089 62 59 -4 -1089 120 11 -104 -1089 -1089 124 -204 48 62 37 -1089 -1 -1089 179 -204 -133 62 -63 -1089 67 62 -221 54 -33 42 137 -1089 -1089 -38 137 -46 -1089 132 -1089 -1089 26 42 59 -4 -233 -138 -121 -4 113 132 -1089 -1089 26 -238 196 -1089 -1089 62 59 -4 -1089 108 -121 -4 -133 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 2 position-specific probability matrix -------------------------------------------------------------------------------- letter-probability matrix: alength= 4 w= 21 nsites= 19 E= 6.0e-002 0.000000 0.263158 0.052632 0.684211 0.157895 0.315789 0.157895 0.368421 0.105263 0.789474 0.105263 0.000000 1.000000 0.000000 0.000000 0.000000 0.315789 0.684211 0.000000 0.000000 0.421053 0.368421 0.210526 0.000000 0.631579 0.263158 0.105263 0.000000 0.000000 0.578947 0.052632 0.368421 0.421053 0.315789 0.000000 0.263158 0.000000 0.842105 0.052632 0.105263 0.421053 0.157895 0.000000 0.421053 0.421053 0.052632 0.315789 0.210526 0.368421 0.631579 0.000000 0.000000 0.210526 0.631579 0.157895 0.000000 0.684211 0.000000 0.000000 0.315789 0.368421 0.368421 0.210526 0.052632 0.105263 0.105263 0.210526 0.578947 0.684211 0.000000 0.000000 0.315789 0.052632 0.947368 0.000000 0.000000 0.421053 0.368421 0.210526 0.000000 0.578947 0.105263 0.210526 0.105263 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 2 regular expression -------------------------------------------------------------------------------- [TC][TC]CA[CA][ACG][AC][CT][ACT]C[AT][AGT][CA][CA][AT][ACG][TG][AT]C[ACG][AG] -------------------------------------------------------------------------------- Time 6.65 secs. ******************************************************************************** ******************************************************************************** MOTIF 3 MEME width = 20 sites = 6 llr = 106 E-value = 5.0e+001 ******************************************************************************** -------------------------------------------------------------------------------- Motif 3 Description -------------------------------------------------------------------------------- Simplified A ::::::2:5:2:855:52a: pos.-specific C ::22:::::::522::32:2 probability G 825::a:5:a83:3:a25:8 matrix T 2838a:855::2::5::2:: bits 2.2 * * * 2.0 ** * * 1.8 ** * * * 1.5 * ** ** * ** Relative 1.3 ** **** ** * * ** Entropy 1.1 ** ***** ** * * ** (25.5 bits) 0.9 ** ******** * ** ** 0.7 ************* ** ** 0.4 ***************** ** 0.2 ******************** 0.0 -------------------- Multilevel GTGTTGTGAGGCAAAGAGAG consensus T TT G GT C sequence -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 3 sites sorted by position p-value -------------------------------------------------------------------------------- Sequence name Start P-value Site ------------- ----- --------- -------------------- 2703 330 4.55e-11 GTTGTGTTGC GTGTTGTGTGGCAAAGGGAG CACGAAGGAA 10034 103 4.55e-11 CCACGTCAAC GTGTTGTTTGGGAGAGCGAG GACCCCTTTA 6620 120 1.09e-08 AAACGTGGCA GTCTTGAGAGGCAGTGAGAC GTCGTCTCTT 6764 314 1.63e-08 ACATATTGGT GGTTTGTTTGGGCATGAAAG ATGGACGCAT 10831 345 1.84e-08 AATCAACCAA TTTCTGTGAGGCAATGCTAG CATCATCAGA 4190 258 1.96e-08 TCCAGGCACC GTGTTGTTAGATACAGACAG TACAAGTTTA -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 3 block diagrams -------------------------------------------------------------------------------- SEQUENCE NAME POSITION P-VALUE MOTIF DIAGRAM ------------- ---------------- ------------- 2703 4.6e-11 329_[+3]_151 10034 4.6e-11 102_[+3]_378 6620 1.1e-08 119_[+3]_361 6764 1.6e-08 313_[+3]_167 10831 1.8e-08 344_[+3]_136 4190 2e-08 257_[+3]_223 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 3 in BLOCKS format -------------------------------------------------------------------------------- BL MOTIF 3 width=20 seqs=6 2703 ( 330) GTGTTGTGTGGCAAAGGGAG 1 10034 ( 103) GTGTTGTTTGGGAGAGCGAG 1 6620 ( 120) GTCTTGAGAGGCAGTGAGAC 1 6764 ( 314) GGTTTGTTTGGGCATGAAAG 1 10831 ( 345) TTTCTGTGAGGCAATGCTAG 1 4190 ( 258) GTGTTGTTAGATACAGACAG 1 // -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 3 position-specific scoring matrix -------------------------------------------------------------------------------- log-odds matrix: alength= 4 w= 20 n= 9139 bayes= 11.672 E= 5.0e+001 -923 -923 194 -66 -923 -923 -38 166 -923 -55 120 33 -923 -55 -923 166 -923 -923 -923 192 -923 -923 220 -923 -72 -923 -923 166 -923 -923 120 92 86 -923 -923 92 -923 -923 220 -923 -72 -923 194 -923 -923 103 62 -66 160 -55 -923 -923 86 -55 62 -923 86 -923 -923 92 -923 -923 220 -923 86 45 -38 -923 -72 -55 120 -66 186 -923 -923 -923 -923 -55 194 -923 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 3 position-specific probability matrix -------------------------------------------------------------------------------- letter-probability matrix: alength= 4 w= 20 nsites= 6 E= 5.0e+001 0.000000 0.000000 0.833333 0.166667 0.000000 0.000000 0.166667 0.833333 0.000000 0.166667 0.500000 0.333333 0.000000 0.166667 0.000000 0.833333 0.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.166667 0.000000 0.000000 0.833333 0.000000 0.000000 0.500000 0.500000 0.500000 0.000000 0.000000 0.500000 0.000000 0.000000 1.000000 0.000000 0.166667 0.000000 0.833333 0.000000 0.000000 0.500000 0.333333 0.166667 0.833333 0.166667 0.000000 0.000000 0.500000 0.166667 0.333333 0.000000 0.500000 0.000000 0.000000 0.500000 0.000000 0.000000 1.000000 0.000000 0.500000 0.333333 0.166667 0.000000 0.166667 0.166667 0.500000 0.166667 1.000000 0.000000 0.000000 0.000000 0.000000 0.166667 0.833333 0.000000 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 3 regular expression -------------------------------------------------------------------------------- GT[GT]TTGT[GT][AT]GG[CG]A[AG][AT]G[AC]GAG -------------------------------------------------------------------------------- Time 9.76 secs. ******************************************************************************** ******************************************************************************** SUMMARY OF MOTIFS ******************************************************************************** -------------------------------------------------------------------------------- Combined block diagrams: non-overlapping sites with p-value < 0.0001 -------------------------------------------------------------------------------- SEQUENCE NAME COMBINED P-VALUE MOTIF DIAGRAM ------------- ---------------- ------------- 10034 1.13e-12 51_[+1(3.61e-08)]_8_[+2(1.02e-05)]_\ 1_[+3(4.55e-11)]_378 1005 1.23e-06 117_[+1(1.35e-07)]_29_\ [+1(3.82e-05)]_255_[+2(4.54e-07)]_36 10831 2.31e-06 157_[+2(6.42e-06)]_166_\ [+3(1.84e-08)]_136 11227 5.98e-07 [+2(1.26e-05)]_150_[+1(3.29e-08)]_\ 308 11795 2.24e-02 477_[+2(1.52e-06)]_2 12131 4.82e-05 97_[+1(8.50e-08)]_140_\ [+2(4.32e-05)]_221 2001 1.04e-07 109_[+2(6.45e-07)]_149_\ [+1(1.67e-08)]_200 23224 4.38e-05 [+1(1.08e-07)]_297_[+2(1.66e-05)]_\ 161 24409 1.34e-05 72_[+2(3.26e-06)]_111_\ [+1(1.56e-07)]_275 261704 7.20e-03 119_[+2(5.12e-07)]_360 263913 1.65e-02 259_[+2(2.72e-06)]_220 2703 7.43e-14 232_[+1(1.35e-08)]_76_\ [+3(4.55e-11)]_14_[+1(4.18e-06)]_93_[+2(1.52e-06)]_2 4190 1.96e-06 257_[+3(1.96e-08)]_137_\ [+2(3.56e-06)]_65 4814 1.90e-06 399_[+2(7.75e-11)]_80 5801 7.62e-02 341_[+2(8.77e-06)]_138 6620 4.65e-06 119_[+3(1.09e-08)]_240_\ [+2(1.10e-05)]_100 6764 1.82e-10 190_[+2(1.12e-06)]_80_\ [+1(2.07e-07)]_1_[+3(1.63e-08)]_167 9475 2.36e-03 300_[+2(7.23e-07)]_117_\ [+2(2.76e-05)]_41 bd887 9.47e-10 73_[+1(1.68e-10)]_246_\ [+2(1.23e-07)]_139 -------------------------------------------------------------------------------- ******************************************************************************** ******************************************************************************** Stopped because nmotifs = 3 reached. ******************************************************************************** CPU: seaotter.hsd1.wa.comcast.net ********************************************************************************