mast.adh.oops.mini-drosoph.s

EM算法的改进

********************************************************************************
MAST - Motif Alignment and Search Tool
********************************************************************************
	MAST version 3.5.1 (Release date: 2006/02/01 02:08:55)

	For further information on how to interpret these results or to get
	a copy of the MAST software please access http://meme.nbcr.net.
********************************************************************************


********************************************************************************
REFERENCE
********************************************************************************
	If you use this program in your research, please cite:

	Timothy L. Bailey and Michael Gribskov,
	"Combining evidence using p-values: application to sequence homology
	searches", Bioinformatics, 14(48-54), 1998.
********************************************************************************


********************************************************************************
DATABASE AND MOTIFS
********************************************************************************
	DATABASE -stdin (nucleotide)
	Database contains 4 sequences, 499297 residues

	Scores for positive and reverse complement strands are combined.

	MOTIFS /home/meme/TEST/tests/meme.adh.oops (peptide)
	MOTIF WIDTH BEST POSSIBLE MATCH
	----- ----- -------------------
	  1    29   YSASKFAVRMLTRSMAHEYAPHGIRVNCI
	  2    29   KVVLITGCSSGIGKATAKHLHKEGAKVVL

	PAIRWISE MOTIF CORRELATIONS:
	MOTIF     1
	----- -----
	   2   0.30
	No overly similar pairs (correlation > 0.60) found.

	Random model letter frequencies (from non-redundant database):
	A 0.070 C 0.024 D 0.040 E 0.052 F 0.040 G 0.074 H 0.029 I 0.041 K 0.052 
	L 0.096 M 0.017 N 0.032 P 0.065 Q 0.042 R 0.067 S 0.084 T 0.052 V 0.059 
	W 0.016 Y 0.022 
********************************************************************************


********************************************************************************
SECTION I: HIGH-SCORING SEQUENCES
********************************************************************************
	- Each of the following 4 sequences has E-value less than 10.
	- The E-value of a sequence is the expected number of sequences
	  in a random database of the same size that would match the motifs as
	  well as the sequence does and is equal to the combined p-value of the
	  sequence times the number of sequences in the database.
	- The combined p-value of a sequence measures the strength of the
	  match of the sequence to all the motifs and is calculated by
	    o finding the score of the single best match of each motif
	      to the sequence (best matches may overlap),
	    o calculating the sequence p-value of each score,
	    o forming the product of the p-values,
	    o taking the p-value of the product.
	- The sequence p-value of a score is defined as the
	  probability of a random sequence of the same length containing
	  some match with as good or better a score.
	- The score for the match of a position in a sequence to a motif
	  is computed by by summing the appropriate entry from each column of
	  the position-dependent scoring matrix that represents the motif.
	- Sequences shorter than one or more of the motifs are skipped.
	- The frame of the (best) motif match(es) is shown.
	  Frames 1, 2, and 3 are labeled a, b c, respectively.
	- The table is sorted by increasing E-value.
********************************************************************************

SEQUENCE NAME                      DESCRIPTION           FRAME   E-VALUE  LENGTH
-------------                      -----------           -----   -------- ------
gi|7296683|gb|AE003602.1|AE003602  Drosophila melanogaster... b    4.6e-23 297266
gi|7295475|gb|AE003567.1|AE003567  Drosophila melanogaster... c    1.3e-12 104808
gi|7289065|gb|AE002569.1|AE002569  Drosophila melanogaster... a        0.3  12850
gi|7289301|gb|AE002567.1|AE002567  Drosophila melanogaster... a        1.1  84373

********************************************************************************



********************************************************************************
SECTION II: MOTIF DIAGRAMS
********************************************************************************
	- The ordering and spacing of all non-overlapping motif occurrences
	  are shown for each high-scoring sequence listed in Section I.
	- A motif occurrence is defined as a position in the sequence whose
	  match to the motif has SEQUENCE p-value less than 0.0001.
	- The SEQUENCE p-value of a match is the probability of
	  some random subsequence in a set of n,
 where n is the sequence length minus the motif width plus 1,
	  scoring at least as well as the observed match.
	- For each sequence, all motif occurrences are shown unless there
	  are overlaps.  In that case, a motif occurrence is shown only if its
	  p-value is less than the product of the p-values of the other
	  (lower-numbered) motif occurrences that it overlaps.
	- The table also shows the E-value of each sequence.
	- Spacers and motif occurences are indicated by
	   o -d-    `d' residues separate the end of the preceding motif 
		    occurrence and the start of the following motif occurrence
	   o [snf]  occurrence of motif `n' with p-value less than 0.0001
		    in frame f.  Frames 1, 2, and 3 are labeled a, b c.
		    A minus sign indicates that the occurrence is on the
		    reverse complement strand.
********************************************************************************

SEQUENCE NAME                      E-VALUE   MOTIF DIAGRAM
-------------                      --------  -------------
gi|7296683|gb|AE003602.1|AE003602   4.6e-23  101665_[+2b]_357_[+1b]_736_[+2c]_
                                             357_[+1c]_752_[-1b]_357_[-2b]_
                                             765_[-1b]_357_[-2b]_191224
gi|7295475|gb|AE003567.1|AE003567   1.3e-12  87158_[-1c]_17563
gi|7289065|gb|AE002569.1|AE002569       0.3  12850
gi|7289301|gb|AE002567.1|AE002567       1.1  84373

********************************************************************************



********************************************************************************
SECTION III: ANNOTATED SEQUENCES
********************************************************************************
	- The positions and p-values of the non-overlapping motif occurrences
	  are shown above the actual sequence for each of the high-scoring
	  sequences from Section I.
	- A motif occurrence is defined as a position in the sequence whose
	  match to the motif has SEQUENCE p-value less than 0.0001 as 
	  defined in Section II.
	- For each sequence, the first line specifies the name of the sequence.
	- The second (and possibly more) lines give a description of the 
	  sequence.
	- Following the description line(s) is a line giving the length, 
	  combined p-value, and E-value of the sequence as defined in Section I.
	- The next line reproduces the motif diagram from Section II.
	- The entire sequence is printed on the following lines.
	- Motif occurrences are indicated directly above their positions in the
	  sequence on lines showing
	   o the motif number/frame of the occurrence (a minus sign indicates that
	  the occurrence is on the reverse complement strand),
	   o the position p-value of the occurrence,
	   o the best possible match to the motif (or its reverse),
	   o columns whose match to the motif has a positive score (indicated 
	     by a plus sign), and
	   o the protein translation of the match (or its reverse).
********************************************************************************


gi|7296683|gb|AE003602.1|AE003602
  Drosophila melanogaster genomic scaffold 142000013386043 section 3 of 8, 
  complete sequence
  LENGTH = 297266  COMBINED P-VALUE = 1.14e-23  E-VALUE =  4.6e-23
  DIAGRAM: 101665_[+2b]_357_[+1b]_736_[+2c]_357_[+1c]_752_[-1b]_357_[-2b]_
           765_[-1b]_357_[-2b]_191224


                                                        [+2b]
                                                        7.5e-13
                                                        K..V..V..L..I..T..G..C..
                                                        +  +  +  +  +  +  +  +
                                                        K..V..V..L..I..T..G..A..
101617 CTTGAGATCCCCTTACATAATTCTGGATCGGATCATGAATTTCGCGGGCAAAGTGGTCCTTATTACGGGAGCA

       
       
       S..S..G..I..G..K..A..T..A..K..H..L..H..K..E..G..A..K..V..V..L..
       +  +  +  +  +     +  +  +  +     +  +  +  +  +  +           +
       S..S..G..I..G..A..A..T..A..I..K..F..A..K..Y..G..A..C..L..A..L..
101690 AGCTCCGGAATCGGAGCTGCAACCGCCATTAAGTTTGCCAAGTACGGCGCCTGTCTGGCTCTCAATGGACGCA

                                                              [+1b]
                                                              1.2e-09
                                                              Y..S..A..S..K..F..
                                                              +        +  +
                                                              Y..N..I..S..K..M..
102055 CATCGTGAACGTGTCCAGTGTCAATGGGATTCGCTCCTTCCCTGGCGTTCTGGCCTACAACATATCCAAAATG

       
       
       A..V..R..M..L..T..R..S..M..A..H..E..Y..A..P..H..G..I..R..V..N..C..I..
       +  +  +     +  +  +  +     +  +  +  +  +     +  +  +  +  +  +  +  +
       G..V..D..Q..F..T..R..C..V..A..L..E..L..A..A..K..G..V..R..V..N..C..V..
102128 GGAGTGGATCAGTTCACCCGCTGTGTGGCGTTGGAGCTGGCTGCCAAGGGTGTGCGCGTGAACTGCGTGAATC

         [+2c]
         3.2e-11
         K..V..V..L..I..T..G..C..S..S..G..I..G..K..A..T..A..K..H..L..H..K..E..G.
         +  +  +  +  +  +  +  +  +  +  +  +  +     +     +  +  +  +     +     +
         K..V..V..L..I..T..G..A..A..S..G..I..G..A..A..A..A..E..M..F..S..K..L..G.
102931 GCAAAGTGGTGCTTATCACGGGCGCAGCCTCCGGGATCGGGGCCGCCGCGGCGGAGATGTTCTCGAAGCTGGG

       
       
       .A..K..V..V..L..
        +           +
       .A..C..L..A..L..
103004 TGCCTGCCTGGCCCTGGTGGATCGGGAGGAGGAGGGCCTCATATGTGTGATGAAACGCTGCATGAAGATGGGC

               [+1c]
               2.6e-13
               Y..S..A..S..K..F..A..V..R..M..L..T..R..S..M..A..H..E..Y..A..P..H.
               +     +  +  +  +  +  +  +     +  +  +  +  +  +  +  +  +  +  +  +
               Y..N..M..S..K..A..A..V..D..Q..F..T..R..S..L..A..L..D..L..G..P..Q.
103369 TGGTGGCCTACAACATGTCCAAGGCGGCGGTGGACCAGTTTACCCGCTCCCTTGCCCTGGATCTGGGTCCCCA

       
       
       .G..I..R..V..N..C..I..
        +  +  +  +  +  +  +
       .G..V..R..V..N..A..V..
103442 GGGTGTTCGGGTGAATGCGGTCAATCCGGGTGTGATTCGCACCAATCTGCAAAAGGCGGGGGGCATGGACGAG

                                                   [-1b]
                                                   1.3e-12
                                                   ..I..C..N..V..R..I..G..H..P..
                                                     +  +  +  +  +  +  +  +  +
                                                   ..V..S..N..V..R..V..G..K..P..
104172 AGTCCACCACGACGCTGCAGCTCGGTGATGATTACGCCGGGATTCACGGAGTTCACACGCACACCCTTGGGAG

       
       
       A..Y..E..H..A..M..S..R..T..L..M..R..V..A..F..K..S..A..S..Y
       +  +  +  +  +     +  +  +  +     +  +  +  +  +  +        +
       A..L..E..L..A..V..C..R..T..F..Q..D..V..A..A..K..S..V..N..Y
104245 CTAGCTCCAGAGCCACGCACCTGGTGAACTGATCCACGGCAGCCTTGGAAACATTGTATGCTAAGACTCCGGG

                                                         [-2b]
                                                         9.6e-09
                                                         ..L..V..V..K..A..G..E..
                                                           +              +
                                                         ..I..T..L..L..G..G..L..
104610 GCCACTATCTGCTCCGCGGTCTCGTTGAGCTTATCCAAATTCCTGCCCACGATGGTGAGCAGGCCTCCCAGTT

       
       
       K..H..L..H..K..A..T..A..K..G..I..G..S..S..C..G..T..I..L..V..V..K
       +  +  +           +        +  +  +  +  +  +  +  +  +  +  +  +  +
       K..A..L..L..V..S..T..G..A..G..I..G..S..S..A..G..T..V..I..I..V..K
104683 TAGCCAAGAGCACCGAAGTACCCGCTCCAATTCCCGAACTGGCTCCGGTCACGATTATAACTTTATCCTTGAA

                                 [-1b]
                                 1.5e-10
                                 ..I..C..N..V..R..I..G..H..P..A..Y..E..H..A..M..
                                   +  +  +  +  +  +  +  +  +  +  +  +  +  +
                                 ..V..A..N..V..R..V..G..K..P..A..L..E..L..A..I..
105486 ATGTCAGTCACAATCACGCCGGGATTCACGGCGTTTACGCGGACTCCTTTGGGGGCCAGTTCCAGGGCTATGC

       
       
       S..R..T..L..M..R..V..A..F..K..S..A..S..Y
       +     +  +     +  +  +  +  +  +        +
       C..A..T..F..Q..D..V..A..A..K..S..V..N..Y
105559 AGGCCGTGAACTGGTCCACCGCTGCCTTAGACACATTGTAGGCCAGGACACCAGGAAAGGCACGCAGTCCACA

                                       [-2b]
                                       4.6e-11
                                       ..L..V..V..K..A..G..E..K..H..L..H..K..A..
                                         +  +           +     +  +  +  +     +
                                       ..I..V..L..L..G..G..L..K..A..L..H..V..A..
105924 GTCTCCTTCAGCTTCTCCTCGTTGCGACCCACGATGACCAGGAGCCCTCCGAGTTTCGCCAAATGGACGGCGG

       
       
       T..A..K..G..I..G..S..S..C..G..T..I..L..V..V..K
          +     +  +  +  +  +  +  +  +  +  +  +  +  +
       A..S..A..G..I..G..S..S..A..G..T..V..I..I..V..K
105997 CACTTGCTCCGATGCCGGAGCTGGCGCCGGTAACAATAATCACTTTGTCCTTGAAGGATGACATCGTGGGGTG


gi|7295475|gb|AE003567.1|AE003567
  Drosophila melanogaster genomic scaffold 142000013386050 section 54 of 54, 
  complete sequence
  LENGTH = 104808  COMBINED P-VALUE = 3.29e-13  E-VALUE =  1.3e-12
  DIAGRAM: 87158_[-1c]_17563


                                                                            [-1c
                                                                            4.6e
                                                                            ..I.
                                                                              +
                                                                            ..L.
87090  AGCCGCTTCATTCGCCGACTCATTCTCGTACAGTGCTTTCGAGAACTTTGTCCTGATGACTCCTGGAGCCAGG

       ]
       -14
       .C..N..V..R..I..G..H..P..A..Y..E..H..A..M..S..R..T..L..M..R..V..A..F..K..
        +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +     +
       .C..N..V..R..I..G..E..P..A..L..D..K..A..A..A..K..T..L..G..I..L..A..T..K..
87163  CAGTTGACGCGAATGCCCTCGGGCGCCAGATCCTTGGCGGCTGCCTTGGTCAAGCCAATCAGCGCGGTCTTGC

       
       
       S..A..S..Y
       +     +  +
       S..V..S..Y
87236  TGACGGAATAGGCTCCCAGTAGCTATGCGATTAAGATAACGGAGATAAGCATTGAACAACTGAACCGCAGATA


gi|7289065|gb|AE002569.1|AE002569
  Drosophila melanogaster genomic scaffold 142000013385354, complete sequence
  LENGTH = 12850  COMBINED P-VALUE = 7.51e-02  E-VALUE =      0.3
  DIAGRAM: 12850



gi|7289301|gb|AE002567.1|AE002567
  Drosophila melanogaster genomic scaffold 142000013385554, complete sequence
  LENGTH = 84373  COMBINED P-VALUE = 2.83e-01  E-VALUE =      1.1
  DIAGRAM: 84373


********************************************************************************


CPU: rocks-155.sdsc.edu
Time 1.550764 secs.

mast meme.adh.oops -stdin -dna -seqp