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.91.1% 170.3 / 7.2 45pseudogene (918/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base . (19/22);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01
*NC_0009132,173,3620C.95.3% 167.1 / 3.7 45pseudogene (917/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base . (19/22);  minor base G (0/2);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGT  >  NC_000913/2173219‑2173496
                                                                                                                                              ||                                                                                                                                      
aGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccac                                                                                                                                         >  1:1081274/1‑143 (MQ=255)
      cgaCGATCCACGGCGACGGCCACCGCCACGACGAAGCCGTTGGTGGCAAGATCGCCAAACGGCAGCACCTGATGACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccgt                                                                                                                                   <  2:1232147/141‑6 (MQ=255)
       cACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccgta                                                                                                                                  >  1:773630/1‑137 (MQ=255)
             gCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTccc                                                                                                                          >  2:411653/1‑143 (MQ=255)
                              gCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGACTGGGATAAAAATCAGCCTTGCCGATAC‑‑CCCCGCCTCATGTCCCAGCAGCAAcccccgccc                                                                                                         >  2:172133/1‑138 (MQ=255)
                               ccACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCa                                                                                                        <  2:1000242/143‑1 (MQ=255)
                                 aCGAAGAAGCCGCGGGTGGCAAGATCGCCAAACGGCAGCACCTGATTCCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCATCAACGCCGGATCAAg                                                                                                      <  1:950179/143‑1 (MQ=255)
                                   gaagaaGGCGATGGGGACAAGGTCGCCAAACGGCAGCACCTGCTTCGCGCGCCCACAATCAGCTAGTAAAAAGGCGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGc                                                                                                    <  2:563277/143‑1 (MQ=255)
                                          ccGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACGGCAAGTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAg                                                                                             <  2:880845/143‑1 (MQ=255)
                                            gATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCTGCAACGCCGGCTCGCGGCCAATcgaga                                                                                           >  1:1150272/1‑139 (MQ=255)
                                                 tggCAAGACCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAACCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTcc                                                                                      <  1:310111/143‑1 (MQ=255)
                                                  ggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTCTGTCCCAGACGCCACGCCCGAGCAAGGCCAATCAGGCCCTccg                                                                                     >  2:1047517/1‑142 (MQ=255)
                                                   gCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCACGGCCCCTCAGGAACTCCTg                                                                                    >  1:280832/1‑143 (MQ=255)
                                                        aTCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc                                                                               <  1:253876/143‑1 (MQ=255)
                                                         tCGCCAAACGGCAGCACCTGATGACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgcc                                                                              <  2:574475/143‑1 (MQ=255)
                                                           gCCAAACGGCAGCACCAGATTACCCGGCACACAAACAGCAATCAAAATGGTGAGTGGGATACAAATCAGGCTCGCCGATCCCA‑‑CCGCCGTATGTCCCAGCCCCCACGCCGCACCAAGGCCCACCAGTCACCCCCGCCCGCCGa                                                                            >  1:246347/1‑143 (MQ=255)
                                                            ccAAACGGCAACACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAACAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGAAGCAACGCCGGATCAAGGCACAACAGGAAATCCTGAcccgcaga                                                                           >  2:308129/1‑137 (MQ=255)
                                                             cAAACGGCAGCACCTGATTACCCGGCACACACACAGAAATAAGAAGGGAGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGCATGTCCCAGCGGCGACCCCGGAACAAGGCGAAGCAGGAACTCCTGACCGCCGacc                                                                          >  2:962019/1‑141 (MQ=255)
                                                             cAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAAc                                                                          <  1:1175526/143‑1 (MQ=255)
                                                              aaaCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGCACt                                                                         >  1:953511/1‑143 (MQ=255)
                                                               aaCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACtt                                                                        <  1:739696/143‑1 (MQ=255)
                                                                aCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAACAATCAGGCTTGCCGCTACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTc                                                                       <  2:95188/143‑1 (MQ=255)
                                                                 cgccagcaACTGCTTACCCAGCACACAGACAGCAATGAGAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGCATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCg                                                                      <  2:160457/140‑1 (MQ=255)
                                                                        aCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTaa                                                               >  1:735270/1‑143 (MQ=255)
                                                                         ccTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTaaa                                                              <  1:201678/143‑1 (MQ=255)
                                                                          cTGATTACCCGGCACACAAACAGCAATTCAGATGGTGAGCGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAc                                                             <  1:231230/143‑1 (MQ=255)
                                                                               tACCCTGCACACAAACAGCAATTAAAATGGGGAGTGAGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTa                                                        <  1:764376/143‑1 (MQ=255)
                                                                                  ccGGCACACAAACAGCAATTAAAATGGTGAGTGAGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGAAAGACCCAGCGGCGACGCCGGATCAAGGCCAACCAGGAACTCCTGACCGCCGAACTGCGCCGGTAAACGACTACgc                                                     >  2:1029936/1‑142 (MQ=255)
                                                                                         acaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACGCCTGACCGCCGAACTTCGCCTGTACACGACTACGACCCTGCt                                              >  1:210073/1‑143 (MQ=255)
                                                                                            aaCCGCACTTAAAATGCTGCGCGGGCCAAAAATCAGGCTTGCCGCTACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGAGCAAGGCCAATCAGGAACTCCTGACCGCCGACCTTCGCCTGTAAACGACTACGAGCCTGCTTAg                                           <  1:26099/143‑1 (MQ=255)
                                                                                               aGCAATTAACACGGCGAGTGGGATAAAAATCAGGCCTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGa                                        <  2:1184251/143‑1 (MQ=255)
                                                                                                  aagtAAAAGGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATgg                                     <  1:84742/140‑1 (MQ=255)
                                                                                                   aTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACACCGAGCGACCTTAGCGATggg                                    >  1:1189676/1‑143 (MQ=255)
                                                                                                     tAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCg                                  >  1:844159/1‑143 (MQ=255)
                                                                                                     tAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCg                                  >  2:185354/1‑143 (MQ=255)
                                                                                                        aaTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTGAGCGATGGGCggtg                               >  1:165652/1‑140 (MQ=255)
                                                                                                         gtGGTGAGGGGGATAAAAAGCAGGCTTGCCGATAC‑‑CGCCGCCGTATGTCCCAGCAGCAGCGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaa                              <  2:641120/142‑1 (MQ=255)
                                                                                                          cggcgAGTGGGATAAAACTCAGGCCTGCCGATACCA‑‑CCGCCGTATGTCCCCGCAGCCCCGCCCGATCAAGGCCAAGCAGGAACTCCTGACCGCCGAACTTCGCCTGGAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaaa                             <  2:962466/139‑1 (MQ=255)
                                                                                                           ggTGCGTGGGACAAACACCAGGCGTGCCGAT‑‑CACACCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAc                            <  2:9485/143‑1 (MQ=255)
                                                                                                                 tGGGATAAAACTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatcca                      <  1:1047517/143‑1 (MQ=255)
                                                                                                                     aTAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTACACGACTACGAGCCTGCTTCGCGTTTGGCGTTAAACCATCCatgat                  >  1:436717/1‑143 (MQ=255)
                                                                                                                      tAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCCATGAGCGGTAAACCATCCatgatg                 >  1:680249/1‑143 (MQ=255)
                                                                                                                       aaaaaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACGCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATgg                <  1:411653/143‑1 (MQ=255)
                                                                                                                          aaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGtt             >  1:1014363/1‑143 (MQ=255)
                                                                                                                             cAGGCTTGCCGATACCA‑‑CCGCCGGATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTaa          <  1:945924/143‑1 (MQ=255)
                                                                                                                              tggCTTGCCGATAC‑‑CGCCGCCGTATGTCCCAGCAGCAAGGCCGGATCAAGGCCAATGAGGAACGCCGGACCGCCGAACTTCGCCTGTAACCGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAt         <  1:696454/142‑1 (MQ=255)
                                                                                                                                gCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAAGCCAATCAGGAACTCCTGACCGCCGCACTTCGCCTGTAAACGACGACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATCATGGGCATAATCa       >  2:928188/1‑143 (MQ=255)
                                                                                                                                     ccGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCACGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGt  <  1:471648/143‑1 (MQ=255)
                                                                                                                                              ||                                                                                                                                      
AGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGT  >  NC_000913/2173219‑2173496

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: