1.1 TAAATT
1.2 TAGAAT
1.3 TGTAAT
The weight matrix calculated
from the frequency table, compiled by Hawley&McClure 1983 of the three
sequences are
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
TAGAAT is 10 2.18 = 151.36 times more likely to be a – 10 promoter of E.coli than a random sequence2. Use the dot matrix method to identify the similarity between these two sequences.
2.1 ATGCTACGGGTAATC
2.2 GGTAATCATGCTACGG
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|||||||||||
|
|
|
|
|
|||||||||||
|
|
|
|
|
|||||||||||
|
|
|
|
|
|||||||||||
|
|
|
|
|
|||||||||||
|
|
|
|||||||||||||
|
|
|
|
|
|||||||||||
|
|
|
|
|
|||||||||||
|
|
|
|
|
|||||||||||
|
|
|
|
||||||||||||
|
|
|
|
||||||||||||
|
|
|
|
|
|||||||||||
|
|
|
|
||||||||||||
|
|
|
|
|
|||||||||||
|
|
|
|
|
There is also a DNA rearrangement (a line perpendicular to the diagonal line at 11th to 14th of the first sequence)
3. Submit the following sequences to Genbank and explain the results.
3.1 agtcgaacgg aaaggtctct
3.2 agtcgaacgg aaaggtctct tcgagtggcg aacgggtgag taacacgtgg
When submit the following sequences to Genbank the results are
3.1 agtcgaacgg aaaggtctct
Query= prasit
(20 letters)
Database: nt
591,775 sequences; 1,603,993,870 total letters
If you have any problems or questions with the results of this search please refer to the BLAST FAQs
Distribution of 100 Blast Hits
on the Query Sequence
Score E
Sequences producing significant alignments: (bits) Value N
emb|X58888.1|ML16SRRN M.leprae
gene for 16S ribosomal RNA 40 0.008 1
gb|M29575.1|MSGRR16SQ M.kansasii
16S ribosomal RNA 40 0.008 1
emb|AJ007315.1|MCA7315
Mycobacterium canettii 16S rRNA ... 40 0.008 1
emb|Z83862.1|MTCY149 Mycobacterium
tuberculosis H37Rv c... 40 0.008 1
gb|AF152559.1|AF152559 Mycobacterium
sp. 'MCRO 33' 16S ... 40 0.008 1
emb|Z13990.1|MU16SRRN M.ulcerans
16S ribosomal RNA 40 0.008 1
emb|X88923.1|MHR16SRNA M.haemophilum
16S rRNA gene 40 0.008 1
emb|X52919.1|MGA16SRN Mycobacterium
gastri 16S rRNA gene 40 0.008 1
gb|U06638.1|MHU06638 Mycobacterium
haemophilum 16S ribo... 40 0.008 1
3.2 agtcgaacgg aaaggtctct tcgagtggcg aacgggtgag taacacgtgg
Query= prasit2
(50 letters)
Database: Non-redundant GenBank+EMBL+DDBJ+PDB sequences527,570 sequences; 1,523,322,523 total letters
Score E
Sequences producing significant alignments: (bits) Value emb|X70960.1|MH16SRRNA M.heidelbergense 16S rRNA 60 4e-08 gb|AF152560.1|AF152560 Mycobacterium malmoense 16S ribosoma... 60 4e-08 gb|AF152559.1|AF152559 Mycobacterium sp. 'MCRO 33' 16S ribo... 60 4e-08 gb|AF115940.1|AF115940 Uncultured Corynebacterium sp. MTcor... 60 4e-08 emb|AJ131120.1|MTU131120 Mycobacterium tuberculosis 16S rRN... 60 4e-08 emb|AJ007315.1|MCA7315 Mycobacterium canettii 16S rRNA gene... 60 4e-08 emb|Z83862|MTCY149 Mycobacterium tuberculosis H37Rv complet... 60 4e-08 gb|AF059853|AF059853 Mycobacterium avium strain ATCC25291 1... 60 4e-08 gb|AF059851|AF059851 Mycobacterium fortuitum strain ATCC684... 60 4e-08
The short sequence result = emb|AJ007315.1|MCA7315 Mycobacterium canettii 16S rRNA ... 40 0.008 1
The long sequence result =emb|AJ007315.1|MCA7315 Mycobacterium canettii 16S rRNA gene... 60 4e-08The result of submitted sequences show the effect of query sequence length. When blast with the same sequence the longer sequence show higher total score and lower probability to be randomly similar (without homologous) to the database sequence. Thus blast result of longer sequence give more probability to be significantly homologous than the shorter one.