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

GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTC  >  NC_000913/357742‑358028
                                                                                                                                                       |                                                                                                                                             
gATAAGCGTAGCGCCTCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑Gcaacaa                                                                                                                                             <  1:803661/151‑1 (MQ=255)
gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑Gcaaca                                                                                                                                              <  1:223995/150‑1 (MQ=255)
 aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACaaa                                                                                                                                            <  2:41241/151‑1 (MQ=255)
   aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACaaaa                                                                                                                                           >  1:669994/1‑150 (MQ=255)
   aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑c                                                                                                                                         <  2:135736/151‑1 (MQ=255)
      cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑ccgc                                                                                                                                      <  2:970472/151‑1 (MQ=255)
        tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑ccgcc                                                                                                                                     <  2:833948/150‑1 (MQ=255)
                    cAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCagt                                                                                                                         >  2:1033389/1‑150 (MQ=255)
                         ccAGCCGCAGCCCTGTGTCAGCGGCTACCGTGTTTCATTCCCGCCAACACCCGCGCATTTACCCACCCCCATGTGAAAAAATGCTTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑ccgccg                                                                                                                                    <  1:165693/134‑1 (MQ=18)
                          cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                  >  2:113506/1‑151 (MQ=255)
                               gCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                                >  2:771791/1‑148 (MQ=255)
                                   aCCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCggggggaga                                                                                                          >  2:889914/1‑145 (MQ=255)
                                   aCCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                               >  2:244458/1‑145 (MQ=255)
                                    ccTGTGTCGGGGGCTCCCGTGTTTCTTTTCCCCCAACACCCGCGTTTTTACCCAACCCCATGGGAAAAAATCCCTCC‑GGGGCGGCTGTCTGCCAGCTATTTTC‑CCGCC‑GCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAag                                                                                                           <  2:888673/148‑1 (MQ=255)
                                        tgtCAGCGGCTCCGGTGTTTCATTCCCCCAAACACCCGCGCTTTTCCCCACCGCCATGTAAAAAAATGCCTTC‑GGGGCGTCTGTCGCCCAGTTTTTTTC‑CCGCC‑GCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga                                                                                                          <  2:880820/145‑1 (MQ=255)
                                         gtCGGCGGCTCCGGTAATTCTTTCCCCCCAACACCCGCGCTTTTACCCACCCCAATGTGAAAAAATCCCTCC‑GCGGGGGCTTTCTGCCATTTATTTTC‑CCGCC‑GAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAag                                                                                                           <  2:447368/143‑1 (MQ=255)
                                             gCGGCTCCCGTTTTTCTTTCCCCCCAACACCCGCGCTTTTACCCACCCCCTTGTAAAAAAATGCCTCC‑GGGGCGGTTTTCTGCCATTTTTTTTC‑CCGCC‑AAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑tt                                                                                              <  2:708350/150‑1 (MQ=255)
                                                   aCCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑gaaaaaaatttatacc                                                                                        >  1:1319091/1‑138 (MQ=25)
                                                   aCCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAG‑aaaaaaattttaaacc                                                                                        >  2:63188/1‑138 (MQ=25)
                                                       tatttCATTCCCCCCAACACCCGCGATTTTACCCACCGCCATGTCAAAAAAGGCCTTC‑GCGGCGGCTGTCTGCCAGGTTTTTTC‑CCGCC‑GCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGcc                                                                                        <  2:151544/143‑1 (MQ=255)
                                                           tCATTCCCCCCAACACCCACGAATTTCTCCACCCCCATGTACAAAAAGGCGTTC‑GCGGGGGTTGTCTGCACGCTTTTTTC‑CCGCC‑ACAACAAAG‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCa                                                                                   <  2:583011/147‑1 (MQ=255)
                                                             ttcccccccAACACCCGCGCTTTTACCCACCGCCTTGTCAAAAAATCCCTTC‑GGGCGGGTTTTCCGCAACCTATTTTC‑CCGCC‑GCAAAAAAG‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCAtgttg                                                                              <  2:15508/143‑1 (MQ=255)
                                                              tcccccccAACACCCGCCCTTTTA‑CCAACCCCATTGTAAAAAAACCCTTCGGGGGGGGTTTTCTCCCAGTTATTTTC‑CCCCC‑AAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCAtgttgt                                                                             <  1:445121/144‑1 (MQ=255)
                                                                 ccccaaacccccGCGCATTTACCCACCGCCAGGAGAAAAAATCCCTCG‑GGGGCGGCTGTCTGCCGTCTTTTTCC‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTtgtgtg                                                                          <  2:1275957/140‑1 (MQ=255)
                                                                 ccccaaacaccCGCTCTTTTACCCACCCCATTGTAAAAAAAGGCCTTC‑GGGGCGGTTTTCTCCCAGTTTTTTTC‑CCGCC‑GAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTtgtg                                                                            <  2:126057/143‑1 (MQ=255)
                                                                  ccccaccaCCCGCGCTTTCATCCCCCGCTGTGTAAAAAAATTCCCGC‑GGGGCGGCGTCCACCCAGCT‑TTTTT‑CCCCC‑GGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTtgtgtg                                                                          <  2:1307605/147‑1 (MQ=255)
                                                                      aacaccCCCGCATTTACCCACCGCCATGTGAAAAATTGCCTTC‑GCGGCGTCTGTCTCCCATTTATTTTC‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATc                                                                    <  2:185310/151‑1 (MQ=255)
                                                                               gcATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑AAAAAAA‑TTTTAAACCCTCCTGTTTTGTGGGAAACgc                                                                  >  2:610625/1‑144 (MQ=25)
                                                                                  tttccccACCGCCATGTGAAAAAATCCCTGC‑GGGGGGGCTTTCTCCCAGTTTTTTTC‑CCGCC‑GAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTaa                                                         <  2:1225706/145‑1 (MQ=255)
                                                                                  tttACCCACCGCCATGTAAAAAAATGCCTCC‑GCGGGGCCTGCCGCCCAGTTTTTTTC‑CCCCC‑GAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTa                                                          <  2:908709/149‑1 (MQ=255)
                                                                                        cAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGAGGCAGAG‑AAAAAAA‑TTTTAGAACCCCACGTTGTGTGGGGATCGGCGCTGGtgaaaaagt                                                   >  2:1018711/1‑142 (MQ=21)
                                                                                         aaCGCCATGTAAAAAAATCCCTGC‑GGGCCGGTTGTCGCCCATTTTTTTTC‑CCGCC‑GAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgtg                                                  <  1:1165073/150‑1 (MQ=255)
                                                                                             ccATGTGAAAAAATCCCTCC‑GCGGCGGCTGTCTCCCAGCT‑TTTTC‑CCCCC‑GCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAgg                                            <  2:1279255/151‑1 (MQ=255)
                                                                                                  tgCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAA‑AAAAAAA‑TTTAAACCCGCCACGGTTGGTGGGGAACCCTGCTGGGGACAAAGAGGGAAggg                                           >  2:486236/1‑148 (MQ=255)
                                                                                                     aaaaaGCCCTTC‑GGGGCGGCTTTCTGCCTGTTTTTTTC‑CCCCC‑GCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGaa                                         <  2:106927/147‑1 (MQ=255)
                                                                                                     aaaaaCGCCTCC‑GTGCGGGCCTTCTGCCA‑TTATTTTT‑CCCCCAACAATAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAg                                        <  2:356179/148‑1 (MQ=255)
                                                                                                           gCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGGG‑AAAAAAA‑TTTTAACCCCTCCTGTTTTGTGGGGAACGCTGCTTGGGACAAAGTGGGAAGGGGAgggggggg                                 >  2:205143/1‑142 (MQ=255)
                                                                                                              ttC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAAAA‑TTTAAACCCGCCA‑GTTTTGTGGGGAAAGCTTCTGGGGACAAAGTGGGAAGGGGAggggggggg                                >  2:420308/1‑139 (MQ=255)
                                                                                                                     gcggtTGTCTCCCAGTTTTTTTC‑CCGCC‑CAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa                      <  2:998501/151‑1 (MQ=255)
                                                                                                                     gcggcTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAAAG‑AAAAAAA‑ATTTAACCCCCCCTGTTGTGTGGGGATCCCTGCTGGTGAAAAAGTGGGAAGGGGAGGGGGGGGGATCAGATcc                       >  2:1186128/1‑150 (MQ=255)
                                                                                                                       ggcTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAA‑AAAAAAC‑ATTAAACCCGCCCTGTTTTGTGGGGATCACCGCCTGGAACAAATTGGGAAGGGGAGGgcgggggaagagctgc                       >  2:224169/1‑136 (MQ=255)
                                                                                                                           gTCTGCAAGTTTTTTTC‑CCCCC‑CAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc                 <  2:549864/150‑1 (MQ=255)
                                                                                                                             cTGCAAGTTTTTTTC‑CCGCC‑GAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAg                     <  2:1233299/144‑1 (MQ=255)
                                                                                                                               tccAGCTTTTTTTCCCGCC‑GCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                  <  1:200379/145‑1 (MQ=255)
                                                                                                                                 cAGCTATTTTT‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAAAAACATTTAAACCC‑CCCCGTTTTGTGGGAATCGCTGCCGGTGAAAAAGTGGAAAAGGGAGGGCGGGGGAACAACTcc                       >  2:394392/1‑138 (MQ=255)
                                                                                                                                           tcccgcc‑gAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc  <  2:401488/148‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                             
GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTC  >  NC_000913/357742‑358028

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