Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,557,6450TG57.9% 17.1 / 21.7 19L782L (CTA→CTCdosPoxygen‑sensing cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base T (1/7);  new base G (11/0);  total (12/7)
Fisher's exact test for biased strand distribution p-value = 1.59e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-04
Rejected as consensus: Frequency below/above cutoff threshold.

CCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCG  >  NZ_CP009273/1557564‑1557723
                                                                                 |                                                                              
cccGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCg                                                                        <  2:1166777/90‑1 (MQ=255)
          cGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGGAGGGGGGGGGGAaaaaaa                                                              >  1:296648/1‑90 (MQ=255)
           gggCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAat                                                             <  1:1003413/90‑1 (MQ=255)
           gggCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAat                                                             <  2:1355994/90‑1 (MQ=255)
           gggCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAat                                                             <  1:1002303/90‑1 (MQ=255)
                agagGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATccc                                                        <  2:48384/90‑1 (MQ=255)
                           tCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAAcgcg                                             <  2:1072376/90‑1 (MQ=255)
                                   tCAGCGGTAACACGCTGCTCATCCAGCCTTGAATTTCTTCGGCGGGGAGGGGGGGGGGAAAAAAAAATCCCCGAATAAAGCGACAggggg                                     >  1:180259/1‑89 (MQ=255)
                                      gCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCt                                  <  1:153734/90‑1 (MQ=255)
                                                  tgctCATCCAGCCTGGAATTTCTTCGGCGGGGAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCAtct                      >  2:472877/1‑90 (MQ=255)
                                                   gctCATCCAGCCTGGAATTTCTTCGGCGGGGAGGGGGGGGGGAAAGAAAAATCCCCCAATAAACCGACAGTGGATCTTGCGTAGCAtctc                     >  1:1240927/1‑90 (MQ=255)
                                                   gctCATCCAGCCTGGAATTTCTTCGGCGGGGAGGGGGGGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGGTCTTGCGTAGCAtctc                     >  1:283279/1‑90 (MQ=255)
                                                   gctCATCCAGCCTGGAATTTCTTCGGCGGGGAGGGGGCGGGGAAAGAAATATCCCTGAATAACGCGACAGGGGGTCTTGCGTAGCAtctc                     >  1:67971/1‑90 (MQ=255)
                                                           aGCCTGGAATTTCTTCGGCGGGTAGGGGGGGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGct             >  2:1418090/1‑90 (MQ=255)
                                                           aGCCTGGAATTTCTTCGGCGGGGAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGct             >  2:267779/1‑90 (MQ=255)
                                                           aGCCTGGAATTTCTTCGGCGGGGAGGGGGCCGGAAAAAAAAAATCCCCCAATAACCCGGCAATGGATCTTGCCTCGCATCCCCAATTGct             >  2:1336154/1‑90 (MQ=255)
                                                             ccTGGAATTTCTTCGGCGGGGAGGGGGCGGGGAAAAAAATATCCCCGAAAAACACGACAGGGGGTCTTGCGGAGCATCTCAAATTTctct           >  1:1190536/1‑90 (MQ=255)
                                                             ccTGGAATTTCTTCGGCGGGGAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGGGGATCTTGCGTAGCATCTCAAATTGctct           >  1:247644/1‑90 (MQ=255)
                                                                      tcttcGGCGGGGAGGGGGCGGGAAAAAAAAAATCCCTGAATAACGCGACAGTGGGTCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCg  >  1:133379/1‑90 (MQ=255)
                                                                                 |                                                                              
CCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCG  >  NZ_CP009273/1557564‑1557723

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