Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,594,1670CT60.5% 4.8 / 48.4 43S64F (TCT→TTT) bglJDNA‑binding transcriptional activator BglJ
Reads supporting (aligned to +/- strand):  ref base C (3/14);  new base T (26/0);  total (29/14)
Fisher's exact test for biased strand distribution p-value = 8.68e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.12e-07
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTT  >  NZ_CP009273/4594088‑4594244
                                                                               |                                                                             
cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttttc                                                                               <  2:604312/80‑1 (MQ=255)
 tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGaa                                                                    >  1:550507/1‑90 (MQ=255)
 tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGaa                                                                    <  1:742298/90‑1 (MQ=255)
  ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCCTTAGAAg                                                                   >  1:473816/1‑90 (MQ=255)
    aTACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTg                                                                 <  2:33134/90‑1 (MQ=255)
     tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGa                                                                <  1:299259/90‑1 (MQ=255)
         tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCCTGAGAAGGGAgggg                                                            >  1:148237/1‑87 (MQ=255)
          ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGGGCa                                                           >  2:56858/1‑90 (MQ=255)
          ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAATGGGCGCa                                                           >  2:732704/1‑90 (MQ=255)
          ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGAGCGCa                                                           >  2:621782/1‑90 (MQ=255)
           taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCag                                                          <  1:357424/90‑1 (MQ=255)
           taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCag                                                          <  2:87910/90‑1 (MQ=255)
             ccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagag                                                        <  2:618764/90‑1 (MQ=255)
              cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGGGGGCagaga                                                       >  2:458718/1‑90 (MQ=255)
               aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTTATTTTTTTTTTTTCTGCCCTGTGAAGTGGGGGCagggag                                                      >  2:430225/1‑90 (MQ=255)
               aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttttc                                                                               <  2:485080/65‑1 (MQ=255)
               aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCagagag                                                      >  2:603485/1‑90 (MQ=255)
               aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCagagag                                                      >  2:299259/1‑90 (MQ=255)
               aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGAGCGCagagag                                                      >  1:114055/1‑90 (MQ=255)
               aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGAGCGCagagag                                                      >  1:75432/1‑90 (MQ=255)
               aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGCGCcgggag                                                      >  1:377/1‑90 (MQ=255)
               aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGGGAAGGGAGCGGagaggg                                                      >  2:15425/1‑88 (MQ=255)
                gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTGCCCTGAGAAGAGGGCGCAGAGAgg                                                     >  1:378093/1‑90 (MQ=255)
                gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTGCCATGAGGAGTGAGGGCAGAGAgg                                                     >  2:472938/1‑90 (MQ=255)
                 tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAggg                                                    >  1:474138/1‑90 (MQ=255)
                    ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGAtt                                                 >  2:99113/1‑90 (MQ=255)
                         cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTg                                                                             >  2:635811/1‑57 (MQ=255)
                         cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCATGAGAAGTGAGCGGAGAGAGGGATTATCtt                                            >  2:714419/1‑90 (MQ=255)
                         cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTg                                                                             <  1:635811/57‑1 (MQ=255)
                         cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCCTGAGAAGGGGGCGCCGAGAGGGGTTATCtt                                            >  1:635536/1‑90 (MQ=255)
                         cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGGGCAGAGAGGGGTTAtttt                                            >  1:527352/1‑90 (MQ=255)
                          aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGAGCGCAGAGAGGGGTTATCTTg                                           >  2:335154/1‑90 (MQ=255)
                                   tGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGAGAGGGATTTTTTTGCctgactga                                  >  1:163956/1‑90 (MQ=255)
                                       aTGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGGAGTGAGGGCGGAGAGGGGGTATTTTTTCTGACTGAACTg                              >  1:659035/1‑90 (MQ=255)
                                          tCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACCGGCg                           <  1:485064/90‑1 (MQ=255)
                                            gcgGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGGGCAGAGGGGGATTATCTTGCCTGACTGAACTGGCGAt                         >  2:611591/1‑90 (MQ=255)
                                            gcgGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGCGCAGAGAGGGGTTATTTTTCCTGACTGAACTGGCGAt                         >  1:60168/1‑90 (MQ=255)
                                              ggATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGGGCAGAGAGGGGTTATCTTTCCTGACTGAAATGGGGATTa                       >  2:783887/1‑90 (MQ=255)
                                              ggATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTa                       <  1:198960/90‑1 (MQ=255)
                                                  ctcGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAAGtt                   >  2:414052/1‑90 (MQ=255)
                                                             cggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAccc        <  1:621782/90‑1 (MQ=255)
                                                                tCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGc     <  1:335154/90‑1 (MQ=255)
                                                                  aTTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGt   <  1:783887/90‑1 (MQ=255)
                                                                   ttttttCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGtt  <  2:635536/90‑1 (MQ=255)
                                                                   ttttttCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGtt  <  1:15425/90‑1 (MQ=255)
                                                                               |                                                                             
CTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTT  >  NZ_CP009273/4594088‑4594244

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: ‑