Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8800TC40.8% 68.9 / 11.8 49K271E (AAA→GAA) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (22/7);  new base C (0/20);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 3.29e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTC  >  NC_000913/357730‑358028
                                                                                                                                                        |                                                                                                                                                         
aCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt                                                                                                                                                           <  1:172983/151‑1 (MQ=255)
   cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTcc                                                                                                                                                         <  1:1125743/150‑1 (MQ=255)
   cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTc                                                                                                                                                          <  2:1588353/149‑1 (MQ=255)
     tAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                                                                                                      >  2:1482480/1‑151 (MQ=255)
         ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgccg                                                                                                                                                     <  2:1011129/148‑1 (MQ=255)
         ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGca                                                                                                                                                   <  2:1416226/150‑1 (MQ=255)
                aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                                          >  1:788458/1‑151 (MQ=255)
                    tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgcc                                                                                                                                       >  2:516436/1‑150 (MQ=255)
                      tcgcATCAGGCAATTCCAGCCGCAGCCCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAACATGCCTTCGCGGCGGCCGTCGCCCAGCTAATTTTCCGCCGCAAGACAACCGCCGt                                                                                                                                     >  2:408778/2‑150 (MQ=18)
                      gcgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCAGt                                                                                                                                     <  1:1587500/150‑1 (MQ=255)
                           tCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctctc                                                                                                                               >  1:837912/1‑151 (MQ=255)
                             aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                              >  1:1130935/1‑150 (MQ=255)
                                 attcccACCCGCAGCCTTGTGCCGGCGCCTCCCGTGTTCCTTTCCCCCCAACACCCGCGCTTTTACCCAACGCCATGTGAAAAAATGCTTTCGCGGCGTCTGTCTCCCATTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCagtag                                                                                                                         <  1:1407055/147‑1 (MQ=255)
                                   ttCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTg                                                                                                                       >  2:565791/1‑151 (MQ=255)
                                     ccAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                                                     >  2:221257/1‑151 (MQ=255)
                                             agccttgtgTCAGCGCCTCCCGTGATTCTTTTCCCCCACCACCCGGTCTTTACCCCACCGCCATGTGCAAAAATGCTTCCGCGGGGGCTGCCTGCCAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCa                                                                                                              <  2:489743/144‑1 (MQ=255)
                                                  tgtgTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGcga                                                                                                             >  1:497668/1‑143 (MQ=255)
                                                   gtgtCAGCGGTTCCCGTGTTTCTTTCCCGCAAACACCCGCACATTTACCCACCCCCATGGCAAAAATGCCCTGCGGGGGGGCTTTCTGCCGTCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaa                                                                                                       <  2:1151513/149‑1 (MQ=255)
                                                    tgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTCTTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggggagg                                                                                                          >  2:838638/1‑140 (MQ=255)
                                                           ggCTCCGGTGATTCATTCCCGCCAACACCCGCGCATTTACCCACCGCCATGTGAAAAAATGCCTTCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTa                                                                                             <  1:370328/150‑1 (MQ=255)
                                                           ggCTACGGGAATTCATTCCCCCAAACACCCGCGCTTTTACCCACCGCCATGTGAAAAAAGCCCTTGGGGGGGGCTGTCGGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCttt                                                                                              <  1:1195636/149‑1 (MQ=255)
                                                              tACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAGAttttaaac                                                                                         >  1:1123604/1‑146 (MQ=25)
                                                                      ttCTTTCCCCCCAACACCCGCGAATTAACCCACCGCAATGGGAAAAAGGCCCTGCGGGGGGGTTGTCTGCCAGTTTTTTTTCCGCGGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTaa                                                                                            <  2:1632082/140‑1 (MQ=255)
                                                                         aTTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                                           >  2:269524/1‑93 (MQ=255)
                                                                         aTTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                                           <  1:269601/93‑1 (MQ=255)
                                                                         aTTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAACGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAACAGATTTaaaaac                                                                                        >  2:1473851/1‑137 (MQ=25)
                                                                          ttcccGCCAACACCCGCGCTTTTACCCAACGCCATGTGAAAAAATGCCTCCGCGGCGGCTGTCGGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttgt                                                                              <  1:661929/147‑1 (MQ=255)
                                                                                     aacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCCAGGGGAAATAACTTTagcaagtcgttct                                                                                 >  2:466029/1‑125 (MQ=25)
                                                                                     aacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAATTTTAACCCGTCCtgttgt                                                                              >  1:749339/1‑139 (MQ=37)
                                                                                        cgcgCATTTACCCACCGCATTGTGAAAAAATGCGTGCGGGGCGCCTGTCGGCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgct                                                               <  2:928454/151‑1 (MQ=255)
                                                                                         gcgcATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCCAAAGAAAAAAATTTAATCTCGT‑‑TGTTGTGTgggggtatcg                                                                  >  1:1473668/1‑137 (MQ=25)
                                                                                                    cAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAAATTTTAAACCCTCATGTTTTGTGGGCATCGCGGTTggggacca                                                       >  2:516924/1‑141 (MQ=25)
                                                                                                     aaCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGTGAGAGGAAAAAATTTATATCCGCCATGATGTGTGGGaa                                                                       >  1:39863/1‑130 (MQ=25)
                                                                                                      ccGCCAGGGGAAAAAATCCCTTCGCGGCGCCTTTCGGGCAGCTTTTTTCCCCCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGa                                                <  2:342882/150‑1 (MQ=255)
                                                                                                         ccATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAg                                             <  1:566042/151‑1 (MQ=255)
                                                                                                            ggtgAAAAAAGGCCTTCGCGGGGCCTGTCTGCCGGTTATTTTCCCGCCGAAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAg                                             <  2:654218/147‑1 (MQ=255)
                                                                                                             tgacAAAAATGCTTCCGCGGCGTCTTCCTGCAAGCTATTTTCCCGCCCAAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGATATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTTGAAAgg                                            <  2:778216/145‑1 (MQ=255)
                                                                                                             gtgAAAAAATGCTTCCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGaa                                         <  1:1525079/151‑1 (MQ=255)
                                                                                                             gtacAAAAATGCCTCCGCGGCGGCTTTCTGCCGGTTATTTTCCCGCCCCAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                          <  2:238843/147‑1 (MQ=255)
                                                                                                             gtaaaaaaaTCCCTTCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg                                           <  1:1045546/145‑1 (MQ=255)
                                                                                                                 aaaaaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAAAATTAAACCCCCCATGTTGTGTGGGGAACGCCGCTGGGAACAAAGTGGAAAaggggggg                                      >  1:1252837/1‑145 (MQ=255)
                                                                                                                 aaaaaTCCCTCCGCGGCGGTTGTCCCCCATTTGTTTTCCCGCCCAAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGt                                       <  2:1374865/149‑1 (MQ=255)
                                                                                                                 aaaaaGCCCTCCGGGGGGGCTGTCGCCCAGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgc                                     <  1:342298/151‑1 (MQ=255)
                                                                                                                   aaaCGCCTTCCGCGGCGGCGGTCTCCAGGCTTTTTCCCCGCCCAAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcg                                    <  2:1163311/151‑1 (MQ=255)
                                                                                                                       gCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAACCTTAAACCCCCCCTGTTTTGTGGGAATCCCCGCTTGGAAAAAAGGGGGAAGGGGAgggggggg                                 >  1:1186139/1‑142 (MQ=255)
                                                                                                                               gcggcgGCTGTCTGCCAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCa                            <  1:691808/148‑1 (MQ=255)
                                                                                                                                cggcggCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGCGGAAGAGAAAAAAATTTTAACCCGTAATGTTTTGGGGGGATCGGTGACTGTAACAATGAGGAAAAGGGAAgggggcgga                                >  1:1286499/1‑136 (MQ=255)
                                                                                                                                 ggggccTGTCTGCCTGTTATTTTCCCGCGAAACAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGaa                               <  2:1123426/138‑1 (MQ=255)
                                                                                                                                             ccAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa           <  1:699849/151‑1 (MQ=255)
                                                                                                                                                 cTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGGAAAAAAATTTAAACCCCTATGGTTTGGGGGGATCGCTGCTGGGGAAAAAGTGGGAAGGGGAGGGGGGGGGAAAAGCTCCaaaaaaggggtaa          >  1:289904/1‑139 (MQ=255)
                                                                                                                                                  ttttttcCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTaa          <  2:124991/145‑1 (MQ=255)
                                                                                                                                                    ttttcCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAtt        <  1:350110/147‑1 (MQ=255)
                                                                                                                                                    ttttcCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg     <  1:633655/150‑1 (MQ=255)
                                                                                                                                                      tttCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCCAGAGAAAAAACTTTTAACCCCGCCTTGTTTGTGGGCAACGCTGCTGGGGACAATGTGGGGAGGGGAGGGGGGGGGAACAACCCCCGAAACGCGCCAATTTAGttt  >  1:617385/1‑150 (MQ=255)
                                                                                                                                                        |                                                                                                                                                         
ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTC  >  NC_000913/357730‑358028

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