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.98.1% 227.2 / ‑3.3 53pseudogene (918/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (21/31);  total (22/31)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
*NC_0009132,173,3620C.100.0% 230.6 / NA 52pseudogene (917/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (21/31);  total (21/31)

TTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCC  >  NC_000913/2173223‑2173494
                                                                                                                                          ||                                                                                                                                      
ttCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACACACAGCAATTAAAATGGTGAGCGGGATAAAAATCAGGCTTGACGAGacccccgcc                                                                                                                                     >  2:47953/1‑140 (MQ=255)
   cACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGCTGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATAAAAATGGAGAGTGAGCTCCAACGCAGGCATGCCGCTACCACcccccca                                                                                                                                  >  2:638226/1‑142 (MQ=255)
      gATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGt                                                                                                                             <  1:393363/143‑1 (MQ=255)
        tGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCATTCAAAAGGTTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTcc                                                                                                                           >  1:551549/1‑143 (MQ=255)
           aCGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGGATGTccccg                                                                                                                        >  2:1007599/1‑141 (MQ=255)
            cGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAACATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagc                                                                                                                       <  1:775311/143‑1 (MQ=255)
             ggCGACCGCCATCGCCACGAAGACGCCCATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagca                                                                                                                      <  1:74583/143‑1 (MQ=255)
                       gtCGCCACGATGCAGCCGATCGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAACCAGCACTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCgg                                                                                                            <  2:744544/142‑1 (MQ=255)
                        tCGCCACGAGGCAGCCGTTGGTGGCAAGACCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGa                                                                                                           <  2:672344/143‑1 (MQ=255)
                        tCGCCACGAAGAAGACGCGGGTCGCAAGACTGCCATACGTCAGCACCTGAGTACCCGGCACACTAACAGCTTTTAAAATGTTCAGTGGGAAAAGAATCAGGCTTGCCGATACCA‑‑CCGCCTTATGTCCCAGCAGCAACGCCGGa                                                                                                           <  2:524534/143‑1 (MQ=255)
                           ccACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGGAACGCCGGGTCa                                                                                                        >  1:63784/1‑143 (MQ=255)
                              cGAAGAAGCCGAGGGTGGCAAGATCGCCAAACGGGAGCACCTGACTACCCGGCACACAAACAGCAATTAAAATGTTGACTGGTATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAgg                                                                                                     <  1:255859/143‑1 (MQ=255)
                               gaaaaaGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACACGGCACACAAACAGCAATTAAAATGGTGAGTGGGAAAAAAATCAAGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCCGCAACGACGGAGCAAGGc                                                                                                    >  2:331611/1‑143 (MQ=255)
                                aagaagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAAATAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGcc                                                                                                   <  2:901297/143‑1 (MQ=255)
                                      ccGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAGAATCAGGCTTGCCGATCCCA‑‑CCGCCGTATGCCCCCGCCGCAACGCCGGATCAAGGCCCATccg                                                                                             >  2:13964/1‑141 (MQ=255)
                                        gATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCACCGCCGGATCAAGGCCAATCAGGa                                                                                           >  1:759981/1‑143 (MQ=255)
                                            gtggCAAGATCGCCAAACGGCAGCACCTGATCACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGGTAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAGGGCCAATCAGGAACTc                                                                                       >  1:101304/1‑143 (MQ=255)
                                               gCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTg                                                                                    >  1:970952/1‑143 (MQ=255)
                                                    aTCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGGAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc                                                                              >  1:31759/1‑143 (MQ=255)
                                                    aTCACCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGTATCAAGGCCAATCAGGAACTCCTGAccgc                                                                              <  2:857564/143‑1 (MQ=255)
                                                      cggcAAACGGCAGCACCTGATTACACGGCACACAAACAGCAATTAAAATGGTGAGTGGGAGAAAAATCAGGCTTGCCGATACCA‑‑TCGCCGTATGTCCCAGCAGCACCGCCGCCACACGGCCACACAGACACCCCCGAccgccg                                                                            >  2:435462/4‑143 (MQ=255)
                                                       gCCAAACGGCAGCACCTGATTACCCGGCACACAAACCGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCCTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGa                                                                           >  2:795390/1‑143 (MQ=255)
                                                       gCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGa                                                                           <  1:215964/143‑1 (MQ=255)
                                                            cgggCAGCACCTGATTACCCGTCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGTCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGATCTCCTGACCGCCGAACTTc                                                                     <  2:115220/141‑1 (MQ=255)
                                                                cagcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCt                                                                 <  1:727073/143‑1 (MQ=255)
                                                                  gcaCCTGATTATCCGGCACACAAACAGCAATGGAGATGGTGAGTGGGATAAAAACCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGt                                                               <  2:685141/143‑1 (MQ=255)
                                                                  gcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGt                                                               <  1:1040930/143‑1 (MQ=255)
                                                                   caCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGcc                                                                                                   >  1:101055/1‑108 (MQ=255)
                                                                   caCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGcc                                                                                                   <  2:101055/108‑1 (MQ=255)
                                                                     ccTGGTTACCCGGCACACAACCAGCAATTAAAGTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTaaa                                                            <  2:672355/143‑1 (MQ=255)
                                                                     ccTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTaaa                                                            >  2:710688/1‑143 (MQ=255)
                                                                      cTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAGATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAc                                                           >  1:483895/1‑143 (MQ=255)
                                                                           tACCCGGCACACAAGCAGCAAGTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGCTACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAACGACAATCCGAAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTa                                                      <  2:202992/143‑1 (MQ=255)
                                                                                ggCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGCACTCCTGACCGCCGACCTTCGCCTGTACCCCACTACGAGc                                                 >  1:128816/1‑143 (MQ=255)
                                                                                      caAACAGCAATTAAACTGGTGGGTGGGAAGAACATCAGGTCTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAGCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCtt                                           <  1:376613/143‑1 (MQ=255)
                                                                                       aaaCAGCAATTAAAATGGTGAGTGGGATAAAAATGAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTa                                          <  1:639880/143‑1 (MQ=255)
                                                                                       aaaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATATCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAACCGACTACGAACCTGCtca                                          >  1:198281/1‑141 (MQ=255)
                                                                                         gcgacAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGc                                        <  2:704053/139‑1 (MQ=255)
                                                                                            gCACTTAAACTGGCGACTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                     <  1:396679/143‑1 (MQ=255)
                                                                                             cAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATg                                    <  1:948656/143‑1 (MQ=255)
                                                                                              aggtAAAAGGGGGAGTGGGATAAAGATCAGGCTTGCCGATACCA‑‑CCGCCGGATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATgg                                   <  2:952664/140‑1 (MQ=255)
                                                                                              aaTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATgg                                   <  1:934521/143‑1 (MQ=255)
                                                                                              aaTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGAGCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGACTGTAAACGACAACACGCCTGCATAGCGATgg                                   >  2:7033/1‑143 (MQ=255)
                                                                                                               ggATAAAAATCAGGCTTGCCGATACCA‑‑ACGCCGTATGTCCCAGCAGCAACGCCGAATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGAAAACGACTACGAGCCTGCATAGCGATGGGCGTTAAAccagccaac                  >  2:614194/1‑141 (MQ=255)
                                                                                                                  caaacgtCAGGCTCGCCGACACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGCCCAATCAGGAACGCCTGGCCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgatg               <  2:896175/137‑1 (MQ=255)
                                                                                                                   aaaacTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATAAGGAGCTCCTGCCCGCGGAACGTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATgg              <  1:1197347/143‑1 (MQ=255)
                                                                                                                   aaaaaTCAGGCTCGCCGATACCA‑‑CCGCCGTATGTACCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATgg              <  1:1068369/143‑1 (MQ=255)
                                                                                                                     aaaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGCTAAACCATCCATAATGGGt            >  1:924117/1‑143 (MQ=255)
                                                                                                                      aagcccgCTTGCCGTTACCA‑‑CCGCCGGATGTCCCCGCAGCAACGCCGGATCAAGGCCAATCAGGAACGCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGtt           <  2:384844/137‑1 (MQ=255)
                                                                                                                       aTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGttt          <  1:590650/143‑1 (MQ=255)
                                                                                                                       aTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTAAGGCAACGTCAACTAACACCGCGAAACATTACGATCCTGCTTAGCGCTGAGCGCTGAACCCTCCATGATGGGttt          >  2:532358/1‑143 (MQ=255)
                                                                                                                         cAGGCTTCCCGATCCCA‑‑CCGCGGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTaa        <  2:810216/143‑1 (MQ=255)
                                                                                                                              ttGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTCAGCGATGGGCAGTACACCATCCATGATGGCTTTCATCAc  >  2:722524/1‑142 (MQ=255)
                                                                                                                               tGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc  <  1:1024475/143‑1 (MQ=255)
                                                                                                                                          ||                                                                                                                                      
TTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCC  >  NC_000913/2173223‑2173494

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: