Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA50.0% 15.8 / 24.3 34L267F (TTG→TTT) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (14/3);  new base A (0/17);  total (14/20)
Fisher's exact test for biased strand distribution p-value = 9.77e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCC  >  NC_000913/357742‑358031
                                                                                                                                                     |                                                                                                                                                
gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaaca                                                                                                                                                 <  2:161219/150‑1 (MQ=255)
   aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaaa                                                                                                                                              <  2:122449/150‑1 (MQ=255)
            gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtct                                                                                                                                    >  1:617630/1‑151 (MQ=255)
                          cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑g                                                                                                                     >  1:788951/1‑151 (MQ=255)
                             ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑ggg                                                                                                                   >  2:389249/1‑150 (MQ=255)
                              cgcgacccTGTGTCGCCGGCTACCTTGTTTCTTTCCCGCAACCACCCGCGCATTAACCAACCCCCAGGTCAAAAAATCCCTT‑CGGGGCGGCTGCCGCCCAGTTTTTTTCCCCCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGc                                                                                                                 <  2:407285/145‑1 (MQ=255)
                                  gACCTGTGTCAGCGGCTCCGGTGTTTATTTCCCGCCAACACCCGCCCTTTTCCCCACCCCAATGTGAAAAAATCCCTC‑CGCGGCGGCTGTCTCCCGGTTTTTTTCCCCCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑gggg                                                                                                                  <  1:196613/146‑1 (MQ=255)
                                    cttttgtCAGGGGCTACCGTGATTCTTTTCCCCCAACACCCGCGCTTTTACCCAACGCCTTGTGAAAAAATGCCTT‑CGCGGGGGCTTTCTGCCAGTTATTTTTCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAaga‑g                                                                                                           <  2:758470/146‑1 (MQ=255)
                                       gtgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGCAGA‑G‑aaa                                                                                                       >  2:407010/1‑150 (MQ=255)
                                         gtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑gggg                                                                                                                  >  1:845445/1‑139 (MQ=255)
                                           cAGCGGCCACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAACGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGGAACAAAACCGCCGTTCTCTCCAGTACTGGC‑GGTCGAAga‑g‑acgaaga                                                                                                   >  1:392705/1‑144 (MQ=21)
                                                     cGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAaga‑gggaaaagct                                                                                                  >  1:364790/1‑133 (MQ=255)
                                                     cGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAGC‑G‑aaaaaaa                                                                                                   >  2:612424/1‑136 (MQ=255)
                                                               ttCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGGG‑G‑AAAAAACTTTAAACCCgcc                                                                                       >  2:259955/1‑141 (MQ=25)
                                                                         aacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCCAGA‑G‑GAAAAAATTTTATCCAGTCAtgttgttttggggac                                                                       >  1:619834/1‑139 (MQ=25)
                                                                           ccGCGCTTTCCTCCACCGCCATGTAAAAAAATCCCTC‑CGCGGCGGCTGTCTGCCAGTTTTTTTTCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgct                                                                    <  1:99916/149‑1 (MQ=255)
                                                                                  tttccccACCGCCATGTGAAAAATTGCTTT‑CGGGGCGGTTGTCTGCGAGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTa                                                             <  2:766787/144‑1 (MQ=255)
                                                                                  tttACCCACCCCCATTTAAAAAAAGCCCTC‑CGGGGCGGCTTTCGCCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGt                                                              <  2:145730/148‑1 (MQ=255)
                                                                                       ccAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGA‑GAAAAAAGCTTTAAACCC‑TCCTGTTTTGTGGGGATCggctgtgggtaaa                                                          >  2:420469/1‑135 (MQ=25)
                                                                                                     aaaaaTCCTTC‑CGCGGCGGTTGTCTCCCACCTTTTTTTCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAg                                           <  2:439233/148‑1 (MQ=255)
                                                                                                     aaaaaTCCCTC‑CGCGGGGCCTGTCTCCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGt                                          <  2:3894/149‑1 (MQ=255)
                                                                                                       aaaaGCCTTCCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcg                                       <  2:243066/151‑1 (MQ=255)
                                                                                                       aaaTGCCTC‑CGCGGCGCCTTTCGCCCACTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcgg                                      <  2:294702/151‑1 (MQ=255)
                                                                                                        aaTGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGACGA‑G‑AAAAAGATTTTAACCCCCCATGTTGTGTGGCCATTGCCGCTGGTGAAAAAGTGggaagggaagggggg                                      >  2:406343/1‑145 (MQ=255)
                                                                                                        aaTGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggc                                     <  2:618075/151‑1 (MQ=255)
                                                                                                                 cgcgGCGGCTGTCTCCCAGTTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGc                             <  1:215249/151‑1 (MQ=255)
                                                                                                                  gcggcgGCTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGc                             <  2:306914/150‑1 (MQ=255)
                                                                                                                   cggcggCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAGA‑G‑AAAAAAATTTTAACCCCTCCTTTTGTGTGGGGATCCCTGCTGGTGAAAAAGTGGGAAGGGGAgggggggggagc                                >  2:215273/1‑135 (MQ=255)
                                                                                                                    gggggcttcctcccAGTTTTTTTCCCGCCCAAAAAAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGc                             <  1:673827/136‑1 (MQ=255)
                                                                                                                    ggcggcTTTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTc                           <  2:458246/150‑1 (MQ=255)
                                                                                                                       ggcTGTCTCCCAGTTTTTTCCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAg                        <  1:660753/150‑1 (MQ=255)
                                                                                                                               cccgttaaTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGaa                      <  2:362261/137‑1 (MQ=255)
                                                                                                                                ccAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                     <  2:756897/144‑1 (MQ=255)
                                                                                                                                     tATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAGG‑AAAAAATTTTAAGCCC‑TCCTGTTGTTTGGGGAACGCTGCTGGGGAAAAAGGGGGAAGGGGAGGGGGGGGGAGCAGCCCCCCAACCCCGCCAAtt           >  2:348568/1‑149 (MQ=255)
                                                                                                                                        tttcCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGt       <  1:844925/150‑1 (MQ=255)
                                                                                                                                            cccccacaaaaaAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGcc  <  2:701197/145‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCC  >  NC_000913/357742‑358031

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