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

CGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTA  >  NC_000913/357733‑358019
                                                                                                                                                    |                                                                                                                                             
cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTccg                                                                                                                                            >  2:229829/1‑151 (MQ=255)
  tAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTTTGCCAGCTATTTTTccgcc                                                                                                                                          >  2:359110/1‑151 (MQ=255)
           tAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑aaa                                                                                                                                 <  2:639455/150‑1 (MQ=255)
            aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAAc                                                                                                                               >  1:352699/1‑151 (MQ=255)
               cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAAccgc                                                                                                                            <  2:28351/151‑1 (MQ=255)
               cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGCGATTCATTTCCGCAAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAAccgc                                                                                                                            >  2:201639/1‑151 (MQ=255)
                     gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTtc                                                                                                                       >  2:73675/1‑150 (MQ=255)
                        tCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGTCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGAtctctc                                                                                                                   >  2:409593/1‑151 (MQ=255)
                            gCAATTCAAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCagt                                                                                                               >  2:407856/1‑151 (MQ=255)
                             cAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCagt                                                                                                               >  2:18290/1‑150 (MQ=255)
                                                  ggcacggtctcccGTGATTCATCCCCCCCAACACCCGCCTATTTTCCCACCGCCAGGTGCAAAACGGCTTC‑CGGGGCGCCTGTCGGCAAGTTTTTTTCCCGCCCAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAg                                                                                        <  2:585305/140‑1 (MQ=255)
                                                    cAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAAAG‑gaaaaaat                                                                                      >  1:626449/1‑146 (MQ=255)
                                                     ggCGGCAACCGTAATTCATTCCCCCCACCACCCCCTCTTTCATCCCCCGCCGTGTTAAAAAATGCGTC‑CGCGGCGCCTGTCTCCCA‑TTTTTTTTCCCCCCAAACAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCtt                                                                                     <  2:458505/150‑1 (MQ=255)
                                                      gCGGCTCCGGTGTTTCTTTCCCGCCAACACCCGCGCATTTACCCACCGCCATGTGAAAAAAGGCCTT‑CGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCt                                                                                      <  1:621502/149‑1 (MQ=255)
                                                        ggCTCCCGTGTTTCTTTCCCCCCAACACCCGCGCATTTACCCACCCCCATGTAAAAAAATCCCTT‑CGGGGCGGCTGTCGGCCAGTTTTTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑aa                                                                                            <  1:313683/141‑1 (MQ=255)
                                                        ggCTCCCGTGATTCTTTCCCGCCAACACCCGCGCATTTACCCACCGCCATGTGAAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCttt                                                                                    <  1:611007/149‑1 (MQ=255)
                                                                  aTTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAAAG‑aaaaaagttttaac                                                                                >  1:799094/1‑132 (MQ=25)
                                                                       ttccccccaccaCCCGCGCTTTTACCAACCCCCTTGTAAAAAAAGGCCTC‑CGCGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAtgtt                                                                      <  1:784427/142‑1 (MQ=255)
                                                                       ttcccGCCAACACCCGCGATTTTCCCCAACGCCATGGAAAAAAATGCCTG‑CGGGGGGGTTTCCGCCCAGTTTTTTTCCCGCCAAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTc                                                                           <  2:526978/140‑1 (MQ=255)
                                                                               aacaccCGCTCTTTTCCCCAACGCCATGTAAAAAAATGCCTC‑GGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCAAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc                                                           <  2:440972/151‑1 (MQ=255)
                                                                                 caacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑GAAAAAATTTTAACCCGTCATGTTTTGTggggat                                                            >  1:541427/1‑144 (MQ=25)
                                                                                    ccccgcTTTTACCCACCGCTTTGTGAAAAAATGCCTG‑CGCGGCGGCTGTCTGCAAGTTTTTTTCCCGCCGAAAC‑AAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg                                                       <  2:771707/147‑1 (MQ=255)
                                                                                    ccccgcTTTTACCCACCCCCAGGT‑CAAAAAAGCCTTCCGGGGGGGCTGCCTGCCGGTTTTTTTCCCGCCACAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgct                                                        <  2:736871/146‑1 (MQ=255)
                                                                                    ccGCGCTTTTCCCCACCCCCTTGTGAAAAAATCCTTC‑CGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgct                                                        <  1:350740/149‑1 (MQ=255)
                                                                                                      cctggtaaaaaaaGGCCTC‑CGGGGGGGCTGCCGGCAATTTTTTTTCCCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAg                                    <  2:795735/143‑1 (MQ=255)
                                                                                                       attggcaaaaaaTGCCTT‑GGGGGGGGTTGTCTGCAGGTTTTTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGa                                       <  1:402524/140‑1 (MQ=255)
                                                                                                         tgtgCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAACTTTAA‑CCCGCCCAGTTTTGTGGGCATTGCCTCTGGTGAAAAAATGGAAAggg                                  >  1:614400/1‑150 (MQ=255)
                                                                                                          ttgaaaaatTGCCTG‑CGGGGCGGCTGTCTGCCAGCTTTTTTTCCCCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                 <  1:233969/146‑1 (MQ=255)
                                                                                                                        tgcgggggggcTTTCTCCCAGCTATTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt                     <  1:622457/141‑1 (MQ=255)
                                                                                                                                         cgtgtttttttcCCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCt   <  2:532343/143‑1 (MQ=255)
                                                                                                                                         ccattttttttcCCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa  <  1:389400/144‑1 (MQ=255)
                                                                                                                                             tttttttcCCCCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc        <  2:640081/138‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
CGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTA  >  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: ‑