Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,173,3640C.79.2% 54.6 / 17.4 24intergenic (‑491/+918)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (1/4);  new base . (11/8);  total (12/12)
Fisher's exact test for biased strand distribution p-value = 3.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACG  >  NC_000913/2173125‑2173601
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             
attggccgccTGGTGGGTGTGAAAACCAAGAGTTTGCGTGGCGATCCACAGGGTGATGCTCATAATGAAGACACCGGAGATTAAGGTGCGGAACAGAGTCCCCCGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGAGCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATAAAAATGGTGAGTGGGATAAAAATCAGGCTCGCCGATACCACC‑‑GCCGTATGTCCca                                                                                                                                                                                                                                  <  2:131755/247‑1 (MQ=255)
            ttGGTTGTGAAAACCAATAATTGGCGTAGCGATCCCCAGTTTGATGCTCAAAAGGCTGACCCCCGACATTAAGGTGCGAAACACATTCCCCCGATGCACGGCGACCGCCATCACCACGAATAAGCCAATGATGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACGGCAATTAAAATGTTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGcc                                                                                                                                                                                                                       <  2:307/250‑1 (MQ=255)
                                      ccgcgATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGAGGCACGGCGCCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGa                                                                                                                                                                                            <  1:227845/250‑1 (MQ=255)
                                             ccACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACACACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgccg                                                                                                                                                                                      >  2:111897/1‑250 (MQ=255)
                                                     tGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGcc                                                                                                                                                                             >  1:39204/1‑251 (MQ=255)
                                                                        ccccGAGATTAAGGTGCGGAACAGATTTCCACGAGGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCt                                                                                                                                                          <  1:4211/250‑1 (MQ=255)
                                                                                  aaggggcggacaataTTTCCACGAGGCACGGGGACCGCCACCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                                                                                                                                <  1:9427/237‑1 (MQ=255)
                                                                                        gggaagagttttCCACGATGCAGGGGGACCGCCACCGCCACGAAGAAGCCTATGGTGGCAAGAGCGCCAAACGGCGGCGCCTGATTACCCGGCACACAAACAGCAATTAAAATGGGGAGGGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCGGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGt                                                                                                                                          <  1:215759/242‑1 (MQ=255)
                                                                                              aGATTCCCACGATGCACGCCGCCCGCCATCGCACCGAAGAAGCCGATGTTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGTCCCACAAACCGCAATTAAAATGGTGATGGGGATAAAACTCAGGCTTGCCGATCCCCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGTCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAcc                                                                                                                                    <  1:393253/251‑1 (MQ=255)
                                                                                                     cccgagccACGCCGACCGCAACCGCCAGGAAAAAGCCGATGGTGGAAAGATCCCCAAACGGCAGCACCTGATTCCCCGGCACACAAACAGAAATTAAAATGGTGAGTGGGATAAAAATCAGCCTTGCCAATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatg                                                                                                                             <  2:347801/244‑1 (MQ=255)
                                                                                                                           agCCACGAAAACGCCGATGGGGGCATGATCCCCAAACGGCAGCACCTGAGTACCCGGCACAAAAAGAGCAATTAAAAGGGTAAGTGGAATAAAAGACATTTGTCACGATCCGACC‑‑GCCGTATGTCCAAGCAGCAACGCCGGGTCAAGGCCAATAAGGAACTCCTGACCGCCGAACTTCGCCTGGAAGCGGCTACTAGCCTGCGTAGCGACGTGCGTTAAACCAGCTAGGATGGGTTTAAGCACCCGTGGCa                                                                                                       <  2:406306/250‑1 (MQ=255)
                                                                                                                                 gaagaaGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGGCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGCACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTACTCACCCGTGGCATTagc                                                                                                  >  2:131558/1‑250 (MQ=255)
                                                                                                                                         ggATGGTGGCAAGATCGCAAAACGGCACCACCTGATTCCCCGGCAAACAAACAACAATTAAAATGTGGAGTGGAAAAAAAACAAGGCTGGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACt                                                                                         <  2:143025/250‑1 (MQ=255)
                                                                                                                                           aTGGTGGCAAGTTCGCCACACGGCAGAACCTGATTACCCGGCACACAAACAGCAATTAAAATGGCGCGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCCTACATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgc                                                                                       >  2:157627/1‑251 (MQ=255)
                                                                                                                                                     gATCGCAAGACGGCGGCACCTGATACCCCGGACCCCAACAAGCAATTAAAATGGTGGGTGGTATAAAAATCAGGTTTGCGAATCCACCCGCCGTATGTCCCAGCAGAAACGCCGTATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGATGACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGTGTTTAATTACCCGTGGCATTAGCAGCATCACTGCCGCCGttttt                                                                             <  2:294674/251‑1 (MQ=255)
                                                                                                                                                           cAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAACTCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCACGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCAACCATGATGGGTTTAATCACCCGTGGCATTAGCAGAATCACTGCCGCCGTTTTTACCGcc                                                                       >  2:96926/1‑251 (MQ=255)
                                                                                                                                                              aCGGCAGCACCTGATTACCCGGCACACACACAGCAATTAAAATGGTGAGTGGGATAACAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTACTCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGc                                                                    >  2:211542/1‑251 (MQ=255)
                                                                                                                                                                         tGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATA‑‑ACACCGCCGTATGTCACAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACAGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTACCGCCGTTTTTCACGCCAGCTGCAATACAcc                                                         >  2:216523/1‑251 (MQ=255)
                                                                                                                                                                          gATTACCCGGCACACAAACAGCAATTAAAATGGAGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCACTCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCt                                                        >  1:244764/1‑251 (MQ=255)
                                                                                                                                                                                           ccAGCAATTAAAAGGGTGAGTGGGATAAAAACCAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCcgc                                       <  1:131558/250‑1 (MQ=255)
                                                                                                                                                                                            cAGCAATTAAATTGGTGCGTGCGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGCTCAAGGCAAATCAGGAACTCCTGACCGCCGAAATTCGCCTGTAAACGACTAAGAGCCTGCTTAGCGATGGGCGTTAAACCATCAATGATGGGTTTAATCACCCGTGGAATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGATGCAATACACCTTTTACATCGTAACAcgcg                                      >  2:287804/1‑251 (MQ=255)
                                                                                                                                                                                                             gAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATc                       >  1:158542/1‑250 (MQ=255)
                                                                                                                                                                                                                         caTCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCCGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATcaca           >  2:79180/2‑250 (MQ=255)
                                                                                                                                                                                                                                 ttGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACg  >  2:107074/1‑251 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             
ATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACG  >  NC_000913/2173125‑2173601

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