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

GTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTC  >  NC_000913/357749‑358028
                                                                                                                                                 |                                                                                                                                           
gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑ccgc                                                                                                                                    >  2:749669/1‑150 (MQ=255)
     gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGtt                                                                                                                                >  2:106096/1‑149 (MQ=255)
              atttcccccggCAGCCCTGTGCCGGCGGCTCCCGTTT‑TTC‑TTTCCCCCAACCAC‑CCGCGATTTTACCCACCCCCTTGTACAAAATTGCTTC‑CGGGGCGGTTGTCTGCCAGTTTTTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCagtag                                                                                                                     <  2:153865/141‑1 (MQ=255)
              attcccccccccAGACCTGTGTCGGCGGCTACGGTGA‑TTC‑TTTCCCCCAACCAC‑CCGCGCTTTTACCCACCGCAGTGTGAAAAAATGCTTT‑CGGGGCGGCTTTCTGCCAGTTTTTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCagtag                                                                                                                     <  2:746321/140‑1 (MQ=255)
              attcccAGCCGCAGCCTTGTGTCGGGGCCCCCCGTGT‑TTCTTTCCCGCCAA‑CAC‑CCGCACTTTTCCCCAACCCCATGTGAAAAAGTCCCTC‑CGCGGCGGCTGTCTCCCAGCTATTTTCCCCCCGCAAAAAAA‑CCGCCGTTCTCTCCagtag                                                                                                                     <  2:613064/147‑1 (MQ=255)
               attccccgccccgACCTGTGGCGGCGGCTTCCGGTA‑TTTATTTCCCCCAA‑AAACCCCCGCTTTTCCCCACCG‑CAAGTTGAAAAATGCCTTCCGGGGCGCCTTTCGCCCAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCagtag                                                                                                                     <  1:591658/139‑1 (MQ=255)
                  cccgCCGCAGCCTTGTTCCAGCGCACACCGGGT‑TTCATTCCCCCCAA‑CAC‑CCGCCCTTTACCCCAACGCCATGTGAAAAATGGCTTT‑CGGGGCGGCTGTCTGCCATTTTTTTTCCCGCCGAAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGc                                                                                                                 <  1:596116/148‑1 (MQ=255)
                           gACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgggggcagg                                                                                                        >  2:229730/1‑146 (MQ=255)
                            cccTGTGCCAGGGGCTACCGGGT‑TTCTTTCCCCCAAA‑CAA‑CCGCCCATTTACCCACCCCCGTGGGAAAAAAGGCCTT‑CGGGGGGGCTGTCGGCCAGCTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAagag                                                                                                       <  2:244679/150‑1 (MQ=255)
                             ccTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGccaaaaa                                                                                                      >  2:430780/1‑145 (MQ=255)
                                       cGGCTACCGTGA‑TTCATTTCCGCCAA‑CAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCaacgagcaaaaacttt                                                                                            >  1:675878/1‑137 (MQ=255)
                                           ttccGTGATTTCTTTCCCCCCAA‑CAC‑CCGCGCTTTTACCAAACGCTATGTGCAAAATTGCCTC‑CG‑GGCCGTTTTCCGCCGGGTTTTTTTCCCCCGGAAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGc                                                                                         <  1:141499/149‑1 (MQ=255)
                                             ccGTGA‑TTCATTTCCGCCAA‑CAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAaaaaaaaattttaaacc                                                                                       >  2:217178/1‑137 (MQ=25)
                                                 gA‑TTCATTTCCGCCAA‑CAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGGAAAAGATTTTAGccccccc                                                                                   >  2:35545/1‑146 (MQ=25)
                                                      taTTCCCGCCAA‑CAC‑CGGCGATTTAACCCACCCCCATGTGCAAAATTGCCTT‑CGGGGGGGCGGTCTCCAAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttg                                                                              <  1:747483/150‑1 (MQ=255)
                                                                                    ccacccccttgtcAAAAATTCCTTCCGGGGCGGTTGTCTGCCACTTTTTTTCCCCCCGCAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgt                                                   <  2:744250/139‑1 (MQ=255)
                                                                                                                     ttGTCTCCCATCT‑TTTTTCCCCCCCAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                  <  2:388495/150‑1 (MQ=255)
                                                                                                                     ctttctcccattttttttcCCGCCGAAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                  <  2:62492/137‑1 (MQ=255)
                                                                                                                     cTGTCTGCCAGCTTTTTTTCCGCCGCAAAAAAA‑CCGCCGTCCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                  <  2:363870/151‑1 (MQ=255)
                                                                                                                                ttgtttccccgcccCAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg       <  2:71594/143‑1 (MQ=255)
                                                                                                                                cTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAATTTTAACCCCCCCTTTTTTTTGGGGCTCGCTGCTGGGGACAAAAAGGGGAGGGAAGGGGGGGGGATCACATCCCGAACCccgccacttt       >  2:116652/1‑145 (MQ=255)
                                                                                                                                     ttcccgccgAAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc  <  2:651123/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                           
GTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTC  >  NC_000913/357749‑358028

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