#=== Query ===#
 
>NM_000032:1-52	sapiens aminolevulinate, delta-, synthase 2 (sideroblastic/hypochromic anemia) (ALAS2), nuclear gene encoding mitochondrial protein, mRNA

    .(((((...((((((.(((((......)))))))).)))...)))))...
1   CACCTGTCATTCGTTCGTCCTCAGTGCAGGGCAACAGGACTTTAGGTTCA 50

    ..
51  AG 52


#=== Hits ===#

Rank Score Query-offset DB Str/seq            Offset            Annotation

__________________________________________________________________________

1     68      1-52     NM_000032:1-52         1-52    sapiens aminolevulinate, delta-, synthase 2 (sider
2     21      10-39    NM_014585:151-250      52-81   sapiens solute carrier family 40 (iron-regulated t
3     19      13-35    NM_003234:3444-3543    41-63   sapiens transferrin receptor (p90, CD71) (TFRC), m
4     17      13-35    NM_003234:3844-3943    44-66   sapiens transferrin receptor (p90, CD71) (TFRC), m
4     17      13-35    NM_003234:3944-4043    56-78   sapiens transferrin receptor (p90, CD71) (TFRC), m
6     16      14-34    NM_000146:1-100        20-40   sapiens ferritin, light polypeptide (FTL), mRNA   
7     13      13-35    NM_003234:3394-3493    41-63   sapiens transferrin receptor (p90, CD71) (TFRC), m

===============================================================

Rank: 1  Score: 68   Query: 52 (ss:20,ds:32)
Identity: str: 100%; seq:100% (ss:100%, ds:100%)
Gap: 0 (ss:0, ds:0)  Mismatch: 0 (ss:0, ds:0)
                   .(((((...((((((.(((((......)))))))).)))...))))).....
                   .(((((...((((((.(((((......)))))))).)))...))))).....
NM_000032:1-52:  1 CACCTGTCATTCGTTCGTCCTCAGTGCAGGGCAACAGGACTTTAGGTTCAAG 52
                   ||||||||||||||||||||||||||||||||||||||||||||||||||||
NM_000032:1-52:  1 CACCTGTCATTCGTTCGTCCTCAGTGCAGGGCAACAGGACTTTAGGTTCAAG 52

===============================================================

Rank: 2  Score: 21   Query: 30 (ss:8,ds:22)
Identity: str: 100%; seq:40% (ss:75%, ds:27%)
Gap: 0 (ss:0, ds:0)  Mismatch: 18 (ss:2, ds:16)
                       ((((((.(((((......)))))))).)))
                       ((((((.(((((......)))))))).)))
NM_000032:1-52:     10 TTCGTTCGTCCTCAGTGCAGGGCAACAGGA 39
                       ::::|||::|::||||| ::|::||: :::
NM_014585:151-250:  52 CAACTTCAGCTACAGTGTTAGCTAAGTTTG 81

===============================================================

Rank: 3  Score: 19   Query: 23 (ss:7,ds:16)
Identity: str: 100%; seq:47% (ss:100%, ds:25%)
Gap: 0 (ss:0, ds:0)  Mismatch: 12 (ss:0, ds:12)
                         (((.(((((......))))))))
                         (((.(((((......))))))))
NM_000032:1-52:       13 GTTCGTCCTCAGTGCAGGGCAAC 35
                         ::|||::::||||||::::||::
NM_003234:3444-3543:  41 TATCGGAAGCAGTGCCTTCCATA 63

===============================================================

Rank: 4  Score: 17   Query: 23 (ss:7,ds:16)
Identity: str: 100%; seq:43% (ss:85%, ds:25%)
Gap: 0 (ss:0, ds:0)  Mismatch: 13 (ss:1, ds:12)
                         (((.(((((......))))))))
                         (((.(((((......))))))))
NM_000032:1-52:       13 GTTCGTCCTCAGTGCAGGGCAAC 35
                         ::|||::::||||| ::::||::
NM_003234:3844-3943:  44 TATCGGGAGCAGTGTCTTCCATA 66

===============================================================

Rank: 4  Score: 17   Query: 23 (ss:7,ds:16)
Identity: str: 100%; seq:43% (ss:85%, ds:25%)
Gap: 0 (ss:0, ds:0)  Mismatch: 13 (ss:1, ds:12)
                         (((.(((((......))))))))
                         (((.(((((......))))))))
NM_000032:1-52:       13 GTTCGTCCTCAGTGCAGGGCAAC 35
                         ::|||::::||||| ::::||::
NM_003234:3944-4043:  56 TATCGGGAACAGTGTTTCCCATA 78

===============================================================

Rank: 6  Score: 16   Query: 21 (ss:7,ds:14)
Identity: str: 100%; seq:47% (ss:85%, ds:28%)
Gap: 0 (ss:0, ds:0)  Mismatch: 11 (ss:1, ds:10)
                     ((.(((((......)))))))
                     ((.(((((......)))))))
NM_000032:1-52:   14 TTCGTCCTCAGTGCAGGGCAA 34
                     ::|:||::||||| ::||:::
NM_000146:1-100:  20 TGCTTCAACAGTGTTTGGACG 40

===============================================================

Rank: 7  Score: 13   Query: 23 (ss:7,ds:16)
Identity: str: 100%; seq:30% (ss:71%, ds:12%)
Gap: 0 (ss:0, ds:0)  Mismatch: 16 (ss:2, ds:14)
                         (((.(((((......))))))))
                         (((.(((((......))))))))
NM_000032:1-52:       13 GTTCGTCCTCAGTGCAGGGCAAC 35
                         ::||:::::||| | :::::|::
NM_003234:3394-3493:  41 TATCAGTGACAGAGTTCACTATA 63



 The score distribution is: 

Max: 68    Min: 13
  score distribution:   one '=' equals 1

< 14.0              |=   [1]
14.0<-->15.0        |
15.0<-->16.0        |=   [1]
16.0<-->17.0        |==   [2]
17.0<-->18.0        |
18.0<-->19.0        |=   [1]
19.0<-->20.0        |
20.0<-->21.0        |=   [1]
21.0<-->22.0        |
22.0<-->23.0        |
23.0<-->24.0        |
24.0<-->25.0        |
25.0<-->26.0        |
26.0<-->27.0        |
27.0<-->28.0        |
28.0<-->29.0        |
29.0<-->30.0        |
30.0<-->31.0        |
31.0<-->32.0        |
32.0<-->33.0        |
33.0<-->34.0        |
34.0<-->35.0        |
35.0<-->36.0        |
36.0<-->37.0        |
37.0<-->38.0        |
38.0<-->39.0        |
39.0<-->40.0        |
40.0<-->41.0        |
41.0<-->42.0        |
42.0<-->43.0        |
43.0<-->44.0        |
44.0<-->45.0        |
45.0<-->46.0        |
46.0<-->47.0        |
47.0<-->48.0        |
48.0<-->49.0        |
49.0<-->50.0        |
50.0<-->51.0        |
51.0<-->52.0        |
52.0<-->53.0        |
53.0<-->54.0        |
54.0<-->55.0        |
55.0<-->56.0        |
56.0<-->57.0        |
57.0<-->58.0        |
58.0<-->59.0        |
59.0<-->60.0        |
60.0<-->61.0        |
61.0<-->62.0        |
62.0<-->63.0        |
63.0<-->64.0        |
64.0<-->65.0        |
65.0<-->66.0        |
66.0<-->67.0        |
>67                 |=   [1]



 The score matrices used are : 

      A  C  G  U 
A     1 -1 -1 -1 
C    -1  1 -1 -1 
G    -1 -1  1 -1 
U    -1 -1 -1  1 


     AU  CG  GC  GU  UA UG 
AU   3   1   1   1   1   1
CG   1   3   1   1   1   1
GC   1   1   3   1   1   1
GU   1   1   1   3   1   1
UA   1   1   1   1   3   1
UG   1   1   1   1   1   3


Gap penalty is set as: -6.0
