Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA44.0% 23.1 / 10.7 25L267F (TTG→TTT) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (11/3);  new base A (0/11);  total (11/14)
Fisher's exact test for biased strand distribution p-value = 1.16e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGC‑ATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGA  >  NC_000913/357741‑358038
                                                                                                                                                       |                                                                                                                                                      
tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGC‑ATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaac                                                                                                                                                        >  2:219232/1‑150 (MQ=255)
 gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGC‑ATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaacaa                                                                                                                                                      <  2:579069/151‑1 (MQ=255)
 gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGC‑ATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaac                                                                                                                                                        <  2:341052/149‑1 (MQ=255)
   tAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGC‑ATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACaaa                                                                                                                                                     >  2:29825/1‑150 (MQ=255)
      gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGC‑ATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAccg                                                                                                                                                 <  1:370891/151‑1 (MQ=255)
            cgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGC‑ATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTtc                                                                                                                                           >  2:570210/1‑151 (MQ=255)
             gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGC‑ATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTtct                                                                                                                                          >  2:199486/1‑151 (MQ=255)
                             gCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGC‑ATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                           >  1:169651/1‑150 (MQ=255)
                              ccGCAGCCTTGTGTCAGCGGCTCCCGTGATTCTTTCCCGCCAACACCCGCGC‑TTTTACCCACCGCCTTGTGAAAAAATCCCTCCGGGGCGGCTGTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                            <  2:254116/148‑1 (MQ=255)
                                   gACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGC‑ATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggg                                                                                                                        >  2:325109/1‑146 (MQ=255)
                                     ccTGTGTCGGCGCCTCCCGTGTTTCATTCCCCCCAACACCCGCCT‑TTTTACCCACCCCCATGTGAAAAAATGCCTGCGGGGCGGCTGTCTGCCAGTTATTTTC‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAag‑ag                                                                                                                  <  2:278963/150‑1 (MQ=255)
                                     ccTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGC‑ATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAag                                                                                                                     >  2:215879/1‑148 (MQ=255)
                                         gtggcgccgcctccGGGATTCCTTCCCCCCAACAACCGCCCAATTTACCCAACCCCATGCAAAAAAACGCCTGCGGGCCGGCTGTCCCCAACCTATTTGC‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAG‑AGAAAt                                                                                                              <  1:29821/139‑1 (MQ=255)
                                            cggCGGCTACCGTTATTCACTCCCCCCAACACCCGCTT‑ATTTACCCACCCCCTTGAGAAAAAACCCCTCCGCGGCGGTTTTCTGCCAGTTTTTTTC‑CCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCa                                                                                                                       <  1:466669/137‑1 (MQ=255)
                                                   tACCGTGATTCATTTCCGCCAACAACCGCGC‑ATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAG‑CGaaaaaaacttaaaa                                                                                                    >  1:6814/1‑139 (MQ=25)
                                                                  ccGCCAACACCCGCGT‑TTTTACCCACCCCTTTGTAAAAAATTCCCTCGGCGGCGGTTTTCTCCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTtgtgtg                                                                                     <  2:454333/150‑1 (MQ=255)
                                                                             cgcgC‑TTTTACCCACCCCCGTGTCAAAAAATGCCTCCGCGGCGGCTGTCTCCCAGTTATTTCC‑CCGCCAAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgc                                                                          <  1:53622/150‑1 (MQ=255)
                                                                               cgc‑ATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAG‑AGAAAAAAATTTTAAACCCCCATGTTGTGTGGGGATCgc                                                                             >  2:355375/1‑145 (MQ=25)
                                                                                             ccccATGTAAAAAAAGGCTTTCGGGGCGGCTGTCGGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATGTGGaa                                                         <  2:214872/149‑1 (MQ=255)
                                                                                                   gaacAAAATTGCCTCCGCGGGGTCTG‑CCCCCAGTTTTTTTTCCCGCCAAAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATGTgg                                                           <  1:509728/140‑1 (MQ=255)
                                                                                                        aaaaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAG‑AGAAAAAAATTTTAAACCCCCCTGTTTTGTGGGGATAGCCGC‑TGGGGAAAAAGTGGGAAGGGGAgaggg                                              >  2:37579/1‑147 (MQ=255)
                                                                                                              ctttcgggcgggcTGTCTCCCTCTTTTTTTC‑CCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATGTGGAAAGGGAAGTgcggc                                            <  2:496042/137‑1 (MQ=255)
                                                                                                                   gcggcgGCTGTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGc                                    <  2:251938/150‑1 (MQ=255)
                                                                                                                              cggCCGGCTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc                         <  2:2382/148‑1 (MQ=255)
                                                                                                                                 ccattttttttc‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc                         <  2:92946/140‑1 (MQ=255)
                                                                                                                                      ttttttc‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCt                      <  1:570091/143‑1 (MQ=255)
                                                                                                                                      tATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGTGGGCG‑AGAAATATCTTTAAAACCCCCCTGTTGTGTGGGTATCTCTGCTTGGTACCAA‑GTGGGAAGGGGAGGGAGGGGAGTCAGCGCCCGaaacccccgcatt                  >  1:286973/1‑138 (MQ=255)
                                                                                                                                            t‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGCGGAAGCGGAAAAAAATTAAACCCCCCGATCTT‑TGTGGGAAAAGGCGT‑TGGTGAGAAAGCGGGAAGGGGAGAGCGGGGAATTAACTCCAGAACcgcgggcgct                  >  1:70982/1‑137 (MQ=255)
                                                                                                                                                     aacaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGa  <  2:180415/151‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                      
TGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGC‑ATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGA  >  NC_000913/357741‑358038

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