Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8800TC44.0% 57.7 / 20.0 50K271E (AAA→GAA) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (19/9);  new base C (0/22);  total (19/31)
Fisher's exact test for biased strand distribution p-value = 2.34e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGT  >  NC_000913/357733‑358026
                                                                                                                                                      |                                                                                                                                                      
cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTcc                                                                                                                                                      <  2:42181/150‑1 (MQ=255)
 gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgc                                                                                                                                                    >  1:501254/1‑151 (MQ=255)
 gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgc                                                                                                                                                    >  1:272748/1‑151 (MQ=255)
   aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgccg                                                                                                                                                  <  2:594949/151‑1 (MQ=255)
   aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                                                                                                   <  2:1017791/150‑1 (MQ=255)
    ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCcgccgc                                                                                                                                                 >  2:1094901/1‑151 (MQ=255)
        tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                                             >  1:518230/1‑150 (MQ=255)
          aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCCGCCCTGTGCCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                                                                                          <  1:277835/151‑1 (MQ=255)
            aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaaa                                                                                                                                         <  2:17119/150‑1 (MQ=255)
            aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                                        >  1:123078/1‑151 (MQ=255)
                    cgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtc                                                                                                                                >  2:1051340/1‑151 (MQ=255)
                       aTCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctc                                                                                                                              >  2:911628/1‑150 (MQ=255)
                          aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                           >  1:1062720/1‑150 (MQ=255)
                                  ccAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                                                  >  1:358491/1‑151 (MQ=255)
                                  ccAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                                                  >  1:66079/1‑151 (MQ=255)
                                          agccttgtgTCAGCGGCTCCCGTGATTCATTCCCGCCAACACCCGCGCTTTTAACAAACCCTTTGTGAAAAAATGCTTTCGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCAAAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCaa                                                                                                          <  2:1081186/145‑1 (MQ=255)
                                           gACCTGTGCCAGCGCCTCCCGTGATTCATTCCCGCAAACACCCGCGCATTTACCCACCGCAAGGTGAAAAAAGGCCTCCGGGGCGGCTGTCCGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCa                                                                                                           <  2:26073/148‑1 (MQ=255)
                                            acttttgtCAGCGCCTCCCGTGTTTCTTTCCCCCCAACACCCGCGCTTTTACCCACCGCCATGTGAAAAAAGGCTTCCGGGGCGGCTGTCTGCCATTTTTTTTTCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                              <  2:391629/140‑1 (MQ=255)
                                             ccTGTGTCAGCGGCTACCGTGATTCTTTCCCCCCAACACCCGCGCATTTACCCACCCCCTTGTAAAAAAATCCCTCCGGGGCGGCTGTCTCCCACCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAgaga                                                                                                    <  2:268722/151‑1 (MQ=21)
                                              cTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggagggt                                                                                                           >  1:997022/1‑141 (MQ=255)
                                               tgtgTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGAAagagga                                                                                                    >  2:309451/1‑148 (MQ=255)
                                                  gtCGGCGGCTCCCGTAATTCTTCCCCCCAAACACCCCCGATTTTACCCACCCCCATGGAAAAAATTCCTTCCGGGGCGGCTGTCTCCCGGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaa                                                                                                   <  1:469523/147‑1 (MQ=255)
                                                    cggCGGCTCCCGTTATTCATCCCCGCCAACCCCCGCGCTTTTCCCCACCGCCATGGAAAAAATGGCCTCCGCGGGGGCTGCCGCCCAGTTATTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAg                                                                                               <  1:385346/147‑1 (MQ=255)
                                                      gCGGCTACGGTAATTCTTTCCCGCAAACACCCCCGTTTTAACCCAACCCCATGTCAAAAAATGCCTTCGGGGCGGCTGTCGGCCTGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATa                                                                                                <  2:939751/146‑1 (MQ=255)
                                                            aCCGTGATTCATTCCCCCCAACACCCCCGCTTTAACCCACCGCGGGGTGAAAAAATCCCTAGGCGGCGCTCGTCTGCCAGCTTTTTCCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTa                                                                                          <  2:505210/146‑1 (MQ=255)
                                                             ccGTAATTCTTTCCCCCAAACACCCCCGCTTTTCCCCAACCCCATGTGAAAAAAGGCTTTGGCGGCGGCTGCCGGCAGGCTTTTTTCCCGCGGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCttt                                                                                           <  1:576002/144‑1 (MQ=255)
                                                                 taattctttCCCCCAACCACCCCCGCTTTTACCCCCCGCCAGGTAAAAAAACCCCTCCGGGGCGCTTGTCGCCCGTTTTTTTTCCCCCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc                                                                                      <  1:792276/140‑1 (MQ=255)
                                                                   ttCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAGGAAAAAGTTTAAAccccccctct                                                                              >  2:456940/1‑144 (MQ=25)
                                                                     cATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggcgggagaggga                                                                                                     >  1:187500/1‑118 (MQ=255)
                                                                       ttcccccaaacaccCGCCCTTTTACCCAACGCCATGGAAAAAAATCCCTCCGCGGCGGTTGTCTCCCAGTTTTTTTCCCGCCAAAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAt                                                                                <  1:197719/137‑1 (MQ=255)
                                                                           cGCCACCACCCGCGCTTTTAGCCAACGCGATTGGAAAAACGGCCTTCGCGGCGGCTGTCTCCCATTTATTTTCCCCCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgtt                                                                             <  2:1047781/144‑1 (MQ=255)
                                                                             ccAACACCCGCGCTTTTACCCACCCCCATGTGAAAAAATGCCTCCGGGGGGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCa                                                                     <  2:368693/150‑1 (MQ=255)
                                                                               aacaccCGCGCATTTACCCAACGCCTTGTGAAAAAATGCCTTCGGGGCGGCTTTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc                                                                  <  2:209819/151‑1 (MQ=255)
                                                                                     ccGCGTATTTACCCACCCCCTTGTAAAAAAGTCCCCCCGCGGCGGCTCCCTGCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgct                                                               <  1:1051245/149‑1 (MQ=255)
                                                                                              tacccAACGCCATGTGCAAAATTGCCTCCGCGGCGGCTGTCTGCCAGTTTTTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACa                                                     <  2:467886/147‑1 (MQ=255)
                                                                                                 ccACCGCCAGGTGAAAAAATCCCTTCGGGGCGGCTTTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg                                                  <  2:501296/150‑1 (MQ=255)
                                                                                                   acccccttgtgAAAAAATGCCTTCGGGGCGGCTGTCTGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgtg                                                <  2:781045/143‑1 (MQ=255)
                                                                                                       ccATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGGAAAAAAATTTAACCCCCCCTGTTTTGTGGGGATCGCTGCTggggaaaaaatgg                                               >  1:1015016/1‑136 (MQ=21)
                                                                                                            tgCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAAGTTAAATCCGCCCAGTTTTGTGGGAAACGCCGTTGGGGACAATGAGGAAAggg                                         >  2:312261/1‑148 (MQ=255)
                                                                                                               aaaaaCGCCTCCGGGGGGGCTGTCGCCAATCTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                        <  1:500222/146‑1 (MQ=255)
                                                                                                                  aaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCCAACGGAAAAACGTTTAATCCCTCCCGTTGTGTGGGGAACGCTGCTTGTGACAATGTGGGAAGGGGAgagggg                                 >  2:625971/1‑145 (MQ=255)
                                                                                                                    aacgctctggcggcgGCTGCCGGCCAGTTTTTTCCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggc                                <  1:363883/142‑1 (MQ=255)
                                                                                                                          tcggcggggGCTGTCTCCCGGTTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGaa                             <  1:434163/145‑1 (MQ=255)
                                                                                                                             gggggggcTGCCGGCCATTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCa                          <  1:524112/143‑1 (MQ=255)
                                                                                                                               ggcggcTGTCTGCCATCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCa                          <  1:873583/146‑1 (MQ=255)
                                                                                                                                   gcTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACAGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAGAAAAAAATTTTAACCCCCCCTGTTTTGTGGGGATCGCTGCTGGTTACAAAATGGGAAGGGGAGGGGGGGGGATCaaatcc                     >  2:489899/1‑142 (MQ=255)
                                                                                                                                        cTGCCGGCTTTTTTTCCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                <  1:994824/147‑1 (MQ=255)
                                                                                                                                         tgcaagtttttttcCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggcg                               <  2:912005/122‑1 (MQ=255)
                                                                                                                                          gCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAAAAACTTTAAACCCCCCTGTTGTGTGGGGAGTGCTTCTGGGAACAATGTGGGGAGGGGAGGGGGGGGGAGCAACTCCCAAAACgcgc           >  2:782822/1‑150 (MQ=255)
                                                                                                                                           ccAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc               <  1:728133/145‑1 (MQ=255)
                                                                                                                                             acttATTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa         <  2:384274/146‑1 (MQ=255)
                                                                                                                                                 tttttcCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTaa        <  2:262791/145‑1 (MQ=255)
                                                                                                                                                  tttttCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAAAAAATTTAAACCCCCCTGTTGTTTGGGGATCCCCGCCGGGGACAAAGTGGGAAGGGGAGGGCGGGGGATCAGCTCCACAACCGCCCCAAtt      >  2:1003184/1‑147 (MQ=255)
                                                                                                                                                   ttttCCGCCGAAACAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGt  <  2:76670/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
CGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGT  >  NC_000913/357733‑358026

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