Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,377,9790GA47.8% 11.8 / 13.6 23F80F (TTC→TTTycjWLacI family DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (9/3);  new base A (0/11);  total (9/14)
Fisher's exact test for biased strand distribution p-value = 3.37e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAG  >  NZ_CP009273/1377907‑1378068
                                                                        |                                                                                         
gtgCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACccgc                                                                          >  2:488677/1‑90 (MQ=255)
  gCGACGTTTGTAATATCACGTCAAAACCTTCTTCTTCAGCTTGGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACccgccg                                                                        <  2:132754/90‑1 (MQ=255)
                atCACGTCAAAACCTTCTTCTTCCGCTTGGGCGGTGATGGAATAAAAAATTAAAAAAAAAAACGGATTCCCCGCCGTAGTTTTGGTCGaa                                                          <  2:1208083/90‑1 (MQ=255)
                accACGTCAAACCTTTCTTCTTCAGTTTGGGGGGTGATGGCATGAAAAACTTAAGAAAAAAACGGATTCCCCGCCGTAGTTTTGGTCGaa                                                          <  2:468849/88‑1 (MQ=255)
                accACGTCAAACCTTTCTTCTTCAGCTTTGGCGGTGATGGCATGAAAAACTTCAGAAAAAAACGGATTACCCGCCGTAGTTTTGGTCGaa                                                          <  1:669087/88‑1 (MQ=255)
                 tCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAAc                                                         >  1:767499/1‑90 (MQ=255)
                  cACGTCAAACCTTTCTTCTTCAGCTTGGGGGGTGATGGCATAAAAACCTTAAAAAAAAAAGGTATTCCCCGCCGTAGTTTTGGTCGAACg                                                        <  2:1057915/90‑1 (MQ=255)
                   aCGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACgg                                                       >  2:809982/1‑90 (MQ=255)
                   aCGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACgg                                                       >  2:785540/1‑90 (MQ=255)
                   aCGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACgg                                                       >  2:576759/1‑90 (MQ=255)
                     gTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGt                                                     >  1:1029663/1‑90 (MQ=255)
                          aaCCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCTGTAGTTTTGGTCGAACGGGTAGaaa                                                <  1:304589/90‑1 (MQ=255)
                            ctttcttcttcAGCTTGGGGGGTGATGGCATGAAAAACTTAAAAAAAAAACGGATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATc                                              <  1:534174/88‑1 (MQ=255)
                            ccTTCTTCTTCATTTTTGGGGGTGATGGCATGAAAAACTTAAAAAAAAAACGGATCCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATc                                              <  1:785540/90‑1 (MQ=255)
                            ccTTCTTCTTCAGTTTGGGGGGTGAGGGAATGAAAAACTTAAAAAAAAAACGGATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATc                                              <  2:468033/90‑1 (MQ=255)
                                 ttcttcAGTTTTGGGGGTGATGGCATGCAAACCTTAAGAAAAAAACGGATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCAt                                         <  1:771244/90‑1 (MQ=255)
                                     tcAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGAt                                                                                 >  2:877540/1‑46 (MQ=255)
                                     tcAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGAt                                                                                 <  1:877540/46‑1 (MQ=255)
                                              ggggTGAGGGCATAAAAAACTTAAAAAAAAAACGGATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCaaaa                            <  1:1163776/88‑1 (MQ=255)
                                                ggTGATGGCATGAAAAACTTCAGAAAAAAACGGATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAAcc                          <  1:576759/90‑1 (MQ=255)
                                                           gCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCa               >  2:502652/1‑90 (MQ=255)
                                                                 ctacAGAAAAAACCGAATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAgcgcgc         <  2:333954/87‑1 (MQ=255)
                                                                   tCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCgg       >  1:333954/1‑90 (MQ=255)
                                                                        gaAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAg  >  1:872647/1‑90 (MQ=255)
                                                                        |                                                                                         
GTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAG  >  NZ_CP009273/1377907‑1378068

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