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

TTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGC  >  NC_000913/2173297‑2173419
                                                                ||                                                         
ttACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccac                                                            >  1:3119649/1‑65 (MQ=255)
  aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                          >  1:1325328/1‑64 (MQ=255)
  aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                          >  1:3322749/1‑64 (MQ=255)
    ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgcc                                                        <  1:4073900/65‑4 (MQ=255)
    ccGGCACACAAACAGCAATTAAAATGGAGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgcc                                                        >  1:2065647/1‑62 (MQ=255)
             aaaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTcct                                             >  1:3131215/1‑64 (MQ=255)
              aaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCca                                            <  1:1409805/65‑1 (MQ=255)
                 aGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagca                                         >  1:1594694/1‑65 (MQ=255)
                 aGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagca                                         <  1:3216981/65‑1 (MQ=255)
                 aGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagca                                         >  1:2804331/1‑65 (MQ=255)
                 aGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGCATGTCCcagca                                         >  1:4064190/1‑65 (MQ=255)
                 aGCAATTAAAATGGTGAGTGAGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagca                                         <  1:3683444/65‑1 (MQ=255)
                         aaaGGGTGAGTGGGATAAAAATCAGGCTTGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGcc                                 <  1:3628709/65‑1 (MQ=255)
                            tGGTGAGTGGTATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGa                              <  1:3964835/65‑1 (MQ=255)
                            tGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGa                              >  1:2665340/1‑65 (MQ=255)
                             ggTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGAt                             >  1:2477949/1‑65 (MQ=255)
                                  gTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAgg                        <  1:965563/65‑1 (MQ=255)
                                    gggATCNAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGcc                      >  1:2578458/1‑65 (MQ=255)
                                    gggATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCGAGCAGCAACGCCGGATCAAGGcc                      >  1:3425801/1‑65 (MQ=255)
                                    gggATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGcc                      >  1:754563/1‑65 (MQ=255)
                                       aTAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAt                   <  1:231196/65‑1 (MQ=255)
                                          aaaaTCAGGCTTGCCGATCCCCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAg                <  1:1552191/65‑1 (MQ=255)
                                            aaTCAGGCTTGCCGATACCA‑‑CCGCCGTACGCCCCAGCAGCAACGCCGGATCAAGGCCAACCACGa              >  1:3585967/1‑65 (MQ=255)
                                             aTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGaa             >  1:1694002/1‑65 (MQ=255)
                                              tCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAAc            <  1:2436728/65‑1 (MQ=255)
                                               cAGGCTTGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACt           <  1:2388088/65‑1 (MQ=255)
                                               cAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACt           <  1:2928453/65‑1 (MQ=255)
                                               cAGGCTGGCCGATAC‑‑CTCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACt           <  1:4046636/65‑1 (MQ=255)
                                                 ggCTTGCCGATGC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTcc         <  1:2153749/65‑1 (MQ=255)
                                                 ggCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTcc         <  1:3941546/65‑1 (MQ=255)
                                                  gCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCt        >  1:1002086/1‑65 (MQ=255)
                                                        cGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc  <  1:1571221/65‑1 (MQ=255)
                                                        cGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc  >  1:4258087/1‑65 (MQ=255)
                                                                ||                                                         
TTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGC  >  NC_000913/2173297‑2173419

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