******************************************************************************** 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/47/47.seqs.fa ALPHABET= ACGT Sequence name Weight Length Sequence name Weight Length ------------- ------ ------ ------------- ------ ------ 42840 1.0000 500 42841 1.0000 500 46564 1.0000 500 47853 1.0000 500 22453 1.0000 500 48746 1.0000 500 43418 1.0000 500 43645 1.0000 500 49068 1.0000 500 55040 1.0000 500 49409 1.0000 500 49660 1.0000 500 43902 1.0000 500 43955 1.0000 500 44247 1.0000 500 34391 1.0000 500 44979 1.0000 500 11994 1.0000 500 35303 1.0000 500 45913 1.0000 500 46193 1.0000 500 42624 1.0000 500 42746 1.0000 500 44288 1.0000 500 32345 1.0000 500 46428 1.0000 500 45249 1.0000 500 36086 1.0000 500 44762 1.0000 500 44138 1.0000 500 49993 1.0000 500 45611 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/47/47.seqs.fa -oc motifs/47 -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= 32 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= 16000 N= 32 strands: + sample: seed= 0 seqfrac= 1 Letter frequencies in dataset: A 0.270 C 0.226 G 0.222 T 0.282 Background letter frequencies (from dataset with add-one prior applied): A 0.270 C 0.226 G 0.222 T 0.282 ******************************************************************************** ******************************************************************************** MOTIF 1 MEME width = 12 sites = 27 llr = 282 E-value = 1.9e-020 ******************************************************************************** -------------------------------------------------------------------------------- Motif 1 Description -------------------------------------------------------------------------------- Simplified A :::914:1:a3: pos.-specific C 3:7:9:::9::: probability G 3:2:::a:1:62 matrix T 4a11:6:9:::8 bits 2.2 2.0 * 1.7 * * 1.5 * ** ** * Relative 1.3 * ** **** Entropy 1.1 * ** **** * (15.1 bits) 0.9 **** **** * 0.7 *********** 0.4 ************ 0.2 ************ 0.0 ------------ Multilevel TTCACTGTCAGT consensus C A AG sequence G -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 sites sorted by position p-value -------------------------------------------------------------------------------- Sequence name Start P-value Site ------------- ----- --------- ------------ 45249 236 1.33e-07 TCTTCTTTCC TTCACTGTCAGT CCCTCCATTC 44288 22 1.33e-07 AAGGACGACG TTCACTGTCAGT ATCTTGGCAC 46564 328 1.33e-07 ATGAAAGGAC TTCACTGTCAGT AGAGCCAGGA 42841 228 1.33e-07 ACAGTCAGGT CTCACTGTCAGT AGAAGTAAAT 44979 383 3.19e-07 TAAGAGACTT TTCACAGTCAGT TGTACGGTAT 45611 419 5.37e-07 GATCGAAGCC TTCACTGTCAAT CTTACTGCAA 44762 486 1.11e-06 GTAACGCGCA GTCACAGTCAAT CCC 42624 373 1.11e-06 CCATTGCCAT TTGACTGTCAGT GAACATCGTA 47853 200 1.31e-06 ATGACGACGA GTCACAGTCAGG GACGTCGTCC 34391 232 1.57e-06 AAAATTTCAC CTCACTGTCAAG GTCATATAAA 42840 404 1.78e-06 TTCCCTGCCA GTGACAGTCAGT CCAAAACACT 43645 151 2.36e-06 CGCTAATGTT GTCAATGTCAGT TACAAATTTC 48746 294 2.54e-06 AGGAATCGCG GTCACTGTGAGT CGAGAGAAGC 45913 1 2.84e-06 . CTCAAAGTCAGT CTCACTCATT 11994 186 3.97e-06 TGAAATAGCT CTCTCTGTCAGT TGCTGTGGTT 35303 179 5.26e-06 CTCCCAATAC TTTACTGTCAAT GTAAATTGTC 44138 85 6.29e-06 TTGTTTACAG CTCACTGTCATT TCGAAACGAT 49660 187 1.03e-05 TGCATGCCTC TTCACAGTCGGT GCATTTTCGT 32345 346 1.81e-05 TGCGACGTTT CTTACAGTCAAG TTGTTTTCGA 49993 177 1.96e-05 GTGATATGTT TTGACTGTGAGG ATGCTGAGCC 49068 266 3.42e-05 TTTTGGCTGT CACAATGTCAGT GTTCGACTCA 43902 267 4.85e-05 AGGTTTGGCG TTCACAGAGAAT TGCATGACGA 43418 175 4.85e-05 GACTGATTGT CTGACTGACAAG TAAATCACAA 55040 261 5.52e-05 TTCAACACTT CTTACAGTTAGT CAACGCAGCA 43955 285 5.81e-05 AGTTCCCTTT TTCTCGGTCAAT GATGTCGCTA 46428 411 6.19e-05 GATTTTCTGG GTCAATGTCACG ACACGAATAG 36086 133 6.87e-05 CGTATGTTGT GTGACATTCAAT GACAGCGCGG -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 block diagrams -------------------------------------------------------------------------------- SEQUENCE NAME POSITION P-VALUE MOTIF DIAGRAM ------------- ---------------- ------------- 45249 1.3e-07 235_[+1]_253 44288 1.3e-07 21_[+1]_467 46564 1.3e-07 327_[+1]_161 42841 1.3e-07 227_[+1]_261 44979 3.2e-07 382_[+1]_106 45611 5.4e-07 418_[+1]_70 44762 1.1e-06 485_[+1]_3 42624 1.1e-06 372_[+1]_116 47853 1.3e-06 199_[+1]_289 34391 1.6e-06 231_[+1]_257 42840 1.8e-06 403_[+1]_85 43645 2.4e-06 150_[+1]_338 48746 2.5e-06 293_[+1]_195 45913 2.8e-06 [+1]_488 11994 4e-06 185_[+1]_303 35303 5.3e-06 178_[+1]_310 44138 6.3e-06 84_[+1]_404 49660 1e-05 186_[+1]_302 32345 1.8e-05 345_[+1]_143 49993 2e-05 176_[+1]_312 49068 3.4e-05 265_[+1]_223 43902 4.8e-05 266_[+1]_222 43418 4.8e-05 174_[+1]_314 55040 5.5e-05 260_[+1]_228 43955 5.8e-05 284_[+1]_204 46428 6.2e-05 410_[+1]_78 36086 6.9e-05 132_[+1]_356 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 in BLOCKS format -------------------------------------------------------------------------------- BL MOTIF 1 width=12 seqs=27 45249 ( 236) TTCACTGTCAGT 1 44288 ( 22) TTCACTGTCAGT 1 46564 ( 328) TTCACTGTCAGT 1 42841 ( 228) CTCACTGTCAGT 1 44979 ( 383) TTCACAGTCAGT 1 45611 ( 419) TTCACTGTCAAT 1 44762 ( 486) GTCACAGTCAAT 1 42624 ( 373) TTGACTGTCAGT 1 47853 ( 200) GTCACAGTCAGG 1 34391 ( 232) CTCACTGTCAAG 1 42840 ( 404) GTGACAGTCAGT 1 43645 ( 151) GTCAATGTCAGT 1 48746 ( 294) GTCACTGTGAGT 1 45913 ( 1) CTCAAAGTCAGT 1 11994 ( 186) CTCTCTGTCAGT 1 35303 ( 179) TTTACTGTCAAT 1 44138 ( 85) CTCACTGTCATT 1 49660 ( 187) TTCACAGTCGGT 1 32345 ( 346) CTTACAGTCAAG 1 49993 ( 177) TTGACTGTGAGG 1 49068 ( 266) CACAATGTCAGT 1 43902 ( 267) TTCACAGAGAAT 1 43418 ( 175) CTGACTGACAAG 1 55040 ( 261) CTTACAGTTAGT 1 43955 ( 285) TTCTCGGTCAAT 1 46428 ( 411) GTCAATGTCACG 1 36086 ( 133) GTGACATTCAAT 1 // -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 position-specific scoring matrix -------------------------------------------------------------------------------- log-odds matrix: alength= 4 w= 12 n= 15648 bayes= 10.8286 E= 1.9e-020 -1140 56 22 53 -286 -1140 -1140 177 -1140 164 -26 -134 178 -1140 -1140 -193 -86 191 -1140 -1140 46 -1140 -258 107 -1140 -1140 212 -292 -186 -1140 -1140 172 -1140 191 -100 -292 183 -1140 -258 -1140 30 -261 142 -292 -1140 -1140 0 147 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 position-specific probability matrix -------------------------------------------------------------------------------- letter-probability matrix: alength= 4 w= 12 nsites= 27 E= 1.9e-020 0.000000 0.333333 0.259259 0.407407 0.037037 0.000000 0.000000 0.962963 0.000000 0.703704 0.185185 0.111111 0.925926 0.000000 0.000000 0.074074 0.148148 0.851852 0.000000 0.000000 0.370370 0.000000 0.037037 0.592593 0.000000 0.000000 0.962963 0.037037 0.074074 0.000000 0.000000 0.925926 0.000000 0.851852 0.111111 0.037037 0.962963 0.000000 0.037037 0.000000 0.333333 0.037037 0.592593 0.037037 0.000000 0.000000 0.222222 0.777778 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 regular expression -------------------------------------------------------------------------------- [TCG]TCAC[TA]GTCA[GA][TG] -------------------------------------------------------------------------------- Time 9.15 secs. ******************************************************************************** ******************************************************************************** MOTIF 2 MEME width = 16 sites = 10 llr = 135 E-value = 7.0e+000 ******************************************************************************** -------------------------------------------------------------------------------- Motif 2 Description -------------------------------------------------------------------------------- Simplified A 954:::28911aa38: pos.-specific C 112::a3::3::::24 probability G :4171:52:19::7:6 matrix T ::339:::15:::::: bits 2.2 * 2.0 * ** 1.7 * *** 1.5 * * * *** Relative 1.3 * ** ** ***** Entropy 1.1 * *** ** ****** (19.5 bits) 0.9 * *** ** ****** 0.7 ** ****** ****** 0.4 ** ****** ****** 0.2 ** ************* 0.0 ---------------- Multilevel AAAGTCGAATGAAGAG consensus GTT CG C ACC sequence C A -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 2 sites sorted by position p-value -------------------------------------------------------------------------------- Sequence name Start P-value Site ------------- ----- --------- ---------------- 43955 170 6.00e-09 AGAGCAAACC AGTGTCCAATGAAGAG GATCAGCAAC 43418 84 1.70e-08 CCTGAAGATA ACAGTCGAATGAAGAG TACTACTCTA 44288 106 2.86e-08 CTGACGACAA AAAGTCGAACGAAAAC CCTCCAGTGG 42841 396 2.86e-08 TGCTCAGCTC AATGTCAAACGAAGAG ACACTGAGAT 45913 396 3.15e-07 CCAAACATCA AGATTCGGATGAAGCG ATTCTAGCTG 42840 167 6.08e-07 TGACGGCGCA AGCTTCAAATGAAAAC CAGTCAGTAC 34391 117 7.45e-07 GATGCAACGT AAAGGCGAAAGAAAAG AGAATCTATT 43902 223 1.11e-06 CAACGACCGA CGCGTCCGACGAAGAC ACAAGGATTC 35303 368 1.26e-06 GTCCAAGGGC AAGGTCCATTGAAGCG CCCGTACTTC 47853 479 1.68e-06 TTGTAAAGTG AATTTCGAAGAAAGAC TGTAAC -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 2 block diagrams -------------------------------------------------------------------------------- SEQUENCE NAME POSITION P-VALUE MOTIF DIAGRAM ------------- ---------------- ------------- 43955 6e-09 169_[+2]_315 43418 1.7e-08 83_[+2]_401 44288 2.9e-08 105_[+2]_379 42841 2.9e-08 395_[+2]_89 45913 3.2e-07 395_[+2]_89 42840 6.1e-07 166_[+2]_318 34391 7.4e-07 116_[+2]_368 43902 1.1e-06 222_[+2]_262 35303 1.3e-06 367_[+2]_117 47853 1.7e-06 478_[+2]_6 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 2 in BLOCKS format -------------------------------------------------------------------------------- BL MOTIF 2 width=16 seqs=10 43955 ( 170) AGTGTCCAATGAAGAG 1 43418 ( 84) ACAGTCGAATGAAGAG 1 44288 ( 106) AAAGTCGAACGAAAAC 1 42841 ( 396) AATGTCAAACGAAGAG 1 45913 ( 396) AGATTCGGATGAAGCG 1 42840 ( 167) AGCTTCAAATGAAAAC 1 34391 ( 117) AAAGGCGAAAGAAAAG 1 43902 ( 223) CGCGTCCGACGAAGAC 1 35303 ( 368) AAGGTCCATTGAAGCG 1 47853 ( 479) AATTTCGAAGAAAGAC 1 // -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 2 position-specific scoring matrix -------------------------------------------------------------------------------- log-odds matrix: alength= 4 w= 16 n= 15520 bayes= 11.5428 E= 7.0e+000 174 -118 -997 -997 89 -118 85 -997 57 -18 -115 9 -997 -997 166 9 -997 -997 -115 168 -997 214 -997 -997 -43 41 117 -997 157 -997 -15 -997 174 -997 -997 -149 -143 41 -115 83 -143 -997 202 -997 189 -997 -997 -997 189 -997 -997 -997 15 -997 166 -997 157 -18 -997 -997 -997 82 143 -997 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 2 position-specific probability matrix -------------------------------------------------------------------------------- letter-probability matrix: alength= 4 w= 16 nsites= 10 E= 7.0e+000 0.900000 0.100000 0.000000 0.000000 0.500000 0.100000 0.400000 0.000000 0.400000 0.200000 0.100000 0.300000 0.000000 0.000000 0.700000 0.300000 0.000000 0.000000 0.100000 0.900000 0.000000 1.000000 0.000000 0.000000 0.200000 0.300000 0.500000 0.000000 0.800000 0.000000 0.200000 0.000000 0.900000 0.000000 0.000000 0.100000 0.100000 0.300000 0.100000 0.500000 0.100000 0.000000 0.900000 0.000000 1.000000 0.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.000000 0.300000 0.000000 0.700000 0.000000 0.800000 0.200000 0.000000 0.000000 0.000000 0.400000 0.600000 0.000000 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 2 regular expression -------------------------------------------------------------------------------- A[AG][ATC][GT]TC[GCA][AG]A[TC]GAA[GA][AC][GC] -------------------------------------------------------------------------------- Time 18.44 secs. ******************************************************************************** ******************************************************************************** MOTIF 3 MEME width = 19 sites = 4 llr = 81 E-value = 1.2e+003 ******************************************************************************** -------------------------------------------------------------------------------- Motif 3 Description -------------------------------------------------------------------------------- Simplified A a:58:8:3:5a5::::::: pos.-specific C :a::3:::83::::::::: probability G ::5383:333::aaa:5aa matrix T ::::::a5:::5:::a5:: bits 2.2 * *** ** 2.0 ** * *** ** 1.7 ** * * **** ** 1.5 ** * * **** ** Relative 1.3 ** * * * * **** ** Entropy 1.1 ******* * * ******* (29.3 bits) 0.9 ******* * ********* 0.7 ******* * ********* 0.4 ******************* 0.2 ******************* 0.0 ------------------- Multilevel ACAAGATTCAAAGGGTGGG consensus GGCG AGC T T sequence G G -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 3 sites sorted by position p-value -------------------------------------------------------------------------------- Sequence name Start P-value Site ------------- ----- --------- ------------------- 44138 469 1.32e-10 AATCCCTGCG ACGAGATTGAAAGGGTGGG AAGGGAAAGA 47853 132 3.75e-10 CCAACTTTGC ACAAGATACGATGGGTTGG TCTGTGATCG 22453 252 7.11e-10 CAGAGGGTGA ACAAGGTGCCATGGGTGGG CGGGCTGGGC 49660 334 8.52e-10 CCATGAGAAG ACGGCATTCAAAGGGTTGG TAACGAGATC -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 3 block diagrams -------------------------------------------------------------------------------- SEQUENCE NAME POSITION P-VALUE MOTIF DIAGRAM ------------- ---------------- ------------- 44138 1.3e-10 468_[+3]_13 47853 3.7e-10 131_[+3]_350 22453 7.1e-10 251_[+3]_230 49660 8.5e-10 333_[+3]_148 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 3 in BLOCKS format -------------------------------------------------------------------------------- BL MOTIF 3 width=19 seqs=4 44138 ( 469) ACGAGATTGAAAGGGTGGG 1 47853 ( 132) ACAAGATACGATGGGTTGG 1 22453 ( 252) ACAAGGTGCCATGGGTGGG 1 49660 ( 334) ACGGCATTCAAAGGGTTGG 1 // -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 3 position-specific scoring matrix -------------------------------------------------------------------------------- log-odds matrix: alength= 4 w= 19 n= 15424 bayes= 11.9125 E= 1.2e+003 189 -865 -865 -865 -865 214 -865 -865 89 -865 117 -865 147 -865 17 -865 -865 14 175 -865 147 -865 17 -865 -865 -865 -865 183 -11 -865 17 83 -865 173 17 -865 89 14 17 -865 189 -865 -865 -865 89 -865 -865 83 -865 -865 217 -865 -865 -865 217 -865 -865 -865 217 -865 -865 -865 -865 183 -865 -865 117 83 -865 -865 217 -865 -865 -865 217 -865 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 3 position-specific probability matrix -------------------------------------------------------------------------------- letter-probability matrix: alength= 4 w= 19 nsites= 4 E= 1.2e+003 1.000000 0.000000 0.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.500000 0.000000 0.500000 0.000000 0.750000 0.000000 0.250000 0.000000 0.000000 0.250000 0.750000 0.000000 0.750000 0.000000 0.250000 0.000000 0.000000 0.000000 0.000000 1.000000 0.250000 0.000000 0.250000 0.500000 0.000000 0.750000 0.250000 0.000000 0.500000 0.250000 0.250000 0.000000 1.000000 0.000000 0.000000 0.000000 0.500000 0.000000 0.000000 0.500000 0.000000 0.000000 1.000000 0.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.500000 0.500000 0.000000 0.000000 1.000000 0.000000 0.000000 0.000000 1.000000 0.000000 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 3 regular expression -------------------------------------------------------------------------------- AC[AG][AG][GC][AG]T[TAG][CG][ACG]A[AT]GGGT[GT]GG -------------------------------------------------------------------------------- Time 28.30 secs. ******************************************************************************** ******************************************************************************** SUMMARY OF MOTIFS ******************************************************************************** -------------------------------------------------------------------------------- Combined block diagrams: non-overlapping sites with p-value < 0.0001 -------------------------------------------------------------------------------- SEQUENCE NAME COMBINED P-VALUE MOTIF DIAGRAM ------------- ---------------- ------------- 42840 2.26e-05 166_[+2(6.08e-07)]_221_\ [+1(1.78e-06)]_85 42841 4.35e-08 214_[+1(1.21e-06)]_1_[+1(1.33e-07)]_\ 156_[+2(2.86e-08)]_89 46564 4.31e-04 327_[+1(1.33e-07)]_161 47853 4.55e-11 131_[+3(3.75e-10)]_49_\ [+1(1.31e-06)]_267_[+2(1.68e-06)]_6 22453 2.00e-05 251_[+3(7.11e-10)]_230 48746 2.25e-03 293_[+1(2.54e-06)]_22_\ [+1(2.64e-05)]_161 43418 1.44e-05 83_[+2(1.70e-08)]_75_[+1(4.85e-05)]_\ 314 43645 1.38e-02 95_[+1(2.13e-05)]_43_[+1(2.36e-06)]_\ 338 49068 6.08e-02 265_[+1(3.42e-05)]_223 55040 2.92e-03 78_[+3(6.32e-05)]_163_\ [+1(5.52e-05)]_228 49409 4.38e-01 500 49660 2.65e-07 186_[+1(1.03e-05)]_135_\ [+3(8.52e-10)]_148 43902 2.98e-04 222_[+2(1.11e-06)]_28_\ [+1(4.85e-05)]_222 43955 1.07e-05 169_[+2(6.00e-09)]_99_\ [+1(5.81e-05)]_204 44247 3.57e-01 500 34391 3.04e-05 116_[+2(7.45e-07)]_99_\ [+1(1.57e-06)]_60_[+1(1.81e-05)]_4_[+1(4.85e-05)]_169 44979 5.90e-04 382_[+1(3.19e-07)]_106 11994 1.14e-02 185_[+1(3.97e-06)]_303 35303 9.22e-05 178_[+1(5.26e-06)]_177_\ [+2(1.26e-06)]_117 45913 1.89e-05 [+1(2.84e-06)]_383_[+2(3.15e-07)]_\ 89 46193 2.72e-01 500 42624 7.18e-03 50_[+1(8.17e-05)]_310_\ [+1(1.11e-06)]_116 42746 7.55e-01 500 44288 7.57e-09 21_[+1(1.33e-07)]_72_[+2(2.86e-08)]_\ 228_[+3(5.41e-05)]_132 32345 1.35e-01 345_[+1(1.81e-05)]_143 46428 2.72e-01 410_[+1(6.19e-05)]_78 45249 7.89e-04 235_[+1(1.33e-07)]_253 36086 3.20e-02 132_[+1(6.87e-05)]_356 44762 2.70e-03 485_[+1(1.11e-06)]_3 44138 7.04e-09 84_[+1(6.29e-06)]_372_\ [+3(1.32e-10)]_13 49993 4.05e-02 176_[+1(1.96e-05)]_312 45611 7.04e-04 418_[+1(5.37e-07)]_70 -------------------------------------------------------------------------------- ******************************************************************************** ******************************************************************************** Stopped because nmotifs = 3 reached. ******************************************************************************** CPU: seaotter.hsd1.wa.comcast.net ********************************************************************************