Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,594,1670CT59.6% 1.5 / 56.3 52S64F (TCT→TTT) bglJDNA‑binding transcriptional activator BglJ
Reads supporting (aligned to +/- strand):  ref base C (2/19);  new base T (31/0);  total (33/19)
Fisher's exact test for biased strand distribution p-value = 2.75e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.03e-08
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGC  >  NZ_CP009273/4594080‑4594254
                                                                                       |                                                                                       
cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTg                                                                                       <  2:390442/90‑1 (MQ=255)
cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTg                                                                                       <  1:753418/90‑1 (MQ=255)
  tataAGCTCCATTCCTATTCCAGTCCGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGcc                                                                                     <  2:386130/90‑1 (MQ=255)
     aaGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATg                                                                                  <  2:186208/90‑1 (MQ=255)
        cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATggga                                                                               >  2:641623/1‑90 (MQ=255)
        cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATgaga                                                                               >  1:403158/1‑90 (MQ=255)
        cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCCTgaga                                                                               >  1:747432/1‑90 (MQ=255)
         tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGaa                                                                              >  2:692534/1‑90 (MQ=255)
          ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCCTGGGAAg                                                                             >  1:15212/1‑90 (MQ=255)
          ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCATgaggag                                                                             >  2:274863/1‑90 (MQ=255)
                  ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGGgcgcc                                                                     >  2:765175/1‑89 (MQ=255)
                  ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGGGCa                                                                     >  2:278483/1‑90 (MQ=255)
                  ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCa                                                                     >  2:81275/1‑90 (MQ=255)
                  ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCAATTTTTCTTTTTTTGCCCTTAGAAGGGAgcgcc                                                                     >  2:680803/1‑89 (MQ=255)
                   taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCag                                                                    <  2:229710/90‑1 (MQ=255)
                     ccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagag                                                                  <  1:278483/90‑1 (MQ=255)
                      cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagaga                                                                 <  2:213125/90‑1 (MQ=255)
                       aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGGGAAGTGAGGGCagagag                                                                >  1:61791/1‑90 (MQ=255)
                       aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGAGGGCagagag                                                                >  2:613650/1‑90 (MQ=255)
                       aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGGGGGCagagag                                                                >  2:421168/1‑90 (MQ=255)
                       aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGCGCagagag                                                                >  1:206952/1‑90 (MQ=255)
                       aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGCGCagagag                                                                >  1:37911/1‑90 (MQ=255)
                       aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTTCCTGAGAAGTGAGGGCcgagag                                                                >  1:636243/1‑90 (MQ=255)
                       aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCATGAGGAGTGAGCGCggagag                                                                >  1:364039/1‑90 (MQ=255)
                        gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTTAGAAGTGGGGGCAGGGAgg                                                               >  2:324289/1‑90 (MQ=255)
                        gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTCCTTTTTTTGCCATGAGAAGGGAGCGCAGAGAgg                                                               >  2:702144/1‑90 (MQ=255)
                        gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCCTGAGAAAGGAGGGCCGAGAgg                                                               >  1:764662/1‑90 (MQ=255)
                         tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTTTTTTTTTTTTTTTGCCCTGAGAAGTGAGGGCAGAGAggg                                                              >  2:187243/1‑90 (MQ=255)
                          cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCCTGAGAAGGGAGGGCAGAGAgggg                                                             >  2:374845/1‑89 (MQ=255)
                          cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAgggg                                                             >  2:613411/1‑89 (MQ=255)
                          cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGa                                                             <  1:448497/90‑1 (MQ=255)
                             gAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGGGCGCAGGGAGGGGTTa                                                          >  2:341795/1‑90 (MQ=255)
                                   tttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGc                                                    <  1:117202/90‑1 (MQ=255)
                                   tttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGc                                                    <  2:187713/90‑1 (MQ=255)
                                      cAGGATGCCATGTCGCGGTTCTCGTTTGGGGCGGGCGTGTTTTTTTTTTTTGCCGTTGGTAGTGGGGGGGGCGAGGGGGTATATTGCctg                                                 >  1:92504/1‑90 (MQ=255)
                                      cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGGGGAGGGAGCGCCGAGGGGGGTTATCTTGCctg                                                 >  1:282599/1‑90 (MQ=255)
                                        ggATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGTGAGGGCAGAGAGGGATTATCTTGCctgac                                               >  2:694204/1‑90 (MQ=255)
                                             ccATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGGGAAGTGAGGGCAGGGAGGGGTTATCTTTCCTGACTGAAc                                          >  2:391331/1‑90 (MQ=255)
                                             ccATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGGGGAGGGAGGGATTTTCTTGCCTGACTGAAc                                          >  2:730284/1‑90 (MQ=255)
                                              cATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACt                                         <  2:764662/90‑1 (MQ=255)
                                              cATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACt                                         <  1:606410/90‑1 (MQ=255)
                                              cATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACt                                         <  1:765175/90‑1 (MQ=255)
                                                  tCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGGGCGCGGGGGGGGATTATTTTTCCTGACTGAACTGGCg                                     >  1:181702/1‑90 (MQ=255)
                                                      ggATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGGTTTTTTTTCCTGACTGAACTGGCGGTTa                                 >  2:530828/1‑90 (MQ=255)
                                                      ggATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGCGCAGAGAGGGGGTTTCTTTCCTGACTGAACTGGCGATTa                                 >  1:214024/1‑90 (MQ=255)
                                                         tctcGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGt                              <  1:757128/90‑1 (MQ=255)
                                                           tcGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCATGAGAAGGGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAAATGGCGGTTAAGttt                            >  2:481718/1‑90 (MQ=255)
                                                                      ggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCg                 <  2:618603/90‑1 (MQ=255)
                                                                      ggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCg                 <  2:665108/90‑1 (MQ=255)
                                                                      ggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCg                 <  1:341795/90‑1 (MQ=255)
                                                                       gTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCgg                <  1:81275/90‑1 (MQ=255)
                                                                          aTTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGt             <  2:181702/90‑1 (MQ=255)
                                                                                 ctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTa      <  2:269108/90‑1 (MQ=255)
                                                                                   ttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTAtt    >  2:84637/1‑90 (MQ=255)
                                                                                    tttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTg   <  2:649363/90‑1 (MQ=255)
                                                                                     ttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGc  <  2:92504/90‑1 (MQ=255)
                                                                                       |                                                                                       
CTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGC  >  NZ_CP009273/4594080‑4594254

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