Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,773,0400G.71.4% 19.5 / 11.3 14coding (760/1914 nt)mtlAmannitol‑specific PTS enzyme IICBA component
Reads supporting (aligned to +/- strand):  ref base G (1/3);  new base . (4/6);  total (5/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGTCTATCTTTGTTGAACCGGCGAAAATCCTGTTCCTCAACAACGCCATTAACCACGGTATCTTCTCGCCGCTGGGTATTCAGCAGTCCCATGAACTGGGTAAATCAATCTTCTTCCTGATTGAAGCTAACCCAGGTCCAGGTATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGGTTCTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGC  >  NC_000913/3772810‑3773272
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                        
cGTCTATCTTTGTTGAACCGGCGAAAATCCTGTTCCTCAACAACGCCATTAACCACGGTATCTTCTCGCCGCTGGGTATTCAGCAGTCCCATGAACTGGGTAGATCAATCTTCTTCCTGATTGAAGCTAACCCAGGTCCAGGTATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGG‑‑TATCCACGAAATCTACTTccc                                                                                                                                                                                                                    <  1:437189/251‑1 (MQ=255)
       cTTTGTTGAACCGGCGAAAATCCTGTTCCTCAACAACGACACTAACCACCGTATCTTATCGCCGCTGGGCCTTCAGCAGCCCAATCAACTGGGTAAATCAATCTTAATCCTGATTGACGATAAACCAGGTCCAGGTATGCGCGTGCTGCTGGCGTACCTGTTATTTGGTCGTGGTAGCGCTAAACAGTCTGCGGCCGATGCGGCACTAATCCACTTC‑ACGGGGGTATCCACGAAATCTACTTCCCGTAtgt                                                                                                                                                                                                              >  1:399548/1‑251 (MQ=255)
                                   cTCAACAACGCCATTAACCACGGTATCTTCTCGCCGCTGGGTATTCAGCAGTCACATGAACTGGGTAAATCAATCTTCTTCCTGATTGAAGCTAACCCAGGTCCAGGTATGGGCGTGCTGCTGGCGTACCTGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGG‑‑TATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGc                                                                                                                                                                                   >  1:220108/1‑249 (MQ=255)
                                            gCCATTAACCACGGTATCTTCTCGCCGCTGGGTATTCAGCAGTCCCATGAACTGGGTAAATCAATCTTCTTCCTGATTGAAGCTAACCCAGGTCCAGGTATGGGCGTGCGGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAAGCATCCACTTCC‑‑GGGGGTAGCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCg                                                                                                                                                                                    <  1:176807/239‑1 (MQ=255)
                                                                                cagcagTCCCATGAACTGGGTAAATCAATCTTCTTCCTGATTGAAGCTAACCCAGGTCCAGGTATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCTGGCTGTGCGGCAATCATCCACTTCC‑‑GGGGGTATCCACGAAATCTACTTCCCTTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCgtgt                                                                                                                                                      <  1:245434/233‑1 (MQ=255)
                                                                                                                                    caggtccaggtATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGG‑‑TATCCACGAAATCTACTGCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCg                                                                                                                                                                                    <  1:7432/151‑1 (MQ=255)
                                                                                                                                                               cGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGAGCGGTGCGGCAATCATCCACTTACTGGGG‑‑TATACACGCAATCTACTTCCCGTATGTGCTGCTGAATCCGCGTCTGATCCTCGCAGTCCACCTCGGCGGTATGACTGGCGTGTTCCCGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGGTTCTATCCTTGCTGTACTGGCGATGACACCAAAAAGTGCTTACTTCGc                                                     >  1:230601/1‑251 (MQ=255)
                                                                                                                                                                  aNATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGG‑‑TATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGGTTCTATCCTTGCTGTACTGGCGATGACACCCAAAGGTGCTTACTTCGCTa                                                   >  1:349315/1‑250 (MQ=255)
                                                                                                                                                                                    gTAGCGCTCAACAGTCTGCGGGCGGTGCGGCAATCATACCCTTCCTGGGG‑‑TATCCCCGAAATCTACTTCGCGCCTGTGCTGCTGACTCCGCGTATGATCATCGccg                                                                                                                                                                                 >  1:446219/1‑104 (MQ=255)
                                                                                                                                                                                                           ggTGCGGCAATCATCCACTTCC‑‑GGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTc                                                                                                                                                    <  1:392151/112‑1 (MQ=255)
                                                                                                                                                                                                            gTGCGGCAATCATCCACTTCCTGGGG‑‑TATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGc                                                                                                              <  1:154032/149‑1 (MQ=255)
                                                                                                                                                                                                               cGGCAATCATCCACTTCCTGGGG‑‑TATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATcc                                                                                                                                       <  1:389176/121‑1 (MQ=255)
                                                                                                                                                                                                                  cAATCAGCCACTGCCTGGGG‑‑TATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGGGTCTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGc  <  1:208009/251‑1 (MQ=255)
                                                                                                                                                                                                                          cACTTCCTGGGG‑‑TATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGGTTCTATCCTTgn                                                                                         <  1:114045/156‑2 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                        
CGTCTATCTTTGTTGAACCGGCGAAAATCCTGTTCCTCAACAACGCCATTAACCACGGTATCTTCTCGCCGCTGGGTATTCAGCAGTCCCATGAACTGGGTAAATCAATCTTCTTCCTGATTGAAGCTAACCCAGGTCCAGGTATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGGTTCTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGC  >  NC_000913/3772810‑3773272

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