Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,319,9290TG55.0% 1.5 / 22.9 20G114G (GGT→GGGatoBacetyl‑CoA acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (1/8);  new base G (11/0);  total (12/8)
Fisher's exact test for biased strand distribution p-value = 7.14e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGA  >  NZ_CP009273/2319851‑2320012
                                                                              |                                                                                   
tgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAatat                                                                          >  1:24729/1‑90 (MQ=255)
    tCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATTAgg                                                                      >  2:292056/1‑89 (MQ=255)
     cGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAAAATGAggt                                                                     >  1:310660/1‑88 (MQ=255)
     cGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAggt                                                                     >  1:173691/1‑88 (MQ=255)
      gggTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGttt                                                                    <  2:706686/90‑1 (MQ=255)
           ttAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGccc                                                               <  2:682771/90‑1 (MQ=255)
           ttAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGccc                                                               <  2:24729/90‑1 (MQ=255)
                 gtgtGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGTTTTGCCCCCcact                                                         >  1:525729/1‑90 (MQ=255)
                   gtgGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGTTTTACCCCCtactta                                                       >  1:324375/1‑90 (MQ=255)
                      gcgcTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTACGGAAAATATGAGTTTAGCCCCCTACTTACTc                                                    <  1:532323/90‑1 (MQ=255)
                      gcgcTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGGTTTGCCCCCCACCTACTc                                                    >  2:145877/1‑90 (MQ=255)
                                         ttCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAATATGAGGTTAGCCCCCTCCTTACTCGATGCAAAAGCACGCTCTg                                 >  2:158143/1‑90 (MQ=255)
                                              gcaggTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGTTTAGCCCCCCACTTACTCGATGCAAAAACACGCCCTGGTTAt                            >  2:214108/1‑90 (MQ=255)
                                                  gTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTc                        >  1:179398/1‑90 (MQ=255)
                                                              gCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGAc            <  1:621144/90‑1 (MQ=255)
                                                                aTTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAg          <  1:292056/90‑1 (MQ=255)
                                                                      gCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTAt    <  2:475070/90‑1 (MQ=255)
                                                                        ggggggGATGGGAAATATGAGTTTAGCCCCCCACTTACTCGATGCAAAAGCACGCCCTGGTTATCGTCTTGGAGGCGGGCAGGGTTATGa  >  1:475070/1‑90 (MQ=255)
                                                                             gTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCaa                                             <  2:391417/42‑1 (MQ=255)
                                                                             gTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCaa                                             >  1:391417/1‑42 (MQ=255)
                                                                              |                                                                                   
TGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGA  >  NZ_CP009273/2319851‑2320012

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