Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,319,9290TG52.4% 17.0 / 31.7 21G114G (GGT→GGGatoBacetyl‑CoA acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/10);  new base G (11/0);  total (11/10)
Fisher's exact test for biased strand distribution p-value = 2.84e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.82e-05
Rejected as consensus: Frequency below/above cutoff threshold.

TAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGT  >  NZ_CP009273/2319842‑2320015
                                                                                       |                                                                                      
tAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTAt                                                                                      <  2:149935/90‑1 (MQ=255)
     tATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGgaa                                                                                 >  2:279884/1‑90 (MQ=255)
     tATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGgaa                                                                                 >  2:226284/1‑90 (MQ=255)
       tgtgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGgaaaa                                                                               >  2:45174/1‑89 (MQ=255)
       tgtgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAat                                                                               >  2:71073/1‑90 (MQ=255)
          ggTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATg                                                                            <  1:572768/90‑1 (MQ=255)
           gTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGa                                                                           >  2:245121/1‑90 (MQ=255)
               gggTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGttt                                                                       <  1:584783/90‑1 (MQ=255)
                           tgtgGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGTTTAGCCCCCtactt                                                           >  2:11195/1‑90 (MQ=255)
                            gtgGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATgg                                                                                    <  1:429610/64‑1 (MQ=255)
                            gtgGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATgg                                                                                    >  2:429610/1‑64 (MQ=255)
                                    tGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGc                                                  >  2:295835/1‑90 (MQ=255)
                                         cccAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAg                                             <  1:679046/90‑1 (MQ=255)
                                            aGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCAc                                          <  1:434581/90‑1 (MQ=255)
                                                cATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGctc                                      <  1:45174/90‑1 (MQ=255)
                                                       gcaggTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTAt                               <  1:295835/90‑1 (MQ=255)
                                                                       gCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGAc               <  1:71073/90‑1 (MQ=255)
                                                                         aTTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCCACTTAATCGATGCAAAAACACGCCCTGGTTATTGTCTTGGGGACGGACAg             >  2:436788/1‑90 (MQ=255)
                                                                                 ggggggGATGGGAAATGTGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGGTATCGTCTTGGAGACGGGCAGGGTTATGa     >  2:155924/1‑90 (MQ=255)
                                                                                 ggggggGATGGGAAAAATGAGTTTAGCCCCCCACTTACTCGATGCAAAAACACGCCCTGGTTATCGTCTTGGGGACGGACAGGTTTATGa     >  1:283908/1‑90 (MQ=255)
                                                                                    gggTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGt  <  2:402727/90‑1 (MQ=255)
                                                                                       |                                                                                      
TAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGT  >  NZ_CP009273/2319842‑2320015

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