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

GTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGA  >  NC_000913/357734‑358024
                                                                                                                                                    |                                                                                                                                                   
gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTccgc                                                                                                                                                 <  1:1162370/151‑1 (MQ=255)
  aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTccgccg                                                                                                                                               >  2:150231/1‑151 (MQ=255)
  aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                                                                                                <  2:499723/150‑1 (MQ=255)
       tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                                           >  2:1029316/1‑150 (MQ=255)
        gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaaca                                                                                                                                          >  2:530222/1‑150 (MQ=255)
           aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑c                                                                                                                                     >  1:546255/1‑151 (MQ=255)
            aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑cc                                                                                                                                    >  1:770908/1‑151 (MQ=255)
              cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑ccgc                                                                                                                                  <  1:23165/151‑1 (MQ=255)
                tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑ccgccg                                                                                                                                >  1:227780/1‑151 (MQ=255)
                   cgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGtt                                                                                                                              >  1:1062702/1‑150 (MQ=255)
                     cATCAGGCATTTCCAGCCGCAGACCTGTGTCAGCGGCTCCCGTGATTCATTTCCGCCAACACCCGCGCATTTACCCAACGCAATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAACAAAA‑CCGCCGTtc                                                                                                                             <  1:424512/149‑1 (MQ=255)
                     cATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTtctc                                                                                                                           >  1:277923/1‑151 (MQ=255)
                        cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTcc                                                                                                                        >  2:1049059/1‑151 (MQ=255)
                             aaTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCagt                                                                                                                     >  1:245153/1‑149 (MQ=255)
                             aaTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCagt                                                                                                                     >  1:1072101/1‑149 (MQ=255)
                              aTTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATATTTTCGCCGGAACAAAA‑CCGCCGctctctaccctag                                                                                                                   >  1:952123/1‑143 (MQ=11)
                              aTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCagtagt                                                                                                                  >  2:830113/1‑151 (MQ=255)
                                  cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGc                                                                                                               >  2:356875/1‑150 (MQ=255)
                                         aGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑gggggg                                                                                                        >  2:678940/1‑147 (MQ=255)
                                            cttgttcctgCGGCTCCTGTTTTTCTTTCCCGCAAACACCCGCGAATTCCCCCACCGCCTTGCGAAAAAATC‑CTTG‑CGCGGCGGCTGTCCGCCAGTTTTTTTCCCCCCCCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑gggg                                                                                                          <  2:771676/135‑1 (MQ=255)
                                               ttgacggcgcccaCCGTG‑TTCCTTTCCCCCAACCACCCTCTCTTTAACCAAGCCCAATTTCAAAAAAG‑CCTTCCGGGGCGGCTGGCTCCCATTTTTTTTCCCTCCCAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAgaga                                                                                                   <  1:1026974/136‑1 (MQ=255)
                                               gtgtCAGCGGCTACCTTGATTCTTTCCCCCCAACACCCGCCCATTTACCCAACGCCATGTGAAAAAATG‑CCTT‑CGGGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAag                                                                                                      <  1:902492/145‑1 (MQ=255)
                                                 gtCGGCGGCTCCCGTTTTTCTTCCCCCCCAACACCCGCCCTTTTACCCACCGCCATGTAAAAAAATG‑CCTC‑CGGGGCGGCTGTCGCCCAGTTTTTTTCCCGCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAg                                                                                              <  2:1127381/151‑1 (MQ=255)
                                                      cGGCTCCCGTGATTCATTCCCGCCACCACCCGCCCTTTCACCCAACCCCATGT‑AAAAAATGCCCTC‑CGCGGCGGCTGTCGGCCGGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCttt                                                                                          <  2:461991/150‑1 (MQ=255)
                                                      cGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAGgaaaaaactt                                                                                          >  2:540759/1‑143 (MQ=25)
                                                            cggTGATTCTTTCCCCCCAACACCCCCGCTTTTACCCACCGCCTTGTGAAAAAATG‑CCTG‑CGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGcc                                                                                    <  1:265176/147‑1 (MQ=255)
                                                            ccgtttttCATTCCCCCCAACACCCGCCCTTTTCCCCAACGCCTTGTGAAAAAATC‑CCTT‑CGGGGGGGCTGTCTCCCATTTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATa                                                                                               <  1:1181730/133‑1 (MQ=255)
                                                              gTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGGGAAAAAGAATT‑aaatc                                                                                    >  2:833654/1‑144 (MQ=25)
                                                                  ttcttcccccccAACACCCGCGCTTTTCCCCACCGCAATGAAAAAAAATG‑CCTG‑CGGGGCGGCTGTCCGCCGGTTTTTTTCCCGCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGc                                                                                     <  2:384/132‑1 (MQ=255)
                                                                        ccccccAACACCCGCGCTTTTCCCCACCCCTTTGTAAAAAAATG‑CCTG‑CGGGGGGGCTGTCTCCCGTTTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGCCCGTCATGTtgtgtg                                                                      <  2:1078359/147‑1 (MQ=255)
                                                                         ccGCCAACACCCGCG‑ATTTTCCCAACCCCTTTTGGAAAAAAGCCCTC‑CGGGGCGGCTGTCCGCCGTCTTTTTTCCCCCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGCCCGTCATGTtgtgtg                                                                      <  1:168647/150‑1 (MQ=255)
                                                                                 acccGCGCATTTACCCAACGCCATGTGAAAAAATC‑CCTT‑CGCGGCGGCTGTCGGCCAGTTTTTTTCCCCCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGCCCGTCATGTTGTGTGGCAATCg                                                               <  1:1057613/147‑1 (MQ=255)
                                                                                   ccgccattttACCCAACCCCATGAAAAAAAAGC‑CCTC‑CGGGGCGGTTGTCTCCCGTTTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGCCCGTCATGTTGTGTGGCAATCgctgc                                                           <  2:289570/144‑1 (MQ=255)
                                                                                   ccGCGCTTTTATCCACCGCCATGTGAAAAAATG‑CCTC‑CGCGGCGGCTGTCGGCCATTTTTTTTCCCCCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGCCCGTCATGTTGTGTGGCAATCgctg                                                            <  1:562994/150‑1 (MQ=255)
                                                                                    cgcgCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCGAGAGGAAAAAATTT‑TATCCCCCCCTGTTTTGTGGGGATCggt                                                             >  1:109606/1‑146 (MQ=25)
                                                                                      cccATTTACCCACCCCCATGTAAAAAAATG‑CCTC‑CGCGGGGGCTGTCTGCCAGTTTTTTTCCCCCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGCCCGTCATGTTGTGTGGCAATCgctgct                                                          <  1:36258/147‑1 (MQ=255)
                                                                                       cctttcaccaaccGCTATGTAAAAAAATC‑CCTC‑CGCGGCGGTTGTCGGCCATTTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGt                                                       <  2:965972/139‑1 (MQ=255)
                                                                                        cATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAGAGAAAAAAATTT‑TAACCCCTCATGTTGTGTTGGTATCGCCGCTGGtg                                                      >  2:598460/1‑150 (MQ=37)
                                                                                                     ccATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGCAAAAAAAAAAAGTTT‑AAACCCG‑CCTTTTTTGTGGGGATCCGTGCTGGGGACAAAGTgg                                             >  1:281691/1‑147 (MQ=19)
                                                                                                          gtaaaaaaaG‑CCTTCCGGGGGGGCTGTCTCCCATTTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                      <  1:1107354/147‑1 (MQ=255)
                                                                                                            cAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAGAAAAAAATTTAAAACCCTCCATGTT‑TGTGGGAATCGCCGCCGGGAAAAAAGTGGGAAggg                                       >  1:2460/1‑146 (MQ=255)
                                                                                                             aaaaaGG‑CCGT‑CGCGGCGGCTGTCTGCCTGTTATTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAg                                    <  2:677445/148‑1 (MQ=255)
                                                                                                              aaaaAG‑CCTTCCGCGGCGGTTGTCGGCAACCTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgc                                 <  1:336790/151‑1 (MQ=255)
                                                                                                               aaaTG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGGCAGAGAAAAAAATTT‑TAATCCCCCCTTTTGTGTGGGGATTCCTGCTGGTGACAAAGTGGGAAGGGGAgggggg                               >  2:825293/1‑146 (MQ=255)
                                                                                                                           gggggggcTGTCCCCCATTTATTTTCCCCCCACAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggcg                             <  1:527057/138‑1 (MQ=255)
                                                                                                                            cggcggCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAGAAAAAACTTT‑TAAACCCCCCTCTTTTTTGGGGATCGCTGCTGGGTACAAATTGGGAAAGGGAGGgcgggggaaaaaccc                    >  1:449998/1‑140 (MQ=255)
                                                                                                                               cggcTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAGAGAAAAAGATTT‑TAACCCGTCCTGTTGTGTGGGAATAGCTGCTGGTTAAAAAGTGGGAAGGGGAGGggggcggaaaaaccc                    >  1:1160813/1‑139 (MQ=255)
                                                                                                                                 gcttccTGCCGGTTTTTTTCCCCCAGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCa                        <  2:560124/137‑1 (MQ=255)
                                                                                                                                     ccgccaagtttttttcCCCCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc             <  2:1072092/138‑1 (MQ=255)
                                                                                                                                        gcccattttttttCCCGCACAAAAAAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgcg          <  2:227779/142‑1 (MQ=255)
                                                                                                                                           agttATTTTCCAACCGCAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTaa      <  1:181568/147‑1 (MQ=255)
                                                                                                                                                ttttcCCGCCCAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGa  <  1:870831/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGA  >  NC_000913/357734‑358024

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: ‑