Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,090,8260TG55.6% 15.6 / 26.7 27G67G (GGT→GGG) hemWradical SAM family heme chaperone HemW
Reads supporting (aligned to +/- strand):  ref base T (3/9);  new base G (15/0);  total (18/9)
Fisher's exact test for biased strand distribution p-value = 4.69e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.28e-03
Rejected as consensus: Frequency below/above cutoff threshold.

TGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATT  >  NZ_CP009273/3090742‑3090913
                                                                                    |                                                                                       
tGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGcc                                                                                    <  2:498276/90‑1 (MQ=255)
     aGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGGCGCCGAggc                                                                               >  2:65678/1‑88 (MQ=255)
      gCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGACGCCGCGccc                                                                              >  2:26454/1‑89 (MQ=255)
      gCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGACGCCGAGccc                                                                              >  2:493287/1‑89 (MQ=255)
         tCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGGGGGGCGCCGGGCCCGCt                                                                           >  1:219458/1‑90 (MQ=255)
                aaCGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTccggc                                                                    <  1:760662/90‑1 (MQ=255)
                aaCGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGACGCCGAGCCTGCTTTcccgc                                                                    >  2:416747/1‑90 (MQ=255)
                      cTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCg                                                              <  1:179195/90‑1 (MQ=255)
                       tGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGACGCCGAGCCTGCTTTTCGGCCCGGCGa                                                             >  2:264393/1‑90 (MQ=255)
                             aCGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCggtggt                                                                                         <  1:62475/56‑1 (MQ=255)
                             aCGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCggtggt                                                                                         >  2:62475/1‑56 (MQ=255)
                             aCGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGACCCCGCGCCCGCTTTCCGGCCCGGCGGAGCaaa                                                       >  2:179195/1‑90 (MQ=255)
                              cGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAAc                                                      <  2:427227/90‑1 (MQ=255)
                                aTGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACgc                                                    <  2:23247/90‑1 (MQ=255)
                                  gtgGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGGACGCCGGGCCTGCTTTCCGGCCCGGCGGTGCAAACgctg                                                  >  2:436002/1‑90 (MQ=255)
                                        tACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGGGCCCGCTTTTCCGCCCCGCGATGGAAACGCTGCTgggg                                            >  1:156596/1‑88 (MQ=255)
                                         aCGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGGGGGGCGCCGCGCCCGCTTTTCCGCCCGGCGCTGCAAACGCTGCTggggg                                           >  1:190300/1‑87 (MQ=255)
                                                   ccGTGAAGTAAAGACAATCTTTATTGGCGGTGGGGCGCCGAGCCCGCTTTTCGGCCCCGCGGAGCAAAAGCCGCTGGGAGGGgcgcgtgc                                 >  1:341406/1‑90 (MQ=255)
                                                     gTGAAGTAAAGACAATCTTTATTGGCGGTGGGACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGGCGGCGGGCGTgcgc                               >  2:128603/1‑90 (MQ=255)
                                                     gTGAAGTAAAGACAATCTTTATTGGCGGGGGGGAACCCAGCCCGCTTTCCGGCCCCGCGGTGCAAAAGCTCCTGGGCGGCGGGGGTgggc                               >  1:414027/1‑90 (MQ=255)
                                                      tGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCCGCCCGGCGATGCAAACGCTGCTGGGCGGCGTGCGTGggcg                              >  2:641403/1‑90 (MQ=255)
                                                          gTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTg                          >  2:318245/1‑90 (MQ=255)
                                                            aaaGACAATCTTTATTGGCGGTGGGGCCCCGCGGCTGCTTTCCGGCCCGGCGACGCAAACCCCTCGGGGCGGCGTGCGCGCGCGTTTGcc                        >  2:583298/1‑90 (MQ=255)
                                                            aaaGACAATCTTTATTGGCGGGGGGGAGCCGAGGCCGCTTTTCGGCCCCGCGGTGCAAAACCCGCTTGGCGGGGGGCGGGGGCGGTTGcc                        >  1:733571/1‑90 (MQ=255)
                                                                    tCTTTATTGGCGGTGGGACGCCGAGCCCGCTTTTCCGCCCCGCGATGCAAACGCTGCTGGGCGGCGGGGGGGCGCGGTTTCCGCCggggg                >  1:448630/1‑87 (MQ=255)
                                                                     cTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGc               <  1:266140/90‑1 (MQ=255)
                                                                      tttATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCg              <  1:318245/90‑1 (MQ=255)
                                                                          ttGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATg          <  2:219458/90‑1 (MQ=255)
                                                                             gCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAg       <  1:660391/90‑1 (MQ=255)
                                                                                  ggtACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAAtt  <  1:648548/90‑1 (MQ=255)
                                                                                    |                                                                                       
TGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGTGGTACGCCGAGCCTGCTTTCCGGCCCGGCGATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATT  >  NZ_CP009273/3090742‑3090913

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