Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273546,6890TG56.2% 12.2 / 18.9 16S265R (AGT→AGGybcFcarbamate kinase
Reads supporting (aligned to +/- strand):  ref base T (0/7);  new base G (9/0);  total (9/7)
Fisher's exact test for biased strand distribution p-value = 8.74e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.80e-04
Rejected as consensus: Frequency below/above cutoff threshold.

ACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACC  >  NZ_CP009273/546602‑546773
                                                                                       |                                                                                    
aCGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGGGAgggg                                                                                    >  1:158033/1‑90 (MQ=255)
  gTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGGGGtt                                                                                  >  1:753622/1‑90 (MQ=255)
  gTGCCATTCGCCACGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGtt                                                                                  <  1:707841/90‑1 (MQ=255)
            ccATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGGGGGGGGGTTTGTTAGAAg                                                                        >  1:464872/1‑90 (MQ=255)
              aTGCCACACCGGATGAGTTAGCGCCGTTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGcc                                                                      <  2:65222/90‑1 (MQ=255)
                        ggATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGGGGGGGGGGGTGTTGGAAGCCGCGGGAAAcc                                                            >  1:781980/1‑90 (MQ=255)
                             aGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGt                                                       <  1:543588/90‑1 (MQ=255)
                                  gcgcCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGAtt                                                  <  1:232272/90‑1 (MQ=255)
                                  gcgcCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGGGGGGGGTTTTTTGAGAAGCCCCGGTAAAACCCCGCGGGtt                                                  >  2:717971/1‑90 (MQ=255)
                                       aTTTGCCAAAGCCGATGGTTCGATGGGGCTGAATGTAACGGCGGTGAGGGGGTATGTCAGAAGCCGTGGTAAAACCGCGTGGGTTggggg                                             >  2:251761/1‑89 (MQ=255)
                                       aTTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGGGGGGTTATGTGAGAAGCCGTGGTGAAACCGCGGGGGGTGGGGc                                             >  2:232272/1‑90 (MQ=255)
                                        tttGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGGGGGGGTTTGTCAGAAACCGGGGTAAAACCCCGGGGGTTtggggg                                            >  2:157178/1‑88 (MQ=255)
                                                    gATGGTTCGATGGGGCCGAATGTAACGGCGGTGGGGGGGTATGTCAGAAGCCGTGGGAAACCCCCGGGGGTTGGGGGGGTTTTGCGaaat                                >  2:352888/1‑88 (MQ=255)
                                                    gATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGGGGGTATGTCAGAAGCCGTGGGAAACCCCCGTGGGTTGGGGGGGTATCGCGAAtt                                >  2:66606/1‑90 (MQ=255)
                                                                               gcggTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCgggg     <  2:753622/90‑1 (MQ=255)
                                                                                  gTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGAcc  <  1:157178/90‑1 (MQ=255)
                                                                                  gTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGAcc  <  1:165927/90‑1 (MQ=255)
                                                                                       |                                                                                    
ACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACC  >  NZ_CP009273/546602‑546773

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