Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA45.5% 44.1 / 19.4 44L267F (TTG→TTT) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (13/11);  new base A (0/20);  total (13/31)
Fisher's exact test for biased strand distribution p-value = 4.96e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCC‑GCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCT  >  NC_000913/357741‑358032
                                                                                                                                                      |                                                                                                                                                
tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCC‑GCCGcaac                                                                                                                                                  >  1:955610/1‑150 (MQ=255)
    aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCC‑GCCGCAACAAAAc                                                                                                                                             <  1:1030269/151‑1 (MQ=255)
     aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCC‑GCCGCAACAAAAcc                                                                                                                                            <  2:633357/151‑1 (MQ=255)
       cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCC‑GCCGCAACAAAAccgc                                                                                                                                          <  2:739455/151‑1 (MQ=255)
           gcgcATCGGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCC‑GCCGCAACAAAAccgccg                                                                                                                                        <  2:300977/149‑1 (MQ=255)
                 cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCC‑GCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                                >  1:942805/1‑151 (MQ=255)
                      attcccgcccGCGGCCCTGTGCCGGCGCCCCCGGTGATTCATTCCCCCAAACACCCCCTCTTTTACCAACCCCCTTGTGAAAAATTCCCTCGGCGGCGGCGTTCCGCCAGTTTTTTTCCC‑GCCGCAACAAAACCGCCGTTCTCTCCa                                                                                                                               <  1:569439/139‑1 (MQ=255)
                          ccACCCGCAGCCCTGTGCCAGGGGCTCCCGTGATTCTTTCCCGCCAACACCCGCGCATTTCCCAACCCCCTTGTCAAAAAATGCTTCCGGGGCGGCTGTCTGCCAGTTATTTTCCC‑CCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                                                       <  2:134436/151‑1 (MQ=255)
                          ccACCCGCAGACTTGTGTCAGCGGCTCCCTTTTTTCTTCCCCCCCACCACCCGCGCTTTTACCCACCCCCATGTCAAAAAATCCCTCCGGGGGGGCTGTCTGCCAGTTTTTTTCCC‑GCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                                                       <  1:713278/151‑1 (MQ=255)
                               cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCC‑GCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg‑g                                                                                                                   >  1:832260/1‑149 (MQ=255)
                                  aGACCTGTGTCAGCGGCTCCCGTGATTCATTCCCCCCAACACCCGCGCTTTTACCCACCCCCATGTGAAAAAATCCCTTCGCGGCGGCTGTCTCCCAGCTTTTTTCCC‑GCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGc                                                                                                                 <  2:221807/148‑1 (MQ=255)
                                   gccgtgtgcaggcgcCTCCCGTGATTCTTTCCCGCAAACACCCGCGCTTTTACCCACCCCAATGTAAAAAATGGCCTCCGCGCCGGCTGTCGCCCAGTTTTTTTCCC‑CCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAag                                                                                                              <  2:898341/139‑1 (MQ=255)
                                    cccTGTGTCAGCGCCTACGGTGATTCATTCCCCCAAACAACCGCGCATTTACCCAACCCCATGTAAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGTTATTTTCCC‑GCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAag‑ag                                                                                                           <  2:467816/150‑1 (MQ=255)
                                    cccTGTGGCGGCGGCTCCCGTAATTCTTTCCCGCCAACACCCGCGCTTTTACCCACCGCCATGTGCAAAAAGCCCTCCGGGGCGGCTGTCTGCCAGTTATTTTCCC‑CCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAag‑ag                                                                                                           <  2:893415/150‑1 (MQ=255)
                                    aCCTGTGTCAGCGGCTACCGTAATTCATTCCCGCCAACAACCGCGCATTTACCAACCGCCATGTAAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGTTATTTTCCC‑GCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAag‑a                                                                                                            <  2:390671/150‑1 (MQ=21)
                                     ccTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCC‑GCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGGGag‑ag                                                                                                           >  1:776423/1‑150 (MQ=255)
                                       tgtgTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCC‑GCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg‑ggggcgag                                                                                                            >  2:871764/1‑142 (MQ=255)
                                          gtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCC‑GCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGGTCG‑AGaaa                                                                                                        >  2:1071892/1‑148 (MQ=255)
                                                 gCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCC‑GCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGGAAG‑AGgaaaaaaatt                                                                                                 >  2:49384/1‑141 (MQ=255)
                                                    aCCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCC‑GCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGGGAG‑AGGAAAAGCatt                                                                                                 >  1:957750/1‑145 (MQ=255)
                                                    aCCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCC‑GCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCGAG‑AGaaaaaaattttatgc                                                                                            >  2:556878/1‑138 (MQ=255)
                                                     ccGTGTTTCATTCCCCCAAACACCCGCGAATTTACCAACCCCCTTGTGAAAAAATCCTTCCGGGGCGGCTGTCTCCCATTTTTTTTCCC‑GCCACAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAG‑AGAAATAGCTTTAAGccc                                                                                           <  2:917175/150‑1 (MQ=255)
                                                        tGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCC‑GCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGTGGAAAG‑GGCAAAAACGTTAAGcc                                                                                            >  1:607173/1‑147 (MQ=25)
                                                        tGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCC‑GCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGGCAG‑AGgaagaacattgttatc                                                                                           >  1:29668/1‑134 (MQ=25)
                                                         gATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCC‑GCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCGAG‑AGaaaaaagttttaaac                                                                                            >  2:531156/1‑133 (MQ=255)
                                                           ttcttccccccaaacaccCGCGCTTTTACCCAACGCATTGTAAAAAAAGGCCTCCGGGGCGGCTTCCGGCCAGTTTTTTTCCC‑GCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAG‑AGAAATAGCTTTAAGCCCGTCAt                                                                                      <  1:341218/137‑1 (MQ=255)
                                                               ttccccccAACACCCGCCCTTTTCCCCACCGCCTTGGGAAAAAATGCCTCCGCGGGGGTTGTCTGCCACTTTTTTTCCC‑GCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAG‑AGAAATAGCTTTAAGCCCGTCAtgt                                                                                    <  1:374898/141‑1 (MQ=255)
                                                                         cacccccccttttACCCACCCCCATGTAAAAAACTGCCTCCGCGCCGGCTGCCTGCCAGTTTTTTTCCC‑GCCAAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc                                                                      <  1:528614/141‑1 (MQ=255)
                                                                                     tATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCC‑GCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCGAG‑AGAAAAAAATTTTAATCCGCCACGTTTTGTggggaatgatcgcgg                                                                >  2:703606/1‑132 (MQ=25)
                                                                                          accaccATGTGAAAAAACCCCTCCGGGGCGGCTTTCTCCAAGTTTTTTTCCC‑GCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg                                                        <  2:475069/144‑1 (MQ=255)
                                                                                          accGCCTTGTGAAAAATTCCCTTCGCGGCGGTTTTCGGCAAGTTTTTTTCCC‑GCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgt                                                       <  1:516524/147‑1 (MQ=255)
                                                                                           accataTGTGAAAAAATGCCTCCGCGGCGCCGGTCCCCCACTTTTTTTTCC‑CCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgtg                                                      <  2:99094/144‑1 (MQ=255)
                                                                                              caatttgCAAAAATGCTTTGGCGGCGGCTGTCTGCAACATTTTTTCCC‑CCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAt                                                         <  2:170755/138‑1 (MQ=255)
                                                                                                  gtgCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCC‑GCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCGAGCAGAAAAAACTTTAATCCC‑CCCTGTTTTGTGGGGATCGCCTCTgg                                                                >  1:258969/1‑132 (MQ=255)
                                                                                                  gggaaaaaaTGCCTCCGCGGCGGCTGTCTGCCGTTTTTTTTCCC‑CCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                              <  1:627731/146‑1 (MQ=255)
                                                                                                      aaaaaTCCCTCCGCGGCGGCTGTCGGCCAGTTTTTTTCCC‑GCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTg                                          <  1:73007/150‑1 (MQ=255)
                                                                                                                 cgcgGCGCTTGTCCCCCTTCTTTTTTTCC‑CCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCa                                <  1:76074/149‑1 (MQ=255)
                                                                                                                  gcggggGTTGTCGCCCAGTTTTTTTCCC‑GCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAg                               <  1:362259/149‑1 (MQ=255)
                                                                                                                    ggcggcTTTCTGCCAGTTTTTTTCCC‑CCCAAAACAAAGCCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggcg                                     <  1:379658/141‑1 (MQ=255)
                                                                                                                       ggcTGTCTGCCAGTTTTTTTCCC‑GCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa                          <  2:717084/149‑1 (MQ=255)
                                                                                                                                  cagtttttttCCGGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa               <  1:312661/145‑1 (MQ=255)
                                                                                                                                     ttttttcCC‑GCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg           <  2:296302/148‑1 (MQ=255)
                                                                                                                                       ttttcCC‑GCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg         <  1:155996/150‑1 (MQ=255)
                                                                                                                                             c‑gccacAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCt  <  2:624139/151‑1 (MQ=255)
                                                                                                                                                ccgaAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCg     <  1:315624/146‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
TGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCC‑GCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCT  >  NC_000913/357741‑358032

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: ‑