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

CTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCA‑TGTTGTGTGGCAATCGCTGCTGGT‑AACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTA  >  NC_000913/357731‑358019
                                                                                                                                                       |                                                                                                                                             
cTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTAttttt                                                                                                                                               <  1:675467/150‑1 (MQ=255)
  cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTccg                                                                                                                                            <  2:597493/151‑1 (MQ=255)
        ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGca                                                                                                                                       >  2:733787/1‑150 (MQ=255)
           gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCACACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACATCCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGcaaca                                                                                                                                    <  1:146131/150‑1 (MQ=255)
              aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGcaacaa                                                                                                                                   <  1:747260/148‑1 (MQ=255)
              aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                                >  1:314963/1‑151 (MQ=255)
               aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                               <  1:555643/151‑1 (MQ=255)
               aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                               <  1:385957/151‑1 (MQ=255)
                 cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCTGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgc                                                                                                                             >  1:811103/1‑151 (MQ=255)
                          tCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctctc                                                                                                                    >  1:759031/1‑151 (MQ=255)
                                aaTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                >  2:485783/1‑149 (MQ=255)
                                 aTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtagt                                                                                                             >  1:594350/1‑151 (MQ=255)
                                 aTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtagt                                                                                                             >  2:329554/1‑151 (MQ=255)
                                   cccccccccaccccTGTGTCGGCGCTCACCTTAATTCTTTCCCCCA‑AACACCCGCGCATTTACCAACCCCCATGTGAAAAAATCCCTCCGGGGCGG‑CTGTCTGCAAGTTTTTTTCCCGCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTg                                                                                                            <  1:735079/138‑1 (MQ=255)
                                        ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                       >  1:214886/1‑150 (MQ=255)
                                         cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggg                                                                                                     >  2:700539/1‑150 (MQ=255)
                                            agcccTGTGTCGGCGCCTACCGTTTTTAATTCCCGCA‑AACACCCGCGCTTTTACCCAACCCCAGGTGAAAAAATGCCTCCGCGGGGC‑CTTTCTCCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCaa                                                                                                   <  2:759126/147‑1 (MQ=255)
                                                     tCAGCGGCTCCGGTG‑TTTATTTCCCCCAAACACCCGCGCTTTTACCAAACGCAGTGTCAAAAAATCCTTTCGCGGGGG‑CTGTCTGCAATGTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaaa                                                                                             <  2:382761/147‑1 (MQ=255)
                                                          gccccccGTGTTTCATCCCCCCC‑AACACCCGCGCTTTTCACAACCGCATTGTGAAAAAATGCCTGCGCGGCGG‑CTGTCGGCAAGTTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAgaga                                                                                               <  1:774401/135‑1 (MQ=255)
                                                               ccGTGATTCATTCCCCCC‑ACCAACCGCGCTTTTACCCACCGCCATGTCAAAAAATCCCTGCGGGGCGG‑CTGTCGCCCAGCTTTTTTCCCGCCCCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAg                                                                                   <  2:83495/147‑1 (MQ=255)
                                                                     ttCATTCCCCCC‑AACACCCGCCATTTTACCAAACGCCATGTACAAAAAGCCCTCCGCGGCGG‑CTTCCTGCCGGTTTTTTTCCCCCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCA‑t                                                                          <  1:225395/149‑1 (MQ=255)
                                                                               cc‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCTGAGAGAAAAAATTTAGTCCGTCA‑TGTtgtgt                                                                   >  1:203953/1‑146 (MQ=25)
                                                                                  accacccgcatttttACCCAGCGCCAGGTGAAAAACTGCTTGCGGGGGGT‑CTGTCGCCCAGTTTTTTTCCCCCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAGAAATAGCTTTAAGCCCGTCA‑TGTTGTGTGGCa                                                               <  1:65712/136‑1 (MQ=255)
                                                                                       ccGCGCTTTTCCCCACCCCCATGCGAAAAACTCCTTCCGGGGCGG‑CTGTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCA‑TGTTGTGTGGCAATc                                                            <  1:228436/146‑1 (MQ=255)
                                                                                              tttACCCAACGCCAGGACAAAAATTTCCTCGGGGGCGG‑TTTTCTGCCGGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCA‑TGTTGTGTGGCAATCgctgct                                                      <  1:785538/145‑1 (MQ=255)
                                                                                                       ccccATGTAAAAAAATCCCTCCGCGGCGG‑CTGTCTCCCAGTTTTTTTCCCGCCACAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCA‑TGTTGTGTGGCAATCGCTGCTGGT‑AACAATGTGGaa                                      <  1:853991/149‑1 (MQ=255)
                                                                                                           aTGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAACTTTTAACCCCTCA‑TTTTTTGTTGGGATAGCTGCCGGTGAACaaatgggaaggg                                   >  2:794437/1‑141 (MQ=255)
                                                                                                                aaaaaaTCCTTCCGCGGGGG‑CTGTCGGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCA‑TGTTGTGTGGCAATCGCTGCTGGT‑AACAATGTGGAAAGGGAAGTg                             <  2:116116/150‑1 (MQ=255)
                                                                                                                 aaaaaTGCTTCGGCGGCGG‑TTGTCTGCCAGCTTTTTTTCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCA‑TGTTGTGTGGCAATCGCTGCTGGT‑AACAATGTGGAAAGGGaa                                <  2:159876/147‑1 (MQ=255)
                                                                                                                  aaaactgctgcgcgGGGGCCTGTCTCCCCGTTTTTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCA‑TGTTGTGTGGCAATCGCTGCTGGT‑AACAATGTGGAAAGGGa                                 <  1:729017/136‑1 (MQ=255)
                                                                                                                          ttCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAATTTTAACCCCTCACTGTT‑TGTTGGCAGCACTGCTGGT‑TACAATGGGGGAAGGGGAgggaggggg                       >  2:410447/1‑139 (MQ=255)
                                                                                                                                        tcggaaggtttttttcCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCA‑TGTTGTGTGGCAATCGCTGCTGGT‑AACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcg       <  1:385731/139‑1 (MQ=255)
                                                                                                                                        gcTGCCGGTTTTTTTCCCGCCAAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCA‑TGTTGTGTGGCAATCGCTGCTGGT‑AACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGa           <  1:444663/145‑1 (MQ=255)
                                                                                                                                            ccAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCA‑TGTTGTGTGGCAATCGCTGCTGGT‑AACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa  <  1:285797/151‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                             
CTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCA‑TGTTGTGTGGCAATCGCTGCTGGT‑AACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTA  >  NC_000913/357731‑358019

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