Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA61.2% 4.0 / 48.1 50L267F (TTG→TTTcynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (14/5);  new base A (0/30);  total (14/36)
Fisher's exact test for biased strand distribution p-value = 1.72e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.73e-03

GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTT  >  NC_000913/357742‑358027
                                                                                                                                                           |                                                                                                                                           
gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa                                                                                                                                           >  1:1145070/1‑151 (MQ=255)
    aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                                       >  2:948607/1‑150 (MQ=255)
     gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccg                                                                                                                                     >  1:545504/1‑151 (MQ=255)
     gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccg                                                                                                                                     >  1:1613410/1‑151 (MQ=255)
     gCGTAGCGAATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCTGCTGTCTGCCAGCTATTTTTTCGCCGCAACAAAAccg                                                                                                                                     >  1:445459/1‑151 (MQ=255)
      cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccg                                                                                                                                     <  1:89342/150‑1 (MQ=255)
                  ggCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCag                                                                                                                        >  1:409512/1‑151 (MQ=255)
                   gCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                       >  2:546493/1‑151 (MQ=255)
                         ccACCCGCAGCCCGGTGTCGGCGGCTCCCGTGTTTCTTTCCCCCCACCACCCGCGCTTTTACCAACCCCCATGTAAAAAAATCCCTCCGGGGCGGCTGTCTCCAGGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                                                 <  1:220148/151‑1 (MQ=255)
                             ccGCTGATCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                            >  1:109193/1‑151 (MQ=255)
                             ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                            >  1:1406907/1‑151 (MQ=255)
                                  gACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggagg                                                                                                       >  1:915892/1‑146 (MQ=255)
                                      tgtgTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCCAGAGaaa                                                                                                   >  2:670461/1‑151 (MQ=255)
                                       gtgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGGAGaaa                                                                                                   >  1:1620230/1‑151 (MQ=255)
                                       gtgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggacgggaaa                                                                                                  >  1:1730209/1‑141 (MQ=255)
                                         tccggCGGCTCCTGTGTTTCTTTCCCCCCACCACCCGCGCTTTTACCAACCGCCATGTGAAAAAATCCCTTCGGGGCGGCTGTCTCCAGGTTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAg                                                                                                <  1:1325786/147‑1 (MQ=255)
                                         gtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGggaaaggaaaaaa                                                                                                >  2:879291/1‑148 (MQ=255)
                                          tcggcgcCTACCTTGTTTCTTTCCCCCCAACACCCCCCCTTTTACCCACCCCATTTTGAAAAAATCCTTTCGGGGCGGCTGTCGGCCAGTTTTTTTCCCGCCACAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAg                                                                                                <  2:1593542/147‑1 (MQ=255)
                                          tCAGGGGCTACCGTGATTATTTCCCCCCAACACCCGCGCTTTTACCCACCGCCATGTGAAAAAATGCCTCCGGGGGGGTTGTCGGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAg                                                                                                <  2:1075527/150‑1 (MQ=255)
                                           cgggggCAACGGTGATTCTTTCCCCCAAACACCCGCGCTTTTACCCACCGCAAGGTGAAAAAATCCTTCCGGGGCGGCTGTCTGCCGGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGc                                                                                               <  2:1256241/146‑1 (MQ=255)
                                                  tcccGTTATTCATTCCCCCCAACACCCGCGTTTTTCCCCACCCCCAGTGAAAAAAATCCTTCGGGGGGTCCTTTCTGCAGTCTTTTTCCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc                                                                                       <  1:987096/149‑1 (MQ=255)
                                                  cacggtgtTTCATCCCCCCCACCACCCGCGATTTTATCCAACGCCATGTGAAAAAGTCCCTCCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCGGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc                                                                                       <  1:1090650/147‑1 (MQ=255)
                                                   aCCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGgaaaaaaatttatac                                                                                       >  2:1136716/1‑138 (MQ=25)
                                                            ctttcccGACAACACCCGCGCTTTACTCCACCCCCATGACAAAAAAGACCTCCGCGGCGGCTGCCGCCCGGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttg                                                                             <  1:1075255/146‑1 (MQ=255)
                                                            catccccccaaacccccGCACATTTCTCCAACGCCAGGTGCAAAAATCCCTTCGCGGGGGCTGTCTGCAAGTTATTTTCCCGCCAAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttg                                                                             <  1:403137/136‑1 (MQ=255)
                                                                 ccccaaacaccCCCGCATTAACCCACCGCCGGGTACAAAAATCCCTTGCGGGGGGCTGTCTGCAAGTTTTTTTCCCCCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTgg                                                                        <  1:1061975/146‑1 (MQ=255)
                                                                  ccccAACACCCCCGCATTCCTCCAACCCCGGGCAAAAAAATCCCTCCGGGGCGGCTGTCGGCCAGTTATTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGc                                                                       <  1:1764011/149‑1 (MQ=255)
                                                                   ccaaacaccCGCGATTTTACCCACCCCCATGTGCAAAAAGCCCTCCGGGGGGGCTGTCGCCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCa                                                                      <  2:1275665/148‑1 (MQ=255)
                                                                    caaacaccCGCGATTTTACCCACCGCCATGTCAAAAATTGCCTTCGGGGGGGCTGTCTGCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCaa                                                                     <  2:826089/149‑1 (MQ=255)
                                                                        cacccGCGCTTTTACCCACCGCAATGGGAAAAAATCCCTCCGGGGGGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc                                                                 <  2:1542018/148‑1 (MQ=255)
                                                                        cacccGCGCATTTCCCCACCCCCATGTGAAAAAATGCTTCCGGGGGGGCTTTCGCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc                                                                 <  2:72773/148‑1 (MQ=255)
                                                                                cattaccccAACGCCATGTGCAAAAATGCCTCCGGGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTa                                                         <  1:152231/144‑1 (MQ=38)
                                                                                  ttcccccgaccccATGGGAAAAAATGCCTCCGCGGCGCTTTCGGGCGAGATTTTTTCCCGCAGAACAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTAGTGTGGCAATCGCTGCTGGTAAc                                                       <  1:916950/140‑1 (MQ=255)
                                                                                     acccaccccaaTGTGAAAAAATCCTTCCGCGGGGGCTGTCGGCAACTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAt                                                    <  2:1715860/141‑1 (MQ=255)
                                                                                      cccAACCCCTTGTACAAAAGTCCCTTCGCGGCGGCTTTCTCCCGGTTTTTTTCCCGCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg                                                   <  1:1012489/150‑1 (MQ=255)
                                                                                      ccaaccccctggcaaaaaaaGCCCTGCGCGGCGCTTGTCTGCCGGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg                                                   <  1:292997/136‑1 (MQ=255)
                                                                                       ccacccctatgtcaaaaatTCCCTGCGCGGCGGCTGTCTGCCAGTTTTTTTACCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgt                                                  <  1:976903/137‑1 (MQ=255)
                                                                                             cccgtttcAAAAAAGCCTTTGGGGGGGTTGTTCTGCCATTTTTTTTTCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaaa                                             <  2:889788/144‑1 (MQ=255)
                                                                                               atgtaaaaaaaTCCTTGCGGGGGGGCTGTCGGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg                                          <  2:1093109/145‑1 (MQ=255)
                                                                                                 tgaaaaaaaTGCCCTGGGGGGCGGCTGTCGGCCAGTTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                         <  1:768995/147‑1 (MQ=255)
                                                                                                 gtcaaaaaaTGCCTCCGCGGCTGTTGTCTCCAACTTTTTTTCCCCCGGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                         <  2:1144492/146‑1 (MQ=255)
                                                                                                                            ctttttcccactttttttcCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                 <  2:672838/137‑1 (MQ=255)
                                                                                                                            cgcatggCCCGCTATTTGCCCGCCAAAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGGGGGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                 <  1:652104/145‑1 (MQ=255)
                                                                                                                               ccggCCAGGTATTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcg               <  2:369335/147‑1 (MQ=255)
                                                                                                                                       tttttttcCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg      <  2:1423202/148‑1 (MQ=255)
                                                                                                                                       cttttttcCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg      <  2:1494916/148‑1 (MQ=255)
                                                                                                                                        ttttttccccgccaAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg      <  1:1397620/143‑1 (MQ=255)
                                                                                                                                           ttttCCCCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt  <  2:24324/151‑1 (MQ=255)
                                                                                                                                           tttcCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt  <  1:813270/151‑1 (MQ=255)
                                                                                                                                            tttCCGCCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt  <  2:638152/151‑1 (MQ=255)
                                                                                                                                             ttcccgcccCAAGAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt  <  2:1507096/148‑1 (MQ=255)
                                                                                                                                                           |                                                                                                                                           
GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTT  >  NC_000913/357742‑358027

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: