Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,245,9730TG52.4% 16.1 / 25.0 21N141T (AAC→ACC) plsBglycerol‑3‑phosphate 1‑O‑acyltransferase PlsB
Reads supporting (aligned to +/- strand):  ref base T (2/8);  new base G (11/0);  total (13/8)
Fisher's exact test for biased strand distribution p-value = 2.21e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.23e-03
Rejected as consensus: Frequency below/above cutoff threshold.

TGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGC  >  NZ_CP009273/4245885‑4246062
                                                                                        |                                                                                         
tGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGtt                                                                                          <  1:549911/90‑1 (MQ=255)
tGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGtt                                                                                          <  1:546030/90‑1 (MQ=255)
        agaAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACCTctc                                                                                  >  2:458107/1‑88 (MQ=255)
                             gACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTTACTTTGCCTTTTTTTCGCCCCCgcgcgc                                                             >  2:140672/1‑90 (MQ=255)
                             gACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTTGCCCTTTTTTCTCCCCCgcgcgc                                                             >  2:401845/1‑90 (MQ=255)
                                cGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGAc                                                          <  2:182197/90‑1 (MQ=255)
                                  aGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTCCCCCTTTTTACGCCCCCGCGCGCGACCa                                                        >  2:442514/1‑90 (MQ=255)
                                            cAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCg                                              <  1:293519/90‑1 (MQ=255)
                                            cAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCg                                              <  1:187134/90‑1 (MQ=255)
                                            cAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCg                                              <  1:674111/90‑1 (MQ=255)
                                            cAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTGACTTCTCCTTTTTTCCCCCCCCGCGCGCGGCCAAACAACAccc                                              >  2:593810/1‑89 (MQ=255)
                                            cAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCCCTTCGCCCTTTTTTCGCCCCCGCGCGCGCCCCAACAACAccc                                              >  2:442944/1‑89 (MQ=255)
                                            cAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTCGCCCTTTTTTCCCCCCCGCGCGCGGCCAAACAACAccc                                              >  2:588779/1‑89 (MQ=255)
                                                gAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGacac                                          <  1:401845/90‑1 (MQ=255)
                                                gAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTTGCCTTGCCCTTTTCCCCCCCCCGCCCGCGACCAAACATCACCCacac                                          >  1:377194/1‑90 (MQ=255)
                                                gAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTCGCCTTTTTTACCCCCCCGCGCGCCGCCAAACAACACCgcccc                                          >  2:433956/1‑86 (MQ=255)
                                                   aaaTTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTCGCCTTTTTTACGCCCCCGCGCGCGGCCAAACATCACCGACACTgg                                       >  1:277090/1‑90 (MQ=255)
                                                     aTTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCa                                     <  1:458107/90‑1 (MQ=255)
                                                             aCGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGc                             >  1:39595/1‑90 (MQ=255)
                                                                   ttAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATcc                       <  1:269772/90‑1 (MQ=255)
                                                                     aaGCATACGCAGCGGCGGGGTCACTTCGCCTTTTTTACCCCCCCGCGCGCGGCCCAACAACACCCACACCGGCACCCCCTGCACAccccg                     >  2:600352/1‑88 (MQ=255)
                                                                          tACGCAGCGGCGGGTTCACTTCGCCTTTTTTACGCCCCGGCGCGCGACCAAACATCACCg                                              >  1:406075/1‑60 (MQ=255)
                                                                          tACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCg                                              <  2:406075/60‑1 (MQ=255)
                                                                                        ttCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGc  <  2:77538/90‑1 (MQ=255)
                                                                                        |                                                                                         
TGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGC  >  NZ_CP009273/4245885‑4246062

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