Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273807,1770TG53.8% 10.3 / 16.4 13V67G (GTG→GGG) bioCmalonyl‑ACP O‑methyltransferase BioC
Reads supporting (aligned to +/- strand):  ref base T (0/6);  new base G (7/0);  total (7/6)
Fisher's exact test for biased strand distribution p-value = 5.83e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-03
Rejected as consensus: Frequency below/above cutoff threshold.

TCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGC  >  NZ_CP009273/807091‑807234
                                                                                      |                                                         
tCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGAc                                                        <  1:446097/90‑1 (MQ=255)
      gCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGACGGCCtt                                                  >  2:61107/1‑90 (MQ=255)
               cacCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGGCACGCGCGGGGTGCGGGCTTTGAtctgtc                                         >  1:288283/1‑90 (MQ=255)
                acCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGAAGTCACGCGCCGGGGACGGCCTTACATCTCTcg                                        >  2:314460/1‑90 (MQ=255)
                           ggccgcgGGTTGTGGACCTGGCTGGATGGGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCt                             <  2:313792/87‑1 (MQ=255)
                           ggACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGCCACGCGCGTGGGGCGGGCCTTGATCTCTCGCCGCCCATgtt                             >  1:313792/1‑88 (MQ=255)
                           ggACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGACGGCCTTAGATCTCTCGCCCCCCATGCt                             >  2:207720/1‑90 (MQ=255)
                                cgGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTc                        <  2:384733/90‑1 (MQ=255)
                                    ttGTGGCCCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGc                    <  1:454045/90‑1 (MQ=255)
                                    ttGTGGCCCGGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGc                    <  1:251524/90‑1 (MQ=255)
                                      gtgGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCAc                  <  2:392331/90‑1 (MQ=255)
                                            ctggctggATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGACGGCCTTAGATCTCTCCCCGCCCATGCTTGTTCAGGCACGCCCGa            >  1:35986/1‑90 (MQ=255)
                                                      gAGCCGCCACTGGCGGGAAAGTCACGCGCAGGGGACCGCCTTAGATCTCTCCCCCCCCATGCTTTTTTAGGCACGCCAGAAGGATGCCGc  >  1:142885/1‑90 (MQ=255)
                                                                                      |                                                         
TCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGC  >  NZ_CP009273/807091‑807234

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