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

ACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGC  >  NC_000913/2173299‑2173419
                                                              ||                                                         
aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                          >  1:3087355/1‑64 (MQ=255)
aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                          >  1:2301282/1‑64 (MQ=255)
aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                          >  1:4069871/1‑64 (MQ=255)
aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                          >  1:2725265/1‑64 (MQ=255)
aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                          >  1:3693522/1‑64 (MQ=255)
aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                          >  1:2869087/1‑64 (MQ=255)
aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                          >  1:4409042/1‑64 (MQ=255)
aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                          >  1:1568051/1‑64 (MQ=255)
aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                          >  1:3009712/1‑64 (MQ=255)
aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                          >  1:1409623/1‑64 (MQ=255)
aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                          >  1:52407/1‑64 (MQ=255)
aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                          >  1:830989/1‑64 (MQ=255)
aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                          >  1:1115803/1‑64 (MQ=255)
 cccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgc                                                         <  1:263342/65‑3 (MQ=255)
  ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgcc                                                        >  1:3337237/1‑62 (MQ=255)
  ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgcc                                                        <  1:379974/65‑4 (MQ=255)
  ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgcc                                                        <  1:376938/65‑4 (MQ=255)
  ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgcc                                                        >  1:4090737/1‑62 (MQ=255)
  ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgcc                                                        <  1:1449782/65‑4 (MQ=255)
       acacaAACAGCAATTAAATTGGTGAGTGGGATAAAAATCAGGCTTGCCGATCCCA‑‑CCGCCGTATg                                                 <  1:273636/65‑1 (MQ=255)
           aaaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTcct                                             >  1:205972/1‑64 (MQ=255)
            aacaacAATTAAAATGGTGAGTGGGATAAAAATCAGTCTTGCCGATACCA‑‑CCGCCGTATGTCCca                                            <  1:1651415/65‑1 (MQ=255)
             aCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcag                                           >  1:1880184/1‑65 (MQ=255)
               aGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagca                                         >  1:2376455/1‑65 (MQ=255)
               aGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagca                                         <  1:848740/65‑1 (MQ=255)
                    ttAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAAc                                    <  1:4492063/65‑1 (MQ=255)
                     tAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACg                                   >  1:2133086/1‑65 (MQ=255)
                      aaaaTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGc                                  <  1:1049651/65‑1 (MQ=255)
                       aaaTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGcc                                 >  1:430140/1‑65 (MQ=255)
                          tGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGa                              >  1:2285739/1‑65 (MQ=255)
                           ggTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGAt                             >  1:3029620/1‑65 (MQ=255)
                           ggTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGAt                             >  1:1801458/1‑65 (MQ=255)
                                gTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAgg                        >  1:1976598/1‑65 (MQ=255)
                                gTGGGATAAAAAGCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAgg                        <  1:1614206/65‑1 (MQ=255)
                                  gggATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGcc                      >  1:43713/1‑65 (MQ=255)
                                     aTAAAAATCAGGCTTCCCGACAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAt                   <  1:1028771/65‑1 (MQ=255)
                                        aaaaTCAGGCTTGCCGGTAC‑‑CGCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAg                <  1:4022437/65‑1 (MQ=255)
                                        aaaaTCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAg                <  1:2797358/65‑1 (MQ=255)
                                         aaaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAgg               >  1:3978804/1‑65 (MQ=255)
                                           aTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGaa             >  1:1920504/1‑65 (MQ=255)
                                            tCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAGCGCCGGATCAAGGCCAATCAGGAAc            <  1:362320/65‑1 (MQ=255)
                                             cAGGCTTGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGaagt           <  1:4155182/65‑3 (MQ=255)
                                             cAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACt           <  1:3831197/65‑1 (MQ=255)
                                               ggCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACttc         >  1:935379/1‑63 (MQ=255)
                                               ggCTTGCCGATACCA‑‑CCGCCGTATGGCCCAGCAGCAACGCCGAATCAAGGCCAATCAGGAACttc         >  1:3621570/1‑63 (MQ=255)
                                               ggCTTGCCGATCCCCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTcc         <  1:1548759/65‑1 (MQ=255)
                                                    gCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGGcc    <  1:3910184/65‑1 (MQ=255)
                                                      cGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc  <  1:13791/65‑1 (MQ=255)
                                                      cGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc  >  1:124334/1‑65 (MQ=255)
                                                              ||                                                         
ACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGC  >  NC_000913/2173299‑2173419

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