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% 70.5 / NA 15pseudogene (918/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (12/3);  total (12/3)
*NC_0009132,173,3620C.100.0% 70.5 / NA 15pseudogene (917/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (12/3);  total (12/3)

ATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGC  >  NC_000913/2173296‑2173419
                                                                 ||                                                         
aTTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                          >  1:2597059/1‑67 (MQ=255)
 ttACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgc                                                         >  1:2241175/1‑66 (MQ=255)
 ttACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATacaaccgc                                                         >  1:2285071/1‑66 (MQ=255)
 ttACCCGGCACACAAACAGCAAGTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgc                                                         >  1:2256236/1‑66 (MQ=255)
                aCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagact                                        >  1:2473353/1‑65 (MQ=255)
                  aGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACTA‑‑CCGCCGTATGTCCCAgcagca                                      >  1:2449195/1‑68 (MQ=255)
                  aGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAgcagca                                      >  1:1348617/1‑68 (MQ=255)
                  aGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCATGCCGATACCA‑‑CCGCCGTATGTCCCAgcagca                                      >  1:2266775/1‑68 (MQ=255)
                     aaTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACg                                   <  1:215158/68‑1 (MQ=255)
                     aaTTAAAATGGTGAGTGGGATAAAAACCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACg                                   <  1:1214478/68‑1 (MQ=255)
                        tAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCg                                >  1:2152244/1‑68 (MQ=255)
                          aaaTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGa                              <  1:60356/68‑1 (MQ=255)
                                  aGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGcc                      >  1:1120808/1‑68 (MQ=255)
                                    tGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCaa                    >  1:3125681/1‑68 (MQ=255)
                                     gggATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAt                   >  1:1949802/1‑68 (MQ=255)
                                       gATAAAAACCAGGCTTGCCCATAC‑‑CGCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCa                 <  1:443285/68‑1 (MQ=255)
                                          aaaaaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGa              >  1:2919949/1‑68 (MQ=255)
                                             aaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACt           >  1:3439297/1‑68 (MQ=255)
                                             aaTCAGCCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAGCt           <  1:2204978/68‑1 (MQ=255)
                                                   gCTTGCCGATACCA‑‑CCGCGGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAc     >  1:3407848/1‑68 (MQ=255)
                                                   gCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCCATCAGGAACTCCTGAc     >  1:1386944/1‑68 (MQ=255)
                                                   gCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAc     >  1:3322329/1‑68 (MQ=255)
                                                      tGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc  <  1:2196784/68‑1 (MQ=255)
                                                                 ||                                                         
ATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGC  >  NC_000913/2173296‑2173419

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: