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

ATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACG  >  NC_000913/2173125‑2173601
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                              
attggccgccTGGTGGGTGTGAAAACCAAGAGTTTGCGTGGCGATCCACAGGGTGATGCTCATAATGAAGACACCGGAGATTAAGGTGCGGAACAGAGTCCCCCGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGAGCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATAAAAATGGTGAGTGGGATAAAAATCAGGCTCGCCGATACCACC‑‑GCCGTATGTCCca                                                                                                                                                                                                                                  <  2:131755/247‑1 (MQ=255)
            ttGGTTGTGAAAACCAATAATTGGCGTAGCGATCCCCAGTTTGATGCTCAAAAGGCTGACCCCCGACATTAAGGTGCGAAACACATTCCCCCGATGCACGGCGACCGCCATCACCACGAATAAGCCAATGATGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACG‑GCAATTAAAATGTTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGcc                                                                                                                                                                                                                       <  2:307/250‑1 (MQ=255)
                                      ccgcgATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGAGGCACGGCGCCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGa                                                                                                                                                                                            <  1:227845/250‑1 (MQ=255)
                                             ccACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACACACA‑GCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgccg                                                                                                                                                                                      >  2:111897/1‑250 (MQ=255)
                                                     tGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGcc                                                                                                                                                                             >  1:39204/1‑251 (MQ=255)
                                                                        ccccGAGATTAAGGTGCGGAACAGATTTCCACGAGGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCt                                                                                                                                                          <  1:4211/250‑1 (MQ=255)
                                                                                  aaggggcggacaataTTTCCACGAGGCACGGGGACCGCCACCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                                                                                                                                <  1:9427/237‑1 (MQ=255)
                                                                                        gggaagagttttCCACGATGCAGGGGGACCGCCACCGCCACGAAGAAGCCTATGGTGGCAAGAGCGCCAAACGGCGGCGCCTGATTACCCGGCACACAAACA‑GCAATTAAAATGGGGAGGGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCGGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGt                                                                                                                                          <  1:215759/242‑1 (MQ=255)
                                                                                              aGATTCCCACGATGCACGCCGCCCGCCATCGCACCGAAGAAGCCGATGTTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGTCCCACAAACC‑GCAATTAAAATGGTGATGGGGATAAAACTCAGGCTTGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGTCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAcc                                                                                                                                    <  1:393253/251‑1 (MQ=255)
                                                                                                     cccgagccACGCCGACCGCAACCGCCAGGAAAAAGCCGATGGTGGAAAGATCCCCAAACGGCAGCACCTGATTCCCCGGCACACAAACA‑GAAATTAAAATGGTGAGTGGGATAAAAATCAGCCTTGCCAATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatg                                                                                                                             <  2:347801/244‑1 (MQ=255)
                                                                                                                           agCCACGAAAACGCCGATGGGGGCATGATCCCCAAACGGCAGCACCTGAGTACCCGGCACAAAAAGA‑GCAATTAAAAGGGTAAGTGGAATAAAAGACATTTGTCACGATCCGACC‑‑GCCGTATGTCCAAGCAGCAACGCCGGGTCAAGGCCAATAAGGAACTCCTGACCGCCGAACTTCGCCTGGAAGCGGCTACTAGCCTGCGTAGCGACGTGCGTTAAACCAGCTAGGATGGGTTTAAGCACCCGTGGCa                                                                                                       <  2:406306/250‑1 (MQ=255)
                                                                                                                                 gaagaaGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGGCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGCACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTACTCACCCGTGGCATTagc                                                                                                  >  2:131558/1‑250 (MQ=255)
                                                                                                                                         ggATGGTGGCAAGATCGCAAAACGGCACCACCTGATTCCCCGGCAAACAAACA‑ACAATTAAAATGTGGAGTGGAAAAAAAACAAGGCTGGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACt                                                                                         <  2:143025/250‑1 (MQ=255)
                                                                                                                                           aTGGTGGCAAGTTCGCCACACGGCAGAACCTGATTACCCGGCACACAAACA‑GCAATTAAAATGGCGCGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCCTACATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgc                                                                                       >  2:157627/1‑251 (MQ=255)
                                                                                                                                                     gATCGCAAGACGGCGGCACCTGAT‑ACCCCGGACCCCAACAAGCAATTAAAATGGTGGGTGGTATAAAAATCAGGTTTGC‑GAATCCACCC‑GCCGTATGTCCCAGCAGAAACGCCGTATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGATGACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGTGTTTAATTACCCGTGGCATTAGCAGCATCACTGCCGCCGttttt                                                                             <  2:294674/251‑1 (MQ=255)
                                                                                                                                                           cAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATGGTGAGTGGGATAAAACTCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCACGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCAACCATGATGGGTTTAATCACCCGTGGCATTAGCAGAATCACTGCCGCCGTTTTTACCGcc                                                                       >  2:96926/1‑251 (MQ=255)
                                                                                                                                                              aCGGCAGCACCTGATTACCCGGCACACACACA‑GCAATTAAAATGGTGAGTGGGATAACAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTACTCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGc                                                                    >  2:211542/1‑251 (MQ=255)
                                                                                                                                                                         tGATTACCCGGCACACAAACA‑GCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATA‑‑ACACCGCCGTATGTCACAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACAGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTACCGCCGTTTTTCACGCCAGCTGCAATACAcc                                                         >  2:216523/1‑251 (MQ=255)
                                                                                                                                                                          gATTACCCGGCACACAAACA‑GCAATTAAAATGGAGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCACTCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCt                                                        >  1:244764/1‑251 (MQ=255)
                                                                                                                                                                                           ccA‑GCAATTAAAAGGGTGAGTGGGATAAAAACCAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCcgc                                       <  1:131558/250‑1 (MQ=255)
                                                                                                                                                                                            cA‑GCAATTAAATTGGTGCGTGCGATAAAAATCAGGCTTGCCGATACCACC‑‑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)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                              
ATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACG  >  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: ‑