Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,173,363 Δ2 bp intergenic (‑490/+918) gatD ← / ← gatB galactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,173,3630C.80.8% 59.8 / 13.6 26intergenic (‑490/+919)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (1/4);  new base . (16/5);  total (17/9)
Fisher's exact test for biased strand distribution p-value = 3.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01
*NC_0009132,173,3640C.84.6% 68.8 / 8.5 26intergenic (‑491/+918)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/4);  new base . (17/5);  total (17/9)
Fisher's exact test for biased strand distribution p-value = 8.43e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.53e-01

AACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCG‑ATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGAC  >  NC_000913/2173148‑2173606
                                                                                                                                                                                                                        ||                                                                                                                                                                                                                                                  
aacaaaTAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGCCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCACACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGCTAAAAATCAGGCTTGCCG‑ATACCACC‑‑GCCGTATGTCCCAGCAGCAACGACGGCTCAAGGca                                                                                                                                                                                                                 >  2:163038/5‑249 (MQ=255)
aaCCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAACTCAGGCTTGCCG‑ATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGcc                                                                                                                                                                                                                 >  1:494356/1‑250 (MQ=255)
  acAATCGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGACGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGAGCGCGGGCTAAAAATCAGGCTTGCCG‑ATACCACC‑‑GCCGTATGCCCCAGCAGCAACGCCGGATCAAGGcccac                                                                                                                                                                                                              >  2:283541/2‑248 (MQ=255)
                        acaGGGTGATGCTCATAATGATGCCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGGGAGTGGGATAAAAATCAGGCTTGCCG‑ATACCACC‑‑GCCGTATGTCCCAGCAGCACCGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAAc                                                                                                                                                                                        >  2:173475/1‑251 (MQ=255)
                                   cTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCG‑ATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTa                                                                                                                                                                              >  2:15644/1‑250 (MQ=255)
                                        aaTGATGACACCCGCGAATAAGGTGCGGAACAGACTTCCACAATGCAAGACGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAACCGGCAGCACCTGATTACCCGGCACACAAAACGCAATTAAAATGGTGAGGGGGATAAAACTCAGGCTTGCCG‑ATACCACC‑‑GCCGTATGTCACAGAAGCAACGCAGGATCACGACCAATCAGGAACTCCTGGCATCCGAACTTCGCCTGTAAACgcc                                                                                                                                                                        >  2:72697/1‑249 (MQ=255)
                                        aaTGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCG‑ATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGAc                                                                                                                                                                        >  1:483161/1‑251 (MQ=255)
                                           gatgaCACCCGAGATTACGGAGCGGAACAGATTTCCACGATGCAAGGCGACCGCCATCGCCACGAAGAAGCCGCTGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCG‑ATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGCACTTCGCCTGTAAACGACTAc                                                                                                                                                                     >  2:476614/1‑251 (MQ=255)
                                                     gagaTTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCG‑ATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAgcttgctt                                                                                                                                                           >  1:391183/1‑251 (MQ=255)
                                                      aaaTTAAGGTGCGGAACAGATTTCCACGGTGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAACCGGCAGCACCTGATTACCCGGCACACACACAGCAATTAACATGGTGAGTGGCCTAAAAATCAGGCTGACCG‑TTACCACC‑‑GCCGTCTGGCCCAGCAGCAACGCCGGATCAATGCCAGTCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGctttcttc                                                                                                                                                          >  2:273699/3‑244 (MQ=255)
                                                               ggCGGAACAGATTCCCACGATGCACGGCGACCTCCATCGCCACGAAGAATCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCATTTAAAATGGTGAGTGGGATAAAAAGCAGGCTTGCCG‑ATACCACC‑‑TCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGc                                                                                                                                                 <  1:85212/250‑1 (MQ=255)
                                                                             ccACGATGCACGGCGACCGCCATCGCCACGCAGAAGCCGATGGTGGCAAGATCGCCCAACCGCAGCACATGATTACCCGGCACCCAATCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCG‑ATCCCACCC‑‑CCGTATGTCCAAGCAGCATCGCCGGATCAAGGCCCATCAGGAACTCCTGACCGCCGAACTTCGCCTGTGACCCTCTTCGAGCTTGCTTAGCGATGGGCGTTAACccatccct                                                                                                                                   >  2:154823/1‑249 (MQ=255)
                                                                               aCGATGCCCGGCGCCCGCCATCGCCACGAAGAAGCCGATGGGGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGCCACACAAACAGCAATTAAAAGGGGGAGTGAGATAAAAATCAGGCTTGCCG‑AT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCatga                                                                                                                                 <  1:495643/251‑1 (MQ=255)
                                                                                                      ccaagaagaaGCCGATGGTGGCAATATCGCCAATCGGCAGCACCTGATTACCCGGCACACAATCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCG‑ATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCAtt                                                                                                          <  1:65048/247‑1 (MQ=255)
                                                                                                               ggCCGATGGTGGCAAGATCGCCAAACGGCGGCACCTGATTACCCGGCACACAACCAGCAATTAAAAGGGTGAGTGGGAAAAAAATCAGGCTTCCCG‑ATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATc                                                                                                 <  2:158948/250‑1 (MQ=255)
                                                                                                                   gATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCG‑ATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTg                                                                                             >  1:441839/1‑251 (MQ=255)
                                                                                                                        gggCAAGATCGCCAAACGGCAGCTCCTGATTACCCGGCACACAAACAGCAATTAAAATGGGGAGTGGGATAAAAATCAGGCTTGCCG‑ATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTGAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgccgcc                                                                                        <  1:321690/250‑1 (MQ=255)
                                                                                                                              gATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCG‑ATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTATCAGCATCACTTCCGCCGttttt                                                                                  >  1:92212/1‑251 (MQ=255)
                                                                                                                                  gccaaactggagcaactgcatACCCGGCCCACAAACAGCAATTAAAATGGTGCGTGGGATAAAAATCAGGCTTGCCG‑ATACCACC‑‑GCCGTATGTCCCCGCAGAAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAAATTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGTCGTTCAACCTTCCATGCTGGGTTTAATCACCCATGGCATTAGAAGCATCACTGCAGCCGTTTTTAccc                                                                              >  2:37779/21‑250 (MQ=255)
                                                                                                                                    cAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCG‑ATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGcc                                                                            >  1:146666/1‑251 (MQ=255)
                                                                                                                                               ccataattccccgccaaacaACCAGAAATAAAATTGGGGAGTGAGATAAAAACCAGCTTTCCCG‑AT‑CC‑CCCCGCCGTATGTCCCATCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGAATTAGCAGCATCACTGCCGCCGTTTTTTCCGCCAGCTGCAAt                                                                   <  2:22812/232‑1 (MQ=255)
                                                                                                                                                    ataccccGGCACAAAAGCAGCAATAAAAATGGGGAGGGTGAAAAAATTCAGGCT‑GCCGGAT‑CC‑CGCCGCCGTATGCCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTGCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCtt                                                            <  1:283541/247‑1 (MQ=255)
                                                                                                                                                                  aaaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCG‑ATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCACGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGATTAGCGATGGGCGTTAAACCATCCTTGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAcccc                                              >  2:165239/1‑251 (MQ=255)
                                                                                                                                                                   aacaacAATAAAAATGGTGAGTGGGAAAAAAACCAGGCTTGCCA‑AT‑CC‑CCCCGCCGTAGGGCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCCCTGCCGCCGTTTTGACCGCCAGCTTCAATACACCTTTGACATCGTAACCcg                                             <  2:494356/251‑1 (MQ=255)
                                                                                                                                                                        cAATTAAAATGGTGATTGGGATAAAAATCAGGCTTGCCG‑ATTCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAgg                                        <  1:304598/251‑1 (MQ=255)
                                                                                                                                                                                                aaaaaacaaggctcacagaaccccccGCCGTAGTTCCCAGCAGAAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCAGCCATGATGGGTTTAAGCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCGTTGACATCGTAACCCGCTAGGAGGCCGATAATCAGCCCCATc                   <  2:56873/229‑1 (MQ=255)
                                                                                                                                                                                                         gctttccagatgccaccgccgTATGCCCCAGCAGAAACCCCGGACCAAGGCCAATCAGGATCTCCTGCCCCCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGTTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCGTTGACATCGTATCCCGCGAGGATGCCGATAATCAGCCCCATCACAAAAccc          <  2:58593/236‑1 (MQ=255)
                                                                                                                                                                                                             cG‑ATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGACAATCAGGAACTCCTGAACGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGGGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATAACAAAACCCACgg      >  1:103308/1‑248 (MQ=255)
                                                                                                                                                                                                                 gataccaccgccgTATGTCCCAGCAGCAAAGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTCAACAATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGAc  >  1:29167/8‑251 (MQ=255)
                                                                                                                                                                                                                        ||                                                                                                                                                                                                                                                  
AACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCG‑ATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGAC  >  NC_000913/2173148‑2173606

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑