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→TTT) cynRtranscriptional 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

GCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATG  >  NC_000913/357747‑358037
                                                                                                                                                |                                                                                                                                                     
gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAAccg                                                                                                                                                >  2:243355/1‑151 (MQ=255)
 cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAAccgc                                                                                                                                               <  1:1209719/151‑1 (MQ=255)
  gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAAccgcc                                                                                                                                              >  1:1014721/1‑151 (MQ=255)
    aGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGt                                                                                                                                            >  2:892951/1‑151 (MQ=255)
           cAGGCAATTCCAGCCGCAGACCTGTGCCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTtctc                                                                                                                                        >  1:726604/1‑148 (MQ=255)
                        ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                                         >  1:841571/1‑150 (MQ=255)
                          gCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                        >  1:975387/1‑149 (MQ=255)
                           cagccttgtgTCAGCGGCTCCCGTGATTCATTCCCGCCAACACCCGCGCTTTTACCCAACGCCATGTGCAAAAATGCTTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCa                                                                                                                      <  1:51395/144‑1 (MQ=21)
                           cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                                         >  1:84397/1‑147 (MQ=255)
                             gcccTGTGTCAGCGGCTCCCGTGATTCTTTCCCCCCACCACCCGCGCATTTACCCACCGCCATGCCAAAAAATGCCTGCGCGGCGGCTGTCTGCCA‑GCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAag                                                                                                                    <  2:410939/148‑1 (MQ=255)
                                 tgtgTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGcgcggg                                                                                                                  >  2:278737/1‑143 (MQ=255)
                                  gtgtCAGCGGCACCCGTAATTCATTCCCGCCAACACCCGCGCTTTTACCCAACGCAATGTGAAAAATTGCCTTCGCGGCGGCTGTCTGCCA‑GTTTTTTTTCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑aa                                                                                                               <  1:571519/149‑1 (MQ=255)
                                     tCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggaaggaaa                                                                                                               >  2:832191/1‑138 (MQ=255)
                                          ggCTCCCGTGTTTCATTCCCCCAAACACCCGCGCTTTTCCCAACCGCCTTGTGCAAAAATCCCTTGGGGGCGGCTGTCTCCCA‑GTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAagag                                                                                                                  <  1:542443/139‑1 (MQ=255)
                                          gcctccgttatttCTTTCCCCCCAACACCCCCCCTTTTACCCACCCCATTGTGAAAAAATGCCTCCGGGGGGTTTGTCTCCCA‑GTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTaa                                                                                                    <  2:231432/140‑1 (MQ=255)
                                               ccGTGATTCATTCCCGCCAACACCCGCGCATTTACCCAACCCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑TTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGccc                                                                                                <  1:118006/150‑1 (MQ=25)
                                                          ttCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAAAG‑AAAAAAA‑TTTaaaaccccccctttttg                                                                                     >  1:29861/1‑135 (MQ=25)
                                                            ccGCCAACACCCGCGCATTTACCCACCCCCATGTGAAAAAATGCTTTCGGGGCGGCTGTCTGCCA‑GCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTtgtgtg                                                                                   <  1:706877/150‑1 (MQ=255)
                                                              gCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAG‑AAAAAGCTTTTAAAACCG‑CCTGTTGTGTGGcg                                                                                >  1:1096841/1‑150 (MQ=25)
                                                                      ccGCCCTTTTACCCAACGCATTTTGAAAAAATGCTTCCGGGGCGGCTTTCTGCCA‑GTTTTTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCgc                                                                           <  1:1146738/148‑1 (MQ=255)
                                                                        gcgcATTTATCCACCGCCATGTGCAAAAATCCTTTCGCGGCGGCTTTTTCCCA‑GCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCgctgct                                                                       <  1:128241/150‑1 (MQ=255)
                                                                            aTTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑AAAAAAA‑TTTTAAACCCCCCTGTTGTGTGGGGATCGCTGCTgg                                                                     >  1:879485/1‑148 (MQ=25)
                                                                             ttactccacccctaTGTGAAAAAAGGCCTTCGCGGCGGTTGTCGGCCA‑GTTTTTTTCCGGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGctgctg                                                                      <  1:258390/133‑1 (MQ=255)
                                                                                  ccACCGCCTTGTAAAAAAATGCCCCCGCGGCGTGTGTCGCCCA‑CCTTTTTTCCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgt                                                            <  2:408038/151‑1 (MQ=255)
                                                                                    acctccttgtgAAAAAAGGCCTCCGCGGGGGCTGTCTCCCG‑GCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgtg                                                           <  1:122420/143‑1 (MQ=255)
                                                                                            gtgCAAAAATTCTTCCGCGGGGGTTGTCTGCCA‑GCTTTTTTTCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg                                                    <  2:192610/149‑1 (MQ=255)
                                                                                             tgaaaaatTCCCTTCGCGGCGGCTGTCGCCCA‑GTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGaa                                                  <  1:1200253/147‑1 (MQ=255)
                                                                                                aaaaaTGCCTTCGCGGCGGCTGTCTGCCA‑GCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgc                                              <  1:775079/151‑1 (MQ=255)
                                                                                                aaaaaTCCCTTCGGGGGGGTTGTCTGCCA‑GTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgc                                              <  1:125454/151‑1 (MQ=255)
                                                                                                aaaaaTCCCTTCGCGGCGGTTGTCGCCCA‑TTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGt                                                <  2:1166190/149‑1 (MQ=255)
                                                                                                 aaaaAGCCCTCCCGGCGGCCTTCCGCCACGCTTTTTCCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGaa                                                  <  2:660123/147‑1 (MQ=255)
                                                                                                       ccTCCGGGGCGGCTGTCTGCCA‑GTTTTTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt                                       <  1:857688/151‑1 (MQ=255)
                                                                                                           cgcgGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAAAGAAAAAAAT‑TTTAATCCCG‑CATTTTTTGGGGGCATAGCTGCTGGGTAAAATGGGGGAAGGGGAgggggggggaacaa                                    >  1:943409/1‑137 (MQ=255)
                                                                                                                     tgccgccgattttttttcCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGa                             <  2:886695/135‑1 (MQ=255)
                                                                                                                          ccA‑GTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcg                         <  1:601815/146‑1 (MQ=255)
                                                                                                                               cTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAA‑GAAAAGA‑GTTTAAACCGCCCCGTTGTGGGGCAATAGCCGCGTGGGACAAAATGGGAAGGGGAGAGGGGCGGATCAGCTCCACAAACgcgcccatt                 >  1:551112/1‑145 (MQ=255)
                                                                                                                                  ttttcCCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg              <  2:1135941/150‑1 (MQ=255)
                                                                                                                                      cccgccgAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc           <  2:1090671/148‑1 (MQ=255)
                                                                                                                                          ccgaAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTc      <  2:541791/150‑1 (MQ=255)
                                                                                                                                             caacaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGGG‑AAAAAAC‑TTTTAACCCCCCATGTTTTGTGGGAATTGCTGCTGGGAAAAAAGTGGGAAGGGGAGAGGGGGGGATCAGCCCCCGAAACGCGCCCATTGAGTTTGCCTCCATg  >  1:430697/1‑151 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
GCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATG  >  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: ‑