Mercurial > repos > iuc > meme_meme
view test-data/meme_output_txt_2.txt @ 0:e416c7c26977 draft
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit 71ac7e12419b8541746ebf8d4ba704cbbd603db1
| author | iuc |
|---|---|
| date | Mon, 21 Dec 2015 06:03:33 -0500 |
| parents | |
| children | 62d1fae3b7d3 |
line wrap: on
line source
******************************************************************************** MEME - Motif discovery tool ******************************************************************************** MEME version 4.11.0 (Release date: Thu Nov 26 17:48:49 2015 +1000) For further information on how to interpret these results or to get a copy of the MEME software please access http://meme-suite.org . 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-suite.org . ******************************************************************************** ******************************************************************************** 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= Galaxy_FASTA_Input ALPHABET= ACGT Sequence name Weight Length Sequence name Weight Length ------------- ------ ------ ------------- ------ ------ chr21_19617074_19617124_ 1.0000 50 chr21_26934381_26934431_ 1.0000 50 chr21_28217753_28217803_ 1.0000 50 chr21_31710037_31710087_ 1.0000 50 chr21_31744582_31744632_ 1.0000 50 chr21_31768316_31768366_ 1.0000 50 chr21_31914206_31914256_ 1.0000 50 chr21_31933633_31933683_ 1.0000 50 chr21_31962741_31962791_ 1.0000 50 chr21_31964683_31964733_ 1.0000 50 chr21_31973364_31973414_ 1.0000 50 chr21_31992870_31992920_ 1.0000 50 chr21_32185595_32185645_ 1.0000 50 chr21_32202076_32202126_ 1.0000 50 chr21_32253899_32253949_ 1.0000 50 chr21_32410820_32410870_ 1.0000 50 chr21_36411748_36411798_ 1.0000 50 chr21_37838750_37838800_ 1.0000 50 chr21_45705687_45705737_ 1.0000 50 chr21_45971413_45971463_ 1.0000 50 chr21_45978668_45978718_ 1.0000 50 chr21_45993530_45993580_ 1.0000 50 chr21_46020421_46020471_ 1.0000 50 chr21_46031920_46031970_ 1.0000 50 chr21_46046964_46047014_ 1.0000 50 chr21_46057197_46057247_ 1.0000 50 chr21_46086869_46086919_ 1.0000 50 chr21_46102103_46102153_ 1.0000 50 chr21_47517957_47518007_ 1.0000 50 chr21_47575506_47575556_ 1.0000 50 ******************************************************************************** ******************************************************************************** COMMAND LINE SUMMARY ******************************************************************************** This information can also be useful in the event you wish to report a problem with the MEME software. command: meme /Users/gvk/work/git_workspace/galaxy/database/files/002/dataset_2490.dat -o /Users/gvk/work/git_workspace/galaxy/database/job_working_directory/001/1929/dataset_2578_files -nostatus -sf Galaxy_FASTA_Input -dna -mod zoops -nmotifs 1 -wnsites 0.8 -minw 8 -maxw 50 -wg 11 -ws 1 -maxiter 50 -distance 0.001 -prior dirichlet -b 0.01 -plib /Users/gvk/work/git_workspace/galaxy/database/files/002/dataset_2577.dat -spmap uni -spfuzz 0.5 model: mod= zoops nmotifs= 1 evt= inf object function= E-value of product of p-values width: minw= 8 maxw= 50 width: wg= 11 ws= 1 endgaps= yes nsites: minsites= 2 maxsites= 30 wnsites= 0.8 theta: spmap= uni spfuzz= 0.5 global: substring= yes branching= no wbranch= no em: prior= dirichlet b= 0.01 maxiter= 50 distance= 0.001 data: n= 1500 N= 30 shuffle= -1 strands: + sample: seed= 0 ctfrac= -1 maxwords= -1 Dirichlet mixture priors file: dataset_2577.dat Letter frequencies in dataset: A 0.294 C 0.231 G 0.257 T 0.217 Background letter frequencies (from dataset with add-one prior applied): A 0.294 C 0.231 G 0.257 T 0.217 ******************************************************************************** ******************************************************************************** MOTIF 1 MEME width = 11 sites = 30 llr = 254 E-value = 5.1e-040 ******************************************************************************** -------------------------------------------------------------------------------- Motif 1 Description -------------------------------------------------------------------------------- Simplified A 3313:9:a798 pos.-specific C 1:3::1:::1: probability G 6756::::::2 matrix T 1:11a1a:3:: bits 2.2 * 2.0 * * 1.8 * * 1.5 * ** * Relative 1.3 * ** * Entropy 1.1 ****** (12.2 bits) 0.9 * ******* 0.7 * ******* 0.4 ** ******** 0.2 *********** 0.0 ----------- Multilevel GGGGTATAAAA consensus AACA T sequence -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 sites sorted by position p-value -------------------------------------------------------------------------------- Sequence name Start P-value Site ------------- ----- --------- ----------- chr21_46046964_46047014_ 13 4.51e-07 AAGGCCAGGA GGGGTATAAAA GCCTGAGAGC chr21_46031920_46031970_ 16 2.22e-06 ATACCCAGGG AGGGTATAAAA CCTCAGCAGC chr21_32202076_32202126_ 14 2.74e-06 CCACCAGCTT GAGGTATAAAA AGCCCTGTAC chr21_46057197_46057247_ 37 4.86e-06 ACAGGCCCTG GGCATATAAAA GCC chr21_45993530_45993580_ 8 4.86e-06 CCAAGGA GGAGTATAAAA GCCCCACAAA chr21_45971413_45971463_ 10 4.86e-06 CAGGCCCTG GGCATATAAAA GCCCCAGCAG chr21_31964683_31964733_ 14 4.86e-06 GATTCACTGA GGCATATAAAA GGCCCTCTGC chr21_47517957_47518007_ 33 6.48e-06 CCGGCGGGGC GGGGTATAAAG GGGGCGG chr21_45978668_45978718_ 5 6.48e-06 CAGA GGGGTATAAAG GTTCCGACCA chr21_32185595_32185645_ 19 6.48e-06 CACCAGAGCT GGGATATATAA AGAAGGTTCT chr21_32410820_32410870_ 22 1.38e-05 AATCACTGAG GATGTATAAAA GTCCCAGGGA chr21_31992870_31992920_ 17 1.38e-05 CACTATTGAA GATGTATAAAA TTTCATTTGC chr21_19617074_19617124_ 40 1.41e-05 CCTCGGGACG TGGGTATATAA chr21_31914206_31914256_ 16 1.61e-05 CCCACTACTT AGAGTATAAAA TCATTCTGAG chr21_46020421_46020471_ 3 1.95e-05 GA GACATATAAAA GCCAACATCC chr21_32253899_32253949_ 18 1.95e-05 CCCACCAGCA AGGATATATAA AAGCTCAGGA chr21_45705687_45705737_ 38 2.16e-05 CGTGGTCGCG GGGGTATAACA GC chr21_47575506_47575556_ 31 3.04e-05 GCTGCCGGTG AGCGTATAAAG GCCCTGGCG chr21_31744582_31744632_ 13 3.04e-05 CAGGTCTAAG AGCATATATAA CTTGGAGTCC chr21_31768316_31768366_ 1 3.67e-05 . AACGTATATAA ATGGTCCTGT chr21_26934381_26934431_ 28 3.93e-05 AGTCACAAGT GAGTTATAAAA GGGTCGCACG chr21_31933633_31933683_ 5 5.65e-05 TCAG AGTATATATAA ATGTTCCTGT chr21_31710037_31710087_ 15 6.24e-05 CCCAGGTTTC TGAGTATATAA TCGCCGCACC chr21_36411748_36411798_ 23 7.15e-05 AGTTTCAGTT GGCATCtaaaa attatataac chr21_46102103_46102153_ 37 1.39e-04 TGCCTGGGTC CAGGTATAAAG GCT chr21_46086869_46086919_ 38 1.39e-04 TGCCTGGGCC CAGGTATAAAG GC chr21_37838750_37838800_ 3 4.81e-04 ga tggttttataa ggggcctcac chr21_31962741_31962791_ 14 8.57e-04 TATAACTCAG GTTGGATAAAA TAATTTGTAC chr21_31973364_31973414_ 8 1.47e-03 aaactta aaactctataa acttaaaact chr21_28217753_28217803_ 27 2.64e-03 GGTGGGGGTG GGGGTTTCACT GGTCCACTAT -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 block diagrams -------------------------------------------------------------------------------- SEQUENCE NAME POSITION P-VALUE MOTIF DIAGRAM ------------- ---------------- ------------- chr21_46046964_46047014_ 4.5e-07 12_[+1]_27 chr21_46031920_46031970_ 2.2e-06 15_[+1]_24 chr21_32202076_32202126_ 2.7e-06 13_[+1]_26 chr21_46057197_46057247_ 4.9e-06 36_[+1]_3 chr21_45993530_45993580_ 4.9e-06 7_[+1]_32 chr21_45971413_45971463_ 4.9e-06 9_[+1]_30 chr21_31964683_31964733_ 4.9e-06 13_[+1]_26 chr21_47517957_47518007_ 6.5e-06 32_[+1]_7 chr21_45978668_45978718_ 6.5e-06 4_[+1]_35 chr21_32185595_32185645_ 6.5e-06 18_[+1]_21 chr21_32410820_32410870_ 1.4e-05 21_[+1]_18 chr21_31992870_31992920_ 1.4e-05 16_[+1]_23 chr21_19617074_19617124_ 1.4e-05 39_[+1] chr21_31914206_31914256_ 1.6e-05 15_[+1]_24 chr21_46020421_46020471_ 1.9e-05 2_[+1]_37 chr21_32253899_32253949_ 1.9e-05 17_[+1]_22 chr21_45705687_45705737_ 2.2e-05 37_[+1]_2 chr21_47575506_47575556_ 3e-05 30_[+1]_9 chr21_31744582_31744632_ 3e-05 12_[+1]_27 chr21_31768316_31768366_ 3.7e-05 [+1]_39 chr21_26934381_26934431_ 3.9e-05 27_[+1]_12 chr21_31933633_31933683_ 5.6e-05 4_[+1]_35 chr21_31710037_31710087_ 6.2e-05 14_[+1]_25 chr21_36411748_36411798_ 7.1e-05 22_[+1]_17 chr21_46102103_46102153_ 0.00014 36_[+1]_3 chr21_46086869_46086919_ 0.00014 37_[+1]_2 chr21_37838750_37838800_ 0.00048 2_[+1]_37 chr21_31962741_31962791_ 0.00086 13_[+1]_26 chr21_31973364_31973414_ 0.0015 7_[+1]_32 chr21_28217753_28217803_ 0.0026 26_[+1]_13 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 in BLOCKS format -------------------------------------------------------------------------------- BL MOTIF 1 width=11 seqs=30 chr21_46046964_46047014_ ( 13) GGGGTATAAAA 1 chr21_46031920_46031970_ ( 16) AGGGTATAAAA 1 chr21_32202076_32202126_ ( 14) GAGGTATAAAA 1 chr21_46057197_46057247_ ( 37) GGCATATAAAA 1 chr21_45993530_45993580_ ( 8) GGAGTATAAAA 1 chr21_45971413_45971463_ ( 10) GGCATATAAAA 1 chr21_31964683_31964733_ ( 14) GGCATATAAAA 1 chr21_47517957_47518007_ ( 33) GGGGTATAAAG 1 chr21_45978668_45978718_ ( 5) GGGGTATAAAG 1 chr21_32185595_32185645_ ( 19) GGGATATATAA 1 chr21_32410820_32410870_ ( 22) GATGTATAAAA 1 chr21_31992870_31992920_ ( 17) GATGTATAAAA 1 chr21_19617074_19617124_ ( 40) TGGGTATATAA 1 chr21_31914206_31914256_ ( 16) AGAGTATAAAA 1 chr21_46020421_46020471_ ( 3) GACATATAAAA 1 chr21_32253899_32253949_ ( 18) AGGATATATAA 1 chr21_45705687_45705737_ ( 38) GGGGTATAACA 1 chr21_47575506_47575556_ ( 31) AGCGTATAAAG 1 chr21_31744582_31744632_ ( 13) AGCATATATAA 1 chr21_31768316_31768366_ ( 1) AACGTATATAA 1 chr21_26934381_26934431_ ( 28) GAGTTATAAAA 1 chr21_31933633_31933683_ ( 5) AGTATATATAA 1 chr21_31710037_31710087_ ( 15) TGAGTATATAA 1 chr21_36411748_36411798_ ( 23) GGCATCTAAAA 1 chr21_46102103_46102153_ ( 37) CAGGTATAAAG 1 chr21_46086869_46086919_ ( 38) CAGGTATAAAG 1 chr21_37838750_37838800_ ( 3) TGGTTTTATAA 1 chr21_31962741_31962791_ ( 14) GTTGGATAAAA 1 chr21_31973364_31973414_ ( 8) AAACTCTATAA 1 chr21_28217753_28217803_ ( 27) GGGGTTTCACT 1 // -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 position-specific scoring matrix -------------------------------------------------------------------------------- log-odds matrix: alength= 4 w= 11 n= 1200 bayes= 5.2854 E= 5.1e-040 -14 -179 114 -112 3 -1155 137 -270 -114 20 86 -71 3 -279 122 -170 -1155 -1155 -295 215 156 -179 -1155 -170 -1155 -1155 -1155 220 172 -279 -1155 -1155 125 -1155 -1155 46 167 -179 -1155 -1155 144 -1155 -63 -270 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 position-specific probability matrix -------------------------------------------------------------------------------- letter-probability matrix: alength= 4 w= 11 nsites= 30 E= 5.1e-040 0.266667 0.066667 0.566667 0.100000 0.300000 0.000000 0.666667 0.033333 0.133333 0.266667 0.466667 0.133333 0.300000 0.033333 0.600000 0.066667 0.000000 0.000000 0.033333 0.966667 0.866667 0.066667 0.000000 0.066667 0.000000 0.000000 0.000000 1.000000 0.966667 0.033333 0.000000 0.000000 0.700000 0.000000 0.000000 0.300000 0.933333 0.066667 0.000000 0.000000 0.800000 0.000000 0.166667 0.033333 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 regular expression -------------------------------------------------------------------------------- [GA][GA][GC][GA]TATA[AT]AA -------------------------------------------------------------------------------- Time 0.17 secs. ******************************************************************************** ******************************************************************************** SUMMARY OF MOTIFS ******************************************************************************** -------------------------------------------------------------------------------- Combined block diagrams: non-overlapping sites with p-value < 0.0001 -------------------------------------------------------------------------------- SEQUENCE NAME COMBINED P-VALUE MOTIF DIAGRAM ------------- ---------------- ------------- chr21_19617074_19617124_ 5.63e-04 39_[+1(1.41e-05)] chr21_26934381_26934431_ 1.57e-03 27_[+1(3.93e-05)]_12 chr21_28217753_28217803_ 1.00e-01 50 chr21_31710037_31710087_ 2.49e-03 14_[+1(6.24e-05)]_25 chr21_31744582_31744632_ 1.22e-03 12_[+1(3.04e-05)]_27 chr21_31768316_31768366_ 1.47e-03 [+1(3.67e-05)]_39 chr21_31914206_31914256_ 6.45e-04 15_[+1(1.61e-05)]_24 chr21_31933633_31933683_ 2.26e-03 4_[+1(5.65e-05)]_35 chr21_31962741_31962791_ 3.37e-02 50 chr21_31964683_31964733_ 1.95e-04 13_[+1(4.86e-06)]_26 chr21_31973364_31973414_ 5.73e-02 50 chr21_31992870_31992920_ 5.52e-04 16_[+1(1.38e-05)]_23 chr21_32185595_32185645_ 2.59e-04 18_[+1(6.48e-06)]_21 chr21_32202076_32202126_ 1.10e-04 13_[+1(2.74e-06)]_26 chr21_32253899_32253949_ 7.78e-04 17_[+1(1.95e-05)]_22 chr21_32410820_32410870_ 5.52e-04 21_[+1(1.38e-05)]_18 chr21_36411748_36411798_ 2.85e-03 22_[+1(7.15e-05)]_17 chr21_37838750_37838800_ 1.90e-02 50 chr21_45705687_45705737_ 8.63e-04 37_[+1(2.16e-05)]_2 chr21_45971413_45971463_ 1.95e-04 9_[+1(4.86e-06)]_30 chr21_45978668_45978718_ 2.59e-04 4_[+1(6.48e-06)]_35 chr21_45993530_45993580_ 1.95e-04 7_[+1(4.86e-06)]_32 chr21_46020421_46020471_ 7.78e-04 2_[+1(1.95e-05)]_37 chr21_46031920_46031970_ 8.89e-05 15_[+1(2.22e-06)]_24 chr21_46046964_46047014_ 1.80e-05 12_[+1(4.51e-07)]_27 chr21_46057197_46057247_ 1.95e-04 36_[+1(4.86e-06)]_3 chr21_46086869_46086919_ 5.54e-03 50 chr21_46102103_46102153_ 5.54e-03 50 chr21_47517957_47518007_ 2.59e-04 32_[+1(6.48e-06)]_7 chr21_47575506_47575556_ 1.22e-03 30_[+1(3.04e-05)]_9 -------------------------------------------------------------------------------- ******************************************************************************** ******************************************************************************** Stopped because requested number of motifs (1) found. ******************************************************************************** CPU: dot1x-cb-51.aset.psu.edu ********************************************************************************