Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,531,8010CG57.1% 15.4 / 27.6 21A353G (GCG→GGG) yhgFTex family protein
Reads supporting (aligned to +/- strand):  ref base C (3/6);  new base G (12/0);  total (15/6)
Fisher's exact test for biased strand distribution p-value = 1.55e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.41e-05
Rejected as consensus: Frequency below/above cutoff threshold.

CCGGACTGCGTGCAACGATGGGCCTCGATCCGGGTCTGCGTACTGGGGTAAAAGTGGCGGTGGTCGATGCCACTGGCAAACTGGTGGCGACCGACACCATTTACCCGCACACCGGACAGGCCGCAAAAGCAGCGATGACCGTTGCTGCGCTGTGTGAAAAACATAACGTTGA  >  NZ_CP009273/3531714‑3531885
                                                                                       |                                                                                    
ccGGACTGCGTGCAACGATGGGCCTCGATCCGGGTCTGCGTACTGGGGTAAAAGTGGCGGTGGTCGATGCCACTGGCAAACTGGTGGcga                                                                                    <  2:62424/90‑1 (MQ=255)
  ggACTGCGTGCAACGATGGGCCTCGATCCGGGTCTGCGTACTGGGGTAAAAGTGGCGGTGGTCGATGCCACTGGCAAACTGGTGGcgacc                                                                                  <  2:516631/90‑1 (MQ=255)
       gCGTGCAACGATGGGCCTCGATCCGGGTCTGCGTACTGGGGTAAAAGTGGCGGTGGGCGATGCCCCCGGCAAACTGGTGGGGGCCGacac                                                                             >  2:585404/1‑90 (MQ=255)
       gCGTGCAACGATGGGCCTCGATCCGGGTCTGCGTACTGGGGTAAAAGTGGCGGGGGGGGGTGCCCCCGGCAAACAGGGGGGGGCCCacac                                                                             >  1:245146/1‑90 (MQ=255)
         gTGCAACGATGGGCCTCGATCCGGGTCTGCGTACTGGGGTAAAAGTGGCGGTGGTCGATGCCCCTGGCAAACTGGTGGGGACCGACACCa                                                                           >  2:520857/1‑90 (MQ=255)
         gTGCAACGATGGGCCTCGATCCGGGTCTGCGTACTGGGGTAAAAGTGGCGGTGGTCGATGCCCCCGGCAAACTGGGGGGGACCGACACCa                                                                           >  2:172707/1‑90 (MQ=255)
          tGCAACGATGGGCCTCGATCCGGGTCTGCGTACTGGGGTAAAAGTGGCGGTGGTCGATGCCACTGGCAAACTGGTGGCGACCGACACCAt                                                                          <  1:28296/90‑1 (MQ=255)
                                     gCGTACTGGGGTAAAAGTGGCGGTGGGGGAAGCCACTGGCAAAATGGTGGGGGCCGACACCCTTTACCCCCCCACCGGGCAGGGCCCaaa                                               >  1:531423/1‑90 (MQ=255)
                                     gCGTACTGGGGTAAAAGTGGCGGTGGGCGATGCCCCTGGCAAAATGGTGGGGACCCACACCCTTTACCCCCACACCCGGCAGGCCGCcaa                                               >  2:419102/1‑90 (MQ=255)
                                     gCGTACTGGGGTAAAAGTGGCGGTGGGCGATGCCACTGGCAAACTGGTGGGGGCCGACACCATTTACCCGCACACCGGGCAGGCCGCcaa                                               >  2:558403/1‑90 (MQ=255)
                                         aCTGGGGTAAAAGTGGCGGTGGGCGATGCCCCTGGCAAAATGGGGGGGGCCCACACCATTTTCCCGCACACCGGACAAGCCGCAAAaaca                                           >  1:3745/1‑90 (MQ=255)
                                           tGGGGTAAAAGTGGCGGTGGTTGATGGCCCTGGCAAACTGGGGGGGGCCGACACCCTTTTCCCCCCCCCCGGGCAGGGCGGAAAagcggg                                         >  2:345045/1‑87 (MQ=255)
                                               gTAAAAGTGGCGGGGGGCGGTGCCACTGGCAAACTGGTGGGGGCCCACAACATTTTCCCCCCCACCCGGCAGGCCGCAAAAACAGCGATg                                     >  1:446741/1‑90 (MQ=255)
                                                  aaaGTGGCGGTGGTCGATGCCACTGGCAAACTGGTGGCGACCGACACCATTTACCCGCACACCGGACAGGCCGCAAAAGCAGCGATGAcc                                  <  2:77907/90‑1 (MQ=255)
                                                  aaaGTGGCGGTGGTCGATGCCACTGGCAAACTGGTGGCGACCGACACCATTTACCCGCACACCGGACAGGCCGCAAAAGCAGCGATGAcc                                  <  1:160443/90‑1 (MQ=255)
                                                  aaaGTGGCGGTGGTCGATGCCACTGGCAAACTGGTGGCGACCGACACCATTTACCCGCACACCGGACAGGCCGCAAAAGCAGCGATGAcc                                  <  1:585404/90‑1 (MQ=255)
                                                    aGTGGCGGTGGTCGATGCCACTGGCAAACTGGTGGCGACCGACACCATTTACCCGCACACCGGACAGGCCGCAAAAGCAGCGATGACCGt                                <  1:172707/90‑1 (MQ=255)
                                                             ggtCGGTGCCACTGGCAAACTGGTGGGGGCCGACACCATTTTCCCCCCCACCCGGCAGGCCGCAAAAACAGCGGTGACCGGTGGTGCGCt                       >  1:241868/1‑90 (MQ=255)
                                                             gggcGATGCCACTGGCAAAATGGGGGGGACCGAAACCAATTTTCCGCACACCGGACCGGCCGCAAAAAGAGGGATGAACGCTGCTGCGCt                       >  1:274810/4‑90 (MQ=255)
                                                                cGATGCCACTGGCAAACTGGTGGCGACCGACACCATTTACCCGCACACCGGACAGGCCGCAAAAGCAGCGATGACCGTTGCTGCGCtgtg                    <  2:241868/90‑1 (MQ=255)
                                                                          ggCAAACTGGTGGCGACCGACACCATTTACCCGCACACCGGACAGGCCGCAAAAGCAGCGATGACCGTTGCTGCGCTGTGTGAAAAACAt          >  1:313861/1‑90 (MQ=255)
                                                                            cAAACTGGTGGCGACCGACACCATTTACCCGCACACCGGACAGGCCGCAAAAGCAGCGATGACCGTTGCTGCGCTGTGTGAAAAACATaa        >  1:234450/1‑90 (MQ=255)
                                                                                  ggtggCGACCGACACCATTTACCCGCACACCGGACAGGCCGCAAAAGCAGCGATGACCGTTGCTGCGCTGTGTGAAAAACATAACGTTGa  >  1:282321/1‑90 (MQ=255)
                                                                                       |                                                                                    
CCGGACTGCGTGCAACGATGGGCCTCGATCCGGGTCTGCGTACTGGGGTAAAAGTGGCGGTGGTCGATGCCACTGGCAAACTGGTGGCGACCGACACCATTTACCCGCACACCGGACAGGCCGCAAAAGCAGCGATGACCGTTGCTGCGCTGTGTGAAAAACATAACGTTGA  >  NZ_CP009273/3531714‑3531885

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: