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.93.1% 113.7 / 3.5 29intergenic (‑490/+919)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base . (16/11);  total (16/13)
Fisher's exact test for biased strand distribution p-value = 1.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01
*NC_0009132,173,3640C.93.1% 114.9 / 2.1 29intergenic (‑491/+918)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base . (16/11);  total (16/13)
Fisher's exact test for biased strand distribution p-value = 1.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

CCAGTTGGGTGTGCAAACCAATAGTTTG‑CGTCGCGATCCACAGGGTG‑ATGCTCATAATGATGACACCCG‑AGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACA  >  NC_000913/2173133‑2173607
                                                                                                                                                                                                                                         ||                                                                                                                                                                                                                                                   
cgaGTTGGGTG‑GGAAACCAAAAATTTGGCTTCGAGATCCACAGGGTG‑AGCCTCATAAGGATGCCACCGG‑AGATAAAGGTGCGAAAAAGTTTCCCACAAGGACCGGCACCCGCCACCGCCACGAAGAAGCCGATGGGGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAAc                                                                                                                                                                                                                                <  2:79650/249‑1 (MQ=255)
  aggggggggggaaacccAATTGTTTG‑GTTCGCAATCCAC‑GGGTGTATGCTCATAATGATGACACCCG‑AGATTAAGGGGCGGAACAGATTCCCACGATGCACGGCGCCCGCCATCGCCCCGAAGAAGCCGATGGTGGCAAGATCGCAAAACGGCAGCCCCTGATTCCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCATATACCACC‑‑GCCGTATGTCCCAGCAGCAACGc                                                                                                                                                                                                                              <  2:346865/236‑1 (MQ=255)
     tGGGTGTGCAAACCAATAGTTTG‑CGTCGCGATCCACAGGGGG‑ATGCTCATAATGATGACACCCG‑AGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACCCAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCACAGCAGCAACGCCgg                                                                                                                                                                                                                           >  1:128704/1‑251 (MQ=255)
              aaaCCAATAGTTTG‑CGTCGCGATCCACAGGGTG‑ATGCTCATAATGATGACACCCG‑AGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGcc                                                                                                                                                                                                                  >  2:240689/1‑251 (MQ=255)
                 ccAATATTTTG‑CGTCGCGATCCACAGGGTG‑ATGCTCATAATGATGACACCCG‑AGATTAAGGTGCGGAACAGCTTTCCACGAGGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAACCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAt                                                                                                                                                                                                               <  1:111161/251‑1 (MQ=255)
                         ttG‑CGTCGCGATCCACAGGGGG‑AGGCTCATAATGATGACCCCCG‑AGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAAc                                                                                                                                                                                                        <  1:26782/250‑1 (MQ=255)
                          gg‑CGTGGCGACCCACAGGGGG‑ATGATCATAATGAGGAC‑CCCGTAGATTAAGGTGGGGACAAGATCTCCAGGATGCCCGGCGCCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTCCCCGGCACACAAACAGCAATTAAAATGGTGAGGGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTc                                                                                                                                                                                                      <  2:124904/250‑1 (MQ=255)
                              gTCGCGATCCACAGGGTG‑ATGCTCATAATGATGACACCCG‑AGATTTAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCACACGGCAGCACCTGATTACCCGGCACACAAACCGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCt                                                                                                                                                                                                    >  1:38075/1‑250 (MQ=255)
                                        acaGGGTG‑ATGCTCATAATGATGACACCCG‑AGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCCGGAACTCCTGACCGCCGAAc                                                                                                                                                                                         >  2:273357/1‑251 (MQ=255)
                                                   gCTCAAAATGATGCCCCCCG‑AGATAAAGGTGCGGAACAGATTCCCAAGAGGCACGGCGCCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGGGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTa                                                                                                                                                                               <  1:141383/251‑1 (MQ=255)
                                                    cTCATAATGATGACACCCG‑AGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTa                                                                                                                                                                               >  1:505585/1‑250 (MQ=255)
                                                      caaaatgaggccccccG‑AGATTAGGGTGCGGAACAGATTCCCAGGATGCACGGCGACCGCCAGCGCCACGAAGAAGCCGAGGGTGGCAAGATCGCAAACCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGTAACTCCTGACCGCCGAACTTCGCCTGTAAAc                                                                                                                                                                            <  1:422972/238‑1 (MQ=255)
                                                                     ccgaaataaaGGTGCGGAACAGATTCCCACGATGCACGGCGACCGCCAGCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCCCCTGATTCCCCGGCACACAAACAGCAATTAAAATGGTGAGGGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGc                                                                                                                                                              <  1:46187/243‑1 (MQ=255)
                                                                          aTTAAGGTGCGGAACAGATTGCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGAACGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATc                                                                                                                                                                                                                        <  1:468056/188‑1 (MQ=255)
                                                                                              gccACGATGCACGGCGAACGCCATCGCCACGACGAAGCCGATGGAGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTACCAGCAGAAACGCCGGATCAAGGCCAATCAGGCACTCCTGACCGCCGAACTCCGCCTGAAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAAccatcc                                                                                                                                      >  2:300664/2‑250 (MQ=255)
                                                                                                                 gCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCg                                                                                                                  >  1:198905/1‑251 (MQ=255)
                                                                                                                           cGAAGAAGCCGAGGGGGGCAAGATCGCCAACCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGATGGGGATAAAAACCAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagc                                                                                                        <  1:240689/251‑1 (MQ=255)
                                                                                                                            aaagaagCCGATGGTGGAAAGATCGCCAAACGGCAGCACCGGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGAAAAAAATCAGGCTGGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagca                                                                                                       <  1:226668/250‑1 (MQ=255)
                                                                                                                               aaagCGGATGGGGGCAAGACCGCCAAACGGCAGCCCCTGATTACCCGGAACACAAACGGCAATTAAAATGGTGAGTGGGATAAAAACCAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAgcagca                                                                                                    <  1:387004/250‑1 (MQ=255)
                                                                                                                                          tggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgccgcc                                                                                         >  1:339950/1‑251 (MQ=255)
                                                                                                                                                   cGCCAAACGGCAGCACCTGTTTACCCGGCACACAAACAGAAATTAAATTGGTGAGTGGGATAAAAATCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTAcc                                                                                <  1:136593/251‑1 (MQ=255)
                                                                                                                                                              agcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCa                                                                      >  1:183961/1‑250 (MQ=255)
                                                                                                                                                                   cTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATacac                                                                >  2:250643/1‑251 (MQ=255)
                                                                                                                                                                   cTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCACCGTTTTTACCGCCAGCTGCAATacac                                                                >  2:338213/1‑251 (MQ=255)
                                                                                                                                                                              gCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCg                                                     >  2:143979/1‑251 (MQ=255)
                                                                                                                                                                                     aaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCcg                                              >  1:275547/1‑251 (MQ=255)
                                                                                                                                                                                                  aTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGAt                                 >  2:196662/1‑251 (MQ=255)
                                                                                                                                                                                                          tGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGcc                         >  2:327553/1‑251 (MQ=255)
                                                                                                                                                                                                                    aTCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACaa                >  1:53438/1‑250 (MQ=255)
                                                                                                                                                                                                                                   ataccaccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACa  <  1:273357/245‑1 (MQ=255)
                                                                                                                                                                                                                                         ||                                                                                                                                                                                                                                                   
CCAGTTGGGTGTGCAAACCAATAGTTTG‑CGTCGCGATCCACAGGGTG‑ATGCTCATAATGATGACACCCG‑AGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACA  >  NC_000913/2173133‑2173607

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