******************************************************************************** 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= temp.fasta ALPHABET= ACDEFGHIKLMNPQRSTVWY Sequence name Weight Length Sequence name Weight Length ------------- ------ ------ ------------- ------ ------ Q08787(my_protein) 1.0000 1275 AACS_DANRE 1.0000 673 FAT1_COCHE 1.0000 643 LCFH_HAEIN 1.0000 607 ORA_PLAFN 1.0000 701 Y5435_MYCS2 1.0000 545 ACSL6_MOUSE 1.0000 697 BCLA_THAAR 1.0000 527 FAA22_MYCMM 1.0000 702 SRFAA_BACSU 1.0000 3587 ******************************************************************************** ******************************************************************************** COMMAND LINE SUMMARY ******************************************************************************** This information can also be useful in the event you wish to report a problem with the MEME software. command: meme temp.fasta -mod zoops -nmotifs 3 -prior dirichlet -nostatus -protein -text model: mod= zoops nmotifs= 3 evt= inf object function= E-value of product of p-values width: minw= 8 maxw= 50 minic= 0.00 width: wg= 11 ws= 1 endgaps= yes nsites: minsites= 2 maxsites= 10 wnsites= 0.8 theta: prob= 1 spmap= pam spfuzz= 120 global: substring= yes branching= no wbranch= no em: prior= dirichlet b= 0.01 maxiter= 50 distance= 1e-05 data: n= 9957 N= 10 sample: seed= 0 seqfrac= 1 Letter frequencies in dataset: A 0.082 C 0.010 D 0.059 E 0.075 F 0.039 G 0.072 H 0.024 I 0.060 K 0.057 L 0.101 M 0.023 N 0.036 P 0.045 Q 0.040 R 0.050 S 0.058 T 0.057 V 0.068 W 0.012 Y 0.034 Background letter frequencies (from dataset with add-one prior applied): A 0.081 C 0.010 D 0.059 E 0.075 F 0.039 G 0.072 H 0.024 I 0.060 K 0.057 L 0.101 M 0.023 N 0.036 P 0.045 Q 0.040 R 0.050 S 0.058 T 0.057 V 0.068 W 0.012 Y 0.034 ******************************************************************************** ******************************************************************************** MOTIF 1 width = 26 sites = 10 llr = 488 E-value = 1.1e-057 ******************************************************************************** -------------------------------------------------------------------------------- Motif 1 Description -------------------------------------------------------------------------------- Simplified A :1215:1::::::::2::22:::13: pos.-specific C ::::1::1::::::::::1:::::1: probability D :55:::::::::::::::::::1::2 matrix E 1::::::::::::::1:::::::::: F 1:1:12::1:::::::::::1::::: G 1::::1:::::a::a:::7:::::31 H ::11:::::::::::::::::3:6:: I :::1:142::::::::::::321::: K :::::::::::::::2:a:::::::: L 2::311211::::::1:::::3:::: M :::1::12:::::::::::13::::: N :2:::1:::::::::1:::2:::1:6 P 4::2:::::::::::1a::::::::: Q :::::::1::::::::::::1::::: R ::::1:::::::::::::::::1:21 S 121:11:::2a:21:1::::::3:1: T :::::1:2:8::89:1:::1:141:: V :::1::11:::::::::::421:::: W ::::::1::::::::::::::::::: Y :::::2::8::::::::::::::1:: bits 6.6 6.0 5.3 4.6 * Relative 4.0 * ** ** ** Entropy 3.3 ******* ** * (70.3 bits) 2.6 * ******* *** * * * 2.0 *** * ********* ********** 1.3 ************************** 0.7 ************************** 0.0 -------------------------- Multilevel PDDLAFIIYTSGTTGAPKGVIHTHAN consensus LNAP YLM S S K AAMLS GD sequence S T NVI R -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 sites sorted by position p-value -------------------------------------------------------------------------------- Sequence name Start P-value Site ------------- ----- --------- -------------------------- Q08787(my_protein) 606 1.61e-28 QASDPATAID PNDPAYIMYTSGTTGKPKGNITTHAN IQGLVKHVDY ACSL6_MOUSE 267 2.51e-24 ENHHAPVPPR PDDLSIVCFTSGTTGNPKGAMLTHGN VVADFSGFLK BCLA_THAAR 178 4.84e-23 SEQFEAAPTR PDDHCFWLYSSGSTGAPKGTVHIHSD LIHTAELYAR SRFAA_BACSU 1647 5.41e-23 DARNLALPLD SAAMANLTYTSGTTGTPKGNIVTHAN ILRTVKETNY ORA_PLAFN 197 5.41e-23 IANVTVQNED PNFIASIVYTSGTSGKPKGVMLSNRN LYNGVIPPCD LCFH_HAEIN 184 7.50e-23 ELTQRLNQKQ LSDLFTIIYTSGTTGEPKGVMLDYAN LAHQLETHDL FAA22_MYCMM 144 1.72e-22 VGPADYEILG GDAAAYATYTSGTTGPPKAAIHRHCD VFAFVEAMCR AACS_DANRE 282 7.74e-22 DPQLEFEQLP FSHPLFIMYSSGTTGAPKCMVHSAGG TLIQHLKEHI Y5435_MYCS2 151 1.93e-21 ADPIDPVETG EDDLALMQLTSGSTGSPKAVQITHRN IHSNAEAMFV FAT1_COCHE 227 1.28e-20 IPDSRTENIE LDSVRGLIYTSGTTGLPKGVFISTGR ELRTDWSISK -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 block diagrams -------------------------------------------------------------------------------- SEQUENCE NAME POSITION P-VALUE MOTIF DIAGRAM ------------- ---------------- ------------- Q08787(my_protein) 1.6e-28 605_[1]_644 ACSL6_MOUSE 2.5e-24 266_[1]_405 BCLA_THAAR 4.8e-23 177_[1]_324 SRFAA_BACSU 5.4e-23 1646_[1]_1915 ORA_PLAFN 5.4e-23 196_[1]_479 LCFH_HAEIN 7.5e-23 183_[1]_398 FAA22_MYCMM 1.7e-22 143_[1]_533 AACS_DANRE 7.7e-22 281_[1]_366 Y5435_MYCS2 1.9e-21 150_[1]_369 FAT1_COCHE 1.3e-20 226_[1]_391 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 in BLOCKS format -------------------------------------------------------------------------------- BL MOTIF 1 width=26 seqs=10 Q08787(my_protein) ( 606) PNDPAYIMYTSGTTGKPKGNITTHAN 1 ACSL6_MOUSE ( 267) PDDLSIVCFTSGTTGNPKGAMLTHGN 1 BCLA_THAAR ( 178) PDDHCFWLYSSGSTGAPKGTVHIHSD 1 SRFAA_BACSU ( 1647) SAAMANLTYTSGTTGTPKGNIVTHAN 1 ORA_PLAFN ( 197) PNFIASIVYTSGTSGKPKGVMLSNRN 1 LCFH_HAEIN ( 184) LSDLFTIIYTSGTTGEPKGVMLDYAN 1 FAA22_MYCMM ( 144) GDAAAYATYTSGTTGPPKAAIHRHCD 1 AACS_DANRE ( 282) FSHPLFIMYSSGTTGAPKCMVHSAGG 1 Y5435_MYCS2 ( 151) EDDLALMQLTSGSTGSPKAVQITHRN 1 FAT1_COCHE ( 227) LDSVRGLIYTSGTTGLPKGVFISTGR 1 // -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 position-specific scoring matrix -------------------------------------------------------------------------------- log-odds matrix: alength= 20 w= 26 n= 9707 bayes= 10.1732 E= 1.1e-057 -997 -997 -997 42 137 48 -997 -997 -997 99 -997 -997 315 -997 -997 79 -997 -997 -997 -997 30 -997 308 -997 -997 -997 -997 -997 -997 -997 -997 247 -997 -997 -997 179 -997 -997 -997 -997 129 -997 308 -997 137 -997 208 -997 -997 -997 -997 -997 -997 -997 -997 79 -997 -997 -997 -997 30 -997 -997 -997 -997 -997 208 75 -997 157 212 -997 215 -997 -997 -997 -997 55 -997 -997 262 329 -997 -997 137 -997 -997 -997 -997 -1 -997 -997 -997 -997 101 79 -997 -997 -997 -997 -997 -997 -997 -997 237 48 -997 75 -997 -1 -997 147 -997 -997 -997 79 82 -997 -997 255 30 -997 -997 -997 -997 -997 -997 274 -997 99 212 -997 -997 -997 -997 -997 -997 55 302 -997 -997 329 -997 -997 -997 -997 -997 174 -997 -1 312 -997 -997 131 -997 -997 182 55 -997 -997 -997 -997 -997 -997 137 -997 -997 -997 -997 -1 -997 -997 -997 -997 -997 -997 -997 -997 -997 455 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 179 382 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 412 -997 -997 -997 -997 -997 -997 -997 -997 -997 380 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 179 382 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 79 399 -997 -997 -997 -997 -997 -997 -997 -997 380 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 129 -997 -997 42 -997 -997 -997 -997 181 -1 -997 147 115 -997 -997 79 82 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 448 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 413 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 129 329 -997 -997 -997 329 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 129 -997 -997 -997 -997 -997 -997 -997 -997 -997 212 247 -997 -997 -997 -997 82 255 -997 -997 -997 -997 -997 -997 137 -997 -997 233 -997 -997 370 -997 -997 131 -997 -997 -997 155 -997 -997 -997 -997 -997 -997 -997 -997 366 174 -997 157 -997 -997 -997 -997 -997 -997 82 55 -997 -997 -997 -997 76 -997 -997 -997 -997 75 -997 -997 -997 -997 -997 -997 101 238 282 -997 -997 -997 30 -997 -997 -997 -997 -997 466 -997 -997 -997 -997 147 -997 -997 -997 -997 82 -997 -997 155 188 329 -997 -997 -997 206 -997 -997 -997 -997 -997 -997 -997 -997 201 79 -997 -997 -997 -997 -997 -997 176 -997 -997 48 -997 -997 -997 -997 -997 405 -997 -997 101 -997 -997 -997 -997 -997 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 position-specific probability matrix -------------------------------------------------------------------------------- letter-probability matrix: alength= 20 w= 26 nsites= 10 E= 1.1e-057 0.000000 0.000000 0.000000 0.100000 0.100000 0.100000 0.000000 0.000000 0.000000 0.200000 0.000000 0.000000 0.400000 0.000000 0.000000 0.100000 0.000000 0.000000 0.000000 0.000000 0.100000 0.000000 0.500000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.200000 0.000000 0.000000 0.000000 0.200000 0.000000 0.000000 0.000000 0.000000 0.200000 0.000000 0.500000 0.000000 0.100000 0.000000 0.100000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.100000 0.000000 0.000000 0.000000 0.000000 0.100000 0.000000 0.000000 0.000000 0.000000 0.000000 0.100000 0.100000 0.000000 0.300000 0.100000 0.000000 0.200000 0.000000 0.000000 0.000000 0.000000 0.100000 0.000000 0.000000 0.500000 0.100000 0.000000 0.000000 0.100000 0.000000 0.000000 0.000000 0.000000 0.100000 0.000000 0.000000 0.000000 0.000000 0.100000 0.100000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.200000 0.100000 0.000000 0.100000 0.000000 0.100000 0.000000 0.100000 0.000000 0.000000 0.000000 0.100000 0.100000 0.000000 0.000000 0.200000 0.100000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.400000 0.000000 0.200000 0.100000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.100000 0.100000 0.000000 0.000000 0.100000 0.000000 0.000000 0.000000 0.000000 0.000000 0.200000 0.000000 0.100000 0.200000 0.000000 0.000000 0.100000 0.000000 0.000000 0.200000 0.100000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.100000 0.000000 0.000000 0.000000 0.000000 0.100000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.800000 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.200000 0.800000 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 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 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.200000 0.800000 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.100000 0.900000 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.200000 0.000000 0.000000 0.100000 0.000000 0.000000 0.000000 0.000000 0.200000 0.100000 0.000000 0.100000 0.100000 0.000000 0.000000 0.100000 0.100000 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.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.200000 0.100000 0.000000 0.000000 0.000000 0.700000 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.200000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.100000 0.200000 0.000000 0.000000 0.000000 0.000000 0.100000 0.400000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.100000 0.000000 0.000000 0.300000 0.000000 0.000000 0.300000 0.000000 0.000000 0.100000 0.000000 0.000000 0.000000 0.200000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.300000 0.200000 0.000000 0.300000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.100000 0.100000 0.000000 0.000000 0.000000 0.000000 0.100000 0.000000 0.000000 0.000000 0.000000 0.100000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.100000 0.300000 0.400000 0.000000 0.000000 0.000000 0.100000 0.000000 0.000000 0.000000 0.000000 0.000000 0.600000 0.000000 0.000000 0.000000 0.000000 0.100000 0.000000 0.000000 0.000000 0.000000 0.100000 0.000000 0.000000 0.100000 0.300000 0.100000 0.000000 0.000000 0.000000 0.300000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.200000 0.100000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.200000 0.000000 0.000000 0.100000 0.000000 0.000000 0.000000 0.000000 0.000000 0.600000 0.000000 0.000000 0.100000 0.000000 0.000000 0.000000 0.000000 0.000000 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 regular expression -------------------------------------------------------------------------------- [PL][DNS][DA][LP]A[FY][IL][IMT]Y[TS]SG[TS]TG[AK]PK[GA][VAN][IMV][HLI][TS]H[AGR][ND] -------------------------------------------------------------------------------- Time 6.78 secs. ******************************************************************************** ******************************************************************************** MOTIF 2 width = 18 sites = 7 llr = 244 E-value = 7.8e-008 ******************************************************************************** -------------------------------------------------------------------------------- Motif 2 Description -------------------------------------------------------------------------------- Simplified A :::::::::3::::::1: pos.-specific C ::4::::::::::::::: probability D 1:::::1::::::::1:3 matrix E 9:::::1:::1::::117 F :::::::::3:::1:::: G :::::913:::7:::::: H :::::::::::::::::: I :3:4::::1::::::::: K ::::3:::::6:::311: L :73:::::1::::11::: M ::::::1111:1:::::: N :::::::1::::::61:: P ::::::4:::::::::6: Q ::1::::::::::::::: R ::::3::::::::::3:: S ::::41:41111:::::: T :::::::::::::::::: V :::6::::411::1:::: W ::1::::::::::3:1:: Y ::::::::::::a3:::: bits 6.6 6.0 5.3 4.6 * Relative 4.0 * Entropy 3.3 * * * * * (50.3 bits) 2.6 ****** * ***** ** 2.0 ****************** 1.3 ****************** 0.7 ****************** 0.0 ------------------ Multilevel ELCVSGPSVAKGYWNRPE consensus ILIK G F YK D sequence R -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 2 sites sorted by position p-value -------------------------------------------------------------------------------- Sequence name Start P-value Site ------------- ----- --------- ------------------ SRFAA_BACSU 1846 5.41e-18 GHVQPAGIAG ELCVSGEGLVKGYYNRPE LTEEKFVPHP BCLA_THAAR 377 1.26e-17 AEITTAGVAG ELQISGPSSAVMYWNNPE KTAATFMGEW FAA22_MYCMM 341 2.43e-17 GAAAGPGVEG DLWVRGPSIAESYWNWPE PLLTDEGWLD Q08787(my_protein) 805 2.76e-17 SQLQPFGAVG ELCISGMGVSKGYVNRAD LTKEKFIENP LCFH_HAEIN 418 9.22e-17 KAEVKIGENN EILVRGGMVMKGYYKKPE ETAQAFTEDG ACSL6_MOUSE 507 1.60e-16 LNYWTCKGEG EICVKGPNVFKGYLKDED RTKEALDSDG ORA_PLAFN 437 1.02e-14 YKATDTIPKG ELLIKSDSMFSGYFLEKE STEHAFTNDG -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 2 block diagrams -------------------------------------------------------------------------------- SEQUENCE NAME POSITION P-VALUE MOTIF DIAGRAM ------------- ---------------- ------------- SRFAA_BACSU 5.4e-18 1845_[2]_1724 BCLA_THAAR 1.3e-17 376_[2]_133 FAA22_MYCMM 2.4e-17 340_[2]_344 Q08787(my_protein) 2.8e-17 804_[2]_453 LCFH_HAEIN 9.2e-17 417_[2]_172 ACSL6_MOUSE 1.6e-16 506_[2]_173 ORA_PLAFN 1e-14 436_[2]_247 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 2 in BLOCKS format -------------------------------------------------------------------------------- BL MOTIF 2 width=18 seqs=7 SRFAA_BACSU ( 1846) ELCVSGEGLVKGYYNRPE 1 BCLA_THAAR ( 377) ELQISGPSSAVMYWNNPE 1 FAA22_MYCMM ( 341) DLWVRGPSIAESYWNWPE 1 Q08787(my_protein) ( 805) ELCISGMGVSKGYVNRAD 1 LCFH_HAEIN ( 418) EILVRGGMVMKGYYKKPE 1 ACSL6_MOUSE ( 507) EICVKGPNVFKGYLKDED 1 ORA_PLAFN ( 437) ELLIKSDSMFSGYFLEKE 1 // -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 2 position-specific scoring matrix -------------------------------------------------------------------------------- log-odds matrix: alength= 20 w= 18 n= 9787 bayes= 11.0543 E= 7.8e-008 -945 -945 127 352 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 226 -945 282 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 539 -945 -945 -945 -945 -945 -945 -945 150 -945 -945 -945 183 -945 -945 -945 -945 353 -945 -945 -945 -945 -945 -945 -945 -945 284 -945 -945 -945 -945 -945 -945 -945 -945 -945 306 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 233 -945 -945 -945 -945 -945 252 289 -945 -945 -945 -945 -945 -945 -945 -945 -945 358 -945 -945 -945 -945 -945 -945 -945 -945 -945 131 -945 -945 -945 -945 -945 -945 127 93 -945 99 -945 -945 -945 -945 263 -945 325 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 199 -945 -945 -945 -945 263 198 -945 -945 -945 289 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 126 -945 50 263 -945 -945 -945 -945 131 -945 265 -945 -945 181 -945 -945 -945 289 -945 -945 -945 -945 -945 263 -945 -945 -945 -945 131 -945 106 -945 -945 -945 -945 -945 93 -945 -945 -945 -945 333 -945 -945 -945 -945 -945 -945 131 -945 106 -945 -945 -945 -945 -945 -945 -945 332 -945 -945 -945 -945 263 -945 -945 -945 -945 131 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 487 -945 -945 -945 -945 189 -945 -945 -945 -945 50 -945 -945 -945 -945 -945 -945 -945 106 453 307 -945 -945 -945 -945 -945 -945 -945 -945 233 50 -945 398 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 127 93 -945 -945 -945 -945 133 -945 -945 198 -945 -945 252 -945 -945 -945 353 -945 81 -945 -945 93 -945 -945 -945 -945 133 -945 -945 -945 367 -945 -945 -945 -945 -945 -945 -945 -945 -945 227 325 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -945 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 2 position-specific probability matrix -------------------------------------------------------------------------------- letter-probability matrix: alength= 20 w= 18 nsites= 7 E= 7.8e-008 0.000000 0.000000 0.142857 0.857143 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 0.000000 0.000000 0.000000 0.000000 0.285714 0.000000 0.714286 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.428571 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.285714 0.000000 0.000000 0.000000 0.142857 0.000000 0.000000 0.000000 0.000000 0.142857 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.428571 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.571429 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.285714 0.000000 0.000000 0.000000 0.000000 0.000000 0.285714 0.428571 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.857143 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.142857 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.142857 0.142857 0.000000 0.142857 0.000000 0.000000 0.000000 0.000000 0.142857 0.000000 0.428571 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.285714 0.000000 0.000000 0.000000 0.000000 0.142857 0.142857 0.000000 0.000000 0.000000 0.428571 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.142857 0.000000 0.142857 0.142857 0.000000 0.000000 0.000000 0.000000 0.142857 0.000000 0.428571 0.000000 0.000000 0.285714 0.000000 0.000000 0.000000 0.285714 0.000000 0.000000 0.000000 0.000000 0.000000 0.142857 0.000000 0.000000 0.000000 0.000000 0.142857 0.000000 0.142857 0.000000 0.000000 0.000000 0.000000 0.000000 0.142857 0.000000 0.000000 0.000000 0.000000 0.571429 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.142857 0.000000 0.142857 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.714286 0.000000 0.000000 0.000000 0.000000 0.142857 0.000000 0.000000 0.000000 0.000000 0.142857 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 0.000000 0.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.000000 0.000000 0.142857 0.000000 0.000000 0.000000 0.000000 0.142857 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.142857 0.285714 0.285714 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.285714 0.142857 0.000000 0.571429 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.142857 0.142857 0.000000 0.000000 0.000000 0.000000 0.142857 0.000000 0.000000 0.142857 0.000000 0.000000 0.285714 0.000000 0.000000 0.000000 0.142857 0.000000 0.142857 0.000000 0.000000 0.142857 0.000000 0.000000 0.000000 0.000000 0.142857 0.000000 0.000000 0.000000 0.571429 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.285714 0.714286 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 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 2 regular expression -------------------------------------------------------------------------------- E[LI][CL][VI][SKR]GP[SG]V[AF]KGY[WY][NK]RP[ED] -------------------------------------------------------------------------------- Time 13.11 secs. ******************************************************************************** ******************************************************************************** MOTIF 3 width = 15 sites = 10 llr = 243 E-value = 3.3e-003 ******************************************************************************** -------------------------------------------------------------------------------- Motif 3 Description -------------------------------------------------------------------------------- Simplified A ::::21::1::11:: pos.-specific C ::::::::::1:::1 probability D :::::::::1::::: matrix E :1::::::::::1:: F :::::::::::2::2 G 2:::1:7:::::::: H :::::::::1::::: I 3:1:::::3::1:11 K ::2::::8:1113:: L 315:::3131::171 M ::::::::11::::: N 11:::2:::1:12:: P :5::5::::::1::1 Q :::::::::1::::: R ::1::::::372::3 S :::::2:::::1:1: T :::a25::::::::: V ::::::::2:::111 W 1::::::1::1:1:: Y :21:::::::::::: bits 6.6 6.0 5.3 4.6 Relative 4.0 * Entropy 3.3 * * * (35.1 bits) 2.6 * ***** * 2.0 ********* * ** 1.3 *************** 0.7 *************** 0.0 --------------- Multilevel IPLTPTGKIRRFKLR consensus LYK ANL L RN F sequence G TS V -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 3 sites sorted by position p-value -------------------------------------------------------------------------------- Sequence name Start P-value Site ------------- ----- --------- --------------- SRFAA_BACSU 943 2.58e-14 VPAYFTEVTE IPLTPSGKVDRRKLF ALEVKAVSGT BCLA_THAAR 511 3.09e-14 YPRWIEFVDD LPKTATGKIQRFKLR SA ORA_PLAFN 609 1.22e-13 YLTSKPWDTT NYLTPTLKIRRFNVF KDFSFFIDEV AACS_DANRE 329 7.07e-13 HGNMTFNDVI IYYTTTGWMMWNWLI SSLAVGASVV Y5435_MYCS2 526 1.95e-12 RNVVVLGPGS IPKTPSGKLRRANSV SLVTQ FAA22_MYCMM 473 1.83e-11 VPHRFAVVER LPRTANGKLLRSALR GQTPAKPIWE Q08787(my_protein) 941 2.17e-11 VPQTFTFLDE LPLTTNGKVNKRLLP KPDQDQLAEE ACSL6_MOUSE 668 3.49e-11 IHCDMFSVQN GLLTPTLKAKRPELR EYFKKQIEEL FAT1_COCHE 157 3.49e-11 AIGAATSFLN WNLTGAGLIHCIKLC ECRFVIADID LCFH_HAEIN 577 1.26e-10 LLSQAFSIKL GEITPTLKLRRKVIL ERYRKQIEAM -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 3 block diagrams -------------------------------------------------------------------------------- SEQUENCE NAME POSITION P-VALUE MOTIF DIAGRAM ------------- ---------------- ------------- SRFAA_BACSU 2.6e-14 942_[3]_2630 BCLA_THAAR 3.1e-14 510_[3]_2 ORA_PLAFN 1.2e-13 608_[3]_78 AACS_DANRE 7.1e-13 328_[3]_330 Y5435_MYCS2 1.9e-12 525_[3]_5 FAA22_MYCMM 1.8e-11 472_[3]_215 Q08787(my_protein) 2.2e-11 940_[3]_320 ACSL6_MOUSE 3.5e-11 667_[3]_15 FAT1_COCHE 3.5e-11 156_[3]_472 LCFH_HAEIN 1.3e-10 576_[3]_16 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 3 in BLOCKS format -------------------------------------------------------------------------------- BL MOTIF 3 width=15 seqs=10 SRFAA_BACSU ( 943) IPLTPSGKVDRRKLF 1 BCLA_THAAR ( 511) LPKTATGKIQRFKLR 1 ORA_PLAFN ( 609) NYLTPTLKIRRFNVF 1 AACS_DANRE ( 329) IYYTTTGWMMWNWLI 1 Y5435_MYCS2 ( 526) IPKTPSGKLRRANSV 1 FAA22_MYCMM ( 473) LPRTANGKLLRSALR 1 Q08787(my_protein) ( 941) LPLTTNGKVNKRLLP 1 ACSL6_MOUSE ( 668) GLLTPTLKAKRPELR 1 FAT1_COCHE ( 157) WNLTGAGLIHCIKLC 1 LCFH_HAEIN ( 577) GEITPTLKLRRKVIL 1 // -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 3 position-specific scoring matrix -------------------------------------------------------------------------------- log-odds matrix: alength= 20 w= 15 n= 9817 bayes= 10.1894 E= 3.3e-003 -997 -997 -997 -997 -997 148 -997 233 -997 157 -997 147 -997 -997 -997 -997 -997 -997 302 -997 -997 -997 -997 42 -997 -997 -997 -997 -997 -1 -997 147 348 -997 -997 -997 -997 -997 -997 255 -997 -997 -997 -997 -997 -997 -997 75 181 231 -997 -997 -997 -997 101 -997 -997 -997 -997 155 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 414 -997 -997 -997 129 -997 -997 -997 -997 48 -997 -997 -997 -997 -997 -997 348 -997 -997 -997 182 -997 -997 -997 30 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 247 -997 -997 -997 179 314 -997 -997 -997 -997 -997 -997 -997 -997 329 -997 -997 -997 157 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 -997 381 -1 -997 -997 -997 -997 -997 -997 -997 -997 302 -997 30 -997 -997 -997 -997 -997 -997 233 -997 157 212 -997 -997 -997 -997 -997 -997 155 -997 -997 -997 -997 76 -997 -997 -997 208 -997 81 -1 212 147 -997 131 259 -997 -997 -997 -997 -997 -997 329 -997 -997 -997 -997 -997 -997 81 -997 -997 -997 -997 -997 381 -997 -997 -997 302 -997 30 -997 -997 -997 237 -997 -997 75 81 -997 -997 147 115 -997 201 79 -997 -997 -997 -997 30 -997 -997 42 -997 -997 -997 -997 240 -1 -997 247 -997 -997 -997 -997 -997 55 302 -997 -997 -997 -997 -997 -997 -997 -997 75 -997 279 -997 -997 -997 -997 -997 79 -997 55 -997 -997 -997 329 -997 -997 237 -997 -997 75 -997 -1 -997 -997 115 -997 259 -997 -997 55 -997 -997 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 3 position-specific probability matrix -------------------------------------------------------------------------------- letter-probability matrix: alength= 20 w= 15 nsites= 10 E= 3.3e-003 0.000000 0.000000 0.000000 0.000000 0.000000 0.200000 0.000000 0.300000 0.000000 0.300000 0.000000 0.100000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.100000 0.000000 0.000000 0.000000 0.000000 0.100000 0.000000 0.000000 0.000000 0.000000 0.000000 0.100000 0.000000 0.100000 0.500000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.200000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.100000 0.200000 0.500000 0.000000 0.000000 0.000000 0.000000 0.100000 0.000000 0.000000 0.000000 0.000000 0.100000 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.200000 0.000000 0.000000 0.000000 0.000000 0.100000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.500000 0.000000 0.000000 0.000000 0.200000 0.000000 0.000000 0.000000 0.100000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.200000 0.000000 0.000000 0.000000 0.200000 0.500000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.700000 0.000000 0.000000 0.000000 0.300000 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.800000 0.100000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.100000 0.000000 0.100000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.300000 0.000000 0.300000 0.100000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.200000 0.000000 0.000000 0.000000 0.000000 0.100000 0.000000 0.000000 0.000000 0.100000 0.000000 0.100000 0.100000 0.100000 0.100000 0.000000 0.100000 0.300000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.100000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.100000 0.000000 0.000000 0.000000 0.000000 0.000000 0.700000 0.000000 0.000000 0.000000 0.100000 0.000000 0.100000 0.000000 0.000000 0.000000 0.200000 0.000000 0.000000 0.100000 0.100000 0.000000 0.000000 0.100000 0.100000 0.000000 0.200000 0.100000 0.000000 0.000000 0.000000 0.000000 0.100000 0.000000 0.000000 0.100000 0.000000 0.000000 0.000000 0.000000 0.300000 0.100000 0.000000 0.200000 0.000000 0.000000 0.000000 0.000000 0.000000 0.100000 0.100000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.100000 0.000000 0.700000 0.000000 0.000000 0.000000 0.000000 0.000000 0.100000 0.000000 0.100000 0.000000 0.000000 0.000000 0.100000 0.000000 0.000000 0.200000 0.000000 0.000000 0.100000 0.000000 0.100000 0.000000 0.000000 0.100000 0.000000 0.300000 0.000000 0.000000 0.100000 0.000000 0.000000 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 3 regular expression -------------------------------------------------------------------------------- [ILG][PY][LK]T[PAT][TNS][GL]K[ILV]RR[FR][KN]L[RF] -------------------------------------------------------------------------------- Time 19.26 secs. ******************************************************************************** ******************************************************************************** SUMMARY OF MOTIFS ******************************************************************************** -------------------------------------------------------------------------------- Combined block diagrams: non-overlapping sites with p-value < 0.0001 -------------------------------------------------------------------------------- SEQUENCE NAME COMBINED P-VALUE MOTIF DIAGRAM ------------- ---------------- ------------- Q08787(my_protein) 1.08e-42 271_[1(3.81e-05)]_308_[1(1.61e-28)]_173_[2(2.76e-17)]_118_[3(2.17e-11)]_320 AACS_DANRE 3.24e-26 3_[2(9.59e-05)]_260_[1(7.74e-22)]_21_[3(7.07e-13)]_330 FAT1_COCHE 1.85e-22 156_[3(3.49e-11)]_55_[1(1.28e-20)]_333_[3(2.43e-06)]_43 LCFH_HAEIN 7.62e-37 183_[1(7.50e-23)]_208_[2(9.22e-17)]_141_[3(1.26e-10)]_16 ORA_PLAFN 9.58e-38 196_[1(5.41e-23)]_214_[2(1.02e-14)]_154_[3(1.22e-13)]_78 Y5435_MYCS2 1.60e-24 150_[1(1.93e-21)]_349_[3(1.95e-12)]_5 ACSL6_MOUSE 2.03e-38 266_[1(2.51e-24)]_214_[2(1.60e-16)]_143_[3(3.49e-11)]_15 BCLA_THAAR 1.33e-41 177_[1(4.84e-23)]_173_[2(1.26e-17)]_116_[3(3.09e-14)]_2 FAA22_MYCMM 1.09e-37 143_[1(1.72e-22)]_171_[2(2.43e-17)]_114_[3(1.83e-11)]_215 SRFAA_BACSU 1.67e-39 444_[2(5.24e-05)]_135_[1(8.08e-13)]_181_[2(7.82e-07)]_120_[3(2.58e-14)]_689_[1(5.41e-23)]_173_[2(5.41e-18)]_118_[3(5.90e-12)]_693_[1(2.41e-11)]_161_[2(1.34e-10)]_116_[3(1.36e-09)]_196_[2(9.43e-05)]_131_[3(5.07e-05)]_202 -------------------------------------------------------------------------------- ******************************************************************************** ******************************************************************************** Stopped because nmotifs = 3 reached. ******************************************************************************** CPU: kodomo ********************************************************************************