Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA41.7% 48.7 / 15.5 36L267F (TTG→TTT) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (17/4);  new base A (0/15);  total (17/19)
Fisher's exact test for biased strand distribution p-value = 7.21e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAACAACCGCGCATTTATCCAACG‑CCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTC  >  NC_000913/357743‑358033
                                                                                                                                                     |                                                                                                                                                  
aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAACAACCGCGCATTTATCCAACG‑CCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa                                                                                                                                                  <  2:494478/150‑1 (MQ=255)
     cGTAGCGCATCAGACAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAACAACCGCGCATTTATCCAACG‑CCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccg                                                                                                                                             <  2:258556/150‑1 (MQ=255)
      gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAACAACCGCGCATTTATCCAACG‑CCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgc                                                                                                                                            <  1:283352/150‑1 (MQ=255)
       tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAACAACCGCGCATTTATCCAACG‑CCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgccg                                                                                                                                          >  1:633252/1‑151 (MQ=255)
                         cgcccGCAGCCCTGTGTGGGC‑GCTTGCGGGATTTCATTCCCCCAACCACCCGCCCTTTTTCCCACCG‑CCTTGTAAAAAAATGCTTTCGGGGCGGTTGTCGGCCACTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑g                                                                                                                       <  1:713337/148‑1 (MQ=255)
                           gCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAACAACCGCGCATTTATCCAACG‑CCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑ggg                                                                                                                     >  1:946137/1‑151 (MQ=255)
                             cGCAGACTTGTTCCAGCGCCTACCGTGT‑TTCATTCCCCCAAACACCCGCGCTTTTCCCCACCG‑CCGCGCGAAAAATTGCCCCCGCGGGGGCTGTCTGCCAGCTATTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑g                                                                                                                       <  2:80725/147‑1 (MQ=255)
                               cAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAACAACCGCGCATTTATCCAACG‑CCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑gggg                                                                                                                    >  2:344979/1‑148 (MQ=255)
                                      gtgtCAGCGCCAACGGGGT‑TTAATTCCCCCCAACACCCGCGATTTTACCCACCC‑CCATGTCAAAAAATCCTTCCGCGGCGGCTTTCGGCCAGCTTTTTTCCCGCCGAAAAAAAGCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑aa                                                                                                           <  1:398476/149‑1 (MQ=255)
                                       ttccggcgcttcccGTGT‑TTCATCCCCGCAACCACCCGCGCTTTCCCCCACCG‑CCATGTGAAAAAAGCCCTGCGGGGCGGTTGTCTGCCAGTTATTTTCCCGCCAAAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAgag‑a                                                                                                            <  2:107420/135‑1 (MQ=255)
                                         ccAGCGTCTCCCGTGA‑TTCTTTCCCGCAAACACCCCCGCTTTTACCCACCC‑CCAGGGAAAAAAATGCCTGCGGGGGGTTTGTCTGCCGGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGc                                                                                                      <  2:677918/150‑1 (MQ=255)
                                          cggCGGCTACCGTTT‑TTCATTCCCCCCAACACCCGCCTTTTTACCCACCC‑CCATGTCAAAAAACGCCTCCGCGGCGGCTGTCGGCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAg                                                                                                       <  2:313211/147‑1 (MQ=255)
                                          cAGCGGCTCCCGTGT‑TTCATTCCCCCCAACACCCGCGCATTTACCCAACG‑CTATGTCAAAAAATGCCTCCGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑aaa                                                                                                          <  1:639819/146‑1 (MQ=255)
                                           ccgccgcctccctttttCCTTTCCCCCCACCACCCCCCCTTTTCCCAACCCCCATGTGAAAAAATCCCTCCGGGGCGGTTGTCGGCCGGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGc                                                                                                      <  2:668496/136‑1 (MQ=255)
                                               gcctcccgtatttCTTTCCCCCCACCAACCGCCCATTACCCCAACC‑CCATGTAAAAAAACCCTTCCGCGGGGGCTGTCGGCCGGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCttt                                                                                                   <  2:365774/138‑1 (MQ=255)
                                               gCTACCGTGA‑TTCATTTCCGCCAACAACCGCGCATTTATCCAACG‑CCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAG‑aaaaaaattt                                                                                                   >  1:14579/1‑141 (MQ=255)
                                                cTACCGTGA‑TTCATTTCCGCCAACAACCGCGCATTTATCCAACG‑CCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGgaaag‑aaaaaaaatt                                                                                                   >  1:873355/1‑140 (MQ=25)
                                                 tACCGTGA‑TTCATTTCCGCCAACAACCGCGCATTTATCCAACG‑CCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAG‑gaaaaaatttttga                                                                                               >  1:723870/1‑139 (MQ=25)
                                                 tACCGTGA‑TTCATTTCCGCCAACAACCGCGCATTTATCCAACG‑CCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAG‑aaaaaaattt                                                                                                   >  1:161260/1‑139 (MQ=255)
                                                  cccGTTA‑TTCTTTCCCCCCAACACCCGCGTTTTTTCCCACCG‑CCATGTACAAAAATGCCTCCGCGGCGGTTTTCGGCCACTTTTTTTCCCGCGAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGccc                                                                                             <  2:797092/150‑1 (MQ=255)
                                                     gTGA‑TTCATTTCCGCCAACAACCGCGCATTTATCCAACG‑CCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGGAG‑AAAAACTTTTAAAcc                                                                                              >  1:893819/1‑148 (MQ=25)
                                                            cATTTCCGCCAACAACCGCGCATTTATCCAACG‑CCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAG‑gaaaaaacttt                                                                                                  >  2:379016/1‑129 (MQ=25)
                                                              ttccccccAACACCCCCGCTTTTACCCACCC‑CCATGTAAAAAATTCCCTCCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGt                                                                                           <  2:895394/136‑1 (MQ=255)
                                                                 ccGCCAACACCCGCGTATTTACCCACCG‑CCATGTGAAAAAATGCTTTCGCGGCGGCTGTCTCCCAGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTgg                                                                               <  1:163057/151‑1 (MQ=255)
                                                                         aacCGCGCATTTATCCAACG‑CCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAG‑AAAAAAATTTTAGCCCCTCCTGTTGTGTggggatagc                                                                        >  1:235069/1‑143 (MQ=25)
                                                                           ccGCGCTTTTACCCACCG‑CCTTGTAAAAAAAGGCTTCCGCGGCGCCTGTCTGCAAGTTTTTTTCCCGCCAAAAAAAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgc                                                                     <  2:763021/151‑1 (MQ=255)
                                                                           ccGCGCTTTTACCCACCC‑CCTTGTCAAAAAATGCCTCCGCGGCGGTTGTCGCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc                                                                        <  1:941688/148‑1 (MQ=255)
                                                                                    tATCCAACG‑CCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAGAAAAAAAATTTAAACCC‑CCCTGTTGTGTGGGGATCGCTGCTggggaca                                                             >  1:12474/1‑145 (MQ=25)
                                                                                                aTGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTTCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGCAGAG‑AAAAAAATTTTAAACCCCCCTTTTGTGTGGGAAACGCCGCGGGTGAA‑AAAATGGGAAggg                                                >  1:422698/1‑151 (MQ=255)
                                                                                                  gtgAAAAAATGCCTCCGGGGCGGCTGTCTGCCAGTTATTTTCCCCCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAggg                                                <  1:26382/149‑1 (MQ=255)
                                                                                                   tgCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAAAG‑AAAAAAAATTAAATACGTCCTGTTGTGTggggctcgccgccg                                                                   >  2:135974/1‑118 (MQ=255)
                                                                                                            gCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAGAGAAAAAAAATTTAAACCC‑CCATGTTGTGTGGGGAACGCTGCTGGGGAC‑AAAATGGGAAGGGGAgggggggg                                      >  1:143206/1‑142 (MQ=255)
                                                                                                             ccTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGGG‑AAAAAAATTTTAACCCCCCCTTTTTTTGGGGAATCGCTGGTGGGGAC‑AAAATGGGAAGGGAAGGggggcg                                      >  2:546068/1‑148 (MQ=255)
                                                                                                                 cgcgGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCGAGAG‑AAAAAACTTTAAAACCGTCCTGTTTTGTGGGGATTGCTGCGGGGGAC‑AAAATGGGAAGGGaggggggggg                                      >  2:519447/1‑135 (MQ=255)
                                                                                                                                   acttttttcCCGCCGCAACAAAGCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa                <  2:150726/144‑1 (MQ=255)
                                                                                                                                          tcccgccgCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt        <  1:165736/147‑1 (MQ=255)
                                                                                                                                              gccgcAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGCAAAG‑AAAAAAATTTTAAACCCCCCTGTTTTGTGGGAATCCCCGCTGGGGACAAATTTGGAAGGGGAAG‑GGGGGGGATCAAATCCAGAACCACGCCAAATTAGTTTGcccc  >  1:269725/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
ATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAACAACCGCGCATTTATCCAACG‑CCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTC  >  NC_000913/357743‑358033

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