Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,594,1670CT54.5% 15.4 / 79.5 67S64F (TCT→TTT) bglJDNA‑binding transcriptional activator BglJ
Reads supporting (aligned to +/- strand):  ref base C (9/21);  new base T (36/0);  total (46/21)
Fisher's exact test for biased strand distribution p-value = 1.61e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.27e-08
Rejected as consensus: Frequency below/above cutoff threshold.

CCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGG  >  NZ_CP009273/4594078‑4594256
                                                                                         |                                                                                         
cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCTCGGATCTCGTTTGCGGCGGTCATtttttctttttc                                                                                           <  1:650006/90‑1 (MQ=255)
  cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTg                                                                                         <  1:538107/90‑1 (MQ=255)
   ttATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttttttttc                                                                                           >  1:1708334/1‑87 (MQ=255)
   ttATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttttc                                                                                           <  2:1708334/87‑1 (MQ=255)
         gCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATgag                                                                                  <  2:984853/90‑1 (MQ=255)
         gCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATgag                                                                                  <  1:187932/90‑1 (MQ=255)
          cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCCTTTTTTCTTTTTTTTCCCTgaga                                                                                 >  2:98723/1‑90 (MQ=255)
           tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGaa                                                                                >  2:1199864/1‑90 (MQ=255)
           tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGaa                                                                                >  1:469095/1‑90 (MQ=255)
            ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCCTGAGAAg                                                                               >  1:241832/1‑90 (MQ=255)
            ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTgaggag                                                                               >  2:180875/1‑90 (MQ=255)
            ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCAt                                                                                     >  1:1512973/1‑84 (MQ=255)
            ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCAt                                                                                     <  2:1512973/84‑1 (MQ=255)
             cATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAgg                                                                              >  2:599065/1‑89 (MQ=255)
              aTACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttttc                                                                                           <  1:1688494/76‑1 (MQ=255)
                 ccTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgc                                                                          <  1:946676/90‑1 (MQ=255)
                   tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGGgcgc                                                                        >  1:1613503/1‑90 (MQ=255)
                   tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCTTGAGGAGGGGgcga                                                                        >  2:275944/1‑89 (MQ=255)
                    ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTg                                                                                         >  1:967917/1‑72 (MQ=255)
                    ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGGGCa                                                                       >  2:1580981/1‑90 (MQ=255)
                    ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGcc                                                                                       >  1:708793/1‑74 (MQ=255)
                    ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGCGCa                                                                       >  2:719700/1‑90 (MQ=255)
                    ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTg                                                                                         <  2:967917/72‑1 (MQ=255)
                    ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGcc                                                                                       <  2:708793/74‑1 (MQ=255)
                    ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCATTGGAAGGGAGGGCa                                                                       >  2:598882/1‑90 (MQ=255)
                      aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGGGCGCaga                                                                     >  1:360955/1‑90 (MQ=255)
                      aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGAGAGCGCaga                                                                     <  1:1570980/90‑1 (MQ=255)
                        cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagaga                                                                   <  2:1344912/90‑1 (MQ=255)
                        cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagaga                                                                   <  1:98723/90‑1 (MQ=255)
                        cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagaga                                                                   <  2:962996/90‑1 (MQ=255)
                         aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTTTTTTTTTTTTTTTTCCATTAGAAGGGGGGGCagaggg                                                                  >  1:1096674/1‑88 (MQ=255)
                         aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCCTGGGAAGTGAGCGCagaggg                                                                  >  2:1716399/1‑88 (MQ=255)
                         aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCCTGAGAAAGGAGCGGagagag                                                                  >  2:253676/1‑90 (MQ=255)
                         aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttttc                                                                                           <  2:438587/65‑1 (MQ=255)
                         aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGTATTTTTTTTTTTTTTGCCCTGAGGAGGGAGGggagggag                                                                  >  2:1230383/1‑90 (MQ=255)
                          gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCCTGAGAAGGGGGCGCCGAGAgg                                                                 >  1:930813/1‑90 (MQ=255)
                          gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGAgggg                                                                 >  2:624214/1‑90 (MQ=255)
                          gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAgg                                                                 >  1:513669/1‑90 (MQ=255)
                           tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGGAGGGGGGGCAGAGgggg                                                                >  1:461942/1‑90 (MQ=255)
                           tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGGGAAGTGAGCGGAGAGgggg                                                                >  1:973348/1‑90 (MQ=255)
                           tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAggg                                                                >  1:129396/1‑90 (MQ=255)
                              ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGTGGTCATTTTTTCTTTTTCTGCCATg                                                                                    <  2:1388584/67‑1 (MQ=255)
                              ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATg                                                                                    >  1:1388584/1‑67 (MQ=255)
                              ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTGTGCCATGAGGAGTGAGCGCAGAGAggggta                                                             >  1:644273/1‑87 (MQ=255)
                                aGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTAt                                                           <  1:719700/90‑1 (MQ=255)
                                 gTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCGGAGAGGGGGTATc                                                          >  1:43126/1‑90 (MQ=255)
                                 gTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCCTGAGAAGTGAGCGCAGAGAGGGGTTATc                                                          >  1:1175561/1‑90 (MQ=255)
                                 gTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTCTGCCCTGAGAAAGGGGGGCAGGGAGGGGGTATc                                                          >  1:687610/1‑90 (MQ=255)
                                   cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGAGAGGGGTTATCtt                                                        >  2:1626683/1‑90 (MQ=255)
                                   cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGGGCAGAGGGGGGTTATCtt                                                        >  2:290004/1‑90 (MQ=255)
                                   cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGcc                                                                                       >  1:324138/1‑59 (MQ=255)
                                   cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGcc                                                                                       >  1:901387/1‑59 (MQ=255)
                                   cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGcc                                                                                       <  2:324138/59‑1 (MQ=255)
                                   cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGcc                                                                                       <  2:901387/59‑1 (MQ=255)
                                     tttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGc                                                      <  1:290004/90‑1 (MQ=255)
                                         aGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCCTGAGAAGGGAGGGCAGAGGGGGGGTATCTTGCctga                                                  >  2:1665240/1‑90 (MQ=255)
                                             tGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttttc                                                                                           >  1:1680726/1‑45 (MQ=255)
                                             tGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttttc                                                                                           <  2:1680726/45‑1 (MQ=255)
                                             tGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGGGAAGGGAGGGCAGGGGGGGGTTATCTTTCcttactga                                              >  2:1440592/1‑90 (MQ=255)
                                              gCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTCTGCCATGAGGAGAGAGCGCAGAGAGGGGTTATCTTTCCTGACTGaa                                             >  1:1720169/1‑90 (MQ=255)
                                                      gcgGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGGGCGCAGAGAGGGGTTATTTTGCCTGACTGAACTGGCggt                                     >  1:97118/1‑88 (MQ=255)
                                                        ggATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGGAGGGGGGGGGGGGGGGGGTTTTTTTGCCTTACTTGACTTGCGAttt                                   >  2:598268/1‑89 (MQ=255)
                                                        ggATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTTCCTGACTGAACTGGCGATTa                                   >  1:1416466/1‑90 (MQ=255)
                                                           tctcGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCCTGAGGAATGAGCTCTGAGGGGGGTTATTTTGCCTGACTTGACTGGCGTTTaaat                                >  1:531558/1‑88 (MQ=255)
                                                             tcGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGGGCGCAGAGGGGGGTTATCTTGGCTGACTTAACTGGCGGTTTAGttt                              >  1:882298/1‑90 (MQ=255)
                                                             tcGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGttt                              <  2:469095/90‑1 (MQ=255)
                                                                 ttGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCcgc                          <  1:459370/90‑1 (MQ=255)
                                                                   gcggcgGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGCGCCGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAAGGTTTcgcgg                        >  1:1595185/1‑89 (MQ=255)
                                                                   gcggcgGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGTTTTTTTGCCTGACTGAACTGGGGATTAAGTTTTcgcgg                        >  1:1645105/1‑89 (MQ=255)
                                                                   gcggcgGTCATTTTTTCTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGt                        >  1:32488/1‑90 (MQ=255)
                                                                      gcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAcc                     <  1:1559829/90‑1 (MQ=255)
                                                                      gcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAcc                     <  1:1230383/90‑1 (MQ=255)
                                                                       cggTCATTTTTTCTTTTTCTGCCATGAAAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAccc                    <  2:937318/90‑1 (MQ=255)
                                                                          tCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGc                 >  1:978515/1‑90 (MQ=255)
                                                                               ttttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTa            <  1:447847/90‑1 (MQ=255)
                                                                                tttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAg           >  1:655414/1‑90 (MQ=255)
                                                                                       ttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGc    <  1:1626683/90‑1 (MQ=255)
                                                                                       ttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGc    <  2:1175561/90‑1 (MQ=255)
                                                                                         ctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCgg  <  1:253676/90‑1 (MQ=255)
                                                                                         ctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCgg  <  2:443529/90‑1 (MQ=255)
                                                                                         ctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCgg  <  1:598882/90‑1 (MQ=255)
                                                                                         |                                                                                         
CCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGG  >  NZ_CP009273/4594078‑4594256

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