#

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

If you use MAST in your research, please cite the following paper:
Timothy L. Bailey and Michael Gribskov, "Combining evidence using p-values: application to sequence homology searches", Bioinformatics, 14(1):48-54, 1998.

Inputs | Search Results | Program information | Documentation

## Inputs

Top

#### Sequence Databases

The following sequence database was supplied to MAST.

Database | Sequence Count | Residue Count | Last Modified
--- | --- | --- | ---
adh.s | 33 | 9996 | Wed Jan 11 17:17:32 2012
Total | 33 | 9996 |

#### Motifs

The following motifs were supplied to MAST from "meme.adh.zoops" last modified on Wed Jan 11 17:17:55 2012.

|  |  | Similarity
--- | --- | --- | ---
Motif | Width | Best possible match | 1 | 2
1 | 21 | G K V V L I T G C S S G I G K A T A K H L | - | 0.29
2 | 29 | S V Y C A S K F A V R M L T R S M A M E Y A P H G I R V N | 0.29 | -

## Search Results

Top

#### Top Scoring Sequences

Each of the following 33 sequences has an E-value less than 10.
The motif matches shown have a position p-value less than 0.0001.
Click on the arrow (↧) next to the E-value to view more information about a sequence.

| Motif 1

| Motif 2

Sequence | E-value |  | Block Diagram
--- | --- | --- | ---
| 0 100 200 300 400 500 600 700 800 900
YRTP_BACSU | 9e-44 | ↧ |
BUDC_KLETE | 1.1e-42 | ↧ |
FIXR_BRAJA | 3.8e-37 | ↧ |
AP27_MOUSE | 2.5e-35 | ↧ |
NODG_RHIME | 3.2e-34 | ↧ |
HDHA_ECOLI | 4.1e-34 | ↧ |
DHGB_BACME | 4.5e-34 | ↧ |
DHB2_HUMAN | 1.2e-32 | ↧ |
HDE_CANTR | 1.3e-32 | ↧ |
DHII_HUMAN | 5e-32 | ↧ |
DHMA_FLAS1 | 5.4e-32 | ↧ |
FVT1_HUMAN | 8e-32 | ↧ |
RIDH_KLEAE | 9.1e-32 | ↧ |
YINL_LISMO | 1.1e-31 | ↧ |
2BHD_STREX | 2.4e-31 | ↧ |
HMTR_LEIMA | 5.5e-31 | ↧ |
ENTA_ECOLI | 3.6e-30 | ↧ |
BDH_HUMAN | 4e-30 | ↧ |
GUTD_ECOLI | 2.2e-29 | ↧ |
DHES_HUMAN | 3.6e-29 | ↧ |
BA72_EUBSP | 1.9e-28 | ↧ |
3BHD_COMTE | 3.1e-28 | ↧ |
DHB3_HUMAN | 2.6e-27 | ↧ |
RFBB_NEIGO | 4.3e-25 | ↧ |
LIGD_PSEPA | 4.3e-25 | ↧ |
BPHB_PSEPS | 1.3e-24 | ↧ |
DHCA_HUMAN | 1.5e-22 | ↧ |
ADH_DROME | 5e-20 | ↧ |
YURA_MYXXA | 8.4e-20 | ↧ |
FABI_ECOLI | 1.8e-19 | ↧ |
PCR_PEA | 2e-14 | ↧ |
CSGA_MYXXA | 3.9e-13 | ↧ |
MAS1_AGRRA | 4.2e-13 | ↧ |

| Motif 1

| Motif 2

Top

##### MAST version

4.8.0 (Release date: Wed Jan 11 15:32:51 EST 2012)

##### Reference

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

##### Command line summary

mast meme.adh.zoops adh.s
Background letter frequencies (from non-redundant database):
A: 0.073 C: 0.018 D: 0.052 E: 0.062 F: 0.040 G: 0.069 H: 0.022 I: 0.056 K: 0.092 L: 0.023 M: 0.046 N: 0.051 P: 0.052 Q: 0.074 R: 0.059 S: 0.000 T: 0.064 V: 0.001 W: 0.033 Y: 0.000
Result calculation took 1.043 seconds
show model parameters...

##### Model parameters

max_correlation = 0.60 remove_correlated = value: "n" strand_handling = value: "protein" translate_dna = value: "n" max_seq_evalue = 10 adj_hit_pvalue = value: "n" max_hit_pvalue = 0.0001 max_weak_pvalue = 0.0001 host = tlb-takumi-lt.imb.uq.edu.au when = Wed Jan 11 17:18:05 2012
hide model parameters...

## Explanation of MAST Results

Top

#### The MAST results consist of

- The inputs to MAST including:
  1. The sequence databases showing the sequence and residue counts. [View]
  2. The motifs showing the name, width, best scoring match and similarity to other motifs. [View]
  3. The nominal order and spacing diagram.
- The search results showing top scoring sequences with tiling of all of the motifs matches shown for each of the sequences. [View]
- The program details including:
  1. The version of MAST and the date it was released. [View]
  2. The reference to cite if you use MAST in your research. [View]
  3. The command line summary detailing the parameters with which you ran MAST. [View]
- This explanation of how to interpret MAST results.

#### Inputs

MAST received the following inputs.

##### Sequence Databases

This table summarises the sequence databases specified to MAST.

Database
  The name of the database file.
Sequence Count
  The number of sequences in the database.
Residue Count
  The number of residues in the database.

##### Motifs

Summary of the motifs specified to MAST.

Name
  The name of the motif. If the motif has been removed or removal is recommended to avoid highly similar motifs then it will be displayed in red text.
Width
  The width of the motif. No gaps are allowed in motifs supplied to MAST as it only works for motifs of a fixed width.
Best possible match
  The sequence that would achieve the best possible match score and its reverse complement for nucleotide motifs.
Similarity
  MAST computes the pairwise correlations between each pair of motifs. The correlation between two motifs is the maximum sum of Pearson's correlation coefficients for aligned columns divided by the width of the shorter motif. The maximum is found by trying all alignments of the two motifs. Motifs with correlations below 0.60 have little effect on the accuracy of the combined scores. Pairs of motifs with higher correlations should be removed from the query. Correlations above the supplied threshold are shown in red text.

##### Nominal Order and Spacing

This diagram shows the normal spacing of the motifs specified to MAST.

#### Search Results

MAST provides the following motif search results.

##### Top Scoring Sequences

This table summarises the top scoring sequences with a Sequence E-value better than the threshold (default 10). The sequences are sorted by the Sequence E-value from most to least significant.

Sequence
  The name of the sequence. This maybe be linked to search a sequence database for the sequence name.
E-value
  The E-value of the sequence. For DNA only; if strands were scored seperately then there will be two E-values for the sequence seperated by a "/". The score for the provided sequence will be first and the score for the reverse-complement will be second.
↧
  Click on this to show additional information about the sequence such as a description, combined p-value and the annotated sequence.
Block Diagram
  The block diagram shows the best non-overlapping tiling of motif matches on the sequence.
  - The length of the line shows the length of a sequence relative to all the other sequences.
  - A block is shown where the positional p-value of a motif is less (more significant) than the significance threshold which is 0.0001 by default.
  - If a significant motif match (as specified above) overlaps other significant motif matches then it is only displayed as a block if its positional p-value is less (more significant) then the product of the positional p-values of the significant matches that it overlaps.
  - The position of a block shows where a motif has matched the sequence.
  - The width of a block shows the width of the motif relative to the length of the sequence.
  - The colour and border of a block identifies the matching motif as in the legend.
  - The height of a block gives an indication of the significance of the match as taller blocks are more significant. The height is calculated to be proportional to the negative logarithm of the positional p-value, truncated at the height for a p-value of 1e-10.
  - Hovering the mouse cursor over the block causes the display of the motif name and other details in the hovering text.
  - DNA only; blocks displayed above the line are a match on the given DNA, whereas blocks displayed below the line are matches to the reverse-complement of the given DNA.
  - DNA only; when strands are scored separately then blocks may overlap on opposing strands.

##### Additional Sequence Information

Clicking on the ↧ link expands a box below the sequence with additional information and adds two dragable buttons below the block diagram.

Description
  The description appearing after the identifier in the fasta file used to specify the sequence.
Combined p-value
  The combined p-value of the sequence. DNA only; if strands were scored seperately then there will be two p-values for the sequence seperated by a "/". The score for the provided sequence will be first and the score for the reverse-complement will be second.
Annotated Sequence
  The annotated sequence shows a portion of the sequence with the matching motif sequences displayed above. The displayed portion of the sequence can be modified by sliding the two buttons below the sequence block diagram so that the portion you want to see is between the two needles attached to the buttons. By default the two buttons move together but you can drag one individually by holding shift before you start the drag. If the strands were scored seperately then they can't be both displayed at once due to overlaps and so a radio button offers the choice of strand to display.

#### Scoring

MAST scores sequences using the following measures.

##### Position score calculation

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.

##### Position p-value

The position p-value of a match is the probability of a single random subsequence of the length of the motif scoring at least as well as the observed match.

##### Sequence p-value

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.

##### Combined p-value

The combined p-value of a sequence measures the strength of the match of the sequence to all the motifs and is calculated by

1. finding the score of the single best match of each motif to the sequence (best matches may overlap),
2. calculating the sequence p-value of each score,
3. forming the product of the p-values,
4. taking the p-value of the product.

##### Sequence E-value

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.

ACGTN