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→AAAcynRtranscriptional 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

TAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACC  >  NC_000913/357735‑358013
                                                                                                                                             |                                                                                                                                             
tAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCCACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                                                                                   >  2:240065/1‑151 (MQ=255)
  ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCcgccgc                                                                                                                                 <  1:706442/151‑1 (MQ=255)
    ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGca                                                                                                                                >  1:665604/1‑150 (MQ=255)
           aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCTTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                        <  2:728787/151‑1 (MQ=255)
            gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                        <  2:786506/150‑1 (MQ=255)
                     aTCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctct                                                                                                              >  1:17788/1‑151 (MQ=255)
                     aTCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctct                                                                                                              >  1:476869/1‑151 (MQ=255)
                       cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctct                                                                                                              >  2:734458/1‑149 (MQ=255)
                            attcccAGCCGCAGCCCTGTGTCAGCGGCTCCCGTGATTCTTTTCCGCCAACACCCGCGCATTTACCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTTTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCagtag                                                                                                       <  1:181352/147‑1 (MQ=21)
                                   gCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                >  1:430604/1‑151 (MQ=255)
                                    ccccAGACCTTTTTCCGCGCCTACCGTGATTCATTTCCCCCAACCCCCGCGCATTTTTCCAACGCCTGGTGCAAAAATCCCTTCCCGGGGTCTGTCGTCCAGCTTTTTTCCCGCCGCAACAAAACCGCCGTtc                                                                                                                 <  2:815484/130‑1 (MQ=255)
                                       cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggtg                                                                                               >  1:678193/1‑146 (MQ=255)
                                       cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                >  2:984825/1‑147 (MQ=255)
                                        aGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggg                                                                                              >  2:832828/1‑147 (MQ=255)
                                          cccTGTGCCAGCGGCTCCGTTTTTTCTTTCCCGCAAACACCCGCCGATTTACCCAACCCCATGCACAAAATGCCCTCCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga                                                                                          <  2:508847/149‑1 (MQ=255)
                                          aCCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggga                                                                                            >  1:309704/1‑145 (MQ=255)
                                           cttgtgTCGGCGGCTCCCGTATTTCATTCCCGCCAACACCCGCGCTTTTATCCACCCCCATGTGAAAAAATCCCTTCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAag                                                                                           <  1:585679/146‑1 (MQ=255)
                                           ccTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                >  2:834207/1‑143 (MQ=255)
                                             tgtgTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggagggaa                                                                                       >  2:959260/1‑142 (MQ=255)
                                             tgtgCCGGCGGCTACCGTATTTCTTTCCCCCCAAAACCCGCGCTTTTACCCAACCCATTGTGAAAAAAGGCTTCCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga                                                                                          <  1:263130/147‑1 (MQ=255)
                                                            ctggATCCATTCCCGCCAACCCCCGCTCATTTCCCCAACCCCGTGTAAAAAAATGCTTTCGGGGCGGCTTTCTGCCAGTTTTTTTCCCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTa                                                                             <  2:890070/141‑1 (MQ=255)
                                                                     ttcccGCAAAAACCCGCGCTTTTACCCACCGCAATGTGAAAAAATGCCTCCGCGGGGGTTGTCTGCAAGTTTTTTTTCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttg                                                               <  1:34897/146‑1 (MQ=255)
                                                                       cccccaaacaccCGCGCTTTTACCCAACCCAATGTGAAAAAATGCTTCCGCGGCGGCTGTCGGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtg                                                             <  1:841971/143‑1 (MQ=255)
                                                                        ccGCAAACACCCGCGCATTTACCCAGCCCCTTGTAAAAAACTGCCTTCGGGGGGGCGTTCCGCCAGCTTTTTTCCCGCCACAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttgt                                                              <  1:529052/147‑1 (MQ=255)
                                                                             aacaccCGCGCATTTATCCAACCCCAGGTGCAAAAATGCCTTCGGGGGGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGc                                                         <  1:218647/147‑1 (MQ=255)
                                                                                  ccGCGCTTTTCTCCACCCCATGGAAAAAAAGTGCTTCGCGGGCGGCTGTCCCCAACTTTTTTTCCCGCAAAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgct                                                  <  2:560931/149‑1 (MQ=255)
                                                                                     cgcATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGGAAAAAATTTTAACCCCTCCTGTTGTGTGGGAGTCGCTGCTgg                                             >  1:69191/1‑151 (MQ=25)
                                                                                               cAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAAAGAAAAAAATTTTAACCCCCCCTGTTTTGTGGGGAACGCTGCTGGGGACaa                                       >  1:156832/1‑147 (MQ=25)
                                                                                                                           ggcggtTGCCTGCCATATTTTTTGCCGCAAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTc         <  2:971473/150‑1 (MQ=255)
                                                                                                                           ggcggcTGTCTGCAGGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCa             <  1:834432/146‑1 (MQ=255)
                                                                                                                           ggcggcTGTCGCCAATTTTTTTCTCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGc           <  1:485496/148‑1 (MQ=255)
                                                                                                                           ggcggcGGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCa             <  1:593768/146‑1 (MQ=255)
                                                                                                                            gcggcTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGGAAAAAATTTTAAACCCCCCCAGTTTGTGGGAAACGCTGCTGGGGAAAAAATGGGAAGGGGAGTGGGGGGGAACAACTcc        >  1:778456/1‑150 (MQ=255)
                                                                                                                              ggcTGTCTGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa       <  1:835601/149‑1 (MQ=255)
                                                                                                                                   tcggCCAGTTTTTTTCCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa       <  2:464625/141‑1 (MQ=255)
                                                                                                                                    cTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAATTTTAAACCCCCCTGTTGTGGGGGAAACGCTGCTGGGTAAAAAATGGGAAGGGGAGGGGGGGGGAACAGCTCCAGaaac  >  1:641567/1‑146 (MQ=255)
                                                                                                                                             |                                                                                                                                             
TAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACC  >  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: