Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8750AT44.1% 41.9 / 14.6 34N272K (AAT→AAA) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base A (13/6);  new base T (0/15);  total (13/21)
Fisher's exact test for biased strand distribution p-value = 3.87e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCA‑TGTTGTGTGGCAATCGCTGCTGGT‑AACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTA  >  NC_000913/357726‑358019
                                                                                                                                                       |                                                                                                                                                  
cTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTa                                                                                                                                                    <  1:584572/150‑1 (MQ=255)
   aaCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTAtttt                                                                                                                                                >  2:785640/1‑151 (MQ=255)
     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:385957/151‑1 (MQ=255)
                    aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                               <  1:555643/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                                                                                                             >  2:329554/1‑151 (MQ=255)
                                      aTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtagt                                                                                                             >  1:594350/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)
                                                                                                                                             gcTGCCGGTTTTTTTCCCGCCAAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCA‑TGTTGTGTGGCAATCGCTGCTGGT‑AACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGa           <  1:444663/145‑1 (MQ=255)
                                                                                                                                                 ccAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCA‑TGTTGTGTGGCAATCGCTGCTGGT‑AACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa  <  1:285797/151‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                  
CTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCA‑TGTTGTGTGGCAATCGCTGCTGGT‑AACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTA  >  NC_000913/357726‑358019

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