Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,113,6490TG61.1% 12.4 / 18.6 18K330N (AAA→AACglcAglycolate permease GlcA
Reads supporting (aligned to +/- strand):  ref base T (2/5);  new base G (3/8);  total (5/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.29e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCTTTCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCTTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGACTATA  >  NZ_CP009273/3113566‑3113735
                                                                                   |                                                                                      
gCTTTCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCTTTAAAc                                                                                  <  2:186658/90‑1 (MQ=255)
        cacTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCAt                                                                          <  1:276092/90‑1 (MQ=255)
        cacTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCAt                                                                          <  1:338751/90‑1 (MQ=255)
        cacTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCAt                                                                          <  2:319629/90‑1 (MQ=255)
         acTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCTTTAAACGGCTTCATg                                                                         <  1:143177/90‑1 (MQ=255)
         acTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCTTTAAACGGCTTCATg                                                                         <  1:28899/90‑1 (MQ=255)
         acTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATg                                                                         <  2:40423/90‑1 (MQ=255)
              ttgATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCa                                                                    <  2:234406/90‑1 (MQ=255)
                           gggATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCTTTAAACGGCTTCATGGTCCAGATGGTGACCAGc                                                       >  2:212030/1‑90 (MQ=255)
                                  agAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTa                                                >  1:328294/2‑90 (MQ=255)
                                                   gAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGAc                               >  2:263344/1‑90 (MQ=255)
                                                                  cccGGAGCAAATAACGCTTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGAtt                <  2:58110/90‑1 (MQ=255)
                                                                  cccGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGAtt                <  2:317934/90‑1 (MQ=255)
                                                                      gAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGa                                             <  2:101410/57‑1 (MQ=255)
                                                                      gAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGa                                             >  1:101410/1‑57 (MQ=255)
                                                                       aGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCCCCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGccc           <  2:313163/90‑1 (MQ=255)
                                                                            aTAACGCTTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGAc      <  1:250812/90‑1 (MQ=255)
                                                                                cGCTTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGACtata  >  2:198704/1‑90 (MQ=255)
                                                                                   |                                                                                      
GCTTTCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCTTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGACTATA  >  NZ_CP009273/3113566‑3113735

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: