Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,726,6400TA52.4% 9.2 / 28.0 22N360K (AAT→AAAxylGD‑xylose ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (1/9);  new base A (11/0);  total (12/10)
Fisher's exact test for biased strand distribution p-value = 3.40e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.06e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CAAGCCATCGCCCAGGGGATTGCGATGGTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGTAAAAATATTACCCTCGCCGCACTCAATAAATTTACCGGTGGCATTAGCCAGCTTGATGACGCGGCAGAGCAAAAATGTATTC  >  NZ_CP009273/3726557‑3726716
                                                                                   |                                                                            
cAAGCCATCGCCCAGGGGATTGCGATGGTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGTAAAAATATTAcc                                                                        <  2:538714/90‑1 (MQ=255)
         gCCCAGGGGATTGCGATGGTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGTAAAAATATTACCCTCGCCGCa                                                               <  1:281200/90‑1 (MQ=255)
          cccAGGGGATTGCGATGGTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGTAAAAATATTACCCTCGCCGCAc                                                              <  1:363250/90‑1 (MQ=255)
                   ttGCGATGGTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGTAAAAATATTACCCTCGCCGCACTCAATAAAt                                                     >  1:289914/1‑90 (MQ=255)
                   ttGCGATGGTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGTAAAAAAATTACCCCCGCCGCACTCAATAAAt                                                     >  1:616950/1‑90 (MQ=255)
                   ttGCGATGGTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGGAAAAAAATTACCCCCGCCGCACTCAATaaaa                                                     >  1:723133/1‑89 (MQ=255)
                   ttGCGATGGTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGGAAAAAAATTACCCCCCCCGCACCCAATaaaa                                                     >  1:37827/1‑89 (MQ=255)
                    tGCGATGGTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGGAAAAAAATTAACCTCGCCCCCCCCaaaaaaat                                                    >  1:752316/1‑88 (MQ=255)
                    tGCGATGGTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGGTGGGGAAAAAAATAACCCCGCCCCCCCCAAAAAAtt                                                    >  2:596382/1‑90 (MQ=255)
                          ggTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTTGTAAAAAAATTACCCCCCCCCCCCCCAATAAAATTACCgg                                              >  1:538683/1‑90 (MQ=255)
                           gTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGTAAAAATATTACCCTCGCCGCACTCAATAAATTTACCGGt                                             <  2:380906/90‑1 (MQ=255)
                             ccccGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGGAAAAAAATTACCCCCCCCCCCCTCCAAAAAATTTCCCGTgg                                           >  1:385334/1‑90 (MQ=255)
                                              gcgACGGCATCGTTCCGGTAATGGCGGTTGGGAAAAAAAATAACCCCCCCCCCCCCAAAAAAATTACCCGGGGGATTAACCAGCTtgatg                          >  2:210309/1‑90 (MQ=255)
                                                gACGGCATCGTTCCGGTAATGGCGGTTGGTAAAAACATTACCCTCGCCGCACTCAATAAATTTACCGGTGGCATTAGCCAGCTTGATGAc                        <  2:538683/90‑1 (MQ=255)
                                                  cGGCATCGTTCCGGTAATGGCGGTTGGTAAAAATATTACCCTCGCCGCACTCAATAAATTTACCGGTGGCATTAGCCAGCTTGATGAcgc                      <  2:200323/90‑1 (MQ=255)
                                                       tCGTTCCGGTAATGGCGGTTGGTAAAAATATTACCCTCGCCGCACTCAATAAATTTACCGGTGGCATTAGCCAGCTTGATGACGCGGCag                 <  1:225610/90‑1 (MQ=255)
                                                            ccGGTAATGGCGGTTGGTAAAAATATTACCCTCGCCGCACTCAATAAATTTACCGGTGGCATTAGCCAGCTTGATGACGCGGCAGAGCaa            <  2:385926/90‑1 (MQ=255)
                                                             cGGTAATGGCGGTTGGGAAAAAAATTACCCCCCCCCCCCTCAATAAATTTACCGGGGGGATTAGCCAGCTTGATGACGCGGCCGAGCaaa           >  2:295464/1‑90 (MQ=255)
                                                              ggTAATGGCGGTTGGGAAAAAAATTACCCCCGCCCCCCTCCATAAATTTTCCGGGGGGATTAGCCAGCTTGTTGACGCGGCAGAGCaaaa          >  2:175980/1‑90 (MQ=255)
                                                                   tGGCGGTTGGTAAAAATATTACCCTCGCCGCACTCAATAAATTTACCGGTGGCATTAGCCAGCTTGATGACGCGGCAGAGCAAAAATGTa     <  2:176495/90‑1 (MQ=255)
                                                                   tGGCGGTTGGTAAAAAAATTACCCCCCCCCCCCTCAATAAAATTACCGGGGGGATTAACCAGCTTGATGACGCGGCAGAGGaaaaaaata     >  2:266339/1‑86 (MQ=255)
                                                                      cGGTTGGTAAAAATATTACCCTCGCCGCACTCAATAAATTTACCGGTGGCATTAGCCAGCTTGATGACGCGGCAGAGCAAAAATGTATTc  <  1:266339/90‑1 (MQ=255)
                                                                                   |                                                                            
CAAGCCATCGCCCAGGGGATTGCGATGGTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGTAAAAATATTACCCTCGCCGCACTCAATAAATTTACCGGTGGCATTAGCCAGCTTGATGACGCGGCAGAGCAAAAATGTATTC  >  NZ_CP009273/3726557‑3726716

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