Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,163,982 G→A M241I (ATG→ATA thiK → thiamine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,163,9820GA100.0% 141.2 / NA 39M241I (ATG→ATAthiKthiamine kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/20);  total (19/20)

TGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTG  >  NC_000913/1163910‑1164059
                                                                        |                                                                             
tGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATACTCAAAGc                                                                       <  2:232635/81‑1 (MQ=255)
  ccACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATACTCAAAGCa                                                                      <  1:232636/80‑1 (MQ=255)
    aCTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATACTCAAAGCAggg                                                                   >  1:232685/1‑81 (MQ=255)
    aCTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATACTCAAAGCAggg                                                                   >  1:232674/1‑81 (MQ=255)
     cTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATACTCAAAGCAGggt                                                                  <  2:232637/81‑1 (MQ=255)
        gcgcgAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATACTCAAAGCAGggtggt                                                               >  1:232669/1‑81 (MQ=255)
            gAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATACTCAAAGCAGGGTGGttt                                                             >  1:232672/1‑79 (MQ=255)
             aaGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATACTCAAAGCAGGGTGGTTTGAg                                                          >  1:232679/1‑81 (MQ=255)
                 tttATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATACTCAAAGCAGGGTGGTTTGAGTAcc                                                      <  2:232638/81‑1 (MQ=255)
                   tATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATACTCAAAGCAGGGTGGTTTGAGTACCg                                                     <  1:232639/80‑1 (MQ=255)
                   tATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATACTCAAAGCAGGGTGGTTTGAGTACCGc                                                    >  1:232677/1‑81 (MQ=255)
                     tCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATACTCAAAGCAGGGTGGTTTGAGTACCGCTg                                                  >  1:232678/1‑81 (MQ=255)
                     tCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATACTCAAAGCAGGGTGGTTTGAGTACCGCTg                                                  <  2:232640/81‑1 (MQ=255)
                       cGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATACTCAAAGCAGGGTGGTTTGAGTACCGCTgg                                                 <  2:232641/80‑1 (MQ=255)
                        ggCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATACTCAAAGCAGGGTGGTTTGAGTACCGCTGGCg                                               <  1:232642/81‑1 (MQ=255)
                            cAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATACTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACa                                            >  1:232693/1‑80 (MQ=255)
                                  tGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATACTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGc                                     >  1:232686/1‑81 (MQ=255)
                                    gCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATACTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCg                                    <  2:232643/80‑1 (MQ=255)
                                     cGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATACTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGAt                                  <  2:232644/81‑1 (MQ=255)
                                      gtcaggtcagGCGATGGTTTCCCTGGCTGCTGATACTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGAt                                  >  1:232682/1‑80 (MQ=255)
                                      gtcaggtcagGCGATGGTTTCCCTGGCTGCTGATACTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATc                                 <  2:232645/81‑1 (MQ=255)
                                         aggtcaggCGATGGTTTCCCTGGCTGCTGATACTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATcaa                               >  2:232676/1‑80 (MQ=255)
                                            tcaggCGATGGTTTCCCTGGCTGCTGATACTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAAt                           <  1:232646/81‑1 (MQ=255)
                                              aggCGATGGTTTCCCTGGCTGCTGATACTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAAtt                          >  2:232690/1‑80 (MQ=255)
                                                gCGATGGTTTCCCTGGCTGCTGATACTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTAt                       <  2:232647/81‑1 (MQ=255)
                                                 cGATGGTTTCCCTGGCTGCTGATACTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATc                      <  2:232648/81‑1 (MQ=255)
                                                   aTGGTTTCCCTGGCTGCTGATACTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAg                    <  1:232649/81‑1 (MQ=255)
                                                    tGGTTTCCCTGGCTGCTGATACTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAgg                   >  2:232683/1‑81 (MQ=255)
                                                       tttCCCTGGCTGCTGATACTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTg                <  1:232650/81‑1 (MQ=255)
                                                          cccTGGCTGCTGATACTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGcc             >  1:232692/1‑81 (MQ=255)
                                                           ccTGGCTGCTGATACTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGcc             >  1:232703/1‑80 (MQ=255)
                                                           ccTGGCTGCTGATACTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGcc             <  2:232651/80‑1 (MQ=255)
                                                            cTGGCTGCTGATACTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGa           <  1:232652/81‑1 (MQ=255)
                                                              ggCTGCTGATACTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGAt          >  2:232673/1‑80 (MQ=255)
                                                              ggCTGCTGATACTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGAt          >  2:232688/1‑80 (MQ=255)
                                                               gctgctGATACTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGAt          <  2:232653/79‑1 (MQ=255)
                                                               gctgctGATACTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATg         <  2:232654/80‑1 (MQ=255)
                                                               gctgctGATACTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGa        >  2:232689/1‑81 (MQ=255)
                                                                     gATACTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTg  >  1:232691/1‑81 (MQ=255)
                                                                        |                                                                             
TGCCACTCGCGCGAAGATTTATCCGGCGCAATTATGGCGTCAGGTCAGGCGATGGTTTCCCTGGCTGCTGATGCTCAAAGCAGGGTGGTTTGAGTACCGCTGGCGACAAACCGGCGATCAACAATTTATCAGGCTGGCCGATGACACCTG  >  NC_000913/1163910‑1164059

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: