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

TAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACC  >  NC_000913/357735‑358013
                                                                                                                                             |                                                                                                                                             
tAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCCACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑ccgcc                                                                                                                                   >  2:240065/1‑151 (MQ=255)
  ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑Ccgccgc                                                                                                                                 <  1:706442/151‑1 (MQ=255)
    ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGca                                                                                                                                >  1:665604/1‑150 (MQ=255)
           aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCTTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAAcc                                                                                                                        <  2:728787/151‑1 (MQ=255)
            gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAAcc                                                                                                                        <  2:786506/150‑1 (MQ=255)
                     aTCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTtctct                                                                                                              >  1:17788/1‑151 (MQ=255)
                     aTCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTtctct                                                                                                              >  1:476869/1‑151 (MQ=255)
                       cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTtctct                                                                                                              >  2:734458/1‑149 (MQ=255)
                            attcccAGCCGCAGCCCTGTGTCAGCGGCTCCCGTGATTCTTTTCCGCCAACACCCGCGCATTTACCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTTTCTGCCAGTTTTTTT‑C‑CCGCCGCAAAAAAACCGCCGTTCTCTCCagtag                                                                                                       <  1:181352/147‑1 (MQ=21)
                                   gCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                >  1:430604/1‑151 (MQ=255)
                                    ccccAGACCTTTTTCCGCGCCTACCGTGATTCATTTCCCCCAACCCCCGCGCATTTTTCCAACGCCTGGTGCAAAAATCCCTTCCC‑GGGGTCTGTCGTCCAGCTTTTTT‑C‑CCGCCGCAACAAAACCGCCGTtc                                                                                                                 <  2:815484/130‑1 (MQ=255)
                                       cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggtg                                                                                               >  1:678193/1‑146 (MQ=255)
                                       cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                >  2:984825/1‑147 (MQ=255)
                                        aGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggg                                                                                              >  2:832828/1‑147 (MQ=255)
                                          cccTGTGCCAGCGGCTCCGTTTTTTCTTTCCCGCAAACACCCGCCGATTTACCCAACCCCATGCACAAAATGCCCTCCGG‑GGCGGCTGTCTGCCAGTTTTTTT‑C‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga                                                                                          <  2:508847/149‑1 (MQ=255)
                                          aCCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggga                                                                                            >  1:309704/1‑145 (MQ=255)
                                           cttgtgTCGGCGGCTCCCGTATTTCATTCCCGCCAACACCCGCGCTTTTATCCACCCCCATGTGAAAAAATCCCTTCGC‑GGCGGCTGTCTGCCAGTTTTTTT‑C‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAag                                                                                           <  1:585679/146‑1 (MQ=255)
                                           ccTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                >  2:834207/1‑143 (MQ=255)
                                             tgtgTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggagggaa                                                                                       >  2:959260/1‑142 (MQ=255)
                                             tgtgCCGGCGGCTACCGTATTTCTTTCCCCCCAAAACCCGCGCTTTTACCCAACCCATTGTGAAAAAAGGCTTCCGG‑GGCGGCTGTCTGCCAGTTTTTTT‑C‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga                                                                                          <  1:263130/147‑1 (MQ=255)
                                                            ctggATCCATTCCCGCCAACCCCCGCTCATTTCCCCAACCCCGTGTAAAAAAATGCTTTCGG‑GGCGGCTTTCTGCCAGTTTTTTT‑C‑CCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTa                                                                             <  2:890070/141‑1 (MQ=255)
                                                                     ttcccGCAAAAACCCGCGCTTTTACCCACCGCAATGTGAAAAAATGCCTCCGC‑GGGGGTTGTCTGCAAGTTTTTTT‑T‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCAtgttg                                                               <  1:34897/146‑1 (MQ=255)
                                                                       cccccaaacaccCGCGCTTTTACCCAACCCAATGTGAAAAAATGCTTCCGC‑GGCGGCTGTCGGCCAGTTTTTTT‑C‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTtgtg                                                             <  1:841971/143‑1 (MQ=255)
                                                                        ccGCAAACACCCGCGCATTTACCCAGCCCCTTGTAAAAAACTGCCTTCGG‑GGGGGCGTTCCGCCAGCTTTTTT‑C‑CCGCCACAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCAtgttgt                                                              <  1:529052/147‑1 (MQ=255)
                                                                             aacaccCGCGCATTTATCCAACCCCAGGTGCAAAAATGCCTTCGG‑GGGGGCTGTCTGCCAGCTTTTTT‑C‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGc                                                         <  1:218647/147‑1 (MQ=255)
                                                                                  ccGCGCTTTTCTCC‑ACCCCATGGAAAAAAAGTGCTTCGCGGGCGGCTGTC‑CCCAACTTTTTT‑TCCCGCAAAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgct                                                  <  2:560931/149‑1 (MQ=255)
                                                                                     cgcATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑GGAAAAAATTTTAACCCCTCCTGTTGTGTGGGAGTCGCTGCTgg                                             >  1:69191/1‑151 (MQ=25)
                                                                                               cAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAAA‑GAAAAAAATTTTAACCCCCCCTGTTTTGTGGGGAACGCTGCTGGGGACaa                                       >  1:156832/1‑147 (MQ=25)
                                                                                                                           ggcggtTGCCTGCCATATTTTTT‑G‑CCGCAAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTc         <  2:971473/150‑1 (MQ=255)
                                                                                                                           ggcggcTGTCTGCAGGTTTTTTT‑C‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCa             <  1:834432/146‑1 (MQ=255)
                                                                                                                           ggcggcTGTC‑GCCAATTTTTTTCT‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGc           <  1:485496/148‑1 (MQ=255)
                                                                                                                           ggcggcGGTCTCCCAGTTTTTTT‑C‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCa             <  1:593768/146‑1 (MQ=255)
                                                                                                                            gcggcTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGGAAAAAATTTTAAACCCCCCCAGTT‑TGTGGGAAACGCTGCTGGGGAAAAAATGGGAAGGGGAGTGGGGGGGAACAACTcc        >  1:778456/1‑150 (MQ=255)
                                                                                                                              ggcTGTCTGCCAGTTATTTT‑C‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa       <  1:835601/149‑1 (MQ=255)
                                                                                                                                   tcggCCAGTTTTTTT‑C‑CCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa       <  2:464625/141‑1 (MQ=255)
                                                                                                                                    cTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑GAAAAAAATTTTAAACCCCCCTGTTGTGGGGGAAACGCTGCTGGGTAAAAAATGGGAAGGGGAGGGGGGGGGAACAGCTCCAGaaac  >  1:641567/1‑146 (MQ=255)
                                                                                                                                             |                                                                                                                                             
TAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACC  >  NC_000913/357735‑358013

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