MAST - Motif Alignment and Search Tool
MAST version 4.1.1 (Release date: Tue Apr 21 15:00:06 PDT 2009)
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 ./Seed23.fasta (peptide)
Last updated on Fri May 28 14:00:08 2010
Database contains 17 sequences, 2499 residues
MOTIFS ./memout.txt (peptide)
MOTIF WIDTH BEST POSSIBLE MATCH
----- ----- -------------------
1 50 GVTRGQIQLADFLRTVGLRFAIVRGNPYDGKKEGEWIAVALYGTIGAPVK
2 50 WKYCTGKVGSMNSQKVVAAMETAAKSNQVIQTDVYRETHALYHAIMEALY
3 29 MLLQKTHRIGRMAMLLALFDENESPQLSI
PAIRWISE MOTIF CORRELATIONS:
MOTIF 1 2
----- ----- -----
2 0.13
3 0.14 0.11
No overly similar pairs (correlation > 0.60) found.
Random model 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.058
L 0.092 M 0.023 N 0.046 P 0.051 Q 0.041 R 0.052 S 0.074 T 0.059 V 0.064
W 0.013 Y 0.033
SECTION I: HIGH-SCORING SEQUENCES
- Each of the following 17 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
- finding the score of the single best match of each motif
to the sequence (best matches may overlap),
- calculating the sequence p-value of each score,
- forming the product of the p-values,
- 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 table is sorted by increasing E-value.
| Links | Sequence Name | Description | E-value | Length
|
|---|
|
| HUTP_BACAN/8-144
|
| 3.5e-87
| 147
|
|
| HUTP_BACSU/8-146
|
| 1.6e-78
| 147
|
|
| HUTP_BACHD/8-145
|
| 4.2e-78
| 147
|
|
| HUTP_GEOTN/8-147
|
| 6.8e-78
| 147
|
|
| Q1MWD4_HELCL/13-145
|
| 6.9e-07
| 147
|
|
| A1HQA4_9FIRM/15-146
|
| 1.9e-06
| 147
|
|
| Q2IUX8_RHOP2/19-148
|
| 5.6e-06
| 147
|
|
| Q2RUS6_RHORT/7-136
|
| 4.6e-05
| 147
|
|
| Q2RIT2_MOOTA/7-134
|
| 9e-05
| 147
|
|
| Q24V87_DESHY/10-137
|
| 0.00015
| 147
|
|
| A4XKN7_CALS8/10-137
|
| 0.00024
| 147
|
|
| A5D1T0_PELTS/8-135
|
| 0.00052
| 147
|
|
| Q3AAT2_CARHZ/8-135
|
| 0.00054
| 147
|
|
| A6LSD2_CLOB8/7-133
|
| 0.0033
| 147
|
|
| B8FZP6_DESHD/12-143
|
| 0.0075
| 147
|
|
| Q895E7_CLOTE/7-133
|
| 0.011
| 147
|
|
| B8CZC6_HALOH/19-148
|
| 0.066
| 147
|
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 POSITION p-value less than 0.0001.
- 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.
- 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
- Spacers and motif occurences are indicated by
- Spacers and motif occurences are indicated by
- occurrence of motif `n' with p-value less than 0.0001.
- Sequences longer than 1000 are not shown to scale and are indicated by thicker lines.
| Links | Name | Expect | Motifs
|
|---|
|
| HUTP_BACAN/8-144
| 3.5e-87
|
|
|
| HUTP_BACSU/8-146
| 1.6e-78
|
|
|
| HUTP_BACHD/8-145
| 4.2e-78
|
|
|
| HUTP_GEOTN/8-147
| 6.8e-78
|
|
|
| Q1MWD4_HELCL/13-145
| 6.9e-07
|
|
|
| A1HQA4_9FIRM/15-146
| 1.9e-06
|
|
|
| Q2IUX8_RHOP2/19-148
| 5.6e-06
|
|
|
| Q2RUS6_RHORT/7-136
| 4.6e-05
|
|
|
| Q2RIT2_MOOTA/7-134
| 9e-05
|
|
|
| Q24V87_DESHY/10-137
| 0.00015
|
|
|
| A4XKN7_CALS8/10-137
| 0.00024
|
|
|
| A5D1T0_PELTS/8-135
| 0.00052
|
|
|
| Q3AAT2_CARHZ/8-135
| 0.00054
|
|
|
| A6LSD2_CLOB8/7-133
| 0.0033
|
|
|
| B8FZP6_DESHD/12-143
| 0.0075
|
|
|
| Q895E7_CLOTE/7-133
| 0.011
|
|
|
| B8CZC6_HALOH/19-148
| 0.066
|
|
| SCALE
|
| | |
| |
| |
| |
| |
| |
| 1 |
25 |
50 |
75 |
100 |
125 |
|
|---|
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 POSITION 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
- the motif number of the occurrence,
- the position p-value of the occurrence,
- the best possible match to the motif, and
- columns whose match to the motif has a positive score (indicated
by a plus sign).
HUTP_BACAN/8-144
LENGTH = 147 COMBINED P-VALUE = 2.04e-88 E-VALUE = 3.5e-87
DIAGRAM: 35-[2]-[1]-12
[2]
4.0e-37
WKYCTGKVGSMNSQKVVAAMETAAKSNQVIQTDVYRETHA
++ +++++++++++++++++++++++++++++
1 RIGRMAMLLALA---DENESPVLSIPK--GWKYCTGKVGSM-----NSQKVVAAMETAAKSNQVIETDVYRETHA
[1]
3.3e-60
LYHAIMEALYGVTRGQIQLADFLRTVGLRFAIVRGNPYDGKKEGEWIAVALYGTIGAPVK
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
76 LYHAIMEALYGVTRGQIQLADVLRTVGLRFAIVRGTPYDGKKEGEWVAVALYGTIGAPVKGSEHEAIGLGIN
HUTP_BACSU/8-146
LENGTH = 147 COMBINED P-VALUE = 9.12e-80 E-VALUE = 1.6e-78
DIAGRAM: 35-[2]-[1]-12
[2]
3.8e-34
WKYCTGKVGSMNSQKVVAAMETAAKSNQVIQTDVYRETHA
++ +++++++++++++++++++++++++++++
1 RIGRLSVLLLLNE--AEESTQVEELERD-GWKVCLGKVGSM-----DAHKVVAAIETASKKSGVIQSEGYRESHA
[1]
4.3e-54
LYHAIMEALYGVTRGQIQLADFLRTVGLRFAIVRGNPYDGKKEGEWIAVALYGTIGAPVK
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
76 LYHATMEALHGVTRGEMLLGSLLRTVGLRFAVLRGNPYESEAEGDWIAVSLYGTIGAPIKGLEHETFGVGIN
HUTP_BACHD/8-145
LENGTH = 147 COMBINED P-VALUE = 2.47e-79 E-VALUE = 4.2e-78
DIAGRAM: 35-[2]-[1]-12
[2]
2.7e-33
WKYCTGKVGSMNSQKVVAAMETAAKSNQVIQTDVYRETHA
++ +++++++++++++++++++++ ++++ ++
1 RIGQLAILLAVS---EEPGIVEASILSS-DYQSCIGRVGSM-----ESNKIVAAIETAAKKNGFIREGVYREAHA
[1]
9.5e-55
LYHAIMEALYGVTRGQIQLADFLRTVGLRFAIVRGNPYDGKKEGEWIAVALYGTIGAPVK
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
76 LYHAIIEALHGVTRGQVQLGSFQRTVGLRFAIVRGKPYKEEAEGEWIAVALYGTIGAPIKGQEHETIGLGIN
HUTP_GEOTN/8-147
LENGTH = 147 COMBINED P-VALUE = 3.98e-79 E-VALUE = 6.8e-78
DIAGRAM: 35-[2]-[1]-12
[2]
1.4e-32
WKYCTGKVGSMNSQKVVAAMETAAKSNQVIQTDVYRETHA
++ +++++++++++++++++++++++++++++
1 RIGRQALLLAMLDEGEEGAILDELRAS--NWRYCQGRVGAM-----EPQKIVAAIETAAKRHEVVDGSLYRDMHA
[1]
7.5e-55
LYHAIMEALYGVTRGQIQLADFLRTVGLRFAIVRGNPYDGKKEGEWIAVALYGTIGAPVK
+++++++++++++++++++++++++ ++++++++++++++++++++++++++++++++++
76 LYHAILEAVHGVTRGQVELGDLLRTAGLRFAVVRGTPYEQPKEGEWIAVALYGTIGAPVRGLEHEAVGLGIN
Q1MWD4_HELCL/13-145
LENGTH = 147 COMBINED P-VALUE = 4.04e-08 E-VALUE = 6.9e-07
DIAGRAM: 35-[2]-[1]-12
[2]
7.7e-10
WKYCTGKVGSMNSQKVVAAMETAAKSNQVIQTDVYRETHA
+ + + ++ +++ + ++++ ++
1 SLERAAMMLAMSETRESEFRLRQYYGETWHLPCVVTEVGGTIGALQTSGKLTNSVVTAALKAKVIESK-PENVHA
[1]
4.1e-05
LYHAIMEALYGVTRGQIQLADFLRTVGLRFAIVRGNPYDGKKEGEWIAVALYGTIGAPVK
+++++++ + + +++ ++ + + + ++ +
76 VIHATIEALKGIFLDYAQNASL----ALKISVIS--------DGKWLAVGIFGK-SSVHPLTEHCRVGLGFM
A1HQA4_9FIRM/15-146
LENGTH = 147 COMBINED P-VALUE = 1.13e-07 E-VALUE = 1.9e-06
DIAGRAM: 35-[2]-[1]-12
[2]
1.2e-06
WKYCTGKVGSMNSQKVVAAMETAAKSNQVIQTDVYRETHA
+ + ++ ++ ++ + + ++
1 SVATAAMLLALTRTISDEETVKKMLADK-GYRFVVTEVGGKSSLTEFQEKTTKAVIGAALNSGIIGKS-ATNYHA
[1]
2.8e-07
LYHAIMEALYGVTRGQIQLADFLRTVGLRFAIVRGNPYDGKKEGEWIAVALYGTIGAPVK
+ +++ ++ + ++ + ++++ ++++++ + +
76 LLHATDEAKRGIMVNVASSVSL----AVKIAIVR--------DDSWIAVALFGD-SALHPLTNHERAGLGVM
Q2IUX8_RHOP2/19-148
LENGTH = 147 COMBINED P-VALUE = 3.29e-07 E-VALUE = 5.6e-06
DIAGRAM: 35-[2]-[1]-12
[2]
5.3e-09
WKYCTGKVGSMNSQKVVAAMETAAKSNQVIQTDVYRETHA
+ + ++ ++ +++ ++ +++
1 RIGKVAIMLAMT-SEEEEASFKAKISKETNYKFSVTWVSGRK--TEINKTFPRSILNAALKHGTIEDR-SGATHA
[1]
3.6e-05
LYHAIMEALYGVTRGQIQLADFLRTVGLRFAIVRGNPYDGKKEGEWIAVALYGTIGAPVK
+ ++ ++++ + ++ + + + ++ ++ ++ +
76 LIHAGLEALSGLIPTVPGDSSL----KVKVVIIA--------DGHWVVVAAYGE-SAFLPETNHERIGLGVM
Q2RUS6_RHORT/7-136
LENGTH = 147 COMBINED P-VALUE = 2.68e-06 E-VALUE = 4.6e-05
DIAGRAM: 35-[2]-[1]-12
[2]
2.6e-05
WKYCTGKVGSMNSQKVVAAMETAAKSNQVIQTDVYRETHA
+ + + ++ + +++ + +
1 RLGKVAIIAAIS-SREEEEEFKAYLAENSKIKLAVTFVSGVR--TDVTKSFVKSIVACALQGHVIRHQ-HNELHG
[1]
6.5e-08
LYHAIMEALYGVTRGQIQLADFLRTVGLRFAIVRGNPYDGKKEGEWIAVALYGTIGAPVK
++ ++ + + +++ ++ +++ + +++++ + +
76 VIHAGLECLKGISADVAAEASL----KLKVAVVN--------DGKWVAVAAHGE-SAFHPETNHERMGFGVM
Q2RIT2_MOOTA/7-134
LENGTH = 147 COMBINED P-VALUE = 5.27e-06 E-VALUE = 9e-05
DIAGRAM: 35-[2]-62
[2]
1.1e-07
WKYCTGKVGSMNSQKVVAAMETAAKSNQVIQTDVYRETHA
+ +++ ++++ +++ + + +
1 KVAAVAIKMALSESREEEARLKKEYAAG-GISTAAVDYGGEF--INSVQKIIERAVVAAKREGVIKEA-HQEEGA
LYHAIMEALY
++ +++
76 IAGATREALSQIMPKAIGLNV-----GGKIGIAR--------WHDHLSVAIFFGIGLLH--LDEVGIGMGHR
Q24V87_DESHY/10-137
LENGTH = 147 COMBINED P-VALUE = 8.89e-06 E-VALUE = 0.00015
DIAGRAM: 35-[2]-62
[2]
2.6e-08
WKYCTGKVGSMNSQKVVAAMETAAKSNQVIQTDVYRETHA
+ +++ ++++ +++++ + +
1 KIAMAAMRMAMAESREEEVHLKDGFGRH-GIRTVAVDYGGEY--LTAIKKIVERAIVAAKREGVIQET-HGDEGA
LYHAIMEALY
++ +++
76 VAGATREALSQIMPKAMGLNI-----GGKIGIAR--------RDEHLSVAVFFGVGLLH--LDEVAVGLGHR
A4XKN7_CALS8/10-137
LENGTH = 147 COMBINED P-VALUE = 1.40e-05 E-VALUE = 0.00024
DIAGRAM: 35-[2]-[1]-12
[2]
2.1e-05
WKYCTGKVGSMNSQKVVAAMETAAKSNQVIQTDVYRETHA
++ ++++ +++ + + +
1 EVSRAAILIALSQSRQEEKKIQDDFGKI-GIRCAAVDFGGEF--ITSVMKIIERAVVAAKREGVIAEV-HQEEGA
[1]
3.4e-05
LYHAIMEALYGVTRGQIQLADFLRTVGLRFAIVRGNPYDGKKEGEWIAVALYGTIGAPVK
++ ++ + + + + + + + ++++ ++
76 VAGATREAISQIMQKAVGLNV-----GGKIGIAR--------YNEHVAVAIFFGIGLLH--LNEVAIGLGHR
A5D1T0_PELTS/8-135
LENGTH = 147 COMBINED P-VALUE = 3.06e-05 E-VALUE = 0.00052
DIAGRAM: 35-[2]-62
[2]
2.4e-07
WKYCTGKVGSMNSQKVVAAMETAAKSNQVIQTDVYRETHA
+ ++ ++++ +++++ + +
1 KVARVAIEMALTDSREQEKEYKQRFARE-GIKTAAVDYGGDF--ISSVGKIIERAIVASKREGVIREV-HADEGA
LYHAIMEALY
++ +++
76 VAGATREALSQIMPKAIGLNI-----GGKIGIAR--------KDDHISVAVFFGVGLLH--LDEVAIGLGHR
Q3AAT2_CARHZ/8-135
LENGTH = 147 COMBINED P-VALUE = 3.19e-05 E-VALUE = 0.00054
DIAGRAM: 35-[2]-62
[2]
1.7e-07
WKYCTGKVGSMNSQKVVAAMETAAKSNQVIQTDVYRETHA
+ ++++ ++++ +++ + + +
1 KAAKAALMMALTETREEEKELKIKLAEE-NIRAAAVDYGGDF--ISAVNKIVERAVVAAKREGVIKET-HADEGA
LYHAIMEALY
++ ++
76 VAGATREAISQLMTKALGLNV-----GGKIGIAR--------FQDHVAVAVFFGVGLLH--LDEVGVGLGHR
A6LSD2_CLOB8/7-133
LENGTH = 147 COMBINED P-VALUE = 1.96e-04 E-VALUE = 0.0033
DIAGRAM: 35-[2]-62
[2]
5.7e-07
WKYCTGKVGSMNSQKVVAAMETAAKSNQVIQTDVYRETHA
+ ++ + ++++++ ++++ +
1 RVAKIATKMAIC-DRHEEEHLKKVYAEK-GIKVTAVNVGGNI--NSSIAKILESALVAAKRNELIREE-HLHEGA
LYHAIMEALY
++ ++
76 VIGATRDAVIQVANRANGQNI-----GGKIGIAR--------GGEHISVCIFLSIGLLH--LDEVAIGIGHR
B8FZP6_DESHD/12-143
LENGTH = 147 COMBINED P-VALUE = 4.42e-04 E-VALUE = 0.0075
DIAGRAM: 35-[2]-62
[2]
7.8e-06
WKYCTGKVGSMNSQKVVAAMETAAKSNQVIQTDVYRETHA
+ ++ ++ + ++ + ++
1 SLERTALALALTETREKENELKKLMLSQ-NVYCAVTELGGTYSALQPTGKLTHSVISAAINTGAIQKE-PKAIHA
LYHAIMEALY
++++++
76 VVHATLEASNGIFVHTNSNASY----ALKVGIAC--------DRDWIAIAIFGR-SSIHALSEHCRVGLGYM
Q895E7_CLOTE/7-133
LENGTH = 147 COMBINED P-VALUE = 6.64e-04 E-VALUE = 0.011
DIAGRAM: 35-[2]-62
[2]
2.2e-05
WKYCTGKVGSMNSQKVVAAMETAAKSNQVIQTDVYRETHA
+ +++ ++++ + ++ +
1 DVAKASVRMAIS-SREEEKLLESELKKE-EILTVAVDVGGNL--NNSIPKIVERALVASKRCGLIQDC-HLHDGA
LYHAIMEALY
++ ++
76 VAGATREAIIQVVSKANGLNV-----GGKIGIAR--------SKEHLSVCLFMSIGLLH--LNEVVIGLGHR
B8CZC6_HALOH/19-148
LENGTH = 147 COMBINED P-VALUE = 3.89e-03 E-VALUE = 0.066
DIAGRAM: 35-[2]-62
[2]
5.7e-05
WKYCTGKVGSMNSQKVVAAMETAAKSNQVIQTDVYRETHA
+ +++ ++ + +++ ++ +
1 SLGKVATILAIVDDQYDDDIIESLKKE--GYKAVITRAGGKG--EELKNKILRNCLGAAIKGGVITES-IQDQRV
LYHAIMEALY
+ + +
76 LTRCVERAMVDFNTPMLSISGA----GIKIGIAS--------KGPHLAVALYGKLGIPGLDVDHEISGMGIH
Debugging Information
CPU: kodomo-count
Time 0.116007 secs.
mast ./memout.txt -d ./Seed23.fasta -ev 10.000000 -mt 0.000100
Button Help
Links to Entrez database at NCBI
Links to sequence scores (section I)
Links to motif diagrams (section II)
Links to sequence/motif annotated alignments (section III)
This information