Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA53.8% 16.2 / 58.1 39L267F (TTG→TTTcynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (14/4);  new base A (0/21);  total (14/25)
Fisher's exact test for biased strand distribution p-value = 2.03e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.59e-05

GCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATG  >  NC_000913/357747‑358037
                                                                                                                                                |                                                                                                                                                     
gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccg                                                                                                                                                >  2:243355/1‑151 (MQ=255)
 cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgc                                                                                                                                               <  1:1209719/151‑1 (MQ=255)
  gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgcc                                                                                                                                              >  1:1014721/1‑151 (MQ=255)
    aGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGt                                                                                                                                            >  2:892951/1‑151 (MQ=255)
           cAGGCAATTCCAGCCGCAGACCTGTGCCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctc                                                                                                                                        >  1:726604/1‑148 (MQ=255)
                        ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                                         >  1:841571/1‑150 (MQ=255)
                          gCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                        >  1:975387/1‑149 (MQ=255)
                           cagccttgtgTCAGCGGCTCCCGTGATTCATTCCCGCCAACACCCGCGCTTTTACCCAACGCCATGTGCAAAAATGCTTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCa                                                                                                                      <  1:51395/144‑1 (MQ=21)
                           cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                                         >  1:84397/1‑147 (MQ=255)
                             gcccTGTGTCAGCGGCTCCCGTGATTCTTTCCCCCCACCACCCGCGCATTTACCCACCGCCATGCCAAAAAATGCCTGCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAag                                                                                                                    <  2:410939/148‑1 (MQ=255)
                                 tgtgTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGcgcggg                                                                                                                  >  2:278737/1‑143 (MQ=255)
                                  gtgtCAGCGGCACCCGTAATTCATTCCCGCCAACACCCGCGCTTTTACCCAACGCAATGTGAAAAATTGCCTTCGCGGCGGCTGTCTGCCAGTTTTTTTTCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaa                                                                                                               <  1:571519/149‑1 (MQ=255)
                                     tCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggaaggaaa                                                                                                               >  2:832191/1‑138 (MQ=255)
                                          ggCTCCCGTGTTTCATTCCCCCAAACACCCGCGCTTTTCCCAACCGCCTTGTGCAAAAATCCCTTGGGGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAagag                                                                                                                  <  1:542443/139‑1 (MQ=255)
                                          gcctccgttatttCTTTCCCCCCAACACCCCCCCTTTTACCCACCCCATTGTGAAAAAATGCCTCCGGGGGGTTTGTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTaa                                                                                                    <  2:231432/140‑1 (MQ=255)
                                               ccGTGATTCATTCCCGCCAACACCCGCGCATTTACCCAACCCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCATTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGccc                                                                                                <  1:118006/150‑1 (MQ=25)
                                                          ttCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAAAGAAAAAAATTTaaaaccccccctttttg                                                                                     >  1:29861/1‑135 (MQ=25)
                                                            ccGCCAACACCCGCGCATTTACCCACCCCCATGTGAAAAAATGCTTTCGGGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtgtg                                                                                   <  1:706877/150‑1 (MQ=255)
                                                              gCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAGCTTTTAAAACCGCCTGTTGTGTGGcg                                                                                >  1:1096841/1‑150 (MQ=25)
                                                                      ccGCCCTTTTACCCAACGCATTTTGAAAAAATGCTTCCGGGGCGGCTTTCTGCCAGTTTTTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc                                                                           <  1:1146738/148‑1 (MQ=255)
                                                                        gcgcATTTATCCACCGCCATGTGCAAAAATCCTTTCGCGGCGGCTTTTTCCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgct                                                                       <  1:128241/150‑1 (MQ=255)
                                                                            aTTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAATTTTAAACCCCCCTGTTGTGTGGGGATCGCTGCTgg                                                                     >  1:879485/1‑148 (MQ=25)
                                                                             ttactccacccctaTGTGAAAAAAGGCCTTCGCGGCGGTTGTCGGCCAGTTTTTTTCCGGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGctgctg                                                                      <  1:258390/133‑1 (MQ=255)
                                                                                  ccACCGCCTTGTAAAAAAATGCCCCCGCGGCGTGTGTCGCCCACCTTTTTTCCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgt                                                            <  2:408038/151‑1 (MQ=255)
                                                                                    acctccttgtgAAAAAAGGCCTCCGCGGGGGCTGTCTCCCGGCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgtg                                                           <  1:122420/143‑1 (MQ=255)
                                                                                            gtgCAAAAATTCTTCCGCGGGGGTTGTCTGCCAGCTTTTTTTCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg                                                    <  2:192610/149‑1 (MQ=255)
                                                                                             tgaaaaatTCCCTTCGCGGCGGCTGTCGCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGaa                                                  <  1:1200253/147‑1 (MQ=255)
                                                                                                aaaaaTGCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgc                                              <  1:775079/151‑1 (MQ=255)
                                                                                                aaaaaTCCCTTCGGGGGGGTTGTCTGCCAGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgc                                              <  1:125454/151‑1 (MQ=255)
                                                                                                aaaaaTCCCTTCGCGGCGGTTGTCGCCCATTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGt                                                <  2:1166190/149‑1 (MQ=255)
                                                                                                 aaaaAGCCCTCCCGGCGGCCTTCCGCCACGCTTTTTCCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGaa                                                  <  2:660123/147‑1 (MQ=255)
                                                                                                       ccTCCGGGGCGGCTGTCTGCCAGTTTTTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt                                       <  1:857688/151‑1 (MQ=255)
                                                                                                           cgcgGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAAAGAAAAAAATTTTAATCCCGCATTTTTTGGGGGCATAGCTGCTGGGTAAAATGGGGGAAGGGGAgggggggggaacaa                                    >  1:943409/1‑137 (MQ=255)
                                                                                                                     tgccgccgattttttttcCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGa                             <  2:886695/135‑1 (MQ=255)
                                                                                                                          ccAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcg                         <  1:601815/146‑1 (MQ=255)
                                                                                                                               cTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAAGAAAAGAGTTTAAACCGCCCCGTTGTGGGGCAATAGCCGCGTGGGACAAAATGGGAAGGGGAGAGGGGCGGATCAGCTCCACAAACgcgcccatt                 >  1:551112/1‑145 (MQ=255)
                                                                                                                                  ttttcCCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg              <  2:1135941/150‑1 (MQ=255)
                                                                                                                                      cccgccgAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc           <  2:1090671/148‑1 (MQ=255)
                                                                                                                                          ccgaAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTc      <  2:541791/150‑1 (MQ=255)
                                                                                                                                             caacaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGGGAAAAAACTTTTAACCCCCCATGTTTTGTGGGAATTGCTGCTGGGAAAAAAGTGGGAAGGGGAGAGGGGGGGATCAGCCCCCGAAACGCGCCCATTGAGTTTGCCTCCATg  >  1:430697/1‑151 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
GCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATG  >  NC_000913/357747‑358037

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: