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

TTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAA  >  NC_000913/2173297‑2173423
                                                                ||                                                             
ttACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccac                                                                >  1:1765183/1‑65 (MQ=255)
ttACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccac                                                                >  1:3638143/1‑65 (MQ=255)
ttACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccac                                                                >  1:2993461/1‑65 (MQ=255)
  aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                              >  1:997730/1‑64 (MQ=255)
  aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                              >  1:1786440/1‑64 (MQ=255)
  aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                              >  1:4287881/1‑64 (MQ=255)
  aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                              >  1:3583333/1‑64 (MQ=255)
  aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                              >  1:2673695/1‑64 (MQ=255)
    ccGGCACCCAAACAGCAATAAAAATGGTGAGTGGGATAAAAACCAGGCTTGCCGATaccaccgcc                                                            <  1:1745076/65‑4 (MQ=255)
    ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgcc                                                            >  1:1626982/1‑62 (MQ=255)
    ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgcc                                                            >  1:3142476/1‑62 (MQ=255)
    ccGGCACACAAACAGCAATTAAAAGGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgcc                                                            <  1:1359612/65‑4 (MQ=255)
             aaaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTccc                                                 >  1:2475048/1‑65 (MQ=255)
                 aGCAATTAAAATGGTGAGTGGGATAAAAATCCAGCATGCTGATAC‑‑CCTGGCCGAATGTCCcagca                                             >  1:4089946/1‑65 (MQ=255)
                 aGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagca                                             >  1:4001342/1‑65 (MQ=255)
                 aGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagca                                             >  1:2533070/1‑65 (MQ=255)
                      ttAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAAc                                        <  1:3663047/65‑1 (MQ=255)
                       tAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACg                                       >  1:42660/1‑65 (MQ=255)
                       tAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACg                                       >  1:3776083/1‑65 (MQ=255)
                       tAAAATGGTGAGTGGGATAAAAACCAGGCTTGCCGATGCCA‑‑CCGCCGTATGTCCCAGCAGCAACg                                       <  1:3761674/65‑1 (MQ=255)
                          aaTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCg                                    >  1:3574032/1‑65 (MQ=255)
                           aTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCgg                                   >  1:375609/1‑65 (MQ=255)
                            tGGTGAGTGGGATAAAAATCAGGCTTGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGa                                  <  1:4078165/65‑1 (MQ=255)
                             ggTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGAt                                 <  1:432885/65‑1 (MQ=255)
                             ggTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGAt                                 >  1:2388818/1‑65 (MQ=255)
                                      gATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCaa                        <  1:135439/65‑1 (MQ=255)
                                       aTAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAt                       >  1:538020/1‑65 (MQ=255)
                                          aaaaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAg                    <  1:2991743/65‑1 (MQ=255)
                                            aaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCGAGCAGCAACGCCGGATCAAGGCCAATCAGGa                  >  1:3805653/1‑65 (MQ=255)
                                            aaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGa                  >  1:3016252/1‑65 (MQ=255)
                                            aaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGa                  >  1:723046/1‑65 (MQ=255)
                                            aaTCAGCCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGa                  <  1:4033185/65‑1 (MQ=255)
                                             aTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGaa                 >  1:3910740/1‑65 (MQ=255)
                                              tCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAAc                <  1:2212658/65‑1 (MQ=255)
                                              tCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAAc                <  1:4354383/65‑1 (MQ=255)
                                               ccggCTTCCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACt               <  1:2476336/63‑1 (MQ=255)
                                               cAGGTTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACt               <  1:2808767/65‑1 (MQ=255)
                                               cAGGCTTGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACt               <  1:2537691/65‑1 (MQ=255)
                                               cAGGCTGGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACt               <  1:3628253/65‑1 (MQ=255)
                                                  gCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCt            >  1:635004/1‑65 (MQ=255)
                                                        cGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc      >  1:715801/1‑65 (MQ=255)
                                                              ccccccGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGACCTCCTGACCGCCGaa  <  1:2673486/63‑1 (MQ=255)
                                                                ||                                                             
TTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAA  >  NC_000913/2173297‑2173423

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