Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8800TC47.6% 41.3 / 16.1 42K271E (AAA→GAA) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (15/7);  new base C (0/20);  total (15/27)
Fisher's exact test for biased strand distribution p-value = 1.89e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCC‑GCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAG  >  NC_000913/357732‑358025
                                                                                                                                                    |                                                                                                                                                    
tCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTcc                                                                                                                                                    >  1:918274/1‑151 (MQ=255)
 cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTcc                                                                                                                                                    <  1:480910/150‑1 (MQ=255)
         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)
                                                                                                                                                    |                                                                                                                                                    
TCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCC‑GCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAG  >  NC_000913/357732‑358025

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