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

CGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGT  >  NC_000913/357733‑358026
                                                                                                                                                     |                                                                                                                                                    
cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccg                                                                                                                                                   >  2:86869/1‑151 (MQ=255)
   aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgccg                                                                                                                                                >  2:942897/1‑151 (MQ=255)
       cTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                                            >  1:133534/1‑151 (MQ=255)
               cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgc                                                                                                                                    >  1:375463/1‑151 (MQ=255)
                  aGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGt                                                                                                                                 >  1:535071/1‑151 (MQ=255)
                    cgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtc                                                                                                                               >  1:349805/1‑151 (MQ=255)
                     gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtct                                                                                                                              >  2:541577/1‑151 (MQ=255)
                              attcccAGCCGCAGCCTTGTGTCAGCGGCTACCGTGATTCTTTCCCCCAAACACCCGCGCTTTTATCCACCGCCTTGTGAAAAAATCCTTTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCagta                                                                                                                      <  2:763287/146‑1 (MQ=18)
                                tcccAGCCGAAGCCCTGTGTCAGCGCCTACCGTAATTCATTCCCCCCAACACCCGCGCTTTTATCCAACCCCTTGTGAAAAAATGCCTTCGGGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCagtagt                                                                                                                    <  2:163983/148‑1 (MQ=255)
                                 tCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTgg                                                                                                                  >  1:45126/1‑151 (MQ=255)
                                  ccAGCCGCAGACTTGTGTCAGCGGCTCCCGTGATTCATTCCCCCCAACACCCGCGCTTTTACCCAACGCCATGTGAAAAAATGCCTCCGGGGCGGCTGTCTGCCAGCTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCagtagt                                                                                                                    <  2:296652/148‑1 (MQ=255)
                                      ccGGAGCCCTGTGTCTGGGGCTCCCGTGATTCATTCCCGCCAACACCCGCGCATTTACCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                               <  2:268957/149‑1 (MQ=18)
                                      ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                             >  2:421789/1‑151 (MQ=255)
                                        gCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                               >  1:715467/1‑147 (MQ=255)
                                           gccttgtgTTAGGGGCAACCGTATTCCATTCCGCCAACTACCGCGGCTTTTACCCAACCCAATGTGAAAAAATCCTTTCGCGGCGGCTGTCGCCCAGCTATTTCACCCCCCAAACAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCa                                                                                                           <  2:402837/144‑1 (MQ=255)
                                              cgtggccggCGACTCCCGTGATTCTTTCCCGCAAACACCCGCCCATTTACCCACCGCATTGTGCAAAAAGCCCTCGGCGGGGCTTGTCCGCCGGTTTTTTTCCCCCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga                                                                                                        <  1:42925/140‑1 (MQ=255)
                                                             ccttttttCGTTCCCGCCACCAACCGCGCATTTACCCAGCCCTATGTGCAAAAATCCTTCGGCGGGGTCTTCCCGCTATCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAg                                                                                        <  2:188540/141‑1 (MQ=255)
                                                             ccGTTATTCATTCCCGCCAACACCCGCGATTTTATCCAACGCCATGTGAAAAAATGCCTTCGGGGCGGCTGTCTCCCAGCTATTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCg                                                                                    <  2:316783/151‑1 (MQ=255)
                                                              cGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGGGaaaaaaactttaaaccc                                                                                    >  1:526064/1‑136 (MQ=25)
                                                                   ttCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAAAAAAACCATTaacacc                                                                                     >  2:433699/1‑140 (MQ=25)
                                                                       ttcccGCAACCACCCGCCCATTTACCCACCCCCATGGCAAAAAATGCCTCCGCGGCGGTTGTCTGCCAGTTTTTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttgt                                                                           <  1:391667/147‑1 (MQ=255)
                                                                         tCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAAAAATTTTAACCCGCCAAGTtttgtg                                                                        >  2:650954/1‑151 (MQ=25)
                                                                             ccAACACCCGCGCTTTTACCCACCCCCATGTGAAAAAATCCCTTCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCaa                                                                    <  1:65122/151‑1 (MQ=255)
                                                                               accacccGCGCTTTTCCCCACCGCTATGGGAAAAAATCCCTTCGGGGCGGCTTCCTGCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc                                                                  <  1:606031/146‑1 (MQ=255)
                                                                                  acccccgcTTTTACCCAACGCCATTTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGTTATTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgct                                                               <  2:919631/146‑1 (MQ=255)
                                                                                  aacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCCGAGGAAAAGCTTTAAAAGCGGCCTGTTGTGTgggga                                                                    >  2:706203/1‑143 (MQ=25)
                                                                                   acCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGGAAAAAGCTTAAACCCCgcacatttt                                                                            >  2:322970/1‑130 (MQ=25)
                                                                                     cgcgCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAAAGAAAAAAATTTTAACCCCCCACGTTTTGTGGGAATc                                                                  >  2:605593/1‑145 (MQ=25)
                                                                                                     cGCCATTTGCAAAAATCCTTTCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGa                                              <  2:28460/150‑1 (MQ=255)
                                                                                                     cGCCATGTGAAAAAAGGCTTTCGCGGCGGCTGCCTGCCGACTTTTTCCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGa                                              <  2:873466/150‑1 (MQ=255)
                                                                                                         aTGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAAATTTTAAACCCGCATGTTTTGGGGGAAACGCTGCCGGGGACaaaatgggaaagg                                         >  2:511680/1‑140 (MQ=255)
                                                                                                          tgtgAAAAAATGCCTGCGGGGCGGCTGTCTGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAgg                                          <  2:375506/149‑1 (MQ=255)
                                                                                                           gtgAAAAAATCCCTCGGGGGCGGCTGTCTGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                        <  1:942787/150‑1 (MQ=255)
                                                                                                            tgaaaaatTGCTTTCGCGGGGGCTGTCTGCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                        <  2:749586/146‑1 (MQ=255)
                                                                                                            tgCAAAAATCCTTTCGCGGGGGCTGTCTGCCAGTTATTTTTCCCCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGaa                                       <  1:779015/150‑1 (MQ=255)
                                                                                                               aaaaaTGCCTCCGCGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAg                                      <  2:535129/148‑1 (MQ=255)
                                                                                                                            gcggcgGCTGTCTCCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAg                         <  1:749937/149‑1 (MQ=255)
                                                                                                                                        tGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAACTTTTAAACCCCCATGTTTTGTGGGGATAGCTGCTGGTGAAAAAGTGGGAAGGGGAGGGGGGCGGATCAGCCCCAGAAACccgc           >  2:15777/1‑151 (MQ=255)
                                                                                                                                         cccAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCt          <  1:63488/150‑1 (MQ=255)
                                                                                                                                          ccattttttttcCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                <  1:312226/137‑1 (MQ=255)
                                                                                                                                                 ttttcCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg   <  2:334764/150‑1 (MQ=255)
                                                                                                                                                    ttCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGTGCCAGGGAAAAAACATTAAACCCCCCCTGTTGTGTGGCGATAGCCGCTGGTGACAATTTGGGAAGGGGAGTGCGGGGGGTCAGCCCCACAACCCCGCCAAATGaat  >  1:502373/1‑146 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGT  >  NC_000913/357733‑358026

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: