Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,726,6480TC55.0% 7.4 / 20.0 20L363P (CTC→CCC) xylGD‑xylose ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (3/6);  new base C (11/0);  total (14/6)
Fisher's exact test for biased strand distribution p-value = 2.17e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.01e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CATCGCCCAGGGGATTGCGATGGTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGTAAAAATATTACCCTCGCCGCACTCAATAAATTTACCGGTGGCATTAGCCAGCTTGATGACGCGGCAGAGCAAAAATGTATTCTGGAATCAATCCAGCA  >  NZ_CP009273/3726562‑3726732
                                                                                       |                                                                                    
cATCGCCCAGGGGATTGCGATGGTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTAGGGAAAAAAATTATCCCcgc                                                                                   >  1:698930/1‑90 (MQ=255)
   cGCCCAGGGGATTGCGATGGTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGTAAAAATATTACCCTcgccgc                                                                                <  1:873707/90‑1 (MQ=255)
    gCCCAGGGGATTGCGATGGTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGTAAAAATATTACCCTCGCCGCa                                                                               <  2:995430/90‑1 (MQ=255)
              ttGCGATGGTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGTAAAAAAATTACCCTCGCCGCACTCAATAAAt                                                                     >  2:348242/1‑90 (MQ=255)
              ttGCGATGGTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGGTGGTAAAAAAAATAACCCCGCCGCCCCCAATAAAt                                                                     >  1:357369/1‑90 (MQ=255)
              ttGCGATGGTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGGTGGGGAAAAAAATACCCCCCCCCCCCCCAATaaaa                                                                     >  1:901984/1‑89 (MQ=255)
                     ggTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGTAAAAAAATTACCCCCCCCCCCCTCAATAAATTTACCgg                                                              >  2:657513/1‑90 (MQ=255)
                     ggTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGGTGGGAAAAAAATTAACCCCGCCCCCCTCAAAAAAATTACCgg                                                              >  1:346492/1‑90 (MQ=255)
                      gTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGTAAAAATATTACCCTCGCCGCACTCAATAAATTTACCGGt                                                             <  2:815196/90‑1 (MQ=255)
                      gTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGTAAAAATATTACCCTCGCCGCACTCAATAAATTTACCGGt                                                             <  1:445203/90‑1 (MQ=255)
                          ccGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGTAAAAAAAATACCCCCCCCCCCCCCAAAAAAATTACCGGTGGCa                                                         >  1:304693/1‑90 (MQ=255)
                                aCAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGTAAAAATATTACCCTCGCCGCACTCAATAAATTTACCGGTGGCATTAGcc                                                   <  1:775816/90‑1 (MQ=255)
                                aCAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGTAAAAAAATTACCCCCCCCCCCCCCAAAAAAATTTCCCGGGGGGTTAGcc                                                   >  1:210499/1‑90 (MQ=255)
                                  aGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGGTGGTAAAAAAATTACCCCCCCCCCCCTCAAAAAAATTTTCGGGGGGATTAGCCAg                                                 >  1:112724/1‑90 (MQ=255)
                                   gAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGGAAAAAAAATAACCCCCCCCCCCCCAAAAAAATTAACCGGGGGCTTAGCCAGc                                                >  1:193386/1‑90 (MQ=255)
                                   gAAAGCGCGACGGCATCGTTCCGGTAATGGCGGGTGGGAAAAAAATTACCCCCCCCCCCCCCCATAAAATTTACGGGGGGATTAGCCAGc                                                >  1:576323/1‑90 (MQ=255)
                                                        cGGTAATGGCGGTTGGTAAAAATATTACCCTCGCCGCACTCAATAAATTTACCGGTGGCATTAGCCAGCTTGATGACGCGGCAGAGCaaa                           <  2:30425/90‑1 (MQ=255)
                                                              tGGCGGTTGGTAAAAATATTACCCTCGCCGCACTCAATAAATTTACCGGGGGGATTAGCCAGCTTGATGACGCGGGAGAGCAAAAAttga                     >  2:162508/1‑87 (MQ=255)
                                                                   gTTGGTAAAAAAATTACCCCCCCCGCACTCAATAAATTTACCGGGGGGATTAGCCAGCTTGATGACGCGGCAGAGCAAAAAAGTATTTTg                >  1:325850/1‑90 (MQ=255)
                                                                                  tACCCTCGCCGCACTCAATAAATTTACCGGTGGCATTAGCCAGCTTGATGACGCGGCAGAGCAAAAATGTATTCTGGAATCAATCCAGCa  >  1:139459/1‑90 (MQ=255)
                                                                                       |                                                                                    
CATCGCCCAGGGGATTGCGATGGTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGTAAAAATATTACCCTCGCCGCACTCAATAAATTTACCGGTGGCATTAGCCAGCTTGATGACGCGGCAGAGCAAAAATGTATTCTGGAATCAATCCAGCA  >  NZ_CP009273/3726562‑3726732

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