Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_hsa_tpiA 4,604,698 C→A A186E (GCA→GAA)  bglJ → bgl operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA4,604,6980CA100.0% 128.6 / NA 41A186E (GCA→GAA) bglJbgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (21/20);  total (21/20)

ACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCG  >  NC_000913_3_hsa_tpiA/4604557‑4604833
                                                                                                                                             |                                                                                                                                       
acagccTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATaaa                                                                                                                                  <  1:543334/144‑1 (MQ=255)
    acTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCCATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTa                                                                                                                              >  2:524528/1‑149 (MQ=255)
    acTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTa                                                                                                                              >  2:37623/1‑149 (MQ=255)
                   gTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGATAGAAAGATTGCGCTTAATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAAc                                                                                                               >  1:486751/1‑149 (MQ=255)
                   gTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGATAc                                                                                                               >  2:288973/1‑149 (MQ=255)
                    tACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATCTAATGTGATGTCGAAACt                                                                                                              >  2:490184/1‑149 (MQ=255)
                       aTTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTggg                                                                                                           >  1:108580/1‑149 (MQ=255)
                       aTTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTggg                                                                                                           >  1:582173/1‑149 (MQ=255)
                                    ggACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGAc                                                                                              >  1:173882/1‑149 (MQ=255)
                                         ttAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAgg                                                                                         <  2:14923/149‑1 (MQ=255)
                                                  aCGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGa                                                                                <  2:108580/149‑1 (MQ=255)
                                                  aCGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGa                                                                                <  1:301877/149‑1 (MQ=255)
                                                  aCGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGa                                                                                <  1:37623/149‑1 (MQ=255)
                                                          agaAATATTGCGCTTTATGTCACGTGGCTACTCAATTACACAAATTGCCGATCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTAGAAACTGGGCGTCAGTTCTGACGCAGGATTGTTTGAGGCCGCAg                                                                        >  2:253213/1‑149 (MQ=255)
                                                               tatTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATAtt                                                                   >  2:377202/1‑149 (MQ=255)
                                                                atTGCGCTTTATGTCGCGTGGCTACTCAATTACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTc                                                                  <  1:264953/149‑1 (MQ=255)
                                                                 ttGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCt                                                                 >  1:241796/1‑149 (MQ=255)
                                                                 ttGCGCTTTATGTCGCGTGGCTACGCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCt                                                                 <  2:34537/149‑1 (MQ=255)
                                                                   gcgcTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGt                                                               <  2:486751/149‑1 (MQ=255)
                                                                   gcgcTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATAGCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGt                                                               <  1:27399/149‑1 (MQ=255)
                                                                       tttATATCGCGTGGCTCCTCAATGACACGAATTGCCGAGCAGCTTAAACGCAATGTCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGCGCGTCAGTCCTGACGCAGGGTCGTTGGAGGCCGCAGATATTCTCTcatg                                                           >  2:435645/1‑149 (MQ=255)
                                                                           tGTCGGGTGGCTACGCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTtatgtatg                                                       <  2:83187/149‑1 (MQ=255)
                                                                            gTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGc                                                      <  1:574595/149‑1 (MQ=255)
                                                                                  tGGCGAATCAATGACACAAAGTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAAGTTAATGTGATGTCGAAACTGGGCGTCAGTTCCGACGCAGGGTGGTTGGTGGCCGCTGATATTCTGTTATGTATGCGGCAtt                                                <  2:275580/149‑1 (MQ=255)
                                                                                      tACTCAATGGCACAAATTGCCGAGCAGCTTAAACGGAAAATGAAAACGATCCGTGAACATAAATTTAAAGTGATGTCGAAACTGGGCGCCAGGTCTAACGCAGCGTGAGTGGAGGCCGCGGATTTTGTGTTATGTATTCGGCAtcggga                                            >  2:72280/1‑144 (MQ=255)
                                                                                        cTCAATGACGCAATGTGACGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATATATTTAGTGTGATGTCGATACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGTATTGCGaaa                                          <  2:363990/149‑1 (MQ=255)
                                                                                           aaTGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACaa                                       <  1:444929/149‑1 (MQ=255)
                                                                                             tGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGt                                     <  1:528687/149‑1 (MQ=255)
                                                                                                 acaAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGttgttg                                                                                  >  1:414344/1‑100 (MQ=255)
                                                                                                 acaAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGttgttg                                                                                  <  2:414344/100‑1 (MQ=255)
                                                                                                   aaaTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAAtgtg                               >  2:537543/1‑149 (MQ=255)
                                                                                                             gcagcTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCt                     <  2:241796/149‑1 (MQ=255)
                                                                                                                      aCGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGCTGTCGAAAATAGGCGTCAGTTGTGACACAGGGTTGTTGGAGGCAGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGGAAGGTGTTGCATCCCTATTaaggcg            >  2:260296/1‑145 (MQ=255)
                                                                                                                      aCGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGGTGTTGGAGGCCGCAGGTATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAGTCCg            >  1:359676/1‑149 (MQ=255)
                                                                                                                      aCGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAAATGGGCGTCAGTTCTGATGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAACAAAGTAATGTGTTGGATCTCTTTTAATcca            >  2:257315/1‑148 (MQ=255)
                                                                                                                       cgtaaTATCAAAACGATGCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGc           <  2:228920/146‑1 (MQ=255)
                                                                                                                         cAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTATTGTGTTGCATCCCATTTAATCCCCat         >  2:240460/1‑149 (MQ=255)
                                                                                                                             atCAAAACGATCCGTGAACATAATTTTAATGTGATGTCGAAAATGAGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCatgatg     >  1:465408/1‑149 (MQ=255)
                                                                                                                                cacaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCg  <  1:234540/146‑1 (MQ=255)
                                                                                                                                aaaaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCg  <  1:377202/149‑1 (MQ=255)
                                                                                                                                aaaaCGATCCGTGAACATAAATTTAATGTGAGGGCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATTTATGCGGCATTGCGAAACAAGTAATGTGTTGCACTCCTATTAATCCGCATGGTGCCg  >  2:215567/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                       
ACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCG  >  NC_000913_3_hsa_tpiA/4604557‑4604833

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑