Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA45.1% 36.5 / 27.9 52L267F (TTG→TTTcynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (18/10);  new base A (0/23);  total (18/34)
Fisher's exact test for biased strand distribution p-value = 5.73e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAA  >  NC_000913/357740‑358035
                                                                                                                                                       |                                                                                                                                                     
cTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                                                       >  2:919363/1‑151 (MQ=255)
    tAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                                                                                                    >  1:951993/1‑150 (MQ=255)
     aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaaa                                                                                                                                                   >  1:132722/1‑150 (MQ=255)
      aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGAGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                                                  >  2:884333/1‑150 (MQ=255)
       gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccg                                                                                                                                                <  1:1008604/151‑1 (MQ=255)
              gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtct                                                                                                                                         >  2:360138/1‑151 (MQ=255)
                           ccAGCCGCAGCCCTGTGTCAGCGGCTCCCGTGATTCATTTCCGCCAACACCCGCGCTTTTACCCACCGCCATGTGCAAAAATGCTTCCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTg                                                                                                                              <  2:1006094/149‑1 (MQ=255)
                             aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                          >  1:758245/1‑151 (MQ=255)
                               ccccgcaCCTTTGACAGCGGTTCCGGTTATTCATTCCCCCAAACACCCCCACTTTTCTCCACCCCCATGTCAAAAAAGGCCTTCGCGGCGGCTGTCGGCCATTTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                          <  1:115071/143‑1 (MQ=255)
                               ccccAGCCCTGTTTCCGCGGCTCCCGTGTTTAATTCCCGCAACCACCCGCGATTTTCTCCACCGCCATGGGAAAAAATCCCTTCGGGGGGGTTGTCTGCCAGTTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                        <  1:57281/148‑1 (MQ=255)
                               ccGCAGCCCTGTGTCAGCGGCTCCGGTGATTCATTCCCGCCAACACCCGCGCTTTTACCCAACGCCATGTGAAAAAATGCCTCCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCagta                                                                                                                                 <  2:40550/142‑1 (MQ=255)
                               ccGCAGCCCTGTGTCAGCGGCTCCCGTGTTTCATTCCCCCCAACACCCGCGCTTTTACCCAACCCCTTGTGAAAAAATGCCTTCGCGGGGGCTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                           <  1:919383/148‑1 (MQ=255)
                               ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                        >  2:258910/1‑151 (MQ=255)
                                  cAGCCCTGTGTCGGGGGCTCCCGTGATTCATTCCCGCAAACACCCCCGCTTTTACCCACCGCCTTGTGAAAAAATGCCTCCGGGGGGGCTGTCTCCCGTCTATTTTTCCGCGGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                          <  1:950404/146‑1 (MQ=255)
                                   aGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggc                                                                                                                      >  2:8990/1‑147 (MQ=255)
                                            tCAGCGGCTCCCTTGATTCTTTCCCCCCAACACCCGCCCTTTTACCCAACCCCATGTGAAAAAATGCCTCCGCGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAg                                                                                                          <  2:33752/150‑1 (MQ=255)
                                            tCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGaaaaa                                                                                                           >  2:302734/1‑147 (MQ=255)
                                                cGGCTCCCGTGATTCATTTCCCCCAACACCCGCGCTTTTACCCACCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCGTTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCttt                                                                                                     <  2:503977/150‑1 (MQ=25)
                                                 ggCTCCCGTGATTCATTCCCGCCAACACCCGCGCATTTACCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTa                                                                                                    <  2:237164/150‑1 (MQ=255)
                                                   ctcccGTGATTCATTCCCCCCAACACCCGCGCTTTTACCCAACGCCATGTGCAAAAATGCCTTCGCGGGGGCTGTCTGCCAGCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTaa                                                                                                   <  1:823704/146‑1 (MQ=25)
                                                   cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAACttt                                                                                                     >  2:628029/1‑147 (MQ=255)
                                                    tcccGTGATTCATTCCCGCAAACACCCGCGCATTTACCCACCGCCATGTGCAAAAAGGCCTCCGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAg                                                                                                  <  1:938377/147‑1 (MQ=255)
                                                      ccGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGaaaaaaaatttaaac                                                                                                 >  1:61807/1‑137 (MQ=25)
                                                        gTGATTCTTTCCCGCCAACACCCGCGCATTTACCCAACGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGt                                                                                             <  2:637840/150‑1 (MQ=255)
                                                                 ttCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCATGACGAAAAAGCTTTaaaacatcatctt                                                                                        >  1:161670/1‑136 (MQ=25)
                                                                  cccGCCAACAACCGCGCATTTATCCAACGCCATGTGAAAAAATGCCTCCGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtgtg                                                                                  <  2:53662/150‑1 (MQ=255)
                                                                   cccccacaacccGCGCATTAACCCAACGCATGGAGAAAAAATGCTTCCGCGGCGGCTGTCGGCCAGTTTTTTTCCCGCCAAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtg                                                                                    <  1:461832/138‑1 (MQ=255)
                                                                       caacaaCCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAATTTAAACCCGGCATGTTGTGTGGGaa                                                                              >  2:360363/1‑150 (MQ=37)
                                                                         acaacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGGGAAAAAAATTTAAAACCCCCCTGTTGTGTgggga                                                                              >  2:92488/1‑145 (MQ=25)
                                                                            acCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCATAAGAAAAAAGCTTTAAACCCCCTTGTTGTGTGGcgaacc                                                                           >  2:47901/1‑143 (MQ=25)
                                                                             ccGCGCTTTTACCCACCGCCTTGTGAAAAAATCCTTCCGCGGCGGCTTTCCGCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc                                                                            <  1:779138/146‑1 (MQ=255)
                                                                             ccGCGCTTTTACCCACCCCCTTGTGAAAAAATGCCTTCGGGGCGGCTTTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgc                                                                       <  1:280207/151‑1 (MQ=255)
                                                                                 gcttttACCCAACCCCATGTCAAAAAATCCTTTCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGt                                                                   <  2:164231/148‑1 (MQ=255)
                                                                                       aTCCACCGCCATGTGCAAAATTCCCTTCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACa                                                               <  2:405279/149‑1 (MQ=255)
                                                                                          caccGCCATGTGAAAAATTCCTTCCGCGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgtg                                                         <  2:698493/148‑1 (MQ=255)
                                                                                           aaCGCCATGTAAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTgg                                                        <  2:892476/151‑1 (MQ=255)
                                                                                                    tgaaaaaaTGCCTTCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                                 <  2:776130/146‑1 (MQ=255)
                                                                                                       aaaaaTGCCTTCGCGGCGGCTGTCTCCCAGCTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                                 <  1:755450/146‑1 (MQ=255)
                                                                                                         aaaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGAAAGAGAAAAAGATTTTAACCCCCCCTGTTTTGTGGGGAACGCTGCATGTAACAATGGGGGAAAGGGagggggg                                           >  2:910811/1‑145 (MQ=255)
                                                                                                             gCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAACTTTTAACCCCCCCTGTTGTGTGGGGATAGCTGCTGGGAACAAAGTGGAAAGGGGAgggggggg                                        >  2:873545/1‑142 (MQ=255)
                                                                                                              ctttCGGGGCGGTTGTTTGCTAGATTTTTTCCCCCCGAAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGaa                                      <  2:629261/148‑1 (MQ=255)
                                                                                                               cttccgcgGCGGCTGTCTGCCGGCTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGa                                       <  2:374370/145‑1 (MQ=255)
                                                                                                                    gcggcgGCTTTCTGCCAGTTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCa                                   <  2:80058/148‑1 (MQ=255)
                                                                                                                    gcggcgGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGCTCTCTCCAGTAGTGGCGGGGGCAAAGAAAAAACATTTAAACCCCCATGTTGTGTGGGGAACTCTGGTGGTAAAAAAGTGGGAAGGGGAgggggg                                          >  2:436789/1‑136 (MQ=255)
                                                                                                                      ggcggcTTTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCt                                <  2:296178/149‑1 (MQ=255)
                                                                                                                      ggcggcTTCCTGCCGGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTc                               <  1:629295/150‑1 (MQ=255)
                                                                                                                              tcggCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc                        <  2:823446/146‑1 (MQ=255)
                                                                                                                                  ccAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgcg                     <  2:686178/148‑1 (MQ=255)
                                                                                                                                  ccAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc                      <  2:1020616/147‑1 (MQ=255)
                                                                                                                                  ccAGCTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc                      <  2:647921/147‑1 (MQ=255)
                                                                                                                                          tttcCCGCCAAAAGAAAACCGCCGTTCTCGCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGt           <  1:884232/150‑1 (MQ=255)
                                                                                                                                                gccgcAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGGAATAACATTAAAACCCTCCTTGTTGGTGGGAAAGGGTGGTGGGAAAAAAGTGGGGAGGGGAATGGGGGGGATGAGCTCCCAAACACCGCTAAATGAggtag        >  2:144759/1‑143 (MQ=255)
                                                                                                                                                 ccgcAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCa   <  2:663544/151‑1 (MQ=255)
                                                                                                                                                    caacaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAAATTAAAACCCCCCCGTTGTTTGGGGATAGCTGCTGGGGAAAAAATGGGAAGGGAGGGGCGGGGGAACAGCTCCAGAACACCGCTAATTTAGTTTGCCCCaa  >  1:484517/1‑149 (MQ=255)
                                                                                                                                                    caacaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAATTTTAACCCCCCCCTTTGGGTGGGGAGCGCTTCTGGGAAAAAAGTGGGAAAGGGAGTGCGGCGGATCAGACCCCCAAACACGCCAAATGAGTTGGCCTcca  >  1:920282/1‑147 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
CTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAA  >  NC_000913/357740‑358035

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: