Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA56.7% 12.4 / 38.6 30L267F (TTG→TTT) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (8/5);  new base A (0/17);  total (8/22)
Fisher's exact test for biased strand distribution p-value = 2.20e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.64e-04

GCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑ATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCT  >  NC_000913/357747‑358032
                                                                                                                                                  |                                                                                                                                               
gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAccg                                                                                                                                          >  1:362792/1‑151 (MQ=255)
  gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAccgc                                                                                                                                         <  1:433418/150‑1 (MQ=255)
      cccATCAGCCAATCCCAGCCCCAGACCTGTTTCAGCGGCTACCGTGATTCATTCCCGCC‑AA‑CACCCGCGCATTTATCCAACGCCATGTGCAAAATTGCCTCGGCGGCGGCTGTCTGCCAGTTATTTTT‑CCGCCGCAACAAAAccgccg                                                                                                                                       <  2:8652/146‑1 (MQ=18)
           cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                               >  1:689215/1‑151 (MQ=255)
                      aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                    >  2:455066/1‑151 (MQ=255)
                       gCCCCACACTTGTGTCTGCGCCTCCCGTAATT‑ATTTCCCCCAAA‑CAGCCGCCCTTTTTCCCACCTCAATGTAAAAAAATCCTTGCGGGGGGGTTGTCGCCCGGTTTTTTCC‑CCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                                                      <  1:131827/148‑1 (MQ=255)
                         cGCAGCCTTGTGTCAGGGCCTACCGTGTTTCATTCCCCCC‑AA‑CACCCGCGCATTTACCCACCGCCTTGTCAAAAAATGCTTGCGGGGGGGCTGTCTGCCGGTTTTTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGc                                                                                                                 <  1:73522/151‑1 (MQ=255)
                           cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                  >  2:923837/1‑148 (MQ=255)
                               cgtgtgTCAGCGGCTCCGGTGATTCTTTCCCCCC‑AA‑CACCCGCGCTTTTACCCACCGCCATGTGAAAAAAGGCTTCCGGGGCGGCTTTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAag                                                                                                              <  1:558397/146‑1 (MQ=255)
                                    gtCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                  >  2:468941/1‑139 (MQ=255)
                                    gtCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAaggg                                                                                                            >  1:270815/1‑143 (MQ=255)
                                                cctttattctttCCCCC‑AACCACCCGCCCATTTACCCACCCCCTGGGGAAAAAATCCCTCGGCGGCGGCTGTCTGCCAGTTTTTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCaa                                                                                                               <  2:837379/122‑1 (MQ=255)
                                                     ttCATTCCCCCC‑AA‑CACCCGCGCTTTTACCCACCGCATTGTGAAAAAATGCTTCCGGGGCGCCTGTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑a                                                                                       <  1:922392/148‑1 (MQ=255)
                                                       ctttcccccc‑AA‑CACCCGCGTTTTTACCCACCCCCATGGGAAAAATTCCCTTCGCGGCGGTTGTCT‑CCAGTTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑At                                                                                      <  1:267735/139‑1 (MQ=255)
                                                           cccccc‑AA‑CACCCGCCATTTTACCCACCGCCTTGTCAAAAAATGCCTCGGGGGCGGCTGTCGGCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑ATGTtgtgtg                                                                              <  1:1278190/147‑1 (MQ=255)
                                                                  aa‑caccCGCCCATTTACCCACCCCCATGTAAAAAAAGGCTTCCGGGGCGGCTGTCGGCCAGCTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑ATGTTGTGTGGCAATc                                                                        <  1:896938/151‑1 (MQ=255)
                                                                          gcgcTTTTACCAAACCCTATGTAAAAAAATCCTTGCGCGGGGGCTGTCTGCCAGTTTTTTTT‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑ATGTTGTGTGGCAATCgctg                                                                    <  1:1196114/148‑1 (MQ=255)
                                                                                    caaacccgtgggcaaaaaaTGCCTCCGGGGCGGCTGTCTGCAAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑ATGTTGTGTGGCAATCGCTGCTGGTAACAAt                                                         <  2:1323813/135‑1 (MQ=255)
                                                                                      accgcttgggacAAAAATGCCTCGGCGGCGGCTGTCTGCCAGTTATTTTC‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑ATGTTGTGTGGCAATCGCTGCTGGTAACAATGTgg                                                     <  2:187722/140‑1 (MQ=255)
                                                                                              gtgAAAAAATCCCTTCGGGGCGGCTGTCTGCCAGTTTTTTTC‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑ATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg                                               <  2:979972/149‑1 (MQ=255)
                                                                                              ggcaaaaaaTGCCTCCGCGGGGGTTGCCGGCGAGTTTTTTCC‑CCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑ATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGaa                                             <  2:310983/147‑1 (MQ=255)
                                                                                                  aaaatTGCCTCCGCGCGGGCTTTCTGCCAGTTATTTTC‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑ATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGaa                                             <  1:79548/147‑1 (MQ=255)
                                                                                                   aaaaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAAGTTTAAAACCCTCAATGTT‑TTTGGGGATAGCTGCTGGTTACAACGTgg                                                     >  2:783339/1‑138 (MQ=255)
                                                                                                    aattGCCTCCGGGGCGCCTTTCCGCCGTTTTTTTTC‑CCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑ATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcgg                                       <  2:24848/148‑1 (MQ=255)
                                                                                                                ggcggtTGTCTCCCGCTTTTTTTT‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑ATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt                                  <  2:1329885/144‑1 (MQ=255)
                                                                                                                     cTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAAATTTTAAACCCCC‑ATGTTTTGTGGGGAACGCTGCCGGGGACAAAATGGGAAGGGGAGGGCGGGGGAACAACTCCCGaa                       >  1:1253634/1‑150 (MQ=255)
                                                                                                                            ccAGCTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑ATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgcg                  <  1:506326/148‑1 (MQ=255)
                                                                                                                                 ttttttc‑CCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑ATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAtt            <  1:1153583/147‑1 (MQ=255)
                                                                                                                                      ttcccgccgCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑ATGTTGTGTGGCAATCGCTGCTGGTAACAACGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt       <  1:810575/147‑1 (MQ=255)
                                                                                                                                          ccccgaAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑ATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCt  <  1:19340/149‑1 (MQ=255)
                                                                                                                                                aacaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAATTTTAAACCCCC‑CTGTTTTTTGGGGATCGCCGCTGGGAAAAAAATGGGAAGGGGAGGGGGGGGGAGCAGCTCCCGAAACGCGCCAAATTAGttt      >  1:1158706/1‑140 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
GCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑ATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCT  >  NC_000913/357747‑358032

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