Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,476,1790CT30.4% 25.4 / 17.8 25D88N (GAT→AAT) yheODNA‑binding transcriptional regulator YheO
Reads supporting (aligned to +/- strand):  ref base C (7/9);  new base T (3/4);  total (11/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.86e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGGTAACCAGATCTTCAACAGAAGAGGCAAAGTTGACGCTTGAACCGACATCCGGGGTTTCTGGCGGCACAAAGGTGCTCATAATCTGCGAGAAGGGAACATCAAGATTCATATTGATGCACAGCAGACCAATTACACGCTGTTCGCGGTTACGAATCGCGATAGTCAGGGACTTCATTAATACGCCGCTTTTGGCGCGAGTAAAGTAGCATTTAGAAACGCTGCTATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAACCAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATCCTGCAAAGAGTGCAAAACGATTTCACAGTGGGAGCCAATAAGCATCGCTAACCCGTCCACCACCGCTTCGTAGGATTTCAGAATATCAAAATCGGTCTGGTCGAAAGGACGTTGATCCAGTAAATCCAACTCACTGGTTTCGTTGGTTAAAAG  >  NC_000913/3475936‑3476428
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                         
gggTAACAAGATCTTCAACAGAAGAGGCAAAGTTGACGCTTGAACCGACATCCGGGGTTTCTGGCGGCACAAAGGTGCTCATAATCTGCGAGAAGGGAACATCAAGATTCATATTGATGCACAGCAGACCAATTACTCGCTGTTCGCGGTTACGAATCGCGATAGTCAGGGACTTCATTAATACGCCGCTTTTGGCGCGAGTAAAGTAGCATTTAGAAACGCTGCTATCCGCTCCCGTCATATTGTgcag                                                                                                                                                                                                                                                     <  1:185430/250‑1 (MQ=255)
              tCAACAGAAGAGGCAAAGTGGACGCTTGAACCGACAGCCGGGGTTTCTGGCGGCACAAAGGTGCTCATAATCTGCGAGAAGGGAACATCAAGATTCATATTGATGCACAGCAGACCAATTACACGCTGTTCGCGGTTACGAATCGCGATAGTCAGGGACTTCATTAATACGCCGCTTTTGGCGCGAGTAAAGTAGCATTTAGAAACGCTGCTATCCGCTCCCGTCATATCGTGCAGCATACGTAgcg                                                                                                                                                                                                                                          <  2:48994/247‑1 (MQ=255)
                  cAGAAGAGGCAAAGTGGACGCTTGAACCGACATCCGGGGTTTCTGGCGGCACAAAGGTGCTCATAATCTGCGAGAATGGAACATCAAGATTCATATTGATGCACAGCAGACCAATTACACGCTGTTCGCGGTTACGAATCGCGATAGTCAGGGACTTCATTAATACGCCGCTTTTGGCGCGAGTAAAGTAGCATTTAGAAACGCTGCTATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCa                                                                                                                                                                                                                                  <  2:168104/251‑1 (MQ=255)
                      agagGCAAAGTTGACGCTTGAACCGACATCCGGGGTTTCTGGCGGCACAAAGGTGCTCATAATCTGCGAGAAGGGAACATCAAGATTCATATTGATGCACAGCAGACCAATTACACGCTGTTCGCGGTTACGAATCGCGATAGTCAGGGACTTCATTAATACGCCGCTTTTGGCGCGAGTAAAGTAGCATTTAGAAACGCTGCTATCCGCTCCCGTCATATTGTGCAGCATACGTAGCGCCAGGTCAGTa                                                                                                                                                                                                                               >  1:248311/1‑250 (MQ=255)
                        agGCAAAGGTGACGCTTGAAGCGACATCCGGCGGTTCTGGCGGCACAAAGGTGCTCATAATCTGCGAGAAGGTAACATCAAGATTCATCGTGATGGACAGCAGACCAATGACACGCTGTTCGCGTTTACGAATCGTGATAGTCATTGAACGCATTCAGATCCCGTTTTTGTGGCGACTAAAGTAGTATTTATAAACGCTGCTGTCCGCTTCCGTCATATGGTgcagctttcggtgcggcattgcattaatt                                                                                                                                                                                                                            <  2:279289/251‑25 (MQ=255)
                                         gAACCGACATCCGGGGTTTCTGGCGGCACAAAGGTGCTCATAATCTGCGAGAAGGGAACATCAAGATTCATATTGATGCACAGCAGACCAATTACACGCTGTTCGCGGTTACGAATCGCGATAGTCAGGGACTTCATTAATACGCCGCTTTTGGCGCGAGTAAAGTAGCATTTAGAAACGCTGCTATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAACCAATCTTCCg                                                                                                                                                                                                           <  1:241532/251‑1 (MQ=255)
                                                     ggggTTTCTGGCGGCACAAAGGTGCTCATAATCTGCGAGAAGGGAACATCAAGATTCATATTGATGCACAGCAGACCAATTACACGCTGTTCGCGGTTACGAATCGCGATAGTCAGGGACTTCATTAATACGCCGCTTTTGGCGCGAGTAAAGTAGCATTTAGAAACGCTGCTATCCGCTCCCGTCATATTGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAACCAATCTTCCGGCCTGTATGTTc                                                                                                                                                                                               <  2:248311/251‑1 (MQ=255)
                                                     ggggTTTCTGGCGACAGCAAGGTGCTCATAATCTGCGAGCCGGGCCCATCACGATGCCTATTGATGCACAGCAGACCAATTACACGCTGATCGCGGCTACGAATCGCGATAGTCAGGGACTTCATTAATACGCCGCTTTTGGCGCGAGTAACGTAGCCTTTCGAAACGCTGCTATACGCTCACGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAACCCATCTTCCGGCCTGTATGTTc                                                                                                                                                                                               >  2:183260/1‑251 (MQ=255)
                                                                               cATAATCTGCTAGAAGGGAACATCAAGATTCATATTGATGCACAGCAGACCAATTACACGCTGTTCGCGGTTACGAATCGCGATAGTCAGGGACTTCATTAATACGTCGCTTTTGGCGCGAGTAAAGTAGCATTTAGAAACGCTGCTATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAACCAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAAc                                                                                                                                                                     <  1:76676/251‑1 (MQ=255)
                                                                                         gagaAGGGAACATCAAGATTCATATTGATGCACAGCAGACCAATTACACGCTGTTCGCGGTTACGAATCGCGATAGTCAGGGACTTCATTAATACGCCGCTTTTGGCGCGAGTAAAGTAGCATTTAGAAACGCTGCTATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAACCAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACAtttt                                                                                                                                                                >  1:126217/1‑246 (MQ=255)
                                                                                          agaAGGGAAAATCCAGATTCATATTGATGCACAGCAGACCAATTACACGCTGTTCGCGGTTACGAATCGCGATAGTCAGGGACTTCATTAATACGCCGCTTTTGGCGCGAGTAAAGTAGCATTTAGAAACGCTGCTATCCGCTCCCGTCATATTGTGCAGCATACGTAGCGCCCGGTCCGAAATTGGCGAACCAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATcc                                                                                                                                                          >  2:291149/1‑251 (MQ=255)
                                                                                                    aTCAAGATTCATATTGATGCACAGCAGACCAATTACACGCTGTTCGCGGTTACGAATCGCGATAGTCAGGGACTTCATTAATACGCCGCTTTTGGCGCGAGTAAAGTAGCATTTAGAAACGCTGCTATCCGCTCCCGGCATATCGTGCAGCATACGTAGCGCCAGGGCAGTAATTGGCGAAGCAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATCCTGCAAAGAGt                                                                                                                                                <  2:147615/251‑1 (MQ=255)
                                                                                                            tCATATTGATGCACAGCAGACCAATTACACGCTGTTCGCGGTTACGACTCGCGATAGTCAGGGACTTCATTAATACGCCGCTTTTGGCGCGAGTAAAGTAGCATTTAGAAACGCTGCTATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAACCACTCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGACCATTTTAGTTCCTGCAAAGAGTGCAAAACg                                                                                                                                        >  2:24889/1‑251 (MQ=255)
                                                                                                             cATATTGATGCACAGCAGCCCAATTACACGCTGTTCGCGGTTACGAATCGCGATAGTCAGGGACTTCATTAATACGCCGCTTTTGGCGCGAGTAAAGTAGCATTTAGAAACGCTGCTATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAACCAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATCCTGCAAAGAGTGCAAAACg                                                                                                                                        >  1:115393/1‑250 (MQ=255)
                                                                                                                   gATGCACAGCAGACCAATTACACGCTGTTCGCGGTTACGAATCGCGATAGTCAGGGACTTCATTAATTCGCCGCTTTTGGCGCGAGTAAAGTAGCATTGAGAAACGCTGCTATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAACCAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATCCTGCAAAGAGTGCAAAACGATTTcac                                                                                                                                 <  1:183260/251‑1 (MQ=255)
                                                                                                                                  aaTTACACGCTGTTCGCGGTTACGAATCGCGATAGTCAGGGACTTCATTAATACGCCGCTTTTGGCGCGAGTAAAGTAGCATTTAGAAACGCTGCTATCCGCTCCCGTCATATTGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAACCAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATCCTGCAAAGAGTGCAAAACGATTTCACAGTGGGAGCCAATa                                                                                                                   >  2:69442/1‑250 (MQ=255)
                                                                                                                                             gTTCGCGGTTACGAATCGCGATAGTCAGGGACTTCATTAATACGCCGCTTTTGGCGCGAGTAAAGTGGCATTTAGAAACGCTGCTATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAACCAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATCCTGCAAAGAGTGCAAAACGATTTCACAGTGGGAGCCAATAAGCGTCGCTAAc                                                                                                       <  2:35486/251‑1 (MQ=255)
                                                                                                                                                       aCGAATCGCGATAGTCAGGAACTTCATTAATACGCCGCTTTTGGCGCGAGTAAAGTAGCATTTAGAAACGCTGCTATCCGCTCACGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAACCAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGACTGGCTGAACATTTTAGATCCTGCAAAGAGTGCAAAACGATTTCACAGTGGGAGCCACTAAGCATCGCTAACCCGTccaca                                                                                              >  2:273124/1‑249 (MQ=255)
                                                                                                                                                                   aGTCAGGGACTTCAGTAATACGCCGCTTTTGGCGCGAGTAAAGTAGCATTTAGAAACGCTGCTATCCGCTCCCGTCATATTGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAACCAATCTTCCGGCCTGTATGGTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATCCTGCAAAGAGTGCAAAACGATTGCACAGTGGGAGCCAATAAGCATCGCTAACCCGTCCACCACCGCTGCGTAgg                                                                                 <  2:14314/251‑1 (MQ=255)
                                                                                                                                                                       aGGGACTGCATTAATACGCCGCTTTTGGCGCGAGTAAAGTAGCATTTAGAAACGCTGCTATCCGCTCCCGTCATATTGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAACCAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATCCTGCAAAGAGTGCAAAACGATTTCACAGTGGGAGCCAATAAGCATCGCTAACCCGTCCACCACCGCTTCGTAGGAttt                                                                             <  1:69442/251‑1 (MQ=255)
                                                                                                                                                                                        gccgcTTTTGGCGCGAGTAAAGTAGCATTTAGAAACGCTGCTATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAACCAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATCCTGCAAAGAGTGCAAAACGATTTCACAGTGGGAGCCAATAAGCATCGCTAACCCGTCCACCACCGCTTCGTAGGATTTCAGAATATCAAAATCgg                                                            >  2:28436/1‑251 (MQ=255)
                                                                                                                                                                                                       aGTAAAGTAGCATTTAGAAACGCTGCGATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAACCAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATCCTGCAAAGAGTGCAAAACGATTTCACAGTGGGAGCCAATAAGCATCGCTAACCCGTCCACCACCGCTTCGTAGGATTTCAGAATATCAAAATCGGTCTGGTCGAAAGGAc                                             <  2:160568/251‑1 (MQ=255)
                                                                                                                                                                                                          aaaGTAGCATTTAGAAACGCTGCTATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAACCAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATCCTGCAAAGAGTGCAAAACGATTTCACAGTGGGAGCCAATAAGCATCGCTAACCCGTCCACCAACGCTTCGTAGGATTTCAGAATATCAAAATCGGTCTGGTCGAAAGGACGtt                                          <  2:256333/251‑1 (MQ=255)
                                                                                                                                                                                                                                      gCTCCCGTAATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAACCAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATCCTGCAAAGAGTGCAAAACGATTTCACAGTGGGAGCCAATAAGCATCGCTAACCCGTCCACCACCGCTTCGTAGGATTTCAGAATATCAAAATCGGTCTGGTCGAAAGGACGTTGATCCAGTAAATCCAACTCACTGGTTTc              >  2:117856/1‑251 (MQ=255)
                                                                                                                                                                                                                                        tCCCGTCCTATAGCGCAGCATACGTCGCGCCAGGTCAGTAATTGGCGAACCAATCTTCCGACCTGTATGTTCTCCGTTAGCAATGCGAATGGCTGAACATTTTAGCTCCTGCAAAGAGTGCAAAACGATTTCACAGTGGGAGCCAATAAGCATCGCTAACCCGTCCACCACCGCTTCGTAGGATTTCAGAATATCAAAATCGGTCTGGTCGAAAGGACGTTGATCCAGTAAATCCACCTCACTGGTTTCGt            >  2:308607/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   cGTGCAGCAAACGTAGCGCCAGGTCAGTAATTGGCGAACCAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATCCTGCAAAGAGTGCAAAACGATTTCACAGTGGGAGCCAATAAGCATCGCTAACCCGTCCACCACCGCTTCGTAGGATTTCAGAATATCAAAATCGGTCTGGTCGAAAGGACGTTGATCCAGTAAATCCAACTCACTGGTTTCGTTGGTTAAAAg  <  1:19841/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                         
GGGTAACCAGATCTTCAACAGAAGAGGCAAAGTTGACGCTTGAACCGACATCCGGGGTTTCTGGCGGCACAAAGGTGCTCATAATCTGCGAGAAGGGAACATCAAGATTCATATTGATGCACAGCAGACCAATTACACGCTGTTCGCGGTTACGAATCGCGATAGTCAGGGACTTCATTAATACGCCGCTTTTGGCGCGAGTAAAGTAGCATTTAGAAACGCTGCTATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAACCAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATCCTGCAAAGAGTGCAAAACGATTTCACAGTGGGAGCCAATAAGCATCGCTAACCCGTCCACCACCGCTTCGTAGGATTTCAGAATATCAAAATCGGTCTGGTCGAAAGGACGTTGATCCAGTAAATCCAACTCACTGGTTTCGTTGGTTAAAAG  >  NC_000913/3475936‑3476428

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