Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,173,363 Δ2 bp intergenic (‑1/+1) gatC ← / ← gatC pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC/pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,173,3610C.94.2% 211.2 / 6.0 52pseudogene (2/442 nt)gatCpseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base . (28/21);  total (29/23)
Fisher's exact test for biased strand distribution p-value = 5.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01
*NC_0009132,173,3620C.98.0% 209.5 / ‑0.4 52pseudogene (1/442 nt)gatCpseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base . (28/21);  minor base T (0/1);  total (29/23)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCG  >  NC_000913/2173223‑2173495
                                                                                                                                          ||                                                                                                                                     
ttCCACGATGCACGGCGACCGCCAACGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTTATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgcc                                                                                                                                    <  1:776287/143‑4 (MQ=255)
    aCGATGCACGGCGACCGCCATCGCCACGAAGGAGCCGATGGGGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAGGCAGGCTTGCCGATaccaccgccgtat                                                                                                                                <  2:897972/143‑8 (MQ=255)
        tGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTCAAATGGTGCGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTcc                                                                                                                          <  1:1062400/143‑1 (MQ=255)
          cACCGCGACGGCTATCTCCCGGAAGAAGCCTATGGTGGCAAGATTGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCca                                                                                                                        <  2:1339923/143‑1 (MQ=255)
                 aCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATCCCA‑‑CCGCCGTATGTCCCagcagcca                                                                                                                 >  1:488439/1‑141 (MQ=255)
                 aCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTAACCGGCACACAAACAGCAATTAAAATGGTGAGGGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTGCCCACAGCaa                                                                                                                 >  1:475955/1‑143 (MQ=255)
                   cGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATGGCCAAACGGCAGCACCTGATTACCCGGCACACCAACAGCAATTAAAATGGTGCGTGGGATAGAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACg                                                                                                               >  1:259612/1‑143 (MQ=255)
                    gCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCGCATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGc                                                                                                              >  1:1273257/1‑143 (MQ=255)
                         cGCCACGAAGAAGCCGCTGGTGGCAAGATCGCCAAACGGCAGCACCTGAGTCCCCGGCACACAAACAGCAATTAAACTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGTAt                                                                                                         <  2:1034274/143‑1 (MQ=255)
                         cGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGAt                                                                                                         <  1:913115/143‑1 (MQ=255)
                            cACGAAGCAGCCGATGGTGGCAAGAGCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATAC‑‑CTCCGCCGTATGTCCCAGCAGCAACGCCGGATCaa                                                                                                      <  1:607346/143‑1 (MQ=255)
                             aCGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCCGCACCTGATTACCCGGCACACAAACAGCAATTCACATGGTGAGTGGGATAAAAATCAGGCTTGCCGAAACCA‑‑CCGCCGCACGTCCCACCACCCACGCCGGATCAAg                                                                                                     >  2:918877/1‑143 (MQ=255)
                             aCGAAGAAGCCGATGGTGGAAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAg                                                                                                     <  2:327830/143‑1 (MQ=255)
                              cGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAACTCAGGCTTGCCGATACCA‑‑ACGCCGTATGACCCAACAGCCACGACGGACCGAgg                                                                                                    >  2:755767/1‑143 (MQ=255)
                              cGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGCATGTCCCAGCGGCAACGCCTCATCAAgg                                                                                                    >  1:452396/1‑143 (MQ=255)
                               gaagaaGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTCCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGc                                                                                                   <  1:1079967/143‑1 (MQ=255)
                               gaagaaGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCAGAACACATATTACAATTTTTATGGTGGGTGGGGTAAAAATCAGGCTTGCCGATAC‑‑CCCCGCCGCATGTCCCAGCCGCGACGCCGGAACAAgac                                                                                                   >  1:693041/1‑141 (MQ=255)
                                aagaagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCATATGACCCAGCAGCAACGCCGGATCAAGGcc                                                                                                  >  1:882438/1‑143 (MQ=255)
                                   gagCCGGTGGTGGCAGCAGCGCGAAAAGTGAGGACCTGATGTCCCGGCAGACAAACAGCAATGAAAAGGGTGAGTGGGATAAAAATCAGGCTTGCCGATAC‑‑CGCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAt                                                                                               <  2:955185/142‑1 (MQ=255)
                                    agCCGATGGTGGCAAGATCGCCAAACTGCAGCACCTGACTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATc                                                                                              <  1:36257/143‑1 (MQ=255)
                                      ccGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CTGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAg                                                                                            <  2:645584/143‑1 (MQ=255)
                                        gATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGa                                                                                          >  1:1008090/1‑143 (MQ=255)
                                         agggtggCAAGATCGCCAAACGACAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGaa                                                                                         <  1:874050/141‑1 (MQ=255)
                                                      cGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgccg                                                                            <  2:410187/143‑1 (MQ=255)
                                                         cAAACGGCAGCACCTGATTACCCGGCACACAAACCGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAAGCCCATCAGGCACTCCTGACCACCGCAc                                                                         >  2:511323/1‑143 (MQ=255)
                                                            aCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGGACTCCTGACCGCCGAACTTc                                                                      >  1:174501/1‑143 (MQ=255)
                                                              ggCAGCACCTGATTACCCGGCACACAAACATCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                    <  2:428756/143‑1 (MQ=255)
                                                                 agcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGGGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTg                                                                 <  2:1170327/143‑1 (MQ=255)
                                                                   caCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTa                                                               >  1:816109/1‑143 (MQ=255)
                                                                   caCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTa                                                               <  1:897781/143‑1 (MQ=255)
                                                                     tcTGATTACCCGGCACACATACAGCAATTAAAATGGTGAGTGGGATAAAAACCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTaaa                                                             <  2:146304/142‑1 (MQ=255)
                                                                     ccTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTaaa                                                             <  1:27355/143‑1 (MQ=255)
                                                                       tGATTACCCGGCACACAAACAGCAATTAAAATGGCGAGTGGGGTAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCCGACCGCCGAACTTCGCCTGTAAACg                                                           >  2:232306/1‑143 (MQ=255)
                                                                        gATTACCCTGCACACACACAGCTATTAGAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGa                                                          <  1:934339/143‑1 (MQ=255)
                                                                        gATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGACCTTCGCCTGTACACGa                                                          >  1:317598/1‑143 (MQ=255)
                                                                                ggCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAATTCAGGCTTGCTTATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                  <  1:1097764/143‑1 (MQ=255)
                                                                                     accaaCAGCAATTAAAATGGTGAGTGGGATAACAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCt                                             <  2:58592/140‑1 (MQ=255)
                                                                                        aaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAg                                          <  1:1278501/143‑1 (MQ=255)
                                                                                        aaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACAGCCGAACGTCGCCTGTAAACACCTACGAGCCTGATTAg                                          >  2:726332/1‑143 (MQ=255)
                                                                                           aGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGGCACAGCAGGAAGTTCGGATCAAGGCCAATAAGGACCTCCTGACCGCCGAACCTCGCCGCCGAACGACCAAGAGCCTGCtcacgga                                       >  2:1248166/1‑137 (MQ=255)
                                                                                               aTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATggg                                   >  2:364148/1‑143 (MQ=255)
                                                                                                   aaaTGGTGAGTGGGATAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTa                              >  1:1003083/1‑143 (MQ=255)
                                                                                                          gggggggATAAAAATCAGGCTCGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccat                        <  2:1273257/139‑1 (MQ=255)
                                                                                                              gggATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTACACGACTACGAGCCTACTTAGCGATGGGCGTTAAAccatccac                    >  1:182976/1‑142 (MQ=255)
                                                                                                                gATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGCTAAACCATCCagga                  >  1:1327144/1‑143 (MQ=255)
                                                                                                                 aTAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGCGCCTGCTTAGCGATGGGCGTTAAACCATCCatgat                 >  1:779480/1‑143 (MQ=255)
                                                                                                                 aTAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatggt                 >  1:1254056/1‑141 (MQ=255)
                                                                                                                  tAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTGCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGAGAAACCGTCCatgatg                >  1:203691/1‑143 (MQ=255)
                                                                                                                    aaaaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGCTCAACCACCCACGCTggg              >  1:1118206/1‑143 (MQ=255)
                                                                                                                      aaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGtt            <  1:1106280/143‑1 (MQ=255)
                                                                                                                            gCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGGCTACGAGCCTGCTTAGCGACGGGCGGTAAACCGTCGTTGCTGGGGTTACTCa      >  1:212352/1‑143 (MQ=255)
                                                                                                                             cTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAc     >  1:900598/1‑143 (MQ=255)
                                                                                                                                gCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACAACGAGCCTGCTTACCGACGCGCGTTAAACCATCGATGAAGAGTTCAATAAcccc  >  2:1281779/1‑142 (MQ=255)
                                                                                                                                gCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCACCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGAGCGATAAACCATCCATGATGGGTTTAATCACCCg  >  2:579181/1‑143 (MQ=255)
                                                                                                                                          ||                                                                                                                                     
TTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCG  >  NC_000913/2173223‑2173495

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: