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

CGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTG  >  NC_000913/357733‑358023
                                                                                                                                                    |                                                                                                                                                
cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTc                                                                                                                                                 <  2:224275/149‑1 (MQ=255)
   aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCTCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                                                                                             >  2:1092402/1‑150 (MQ=255)
      ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaa                                                                                                                                         <  1:251860/151‑1 (MQ=255)
         gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaaca                                                                                                                                       >  2:554344/1‑150 (MQ=255)
            aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                                   >  2:996983/1‑151 (MQ=255)
                            acattcccAGCCGCAGACTTGTTTCAGGGGCCCCCGTGATTCATTCCCGCAAAAAACCGCCCTTTTACCCAACCCCAGGTCAAAAAATCCCTCCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTtctct                                                                                                                        <  2:1083644/140‑1 (MQ=255)
                             cAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                   >  2:138569/1‑150 (MQ=255)
                              aaTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtagt                                                                                                                >  2:663064/1‑151 (MQ=255)
                                  ccAGCCTCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTg                                                                                                               >  2:109120/1‑149 (MQ=255)
                                  ccAGCCGCGGCCTTGTGTCAGGGGCTCCCTTGTTTCATTTCCCCCAACACCCGCGCTTTTACCCACCGCCATGTAAAAAAATCCTTTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                                             <  1:134995/151‑1 (MQ=255)
                                    aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                           >  2:839109/1‑151 (MQ=255)
                                    aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                           >  1:901648/1‑151 (MQ=255)
                                    aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                           >  2:256701/1‑151 (MQ=255)
                                      ccGCAGCCTTGTGCCAGGGCCTCCCGTGATTCTTTCCCCCCAACACCCGCGCATTTACCCAACCCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTCCCACTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                         <  1:9772/151‑1 (MQ=255)
                                       cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                         >  1:793438/1‑150 (MQ=255)
                                        gCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggt                                                                                                        >  1:386719/1‑149 (MQ=255)
                                        gCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                           >  2:488940/1‑147 (MQ=255)
                                           gcccTTTGTCGGCGGCTACCGTGTTTCATCCCCGCCAACACCCGCGCATTTACCCAACGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAag                                                                                                     <  1:929375/148‑1 (MQ=255)
                                              cTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCACGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggcg                                                                                                     >  1:62353/1‑143 (MQ=255)
                                                        gtccccggtgtTTCATTCCCGCCACCACCCGCGTTTTAACCCACCCCATTGGGAAAAAATGCCTTCGGGGGGGCTGTCGCCAAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAgaga                                                                                                  <  2:759033/133‑1 (MQ=255)
                                                          caaCCTTGCTTCATTCCCGCCAACACCCGCGCTTTTAACCAACGCCTTGGACAAAAATTCCCCCGCGCCGTCTGTCTGCCATCTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTaa                                                                                       <  2:70570/147‑1 (MQ=255)
                                                          cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGGGAAAAAgatgtt                                                                                        >  2:407911/1‑143 (MQ=25)
                                                             cgttaaTTCTTTCCCCCAAACACCCGCGTTTTAACCCACCCCCATGTGAAAAAATGCCTCCGGGGCGCCTGTCTGCCTGTTATTTTCCCCCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCt                                                                                           <  2:675674/137‑1 (MQ=255)
                                                             ccGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGGGAGAAATaacatttaaac                                                                                    >  2:75430/1‑139 (MQ=25)
                                                             ccGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAAAGAAAAAACtttttatcc                                                                                   >  1:956624/1‑144 (MQ=25)
                                                                tGATTCATCCCCCCCAACACCCGCGCTTTTACCCAACGCCTTGTAAAAAAATGCTTCCGCGGGGGCTGTCTCCCAGCTATTTTCCCGCCAAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCg                                                                                  <  2:1070312/148‑1 (MQ=255)
                                                                               accacccGCCCATTTACCAACCGCCTGGTAAAAAATGGCCTGCGCGGCGGCTGTCGGCCATTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAAt                                                                <  1:506855/145‑1 (MQ=255)
                                                                               aacaacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGCGAGAAAAAAAATTTAACCCCCCCTGTTGTGTgggg                                                                  >  2:609613/1‑146 (MQ=21)
                                                                                acaccCGCCCATTCACCCAACCCCTGGTAAAAAAATCCCTTCGCGGCGGTTGTCTGCCAGCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCa                                                                  <  1:612365/147‑1 (MQ=255)
                                                                                        gcattcccccAACTCCATGAACAAAAAGGCCTCCGCGGCGCCTTCCTCCCAGCTTTTTTCCCGCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTgg                                                       <  2:124077/142‑1 (MQ=255)
                                                                                           tttccccacccccATGAGAAAAAATCCCTCCGGGGCGGCTTTCGCCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTa                                                     <  2:486361/138‑1 (MQ=255)
                                                                                           tttACCCAACGCCTTGTGAAAAAATGCTTCCGCGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTa                                                     <  1:627872/149‑1 (MQ=255)
                                                                                               tCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAACTTTTAAACCCTCCTGTTGTGTGGGGATTGCTGCTGGtgaaaa                                                 >  1:423473/1‑144 (MQ=25)
                                                                                                   accgccagtggaaaaaaTCCCTTCGGGGCGGCTGTCGCCCATTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACaa                                                 <  2:133230/134‑1 (MQ=255)
                                                                                                          tgtgCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAAATTTTAACCCCCCCTGTTTTGTGGGAAACGCCGCTGGGGAAAAAGTGGAAAggg                                      >  1:70883/1‑150 (MQ=255)
                                                                                                           gtgAAAAATTCCCTTCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAgg                                       <  1:937868/148‑1 (MQ=255)
                                                                                                              cAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAATTTTAACCCCCCCTGTTGTGTGGGGATCGCGTCTGGGAAAAAAGTGGGAAGGGGAggg                                 >  1:818723/1‑149 (MQ=255)
                                                                                                               aaaagtccattCGCGGCGGCTGTCTCCGCGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAg                                   <  2:48734/139‑1 (MQ=255)
                                                                                                                aaatTGCCTCCGGGGCGGCTGTCCCCCATTTTTTTTCCCCCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcg                               <  1:442887/151‑1 (MQ=255)
                                                                                                                    tGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGCGAGAAAAAAATTTTAAAACCCCCTCTTGTGGGGGAAACGCTGGTGGGGAAAATGTGGGAAGGGGAGGggggcg                            >  2:213372/1‑150 (MQ=255)
                                                                                                                       cTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGGGCAAAAAATTTAAATCCTTCCACTTTGTGGGGGATTGCTGCTGGGAACATTATGGAAAaggggagag                                >  2:813350/1‑137 (MQ=255)
                                                                                                                           gcggcgGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt                         <  1:296838/146‑1 (MQ=255)
                                                                                                                             ggcggcTTTCTCCCAGTTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCt                    <  1:225570/149‑1 (MQ=255)
                                                                                                                                         ccggtttttttcCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTAGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCt       <  1:516877/143‑1 (MQ=255)
                                                                                                                                         ccAGCTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc            <  2:896195/145‑1 (MQ=255)
                                                                                                                                              ttttttcCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg  <  1:841908/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTG  >  NC_000913/357733‑358023

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: