Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,213,5210TG55.6% 30.3 / 37.1 27K12N (AAA→AACiclRglyoxylate bypass operon transcriptional repressor IclR
Reads supporting (aligned to +/- strand):  ref base T (0/12);  new base G (15/0);  total (15/12)
Fisher's exact test for biased strand distribution p-value = 5.75e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.62e-06
Rejected as consensus: Frequency below/above cutoff threshold.

TCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGT  >  NZ_CP009273/4213448‑4213604
                                                                         |                                                                                   
tCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTcg                                                                     <  2:246078/90‑1 (MQ=255)
            aTTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTg                                                         <  1:78778/90‑1 (MQ=255)
              ttCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCGGTTTGCGGGAATggggggg                                                       >  1:87120/1‑86 (MQ=255)
                   gCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTCGCGGGGATGGGGGCGGCCAt                                                  >  1:852998/1‑90 (MQ=255)
                     cACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCCCGTTTCGCGGGAATGGGGGCGAccccga                                                >  1:253436/1‑86 (MQ=255)
                       cgcgTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCGTTTTGCGGGGATGTGGGCGACCCTGTCa                                              >  1:177953/1‑90 (MQ=255)
                         cgTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGt                                            <  1:317944/90‑1 (MQ=255)
                         cgTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGt                                            <  1:95822/90‑1 (MQ=255)
                           ttAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGtct                                          <  2:542126/90‑1 (MQ=255)
                              aaGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCtt                                       <  1:579643/90‑1 (MQ=255)
                                 aCTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTTGCGGGGATGGGGGGGGCCACGACAGGCTCttttt                                    >  1:538011/1‑90 (MQ=255)
                                   tGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGGTTTGCGGGGATGGGGGGGGCCATGACAGTCTCttttttt                                  >  2:634382/1‑89 (MQ=255)
                                   tGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGGTTTGCGGGGAAGGGGGGGGCCACGACAGGCTCttttttt                                  >  2:91955/1‑89 (MQ=255)
                                   tGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGGTTCGCGGGGAAGGGGGGGGCCCTGACAGTCTCttttttt                                  >  2:833223/1‑89 (MQ=255)
                                    gAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCGTTTTGCGGGGGTGGGGGGGGCCAAGACAGTCTCtttttttt                                 >  2:138852/1‑88 (MQ=255)
                                    gAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCGGTTTGCGGGGAAGGGGGGGGCCCTGAAAGGCTCtttttttt                                 >  1:577197/1‑88 (MQ=255)
                                          gTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTTGCGGGGATGGGGGGGGCCACGACAGACTCTTTTTTTTTTAttg                           >  1:43188/1‑88 (MQ=255)
                                          gTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTCGCGGGAATGGGGGGGGCCATGACAGGCTCTTTTTTTTGTATCg                           >  1:854825/1‑90 (MQ=255)
                                              aGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTTGCGGGGATGGGGGGGACCCTGAGAGGCTCTTTTTTTTTTATCGTggg                       >  1:688341/1‑89 (MQ=255)
                                               gTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGaa                      <  1:91955/90‑1 (MQ=255)
                                                        gCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAAtcattttc             <  2:43188/90‑1 (MQ=255)
                                                        gCGGTGGCAACGGCGGGGTTTTTTCCCCGGTTCCCGGGGAAGGGGGGGGCCATGGCCGGCTCTTTTTTTTGTTTTGTGGGAAtacttttt             >  2:257371/1‑89 (MQ=255)
                                                             ggCAACGGCGGGGTTTTTGCCGCGTTTCGCGGGAATGGGGGCGACCATGACAGTCTCTTTTTTTTGTATCGTGGGAAAcatttttatttt        >  1:749135/1‑90 (MQ=255)
                                                              gCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCAttttt       <  2:177953/90‑1 (MQ=255)
                                                              gCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCAttttt       <  2:570494/90‑1 (MQ=255)
                                                                   ggcggGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGt  <  1:404147/90‑1 (MQ=255)
                                                                   ggcggGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGt  <  2:749135/90‑1 (MQ=255)
                                                                         |                                                                                   
TCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGT  >  NZ_CP009273/4213448‑4213604

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