Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,594,1610CT58.2% 0.9 / 75.1 67S62F (TCT→TTT) bglJDNA‑binding transcriptional activator BglJ
Reads supporting (aligned to +/- strand):  ref base C (7/21);  new base T (39/0);  total (46/21)
Fisher's exact test for biased strand distribution p-value = 9.13e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.31e-08
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAG  >  NZ_CP009273/4594076‑4594247
                                                                                     |                                                                                      
aaCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttttttt                                                                                    >  2:15937/1‑90 (MQ=255)
aaCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttttttt                                                                                    >  2:1456022/1‑90 (MQ=255)
  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)
                aTACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAttttttttttttt                                                                                  >  2:1688494/1‑75 (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:962996/90‑1 (MQ=255)
                          cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagaga                                                          <  2:1344912/90‑1 (MQ=255)
                          cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagaga                                                          <  1:98723/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)
                           aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAttttttttttttt                                                                                  >  1:438587/1‑64 (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                                                                              <  2:324138/59‑1 (MQ=255)
                                     cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGcc                                                                              <  2:901387/59‑1 (MQ=255)
                                     cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGcc                                                                              >  1:901387/1‑59 (MQ=255)
                                       tttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGc                                             <  1:290004/90‑1 (MQ=255)
                                           aGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCCTGAGAAGGGAGGGCAGAGGGGGGGTATCTTGCctga                                         >  2:1665240/1‑90 (MQ=255)
                                               tGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttttc                                                                                  <  2:1680726/45‑1 (MQ=255)
                                               tGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttttc                                                                                  >  1:1680726/1‑45 (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)
                                                                                     |                                                                                      
AACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAG  >  NZ_CP009273/4594076‑4594247

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: