Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA56.0% 14.2 / 72.6 50L267F (TTG→TTTcynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (16/6);  new base A (0/28);  total (16/34)
Fisher's exact test for biased strand distribution p-value = 1.52e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.25e-04

TGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGA  >  NC_000913/357741‑358038
                                                                                                                                                         |                                                                                                                                                      
tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                                                        >  2:26463/1‑150 (MQ=255)
 gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa                                                                                                                                                      <  2:435296/151‑1 (MQ=255)
    aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaaa                                                                                                                                                    <  2:37153/150‑1 (MQ=255)
    aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaaa                                                                                                                                                    <  1:352/150‑1 (MQ=255)
     aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                                                  >  1:313012/1‑151 (MQ=255)
                  aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                                      >  1:713679/1‑150 (MQ=255)
                    gCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                                   >  2:41566/1‑151 (MQ=255)
                    gCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                                   >  1:1099945/1‑151 (MQ=255)
                     cAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                                   >  2:1089009/1‑150 (MQ=255)
                         tCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTgg                                                                                                                              >  1:1328342/1‑151 (MQ=255)
                          ccacccccaccctggtgtCAGCGTCTCCGGTTATTCTTTCCCGCCAACACCCGCGCTTTTATCCACCGCTATGTGAAAAAATGCCTTCGCGGCGGCTGTTTGCCAGTTTTTTTTCCCCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTg                                                                                                                               <  2:1104056/135‑1 (MQ=255)
                          ccAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                                                             >  1:929754/1‑151 (MQ=255)
                              ccGCAGCCCTGTGCCAGCGGCTCCCGTGTTTCTTTCCCGCCAACACCCGCGCTTTTACCCACCCCCTTTTCAAAAATGCCCTTGGGGGCGGCTTTCGGCCAGCTTTTTTCCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                            <  2:688148/148‑1 (MQ=255)
                               cgcagcctttgctCAGCGCCTACCGGGTTTCATTCCCGCCAACACCCGCCCTTTTATCCACCGCATTGTGCAAAAATGCCTTCGGGGCGGCTGTCTGCAAGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                            <  2:398269/135‑1 (MQ=255)
                                       tgtgTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAagcgga                                                                                                                 >  2:443109/1‑146 (MQ=255)
                                        gtgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGaaa                                                                                                               >  2:275825/1‑150 (MQ=255)
                                        gtgccggcgcCTCCCGTATTTCTTTCCCCCCAACACCCGCTCTTTTACCCACCGCCAGGTAAAAAAATGCCTCGGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAagag                                                                                                                   <  1:710642/141‑1 (MQ=255)
                                          gtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggtcacaatgacaaa                                                                                                              >  1:326928/1‑137 (MQ=255)
                                           tCAGCGGCTCCCGTGTTTCATTCCCCCCAACACCCGCGCTTTTACCCACCGCCATGGGAAAAATGGCCTTCGCGGCGCCTGTCTCCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAg                                                                                                            <  2:31090/150‑1 (MQ=255)
                                              gcgcCTACCGTTTTTCATTCCCGACAACACCCGCGATTTTACCAAACCCAGTGTGAAAAAATCCCTCCGCGGCGGCTTCCGGCAAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaaa                                                                                                               <  1:616661/144‑1 (MQ=255)
                                                gcctcccGTGATTCATTCCCCCCAACACCCGCCCTTTTACCCACCGCAATGTAAAAAAATGCCTGCGGGGCGCCTGTCTGCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCtt                                                                                                         <  1:191507/143‑1 (MQ=255)
                                                gcctcccGTGATTCATTCCCCCAAACACCCGCGCTTTTCCCCACCGCCATGTGAAAAAAGGCTTCCGGGGCGGTTGTCGGCCGGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAag                                                                                                                     <  1:199946/132‑1 (MQ=255)
                                                gccTACCGTAATTCATTCCCCCCACCACCCGCGCATTTACCCACCGCAAGGTGAAAAAATGCCTTCGGGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATa                                                                                                             <  2:1269329/142‑1 (MQ=255)
                                                  ctcccGTGATTCATTCCCCCCAACACCCGCGCTTTTACCAACCCCAATGTGAAAAAAGCCCTGCGGGGGGGCTGTCTGCCAGTTTTTTTCCCGCGCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTa                                                                                                       <  1:980603/145‑1 (MQ=255)
                                                  caaCCGTAATTATTTCCCGCAAACACCCGCAAATTATTCCACCGCCGTGAAAAAAATCCCTTCGGCGGCGGGTGTCGGCCAGATATTTTCCCCCCTCACAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGc                                                                                                    <  2:436222/149‑1 (MQ=255)
                                                             cATTCCCCCAAACACCCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGGGGCGGTTGTCTGCCAGCTTTTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgt                                                                                           <  2:537106/149‑1 (MQ=37)
                                                               ttccccccAACACCGGCGTATTTACCCACCCCCATGTAAAAAATTGCCTCCGCGGCGGCTTTCGCCCAGTTTTTTTCCCCCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTc                                                                                               <  2:714756/137‑1 (MQ=255)
                                                               ttcccGCCAACACCCGCGCATTTACCCAACGCCATGTGAAAAAATGCCTTCGGGGCGGCTGTCTGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgtt                                                                                          <  2:1073043/145‑1 (MQ=255)
                                                               tccccGCCAAAACCCGCCAATTTACCAACCCCCCTGTGCAAAAACGCTTGCGGGGCGGTTGTCCGCCATCTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtg                                                                                            <  2:234688/143‑1 (MQ=255)
                                                                             cgcacttttACCCACCCCCATGTAAAAAAATGCCTCGGCGGCGGCTGTCTCCCAGCTTTTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgct                                                                         <  2:396197/145‑1 (MQ=255)
                                                                                   tttACCCACCGCCATGTGAAAAAAGCCTTGCGGGGCGGCTGTCTGCCGTTTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGt                                                                     <  1:753979/148‑1 (MQ=255)
                                                                                          aaCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGGGAAAAAAGATTTAATCCCTCATTTTTTGTGGGGATCgctccgggggacaa                                                                 >  1:1290858/1‑135 (MQ=25)
                                                                                           aaccctttgtaaaaaaaTGCCTCCGCGGGGGCTTTCGGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg                                                              <  2:834062/136‑1 (MQ=255)
                                                                                                caggtcgaaaaaTGCCGGCGGGGCGGCTGTCTGCCGGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAgg                                                      <  2:1071704/144‑1 (MQ=255)
                                                                                                   gtgAAAAAGTCCCTCGGCGGGGGTTGTCTCCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                                    <  2:1199636/150‑1 (MQ=255)
                                                                                                       aaaaaTGCCTCCGCGGCGGCTTTCTGCCATTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAg                                                  <  2:249186/148‑1 (MQ=255)
                                                                                                       aaaaaTCCCTTCGCGGCGGCTTTCTCCAAGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGt                                                 <  2:571262/149‑1 (MQ=255)
                                                                                                       aaaaaTCCCTCCGGGGGGGCTGTCTCCCATTTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgc                                               <  1:1202249/151‑1 (MQ=255)
                                                                                                       aaaaaGGCCTGCGGGGCGGTTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaa                                                         <  1:1300540/141‑1 (MQ=255)
                                                                                                        aaaaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGGAAAAAATTTAAAGACCCCCTTTTTTGTGGCGATCGCTGGTGGGGACCAATTGGAAAGGGAAggg                                                >  2:396467/1‑147 (MQ=255)
                                                                                                           aatcccgtgcgcgCCGCCTGTCTCCCTCCTTTTTTCCCCCCGCAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGTGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcgg                                             <  1:474821/142‑1 (MQ=255)
                                                                                                           aacgctttCGGGGGGGCTGTCTGCCGTCTTTTTTCCCGCCAAAACAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggc                                            <  2:807196/145‑1 (MQ=255)
                                                                                                                        ggcggcTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAGAAAAAAATTTTAGACCGGCCTTTTTTGTGGGAATCGCTGCTGGGTAAAAAGTGGGAAAGGGAAGGGGGGGGaaaaacccc                                  >  2:1129257/1‑144 (MQ=255)
                                                                                                                              tGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAACTTTTAACCCCTCCTGTTGTGTGGGGATCGCTGCTGGGAAAAAAGTGGGAAAGGGAGGGGGGGGGATCAACTCCaaaaac                           >  2:477572/1‑149 (MQ=255)
                                                                                                                                    ccagttttttccccgccgCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc                           <  1:1218524/133‑1 (MQ=255)
                                                                                                                                      gtttttttgcccgccaAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa                     <  1:393044/139‑1 (MQ=255)
                                                                                                                                           ttttcCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg               <  2:832256/150‑1 (MQ=255)
                                                                                                                                           ttttcCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCa                                      <  2:671595/127‑1 (MQ=255)
                                                                                                                                               tCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGGAAAAAATTTTAACCCGGCATGTTTTGTGGGGAACGCCGCTGGGTACAAAGTGGAAAGGGGAGAGGGGGGGATCAACCCCCCAAACGCGCTCATTGAGtttg           >  2:643437/1‑148 (MQ=255)
                                                                                                                                               tCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAAATTAAACCCCTCATTTTTTGTGGGCATCGCTGCCGGGAAAAAAGTGGGAAGGGGAGGGGGGGGGAGCAGCCCCACAACCGCGCCCATTGAGt              >  1:851969/1‑147 (MQ=255)
                                                                                                                                                ccgccaAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCg           <  1:60283/149‑1 (MQ=255)
                                                                                                                                                ccgccaAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCg           <  2:226670/149‑1 (MQ=255)
                                                                                                                                                ccccctcaaaaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc            <  2:242294/142‑1 (MQ=255)
                                                                                                                                                       aacaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGa  <  2:1179352/151‑1 (MQ=255)
                                                                                                                                                         |                                                                                                                                                      
TGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGA  >  NC_000913/357741‑358038

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