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

ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑ATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATT  >  NC_000913/357730‑358022
                                                                                                                                                        |                                                                                                                                                
aCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt                                                                                                                                                  >  1:347186/1‑151 (MQ=255)
aCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt                                                                                                                                                  >  2:52385/1‑151 (MQ=255)
  tCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt                                                                                                                                                  >  2:519106/1‑149 (MQ=255)
     tAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑ccgc                                                                                                                                             <  2:109404/150‑1 (MQ=255)
      aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑ccgcc                                                                                                                                            <  2:139544/150‑1 (MQ=255)
       ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑Ccgccgc                                                                                                                                          >  2:116498/1‑151 (MQ=255)
                 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)
                                                                                                                                                        |                                                                                                                                                
ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑ATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATT  >  NC_000913/357730‑358022

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