******************************************************************************** MEME - Motif discovery tool ******************************************************************************** MEME version 4.3.0 (Release date: Sat Sep 26 01:51:56 PDT 2009) 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= seqs.fasta ALPHABET= ACDEFGHIKLMNPQRSTVWY Sequence name Weight Length Sequence name Weight Length ------------- ------ ------ ------------- ------ ------ 1 1.0000 35 2 1.0000 35 3 1.0000 35 4 1.0000 35 5 1.0000 35 6 1.0000 35 7 1.0000 35 8 1.0000 35 9 1.0000 35 10 1.0000 35 11 1.0000 35 12 1.0000 35 13 1.0000 35 14 1.0000 35 15 1.0000 35 16 1.0000 35 17 1.0000 35 18 1.0000 35 19 1.0000 35 20 1.0000 35 21 1.0000 35 22 1.0000 35 23 1.0000 35 24 1.0000 35 25 1.0000 35 26 1.0000 35 27 1.0000 35 28 1.0000 35 29 1.0000 35 30 1.0000 35 31 1.0000 35 32 1.0000 35 33 1.0000 35 34 1.0000 35 35 1.0000 35 36 1.0000 35 37 1.0000 35 38 1.0000 35 39 1.0000 35 40 1.0000 35 41 1.0000 35 42 1.0000 35 43 1.0000 35 44 1.0000 35 45 1.0000 35 46 1.0000 35 47 1.0000 35 48 1.0000 35 49 1.0000 35 50 1.0000 35 51 1.0000 35 52 1.0000 35 53 1.0000 35 54 1.0000 35 55 1.0000 35 56 1.0000 35 57 1.0000 35 58 1.0000 35 59 1.0000 35 60 1.0000 35 61 1.0000 35 62 1.0000 35 63 1.0000 35 64 1.0000 35 65 1.0000 35 66 1.0000 35 67 1.0000 35 68 1.0000 35 69 1.0000 35 70 1.0000 35 71 1.0000 35 72 1.0000 35 73 1.0000 35 74 1.0000 35 75 1.0000 35 76 1.0000 35 77 1.0000 35 78 1.0000 35 79 1.0000 35 80 1.0000 35 81 1.0000 35 82 1.0000 35 83 1.0000 35 84 1.0000 35 85 1.0000 35 86 1.0000 35 87 1.0000 35 88 1.0000 35 89 1.0000 35 90 1.0000 35 91 1.0000 35 92 1.0000 35 93 1.0000 35 94 1.0000 35 95 1.0000 35 96 1.0000 35 97 1.0000 35 98 1.0000 35 99 1.0000 35 100 1.0000 35 101 1.0000 35 102 1.0000 35 103 1.0000 35 104 1.0000 35 105 1.0000 35 106 1.0000 35 107 1.0000 35 108 1.0000 35 109 1.0000 35 110 1.0000 35 111 1.0000 35 112 1.0000 35 113 1.0000 35 114 1.0000 35 115 1.0000 35 116 1.0000 35 117 1.0000 35 118 1.0000 35 119 1.0000 35 120 1.0000 35 121 1.0000 35 122 1.0000 35 123 1.0000 35 124 1.0000 35 125 1.0000 35 126 1.0000 35 127 1.0000 35 128 1.0000 35 129 1.0000 35 130 1.0000 35 131 1.0000 35 132 1.0000 35 133 1.0000 35 134 1.0000 35 135 1.0000 35 136 1.0000 35 137 1.0000 35 138 1.0000 35 139 1.0000 35 140 1.0000 35 141 1.0000 35 142 1.0000 35 143 1.0000 35 144 1.0000 35 145 1.0000 35 146 1.0000 35 147 1.0000 35 148 1.0000 35 149 1.0000 35 150 1.0000 35 151 1.0000 35 152 1.0000 35 153 1.0000 35 154 1.0000 35 155 1.0000 35 156 1.0000 35 157 1.0000 35 158 1.0000 35 159 1.0000 35 160 1.0000 35 161 1.0000 35 162 1.0000 35 163 1.0000 35 164 1.0000 35 165 1.0000 35 166 1.0000 35 167 1.0000 35 168 1.0000 35 169 1.0000 35 170 1.0000 35 171 1.0000 35 172 1.0000 35 173 1.0000 35 174 1.0000 35 175 1.0000 35 176 1.0000 35 177 1.0000 35 178 1.0000 35 179 1.0000 35 180 1.0000 35 181 1.0000 35 182 1.0000 35 183 1.0000 35 184 1.0000 35 185 1.0000 35 186 1.0000 35 187 1.0000 35 188 1.0000 35 189 1.0000 35 190 1.0000 35 191 1.0000 35 192 1.0000 35 193 1.0000 35 194 1.0000 35 195 1.0000 35 196 1.0000 35 197 1.0000 35 198 1.0000 35 199 1.0000 35 200 1.0000 35 201 1.0000 35 202 1.0000 35 203 1.0000 35 204 1.0000 35 205 1.0000 35 206 1.0000 35 207 1.0000 35 208 1.0000 35 209 1.0000 35 210 1.0000 35 211 1.0000 35 212 1.0000 35 213 1.0000 35 214 1.0000 35 215 1.0000 35 216 1.0000 35 217 1.0000 35 218 1.0000 35 219 1.0000 35 220 1.0000 35 221 1.0000 35 222 1.0000 35 223 1.0000 35 224 1.0000 35 225 1.0000 35 226 1.0000 35 227 1.0000 35 228 1.0000 35 229 1.0000 35 230 1.0000 35 231 1.0000 35 232 1.0000 35 233 1.0000 35 234 1.0000 35 235 1.0000 35 236 1.0000 35 237 1.0000 35 238 1.0000 35 239 1.0000 35 240 1.0000 35 241 1.0000 35 242 1.0000 35 243 1.0000 35 244 1.0000 35 245 1.0000 35 246 1.0000 35 247 1.0000 35 248 1.0000 35 249 1.0000 35 250 1.0000 35 251 1.0000 35 252 1.0000 35 253 1.0000 35 254 1.0000 35 255 1.0000 35 256 1.0000 35 257 1.0000 35 258 1.0000 35 259 1.0000 35 260 1.0000 35 261 1.0000 35 262 1.0000 35 263 1.0000 35 264 1.0000 35 265 1.0000 35 266 1.0000 35 267 1.0000 35 268 1.0000 35 269 1.0000 35 270 1.0000 35 271 1.0000 35 272 1.0000 35 273 1.0000 35 274 1.0000 35 275 1.0000 35 276 1.0000 35 277 1.0000 35 278 1.0000 35 279 1.0000 35 280 1.0000 35 281 1.0000 35 282 1.0000 35 283 1.0000 35 284 1.0000 35 285 1.0000 35 286 1.0000 35 287 1.0000 35 288 1.0000 35 289 1.0000 35 290 1.0000 35 291 1.0000 35 292 1.0000 35 293 1.0000 35 294 1.0000 35 295 1.0000 35 296 1.0000 35 297 1.0000 35 298 1.0000 35 299 1.0000 35 300 1.0000 35 ******************************************************************************** ******************************************************************************** COMMAND LINE SUMMARY ******************************************************************************** This information can also be useful in the event you wish to report a problem with the MEME software. command: meme -plib /etc/meme/prior30.plib seqs.fasta model: mod= zoops nmotifs= 1 evt= inf object function= E-value of product of p-values width: minw= 8 maxw= 35 minic= 0.00 width: wg= 11 ws= 1 endgaps= yes nsites: minsites= 2 maxsites= 300 wnsites= 0.8 theta: prob= 1 spmap= pam spfuzz= 120 global: substring= yes branching= no wbranch= no em: prior= megap b= 52500 maxiter= 50 distance= 1e-05 data: n= 10500 N= 300 sample: seed= 0 seqfrac= 1 Dirichlet mixture priors file: prior30.plib Letter frequencies in dataset: A 0.423 C 0.060 D 0.000 E 0.000 F 0.000 G 0.193 H 0.000 I 0.000 K 0.000 L 0.000 M 0.000 N 0.000 P 0.000 Q 0.000 R 0.000 S 0.000 T 0.323 V 0.000 W 0.000 Y 0.000 Background letter frequencies (from dataset with add-one prior applied): A 0.423 C 0.060 D 0.000 E 0.000 F 0.000 G 0.193 H 0.000 I 0.000 K 0.000 L 0.000 M 0.000 N 0.000 P 0.000 Q 0.000 R 0.000 S 0.000 T 0.323 V 0.000 W 0.000 Y 0.000 ******************************************************************************** ******************************************************************************** MOTIF 1 width = 8 sites = 106 llr = 765 E-value = 2.6e-025 ******************************************************************************** -------------------------------------------------------------------------------- Motif 1 Description -------------------------------------------------------------------------------- Simplified A 7::7::66 pos.-specific C ::::::1: probability D :::::::: matrix E :::::::: F :::::::: G 2aa2aa12 H :::::::: I :::::::: K :::::::: L :::::::: M :::::::: N :::::::: P :::::::: Q :::::::: R :::::::: S :::::::: T 1:::::22 V :::::::: W :::::::: Y :::::::: bits 13.4 12.0 10.7 9.4 Relative 8.0 Entropy 6.7 (10.4 bits) 5.3 4.0 2.7 ** ** 1.3 ** ** 0.0 -------- Multilevel AGGAGGAA consensus G TT sequence -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 sites sorted by position p-value -------------------------------------------------------------------------------- Sequence name Start P-value Site ------------- ----- --------- -------- 293 20 4.59e-05 ATAACTAAAG AGGAGGAA AACTACAT 254 21 4.59e-05 AATTGCTAAA AGGAGGAA ATGCTTT 253 23 4.59e-05 AGAATAATTT AGGAGGAA CAAAA 243 26 4.59e-05 CCTAATGATA AGGAGGAA AT 196 23 4.59e-05 AATTTAATTT AGGAGGAA TAAAG 171 18 4.59e-05 GTAATATAAT AGGAGGAA TTAGCAGGAT 142 23 4.59e-05 AACCATAGAT AGGAGGAA AAATT 138 23 4.59e-05 AGTAGATGCA AGGAGGAA GAGAA 137 23 4.59e-05 ATAGACAAAT AGGAGGAA GCAAC 136 23 4.59e-05 GTTGAGTCCT AGGAGGAA TAAAA 127 22 4.59e-05 TTTAAATTTT AGGAGGAA CTAATT 77 22 4.59e-05 TTGAAGGGCT AGGAGGAA CTTACA 73 22 4.59e-05 ATGGTTATTT AGGAGGAA TTTCAA 43 23 4.59e-05 CTATAGCGCA AGGAGGAA TGAAA 42 22 4.59e-05 AATAAAAAGC AGGAGGAA ACAAAG 7 21 4.59e-05 ATGACATTAA AGGAGGAA CAGGCCT 215 24 5.29e-05 ATTGTAAAAT AGGAGGCA TTAT 197 25 5.29e-05 AATTTAGATT AGGAGGCA CTA 277 19 9.37e-05 AATGTATTAA AGGGGGAA AAATTAAAA 201 23 9.37e-05 TAAAAATTAA AGGGGGAA TTTTC 174 23 9.37e-05 TATATAAAAA AGGGGGAA TTATT 129 23 9.37e-05 ATATTATAAG AGGGGGAA ATGTT 101 22 9.37e-05 TTTTCAATAT AGGAGGAG AAGAAC 6 22 9.37e-05 TTATAATATG AGGAGGAG TTTTAA 298 23 1.48e-04 TATGGTCTAA AGGAGGAT AAAAT 264 25 1.48e-04 CTTAATAAAA AGGAGGAT AAA 185 24 1.48e-04 GCATAAAGAT AGGAGGAT TAGT 172 22 1.48e-04 TAATTTTTAA AGGAGGAT ACTATT 162 24 1.48e-04 TGTAAAAGTT AGGAGGAT AAAC 128 22 1.48e-04 AACTAAATAT AGGAGGGA ATACCT 116 25 1.48e-04 TTATAAAGAT AGGAGGAT ATT 99 22 1.48e-04 CAAAAAAATT AGGAGGAT TAAAAG 97 19 1.48e-04 TAATTTGTAG AGGAGGGA GAATAGTTT 83 23 1.48e-04 AGTTAATTTT AGGAGGAT GAAAG 80 24 1.48e-04 TATTTTTAAA AGGAGGAT AAAA 70 24 1.48e-04 CTTTAATTAT AGGAGGGA TTTT 36 24 1.48e-04 GCATTACATT AGGAGGAT ACAC 272 24 1.82e-04 ATATCTCTGA AGGAGGTA AAAT 118 23 1.82e-04 ATATCCTTTA AGGAGGTA ATAAT 112 23 1.82e-04 TTAATAAATC AGGAGGTA TTATT 106 23 1.82e-04 TTAAAAAATA AGGAGGTA AAGTT 95 25 1.82e-04 AATAATAGCC AGGAGGTA AAT 72 23 1.82e-04 ATTAAAATTT AGGAGGTA TTTTT 294 4 2.09e-04 AAT GGGAGGAA ATTAAGAAGG 258 23 2.09e-04 TTAAATTTTA GGGAGGAA TTTAT 247 20 2.09e-04 AAAGGGGAAT GGGAGGAA GATAAATA 86 25 2.09e-04 TTCAAAACTA GGGAGGAA TAA 79 21 2.09e-04 ATATACTAGA GGGAGGAA ACTTAAT 74 21 2.09e-04 TTATTGAAAA GGGAGGAA TTTTATT 65 23 2.09e-04 ATAAGTATTT GGGAGGAA TTTGT 5 23 2.09e-04 CGTTTTTTAT GGGAGGAA TTTTA 232 22 2.23e-04 GATTTATTTT AGGAGGCT AAAAAT 195 22 2.23e-04 CATCTATGCT AGGGGGAG AAAATA 28 22 2.23e-04 GACTTAAAGA AGGGGGAG TATAGT 261 22 2.57e-04 ATAACTAATA AGGGGGAT TAATAA 251 20 2.57e-04 GATTTATAAG AGGGGGGA AGTATGGA 246 22 2.57e-04 TTAAAAATAA AGGGGGAT TTTATT 191 22 2.57e-04 ATTATATCAT AGGAGGGG AATTTA 189 22 2.57e-04 GTATTGTATA AGGGGGAT ATAGTA 153 24 2.57e-04 CCAATTCATA AGGGGGAT TTAT 98 24 2.57e-04 ACTATTAAAG AGGGGGAT ATTA 61 22 2.57e-04 ATCACCTATA AGGGGGGA GCTTTG 45 21 2.57e-04 TATTAATTGA AGGGGGAT TAATACA 27 21 2.57e-04 TTAAAAGGGT AGGGGGAT GGGTAAC 22 23 2.57e-04 TCGACATAAT AGGGGGGA TTGAC 18 22 2.57e-04 TAACTATATT AGGGGGGA GTATTA 2 24 2.57e-04 TATTAATCCT AGGGGGAT TTCA 299 22 3.07e-04 TGAAATATTT AGGGGGTA TACAAT 178 23 3.07e-04 ACATTCTTTA AGGGGGTA TAGAT 168 23 3.07e-04 ATATTGAAAG AGGAGGGT ATGAA 125 24 3.07e-04 ATATTTTATG AGGGGGTA AAAA 59 24 3.07e-04 ATAATTAAAA AGGAGGGT TTTA 41 21 3.07e-04 ATAATAAGGA AGGAGGTG ATAATCC 34 21 3.07e-04 GCTCAAAGGT AGGAGGTG TACTATA 280 24 3.59e-04 ATACATGCGG AGGCGGAA AATA 102 19 3.59e-04 ACAAAATAAA GGGAGGAG AATACAAAG 24 19 3.59e-04 TATACCAAAT GGGAGGAG TGAATTTCT 300 21 4.27e-04 CTGAACTATT TGGAGGAA ATTAGTT 230 22 4.46e-04 ATATTCATTT GGGAGGTA TATGTT 38 25 4.64e-04 TTATAGCTTA AGGGGGTG GTA 203 20 4.92e-04 TAAGGCCGAA GGGAGGCT TCCTATAT 57 23 4.92e-04 TATTGGGATT TGGAGGCA ATTCT 51 25 4.92e-04 AATAAGTATT TGGAGGCA AAT 289 24 5.49e-04 TAACTTTGTT TGGAGGAG TTAT 235 21 5.49e-04 ATAAATTTAA TGGAGGAG TACAAGA 226 24 5.49e-04 TTATATTTTT AGGCGGAT TAAT 10 20 5.49e-04 TATGTGAAAT GGGAGGGG ATAATAAT 285 23 6.20e-04 GAATATAATT TGGAGGGA TAATT 218 23 6.20e-04 ATTTATTTTA GGGAGGTG TATAA 176 23 6.20e-04 AAAATCCATC GGGAGGTG TTTTT 164 24 6.20e-04 TAGATTAAGA TGGAGGAT AAAA 117 22 6.20e-04 AACTAAGTAT GGGAGGGT ATGATT 268 22 6.66e-04 CTATTTTAAA TGGAGGTA ACATAA 257 21 6.66e-04 TAGAATCGTA TGGAGGTA TAATTAA 208 25 6.66e-04 TAATGAGATT TGGAGGTA AAA 152 22 6.66e-04 ACCTATTAAA AGGAGGTC TAGATA 100 24 6.66e-04 TACATATTTA TGGAGGTA AAAA 92 24 6.66e-04 ATATATATAA AGGAGGTC AGAA 20 23 6.66e-04 ATAGTTGTTG TGGAGGTA GTAGT 4 23 6.66e-04 TATAAAATTA TGGAGGTA ATAGT 156 21 6.97e-04 AAATTAAAAT CGGAGGAA ATTTAAT 209 25 7.36e-04 TAGAGTAAAA TGGAGGGG ATT 139 23 7.36e-04 TTATTCCGTT AGGCGGTG ATTTC 221 26 8.84e-04 AATATAAAAG AGGTGGAA AC 68 25 9.66e-04 TATATTTGTG AGGTGGAG TAA 222 25 1.02e-03 TATTAGGAAT AGGTGGGA AAT -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 block diagrams -------------------------------------------------------------------------------- SEQUENCE NAME POSITION P-VALUE MOTIF DIAGRAM ------------- ---------------- ------------- 293 4.6e-05 19_[1]_8 254 4.6e-05 20_[1]_7 253 4.6e-05 22_[1]_5 243 4.6e-05 25_[1]_2 196 4.6e-05 22_[1]_5 171 4.6e-05 17_[1]_10 142 4.6e-05 22_[1]_5 138 4.6e-05 22_[1]_5 137 4.6e-05 22_[1]_5 136 4.6e-05 22_[1]_5 127 4.6e-05 21_[1]_6 77 4.6e-05 21_[1]_6 73 4.6e-05 21_[1]_6 43 4.6e-05 22_[1]_5 42 4.6e-05 21_[1]_6 7 4.6e-05 20_[1]_7 215 5.3e-05 23_[1]_4 197 5.3e-05 24_[1]_3 277 9.4e-05 18_[1]_9 201 9.4e-05 22_[1]_5 174 9.4e-05 22_[1]_5 129 9.4e-05 22_[1]_5 101 9.4e-05 21_[1]_6 6 9.4e-05 21_[1]_6 298 0.00015 22_[1]_5 264 0.00015 24_[1]_3 185 0.00015 23_[1]_4 172 0.00015 21_[1]_6 162 0.00015 23_[1]_4 128 0.00015 21_[1]_6 116 0.00015 24_[1]_3 99 0.00015 21_[1]_6 97 0.00015 18_[1]_9 83 0.00015 22_[1]_5 80 0.00015 23_[1]_4 70 0.00015 23_[1]_4 36 0.00015 23_[1]_4 272 0.00018 23_[1]_4 118 0.00018 22_[1]_5 112 0.00018 22_[1]_5 106 0.00018 22_[1]_5 95 0.00018 24_[1]_3 72 0.00018 22_[1]_5 294 0.00021 3_[1]_24 258 0.00021 22_[1]_5 247 0.00021 19_[1]_8 86 0.00021 24_[1]_3 79 0.00021 20_[1]_7 74 0.00021 20_[1]_7 65 0.00021 22_[1]_5 5 0.00021 22_[1]_5 232 0.00022 21_[1]_6 195 0.00022 21_[1]_6 28 0.00022 21_[1]_6 261 0.00026 21_[1]_6 251 0.00026 19_[1]_8 246 0.00026 21_[1]_6 191 0.00026 21_[1]_6 189 0.00026 21_[1]_6 153 0.00026 23_[1]_4 98 0.00026 23_[1]_4 61 0.00026 21_[1]_6 45 0.00026 20_[1]_7 27 0.00026 20_[1]_7 22 0.00026 22_[1]_5 18 0.00026 21_[1]_6 2 0.00026 23_[1]_4 299 0.00031 21_[1]_6 178 0.00031 22_[1]_5 168 0.00031 22_[1]_5 125 0.00031 23_[1]_4 59 0.00031 23_[1]_4 41 0.00031 20_[1]_7 34 0.00031 20_[1]_7 280 0.00036 23_[1]_4 102 0.00036 18_[1]_9 24 0.00036 18_[1]_9 300 0.00043 20_[1]_7 230 0.00045 21_[1]_6 38 0.00046 24_[1]_3 203 0.00049 19_[1]_8 57 0.00049 22_[1]_5 51 0.00049 24_[1]_3 289 0.00055 23_[1]_4 235 0.00055 20_[1]_7 226 0.00055 23_[1]_4 10 0.00055 19_[1]_8 285 0.00062 22_[1]_5 218 0.00062 22_[1]_5 176 0.00062 22_[1]_5 164 0.00062 23_[1]_4 117 0.00062 21_[1]_6 268 0.00067 21_[1]_6 257 0.00067 20_[1]_7 208 0.00067 24_[1]_3 152 0.00067 21_[1]_6 100 0.00067 23_[1]_4 92 0.00067 23_[1]_4 20 0.00067 22_[1]_5 4 0.00067 22_[1]_5 156 0.0007 20_[1]_7 209 0.00074 24_[1]_3 139 0.00074 22_[1]_5 221 0.00088 25_[1]_2 68 0.00097 24_[1]_3 222 0.001 24_[1]_3 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 in BLOCKS format -------------------------------------------------------------------------------- BL MOTIF 1 width=8 seqs=106 293 ( 20) AGGAGGAA 1 254 ( 21) AGGAGGAA 1 253 ( 23) AGGAGGAA 1 243 ( 26) AGGAGGAA 1 196 ( 23) AGGAGGAA 1 171 ( 18) AGGAGGAA 1 142 ( 23) AGGAGGAA 1 138 ( 23) AGGAGGAA 1 137 ( 23) AGGAGGAA 1 136 ( 23) AGGAGGAA 1 127 ( 22) AGGAGGAA 1 77 ( 22) AGGAGGAA 1 73 ( 22) AGGAGGAA 1 43 ( 23) AGGAGGAA 1 42 ( 22) AGGAGGAA 1 7 ( 21) AGGAGGAA 1 215 ( 24) AGGAGGCA 1 197 ( 25) AGGAGGCA 1 277 ( 19) AGGGGGAA 1 201 ( 23) AGGGGGAA 1 174 ( 23) AGGGGGAA 1 129 ( 23) AGGGGGAA 1 101 ( 22) AGGAGGAG 1 6 ( 22) AGGAGGAG 1 298 ( 23) AGGAGGAT 1 264 ( 25) AGGAGGAT 1 185 ( 24) AGGAGGAT 1 172 ( 22) AGGAGGAT 1 162 ( 24) AGGAGGAT 1 128 ( 22) AGGAGGGA 1 116 ( 25) AGGAGGAT 1 99 ( 22) AGGAGGAT 1 97 ( 19) AGGAGGGA 1 83 ( 23) AGGAGGAT 1 80 ( 24) AGGAGGAT 1 70 ( 24) AGGAGGGA 1 36 ( 24) AGGAGGAT 1 272 ( 24) AGGAGGTA 1 118 ( 23) AGGAGGTA 1 112 ( 23) AGGAGGTA 1 106 ( 23) AGGAGGTA 1 95 ( 25) AGGAGGTA 1 72 ( 23) AGGAGGTA 1 294 ( 4) GGGAGGAA 1 258 ( 23) GGGAGGAA 1 247 ( 20) GGGAGGAA 1 86 ( 25) GGGAGGAA 1 79 ( 21) GGGAGGAA 1 74 ( 21) GGGAGGAA 1 65 ( 23) GGGAGGAA 1 5 ( 23) GGGAGGAA 1 232 ( 22) AGGAGGCT 1 195 ( 22) AGGGGGAG 1 28 ( 22) AGGGGGAG 1 261 ( 22) AGGGGGAT 1 251 ( 20) AGGGGGGA 1 246 ( 22) AGGGGGAT 1 191 ( 22) AGGAGGGG 1 189 ( 22) AGGGGGAT 1 153 ( 24) AGGGGGAT 1 98 ( 24) AGGGGGAT 1 61 ( 22) AGGGGGGA 1 45 ( 21) AGGGGGAT 1 27 ( 21) AGGGGGAT 1 22 ( 23) AGGGGGGA 1 18 ( 22) AGGGGGGA 1 2 ( 24) AGGGGGAT 1 299 ( 22) AGGGGGTA 1 178 ( 23) AGGGGGTA 1 168 ( 23) AGGAGGGT 1 125 ( 24) AGGGGGTA 1 59 ( 24) AGGAGGGT 1 41 ( 21) AGGAGGTG 1 34 ( 21) AGGAGGTG 1 280 ( 24) AGGCGGAA 1 102 ( 19) GGGAGGAG 1 24 ( 19) GGGAGGAG 1 300 ( 21) TGGAGGAA 1 230 ( 22) GGGAGGTA 1 38 ( 25) AGGGGGTG 1 203 ( 20) GGGAGGCT 1 57 ( 23) TGGAGGCA 1 51 ( 25) TGGAGGCA 1 289 ( 24) TGGAGGAG 1 235 ( 21) TGGAGGAG 1 226 ( 24) AGGCGGAT 1 10 ( 20) GGGAGGGG 1 285 ( 23) TGGAGGGA 1 218 ( 23) GGGAGGTG 1 176 ( 23) GGGAGGTG 1 164 ( 24) TGGAGGAT 1 117 ( 22) GGGAGGGT 1 268 ( 22) TGGAGGTA 1 257 ( 21) TGGAGGTA 1 208 ( 25) TGGAGGTA 1 152 ( 22) AGGAGGTC 1 100 ( 24) TGGAGGTA 1 92 ( 24) AGGAGGTC 1 20 ( 23) TGGAGGTA 1 4 ( 23) TGGAGGTA 1 156 ( 21) CGGAGGAA 1 209 ( 25) TGGAGGGG 1 139 ( 23) AGGCGGTG 1 221 ( 26) AGGTGGAA 1 68 ( 25) AGGTGGAG 1 222 ( 25) AGGTGGGA 1 // -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 position-specific scoring matrix -------------------------------------------------------------------------------- log-odds matrix: alength= 20 w= 8 n= 8400 bayes= 5.17665 E= 2.6e-025 73 -262 76 96 69 -36 13 117 70 175 76 173 215 128 119 476 -130 299 -97 8 -1240 -1084 -39 -118 -129 237 -117 -200 -100 -56 -287 -88 -25 -125 -60 -41 -1253 -145 -205 -140 -1240 -1084 -39 -118 -129 237 -117 -200 -100 -56 -287 -88 -25 -125 -60 -41 -1253 -145 -205 -140 79 -108 76 96 69 10 13 117 70 175 76 173 215 128 119 476 -347 299 -97 8 -1240 -1084 -39 -118 -129 237 -117 -200 -100 -56 -287 -88 -25 -125 -60 -41 -1253 -145 -205 -140 -1240 -1084 -39 -118 -129 237 -117 -200 -100 -56 -287 -88 -25 -125 -60 -41 -1253 -145 -205 -140 44 -10 76 96 69 -45 13 117 70 175 76 173 215 128 118 476 -52 299 -97 8 44 -165 76 96 69 -19 13 117 70 175 76 173 215 128 118 476 -47 299 -97 8 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 position-specific probability matrix -------------------------------------------------------------------------------- letter-probability matrix: alength= 20 w= 8 nsites= 106 E= 2.6e-025 0.707547 0.009434 0.000000 0.000000 0.000000 0.150943 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.132075 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.735849 0.028302 0.000000 0.000000 0.000000 0.207547 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.028302 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.575472 0.056604 0.000000 0.000000 0.000000 0.141509 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.226415 0.000000 0.000000 0.000000 0.575472 0.018868 0.000000 0.000000 0.000000 0.169811 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.235849 0.000000 0.000000 0.000000 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 regular expression -------------------------------------------------------------------------------- AGG[AG]GG[AT][AT] -------------------------------------------------------------------------------- Time 104.69 secs. ******************************************************************************** ******************************************************************************** SUMMARY OF MOTIFS ******************************************************************************** -------------------------------------------------------------------------------- Combined block diagrams: non-overlapping sites with p-value < 0.0001 -------------------------------------------------------------------------------- SEQUENCE NAME COMBINED P-VALUE MOTIF DIAGRAM ------------- ---------------- ------------- 1 1.55e-01 35 2 7.18e-03 35 3 3.37e-01 35 4 1.85e-02 35 5 5.82e-03 35 6 2.62e-03 21_[1(9.37e-05)]_6 7 1.28e-03 20_[1(4.59e-05)]_7 8 5.90e-02 35 9 4.04e-01 35 10 1.53e-02 35 11 1.71e-01 35 12 5.41e-01 35 13 1.05e-01 35 14 1.05e-01 35 15 1.71e-01 35 16 1.00e+00 35 17 1.22e-01 35 18 7.18e-03 35 19 1.53e-02 35 20 1.85e-02 35 21 1.55e-01 35 22 7.18e-03 35 23 1.22e-01 35 24 1.00e-02 35 25 3.11e-02 35 26 1.90e-01 35 27 7.18e-03 35 28 6.22e-03 35 29 1.45e-01 35 30 1.90e-01 35 31 1.71e-01 35 32 1.05e-01 35 33 2.27e-02 35 34 8.56e-03 35 35 2.35e-01 35 36 4.14e-03 35 37 6.96e-01 35 38 1.29e-02 35 39 2.43e-01 35 40 3.11e-02 35 41 8.56e-03 35 42 1.28e-03 21_[1(4.59e-05)]_6 43 1.28e-03 22_[1(4.59e-05)]_5 44 1.90e-01 35 45 7.18e-03 35 46 3.66e-01 35 47 1.45e-01 35 48 3.37e-01 35 49 1.05e-01 35 50 1.45e-01 35 51 1.37e-02 35 52 5.06e-01 35 53 3.49e-01 35 54 1.22e-01 35 55 8.67e-01 35 56 5.90e-02 35 57 1.37e-02 35 58 5.00e-01 35 59 8.56e-03 35 60 6.85e-02 35 61 7.18e-03 35 62 1.22e-01 35 63 2.18e-02 35 64 4.37e-01 35 65 5.82e-03 35 66 8.51e-02 35 67 5.90e-02 35 68 2.67e-02 35 69 1.05e-01 35 70 4.14e-03 35 71 1.90e-01 35 72 5.09e-03 35 73 1.28e-03 21_[1(4.59e-05)]_6 74 5.82e-03 35 75 1.71e-01 35 76 8.51e-02 35 77 1.28e-03 21_[1(4.59e-05)]_6 78 1.05e-01 35 79 5.82e-03 35 80 4.14e-03 35 81 3.45e-02 35 82 1.53e-02 35 83 4.14e-03 35 84 4.11e-02 35 85 5.41e-01 35 86 5.82e-03 35 87 9.70e-01 35 88 1.55e-01 35 89 1.71e-01 35 90 1.71e-01 35 91 1.90e-01 35 92 1.85e-02 35 93 1.90e-01 35 94 1.71e-01 35 95 5.09e-03 35 96 1.45e-01 35 97 4.14e-03 35 98 7.18e-03 35 99 4.14e-03 35 100 1.85e-02 35 101 2.62e-03 21_[1(9.37e-05)]_6 102 1.00e-02 35 103 6.11e-01 35 104 6.11e-01 35 105 1.90e-01 35 106 5.09e-03 35 107 6.28e-01 35 108 6.28e-01 35 109 2.98e-01 35 110 3.37e-01 35 111 2.18e-02 35 112 5.09e-03 35 113 1.55e-01 35 114 2.28e-01 35 115 6.11e-01 35 116 4.14e-03 35 117 1.72e-02 35 118 5.09e-03 35 119 1.90e-01 35 120 1.71e-01 35 121 8.67e-01 35 122 6.85e-02 35 123 1.71e-01 35 124 3.49e-01 35 125 8.56e-03 35 126 2.18e-02 35 127 1.28e-03 21_[1(4.59e-05)]_6 128 4.14e-03 35 129 2.62e-03 22_[1(9.37e-05)]_5 130 1.05e-01 35 131 1.71e-01 35 132 5.90e-02 35 133 1.22e-01 35 134 2.18e-02 35 135 1.90e-01 35 136 1.28e-03 22_[1(4.59e-05)]_5 137 1.28e-03 22_[1(4.59e-05)]_5 138 1.28e-03 22_[1(4.59e-05)]_5 139 2.04e-02 35 140 9.38e-01 35 141 1.05e-01 35 142 1.28e-03 22_[1(4.59e-05)]_5 143 1.71e-01 35 144 7.40e-02 35 145 5.90e-02 35 146 1.05e-01 35 147 1.05e-01 35 148 1.45e-01 35 149 1.90e-01 35 150 1.00e+00 35 151 5.90e-02 35 152 1.85e-02 35 153 7.18e-03 35 154 2.31e-02 35 155 1.35e-01 35 156 1.93e-02 35 157 1.22e-01 35 158 1.45e-01 35 159 6.11e-01 35 160 1.90e-01 35 161 3.77e-02 35 162 4.14e-03 35 163 8.67e-01 35 164 1.72e-02 35 165 1.90e-01 35 166 9.91e-01 35 167 1.90e-01 35 168 8.56e-03 35 169 2.43e-01 35 170 1.00e+00 35 171 1.28e-03 17_[1(4.59e-05)]_10 172 4.14e-03 35 173 4.27e-01 35 174 2.62e-03 22_[1(9.37e-05)]_5 175 1.55e-01 35 176 1.72e-02 35 177 1.55e-01 35 178 8.56e-03 35 179 1.90e-01 35 180 2.27e-02 35 181 1.05e-01 35 182 1.71e-01 35 183 7.82e-01 35 184 3.37e-01 35 185 4.14e-03 35 186 4.59e-02 35 187 9.86e-01 35 188 1.00e+00 35 189 7.18e-03 35 190 2.28e-01 35 191 7.18e-03 35 192 1.93e-02 35 193 6.96e-01 35 194 1.72e-02 35 195 6.22e-03 35 196 1.28e-03 22_[1(4.59e-05)]_5 197 1.48e-03 24_[1(5.29e-05)]_3 198 1.22e-01 35 199 8.51e-02 35 200 1.05e-01 35 201 2.62e-03 22_[1(9.37e-05)]_5 202 1.71e-01 35 203 1.37e-02 35 204 1.93e-02 35 205 3.49e-01 35 206 9.80e-01 35 207 4.22e-01 35 208 1.85e-02 35 209 2.04e-02 35 210 6.11e-01 35 211 1.55e-01 35 212 1.53e-02 35 213 2.28e-01 35 214 1.72e-02 35 215 1.48e-03 23_[1(5.29e-05)]_4 216 9.95e-01 35 217 1.22e-01 35 218 1.72e-02 35 219 9.67e-01 35 220 8.67e-01 35 221 2.44e-02 35 222 2.82e-02 35 223 1.90e-01 35 224 4.44e-01 35 225 2.13e-01 35 226 1.53e-02 35 227 1.35e-01 35 228 1.05e-01 35 229 6.22e-03 35 230 1.24e-02 35 231 8.43e-01 35 232 6.22e-03 35 233 1.05e-01 35 234 1.22e-01 35 235 1.53e-02 35 236 4.15e-01 35 237 1.90e-01 35 238 6.85e-02 35 239 9.95e-01 35 240 4.15e-01 35 241 9.70e-01 35 242 2.98e-01 35 243 1.28e-03 25_[1(4.59e-05)]_2 244 6.85e-02 35 245 1.22e-01 35 246 7.18e-03 35 247 5.82e-03 35 248 8.67e-01 35 249 1.05e-01 35 250 3.37e-01 35 251 7.18e-03 35 252 2.13e-01 35 253 1.28e-03 22_[1(4.59e-05)]_5 254 1.28e-03 20_[1(4.59e-05)]_7 255 6.61e-01 35 256 5.81e-01 35 257 1.85e-02 35 258 5.82e-03 35 259 1.90e-01 35 260 1.71e-01 35 261 7.18e-03 35 262 2.28e-01 35 263 2.59e-01 35 264 4.14e-03 35 265 9.97e-01 35 266 4.59e-02 35 267 1.71e-01 35 268 1.85e-02 35 269 6.96e-01 35 270 8.51e-02 35 271 1.71e-01 35 272 5.09e-03 35 273 8.86e-01 35 274 7.33e-01 35 275 1.22e-01 35 276 1.05e-01 35 277 2.62e-03 18_[1(9.37e-05)]_9 278 1.90e-01 35 279 6.11e-01 35 280 1.00e-02 35 281 1.22e-01 35 282 1.55e-01 35 283 3.19e-02 35 284 6.85e-02 35 285 1.72e-02 35 286 9.38e-01 35 287 8.86e-01 35 288 1.71e-01 35 289 1.53e-02 35 290 7.33e-01 35 291 5.81e-01 35 292 2.43e-01 35 293 1.28e-03 19_[1(4.59e-05)]_8 294 5.82e-03 35 295 1.22e-01 35 296 1.45e-01 35 297 3.11e-02 35 298 4.14e-03 35 299 8.56e-03 35 300 1.19e-02 35 -------------------------------------------------------------------------------- ******************************************************************************** ******************************************************************************** Stopped because nmotifs = 1 reached. ******************************************************************************** CPU: kodomo.fbb.msu.ru ********************************************************************************