Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,173,363 Δ2 bp pseudogene (915‑916/1358 nt) gatC ← galactitol‑specific PTS enzyme IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,173,3610C.100.0% 90.1 / NA 19pseudogene (918/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (10/9);  total (10/9)
*NC_0009132,173,3620C.100.0% 90.1 / NA 19pseudogene (917/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (10/9);  total (10/9)

TGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCG  >  NC_000913/2173294‑2173418
                                                                   ||                                                        
tGATTACCCGGCACACAAACAGCAATTAAAATGGAGAGTGGGATAAAAATCAGGCTTGCCGATaccac                                                           <  1:3429998/68‑1 (MQ=255)
    tACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgcc                                                       <  1:424496/68‑4 (MQ=255)
           cacgcaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTc                                              <  1:3109674/68‑1 (MQ=255)
            acacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTcc                                             <  1:3124938/68‑1 (MQ=255)
              acaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCca                                           >  1:3505983/1‑68 (MQ=255)
                aaaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagc                                         >  1:1363767/1‑68 (MQ=255)
                 aaCAGCAATAAAAATGGTGATGGGAATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagca                                        <  1:2066742/68‑1 (MQ=255)
                    aGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAgcagca                                     >  1:2291966/1‑68 (MQ=255)
                    aGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAgcagca                                     >  1:2358397/1‑68 (MQ=255)
                    aGCAATTAAAACGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAgcagca                                     >  1:2590063/1‑68 (MQ=255)
                          tAAAATGGTGAGTGGGACAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCg                               >  1:158151/1‑68 (MQ=255)
                             aaTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGAt                            <  1:667178/68‑1 (MQ=255)
                                ggagAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGGATGTCCCAGCAGCAACGCCGGATCaa                         >  1:3366469/4‑68 (MQ=255)
                                  tGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAgg                       <  1:1941834/68‑1 (MQ=255)
                                  tGAGTGGAATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAgg                       <  1:1975534/68‑1 (MQ=255)
                                          aTAAAAATCGGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAg               <  1:1308663/68‑1 (MQ=255)
                                            aaaaaTCAGGCTTGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGa             <  1:3494300/68‑1 (MQ=255)
                                               aaTCAGGCTTGCCGATAGCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATTAAGGCCAATCAGGAACt          <  1:276945/68‑1 (MQ=255)
                                               aaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACt          <  1:1380741/68‑1 (MQ=255)
                                               aaTCAGGCTTGCCGATACCA‑‑CCGCCGGATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACt          >  1:1793048/1‑68 (MQ=255)
                                                 tCAGGCTTCCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTcc        <  1:1463579/68‑1 (MQ=255)
                                                    ggCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGa     >  1:185532/1‑68 (MQ=255)
                                                     gCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAc    >  1:1798816/1‑68 (MQ=255)
                                                       ttGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccg  <  1:2051169/68‑1 (MQ=255)
                                                                   ||                                                        
TGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCG  >  NC_000913/2173294‑2173418

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