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.87.5% 106.3 / 8.7 32pseudogene (918/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/4);  new base . (12/16);  total (12/20)
Fisher's exact test for biased strand distribution p-value = 2.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009132,173,3620C.90.3% 107.2 / 7.4 32pseudogene (917/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base . (12/16);  total (12/20)
Fisher's exact test for biased strand distribution p-value = 2.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAT  >  NC_000913/2173227‑2173489
                                                                                                                                      ||                                                                                                                               
aCGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccgtat                                                                                                                          >  2:43743/1‑136 (MQ=255)
  gATGCACGGCGACCGCCATGGCCACGAAGACGCCGAGGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGt                                                                                                                      <  2:334928/143‑1 (MQ=255)
  gATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGGGAGGGGGGTAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGt                                                                                                                      >  1:281751/1‑143 (MQ=255)
          gCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGCATGTCCcagcag                                                                                                              >  2:39067/1‑143 (MQ=255)
                gCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCCGCACCGc                                                                                                        >  1:723200/1‑143 (MQ=255)
                    ttgCCACGAACGAGCCGATGGGGGCAAGGTCGCCAACCGGCAGCCGCTGATTACCCGGCACACAGACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGa                                                                                                    <  2:705178/141‑1 (MQ=255)
                    tCGCCACGAAGAAGCCGATGGTGGCAAGATCGGCAGGCGCCAACACCGGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGa                                                                                                    <  2:336289/143‑1 (MQ=255)
                    gcGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCCATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGa                                                                                                    <  2:553489/142‑1 (MQ=255)
                        cACGAAGAAGCCGTTGGTGGCAAGGTCGCCAAACGGCAGCACCTGATTACCCGGAACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCaa                                                                                                <  2:48603/143‑1 (MQ=255)
                         aCGAGGAAGCTGATGGTGGCAAGATCGCGAAACGGCAGCACCTGGTTACCCGGCACACAAACAGCAAGTAAAAGGGTGAGTGGGATAAAAATCAGGCTTGCGGAGACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAg                                                                                               <  2:714276/143‑1 (MQ=255)
                         aCGAAGCAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAg                                                                                               <  2:734509/143‑1 (MQ=255)
                         aCGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCCGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAg                                                                                               <  1:65857/143‑1 (MQ=255)
                                  ccGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAg                                                                                      <  1:43743/143‑1 (MQ=255)
                                                agcGCCAAACGGCGGCGCCTGATTACCCGGCACACAAACAGCAATTAATATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc                                                                        <  2:411563/141‑1 (MQ=255)
                                                   gCCAAACGGCAGCACCTGATTACCCGGCACACAAACCGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGGACTCCTGAccgccgg                                                                     >  1:492870/1‑142 (MQ=255)
                                                          ggCAGCACCTGATTGCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                              <  1:96823/143‑1 (MQ=255)
                                                                aCCCGATTACCCGGCACACAAACAGCAAGTAAAATGGTGAGTGCGATAAAAATCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCTGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTaa                                                        <  1:730799/143‑1 (MQ=255)
                                                                               acacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTg                                         >  2:512773/1‑143 (MQ=255)
                                                                                 acaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGCACTTCGCCTGTAAACGACTACGAGCCGGCt                                       >  1:633894/1‑143 (MQ=255)
                                                                                 acaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCt                                       <  1:124878/143‑1 (MQ=255)
                                                                                   aaaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTa                                     >  1:455692/1‑143 (MQ=255)
                                                                                              aaaaTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTACGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGt                          <  1:483385/143‑1 (MQ=255)
                                                                                                   ggTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAc                     >  1:713074/1‑143 (MQ=255)
                                                                                                       agcgggATAAAAAGCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatc                 <  1:258139/140‑1 (MQ=255)
                                                                                                         tGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccg                                                                         >  2:263122/1‑85 (MQ=255)
                                                                                                         tGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccg                                                                         <  1:263122/85‑1 (MQ=255)
                                                                                                           ggATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAATGCCAATAAGGAACTACTGACCGACGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGTGCGTTAAAccatgcacg             >  2:671535/1‑141 (MQ=255)
                                                                                                             aTAAAGATCGGGCTTGCCGAGAC‑‑CCCCGCCGTATGTCCCGGCACCAACGCCGGATCAAGGCCAATCAGGAACCCCTGACCGCCGAACGTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgat           <  2:340907/143‑1 (MQ=255)
                                                                                                               gaaaaTCACGCATGCCGCTAC‑‑CGCCGCCGTATGGCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATgg         <  1:202058/142‑1 (MQ=255)
                                                                                                               aaaaaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGGTAAACCATCCATGATgg         >  1:427696/1‑143 (MQ=255)
                                                                                                                  aaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGACCTCCTGACCGCCGACCTTCGCCTGTAAACGACTACCAGCCTCCTGCGCGCTGCGCGTTAAACCACCCATGATGGGta      >  2:411911/1‑142 (MQ=255)
                                                                                                                  aaTCAGGCTGGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGtt      <  2:664915/143‑1 (MQ=255)
                                                                                                                      aGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAt  <  1:674194/143‑1 (MQ=255)
                                                                                                                                      ||                                                                                                                               
ACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAT  >  NC_000913/2173227‑2173489

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: