#=== Query ===# >NM_000032:1-52 Homo 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 Homo sapiens aminolevulinate, delta-, synthase 2 ( 2 47 2-52 NM_000146:1-100 1-76 Homo sapiens ferritin, light polypeptide (FTL), mR 3 45 1-52 NM_014585:151-250 2-97 Homo sapiens solute carrier family 40 (iron-regula 4 41 1-51 NM_003234:3844-3943 1-81 Homo sapiens transferrin receptor (p90, CD71) (TFR 5 37 2-47 NM_003234:3394-3493 1-100 Homo sapiens transferrin receptor (p90, CD71) (TFR 5 37 2-51 NM_003234:3444-3543 26-99 Homo sapiens transferrin receptor (p90, CD71) (TFR 5 37 2-49 NM_003234:3944-4043 2-100 Homo sapiens transferrin receptor (p90, CD71) (TFR =============================================================== 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: 47 Query: 43 (ss:13,ds:30) Identity: str: 45%; seq:74% (ss:76%, ds:73%) Gap: 36 (ss:20, ds:16) Mismatch: 8 (ss:0, ds:8) (( ( ((...((( (( (. ( ((((...... )))) ) ))). )))...)) ))) .. . . . (((((.(((( ( ((((.( (.((.(((((..... .))))))))))) )) ) ))))))))).................... NM_000032:1-52: 2 ---AC-C-TGTCATTC-GT-TC---G--TCCTCAGTGC-AGGG-C--AACA-GGACTTTA-GGT---TC------A-------A--G 52 :| | :| | | || | | ||::||||| ::|| | ||| | | |: ||: || | | | NM_000146:1-100: 1 GTCCCGCGGG---T-CTGTCT-CTTGCTTCAACAGTG-TTTGGACGGAAC-AG-A---TCCGGGGACTCTCTTCCAGCCTCCGACCG 76 =============================================================== Rank: 3 Score: 45 Query: 45 (ss:18,ds:27) Identity: str: 33%; seq:73% (ss:72%, ds:74%) Gap: 59 (ss:29, ds:30) Mismatch: 4 (ss:0, ds:4) .( (((( . . . ((( (((. ((((( ...... ))))) ))). )))...) )))).... . . ......((( ((.((((((.................)))))).))(((((((((((.(( (((..... .))) ))))) .))))) ))).))) ....... NM_000032:1-52: 1 CA------CCTG-------------------T-C-A--------------TTC--GTTC-GTCCT-CAGTGC--AGGGC-AACA---GGACTTT----AGGTTCAA----G 52 | ||| | | | |:| :||| | || ||||| || | ||: |:| ||| || | NM_014585:151-250: 2 C-GGACGACCT-GCTGAGCCTCCCAAACCGCTTCCATAAGGCTTTGCCTTTCCAACTTCAG--CTACAGTG-TTAG--CTAAG-TTTGGA----AAGAAGG---AAAAAAG 97 =============================================================== Rank: 4 Score: 41 Query: 51 (ss:19,ds:32) Identity: str: 45%; seq:52% (ss:68%, ds:43%) Gap: 40 (ss:32, ds:8) Mismatch: 14 (ss:0, ds:14) .( ( ((( .. .( (( (( (.(((((...... )))))) )). )))... ))) )). ... . ....................((((((.......... ..(((((((.(((((..... .))))))))) ))) )))))) ..... NM_000032:1-52: 1 CA--------------------C-CTG------TC--AT--TC-GT-TCGTCCTCAGTGC-AGGGCA-ACA-GGACTT-TAG-GTT--CAA 51 | | :|| || | :| :| |||::::||||| ::::|| |: |: ||: | ||| NM_003234:3844-3943: 1 C-ATTTGCAAAATTTCCAGTACCTTTGTCACAATCCTA-ACACATTATCGGGAGCAGTG-TCTTCCATAA-TGT----ATAAAG--AACAA 81 =============================================================== Rank: 5 Score: 37 Query: 39 (ss:13,ds:26) Identity: str: 30%; seq:61% (ss:61%, ds:61%) Gap: 66 (ss:26, ds:40) Mismatch: 10 (ss:0, ds:10) (((( ( . . . ( ( ((( (. ((((( ... ... ))))) ) )).) ) )... ) )))) ((((((.((((....)).))...((((.........(((((((.(( (((...... ))) )))))) )))(((((((......)))))))...)) )))))))) NM_000032:1-52: 2 ACCT-G-----T-----C----A--T----------T-CGT-TC-GTCCT--CAG--TGC--AGGGC-A-ACAG-G------------------------ACTT--T-AGGT 47 :|:| : | | | | | ::| || | | ||| | | | | |: : | | : |:|: NM_003234:3394-3493: 1 GCTTTCTGTCCTTTTGGCACTGAGATATTTATTGTTTATTTATCAG---TGACAGAGT--TCA---CTATAA-ATGGTGTTTTTTTAATAGAATATAATTA---TCGGAAGC 100 =============================================================== Rank: 5 Score: 37 Query: 46 (ss:16,ds:30) Identity: str: 49%; seq:50% (ss:81%, ds:33%) Gap: 33 (ss:7, ds:26) Mismatch: 18 (ss:0, ds:18) (((( (... ((((( (.(((((......)))))) )).))). ..) )))) .... ((((. (((((((((((((.(((((......))))))))) ))) )(((((.....))))).)))))..)))) .. .. NM_000032:1-52: 2 ACCT-GTCA------TTCGT-TCGTCCTCAGTGCAGGGCA-ACAGGAC----------------------TTT-AGGT-TCAA 51 :| :| |:::| |||::::||||||::::|| |: ::| ||: |: | || NM_003234:3444-3543: 26 --TTAAT--AGAATATAATTATCGGAAGCAGTGCCTTCCATAA-TTA-TGACAGTTATACTGTCGGTTTTTTTTAA--AT-AA 99 =============================================================== Rank: 5 Score: 37 Query: 47 (ss:15,ds:32) Identity: str: 40%; seq:46% (ss:80%, ds:31%) Gap: 57 (ss:23, ds:34) Mismatch: 20 (ss:0, ds:20) (((( (. ..( (((( (.(((((...... )))))) )).))) ... ) ) ))). . (((((((((( .(((((.. ....))))).....((((((.....))))))..((((((.(((((..... .))))))))) )) ......)))))..))) )). NM_000032:1-52: 2 --ACCT---GT------CAT---------------------------------TCGT-TCGTCCTCAGTGC-AGGGCA-ACAGGA--CTT-T---A--GGTT--C 49 ::|| : || :::| |||::::||||| ::::|| |: :: ||| : | |:: | NM_003234:3944-4043: 2 GTGTCTGTAT-CGGAGACA-GTGATCTCCATATGTTACACTAAGGGTGTAAGTAATTATCGGGAACAGTG-TTTCCCATAA-TT-TTCTTCATGCAATGAC-ATC 100 The score distribution is: Max: 68 Min: 37 score distribution: one '=' equals 1 < 38.0 |=== [3] 38.0<-->39.0 | 39.0<-->40.0 | 40.0<-->41.0 |= [1] 41.0<-->42.0 | 42.0<-->43.0 | 43.0<-->44.0 | 44.0<-->45.0 |= [1] 45.0<-->46.0 | 46.0<-->47.0 |= [1] 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: -2.0