Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8800TC44.8% 27.0 / 10.3 29K271E (AAA→GAA) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (12/4);  new base C (0/13);  total (12/17)
Fisher's exact test for biased strand distribution p-value = 5.94e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTA  >  NC_000913/357733‑358019
                                                                                                                                                    |                                                                                                                                             
cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccg                                                                                                                                            >  2:229829/1‑151 (MQ=255)
  tAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTTTGCCAGCTATTTTTccgcc                                                                                                                                          >  2:359110/1‑151 (MQ=255)
           tAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                                                                                 <  2:639455/150‑1 (MQ=255)
            aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                               >  1:352699/1‑151 (MQ=255)
               cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgc                                                                                                                            <  2:28351/151‑1 (MQ=255)
               cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGCGATTCATTTCCGCAAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgc                                                                                                                            >  2:201639/1‑151 (MQ=255)
                     gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtc                                                                                                                       >  2:73675/1‑150 (MQ=255)
                        tCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGTCAGCTATTTTTCCGCCGCAACAAAACCGCCGAtctctc                                                                                                                   >  2:409593/1‑151 (MQ=255)
                            gCAATTCAAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                               >  2:407856/1‑151 (MQ=255)
                             cAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                               >  2:18290/1‑150 (MQ=255)
                                                  ggcacggtctcccGTGATTCATCCCCCCCAACACCCGCCTATTTTCCCACCGCCAGGTGCAAAACGGCTTCCGGGGCGCCTGTCGGCAAGTTTTTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAg                                                                                        <  2:585305/140‑1 (MQ=255)
                                                    cAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAAAGgaaaaaat                                                                                      >  1:626449/1‑146 (MQ=255)
                                                     ggCGGCAACCGTAATTCATTCCCCCCACCACCCCCTCTTTCATCCCCCGCCGTGTTAAAAAATGCGTCCGCGGCGCCTGTCTCCCATTTTTTTTCCCCCCAAACAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCtt                                                                                     <  2:458505/150‑1 (MQ=255)
                                                      gCGGCTCCGGTGTTTCTTTCCCGCCAACACCCGCGCATTTACCCACCGCCATGTGAAAAAAGGCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCt                                                                                      <  1:621502/149‑1 (MQ=255)
                                                        ggCTCCCGTGTTTCTTTCCCCCCAACACCCGCGCATTTACCCACCCCCATGTAAAAAAATCCCTTCGGGGCGGCTGTCGGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaa                                                                                            <  1:313683/141‑1 (MQ=255)
                                                        ggCTCCCGTGATTCTTTCCCGCCAACACCCGCGCATTTACCCACCGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCttt                                                                                    <  1:611007/149‑1 (MQ=255)
                                                                  aTTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAAAGaaaaaagttttaac                                                                                >  1:799094/1‑132 (MQ=25)
                                                                       ttccccccaccaCCCGCGCTTTTACCAACCCCCTTGTAAAAAAAGGCCTCCGCGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgtt                                                                      <  1:784427/142‑1 (MQ=255)
                                                                       ttcccGCCAACACCCGCGATTTTCCCCAACGCCATGGAAAAAAATGCCTGCGGGGGGGTTTCCGCCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTc                                                                           <  2:526978/140‑1 (MQ=255)
                                                                               aacaccCGCTCTTTTCCCCAACGCCATGTAAAAAAATGCCTCGGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc                                                           <  2:440972/151‑1 (MQ=255)
                                                                                 caacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGGAAAAAATTTTAACCCGTCATGTTTTGTggggat                                                            >  1:541427/1‑144 (MQ=25)
                                                                                    ccccgcTTTTACCCACCGCTTTGTGAAAAAATGCCTGCGCGGCGGCTGTCTGCAAGTTTTTTTCCCGCCGAAACAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg                                                       <  2:771707/147‑1 (MQ=255)
                                                                                    ccccgcTTTTACCCACCCCCAGGTCAAAAAAGCCTTCCGGGGGGGCTGCCTGCCGGTTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgct                                                        <  2:736871/146‑1 (MQ=255)
                                                                                    ccGCGCTTTTCCCCACCCCCTTGTGAAAAAATCCTTCCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgct                                                        <  1:350740/149‑1 (MQ=255)
                                                                                                      cctggtaaaaaaaGGCCTCCGGGGGGGCTGCCGGCAATTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAg                                    <  2:795735/143‑1 (MQ=255)
                                                                                                       attggcaaaaaaTGCCTTGGGGGGGGTTGTCTGCAGGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGa                                       <  1:402524/140‑1 (MQ=255)
                                                                                                         tgtgCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAACTTTAACCCGCCCAGTTTTGTGGGCATTGCCTCTGGTGAAAAAATGGAAAggg                                  >  1:614400/1‑150 (MQ=255)
                                                                                                          ttgaaaaatTGCCTGCGGGGCGGCTGTCTGCCAGCTTTTTTTCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                 <  1:233969/146‑1 (MQ=255)
                                                                                                                        tgcgggggggcTTTCTCCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt                     <  1:622457/141‑1 (MQ=255)
                                                                                                                                         cgtgtttttttcCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCt   <  2:532343/143‑1 (MQ=255)
                                                                                                                                         ccattttttttcCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa  <  1:389400/144‑1 (MQ=255)
                                                                                                                                             tttttttcCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc        <  2:640081/138‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
CGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTA  >  NC_000913/357733‑358019

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: