Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913328,8030TG56.0% ‑2.3 / 29.3 25T178P (ACC→CCC) betIcholine‑inducible betIBA‑betT divergent operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (5/6);  new base G (14/0);  total (19/6)
Fisher's exact test for biased strand distribution p-value = 2.61e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.56e-02

CTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGA  >  NC_000913/328659‑328952
                                                                                                                                                   |                                                                                                                                                     
ctcGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATc                                                                                                                                                 <  1:441316/151‑1 (MQ=255)
                            gTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCCGGGGATTAGCGCGGGGGTTTTCCAGCGGGTTGCCGCTCCCAGCcccgcg                                                                                                                      >  1:295890/1‑150 (MQ=255)
                               taACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGc                                                                                                                  <  1:26259/151‑1 (MQ=255)
                                    accaTGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGGAATTGCGCGGGGGTTTTTCAGCGGTGTTGCGCCCCGCGCCGCGCgccgccccca                                                                                                             >  1:417656/1‑147 (MQ=255)
                                    accaTGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATaa                                                                                                             <  1:34920/151‑1 (MQ=255)
                                     ccaTGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCAGGGGAATATAGCGGGGTTTTTTCCACCGGTTTTCCCCCCCAGCCGCGCgccgcccaaaa                                                                                                            >  2:366395/1‑146 (MQ=255)
                                       aTGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATcc                                                                                                          <  2:24938/151‑1 (MQ=255)
                                        tGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGGCGGGGATTAGCGGGGGGTTTTTTTCCCGGGTTGCCGCCCCCAGCCccgcgcgggcccaaatc                                                                                                          >  1:183770/1‑139 (MQ=255)
                                           atataAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCGGGGGATTAGTGGGGGGTTTTTTCAACGGGTTGTCGCCCCGAGCCGCGCGCAGCCCTAATccacca                                                                                                      >  1:245014/1‑151 (MQ=255)
                                           atataAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGGATTAGCGCGGGGTTTTTCCAGCGGGTTGCCGCTCAGAGCCGCGCGCCGCCCTAATCCatca                                                                                                      >  2:156961/1‑151 (MQ=255)
                                           atataAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGGAGGGGAATTGGGCGGGGGTTTTTCCCCGGGTTGCCGCTCAGAGGCGCGCGCCGCCCAAAACCatc                                                                                                       >  2:51444/1‑150 (MQ=255)
                                                   ctgctgTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAgcgc                                                                                              <  1:133340/151‑1 (MQ=255)
                                                        gTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGGATTAGGGCGGGGTTTTTTCAGCGGGTTTTCGCTCCGAGCCGCGCGGAGCCCTAATACATCAATAAACGCGGcc                                                                                          >  2:139511/1‑150 (MQ=255)
                                                        gTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGGAATAGCGCGGGGGTTTTTCCACGGGTTGTCGCCCCGAGCCCCGCGCAGGCCCAATCCATCAATCAGCgccgcc                                                                                          >  1:452455/1‑150 (MQ=255)
                                                        gTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCAGGGGAATAGCGGGGGGTTTTTTCAGCGGGTTGCCGCCCCGAGCCCCCCGCAGCCCCAATACATCCAACAGCGCGGcc                                                                                          >  1:578822/1‑150 (MQ=255)
                                                                   cGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCAGGGGAATAGCGGGGGGGTTTTTCAGCGGGTTGCCGCCCCGAGCCCCGCGCGGCCCCAAACCATCAAACAACGcggcccggccccggcc                                                                               >  1:274856/1‑144 (MQ=255)
                                                                          tCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCCGGGGATTAACGCGGGGTTTTTTCAACGGGTTTTCGCCCCCAGCCGCGCGCGGCCCTAATTCATCAATCAGCGGGGCCCGGCCGGCGcccccctc                                                                         >  2:148404/1‑144 (MQ=255)
                                                                                               gTGGGTAGATGCTGAGTGATAAAGTGGCGGGGGAGGGGAATAGCGGGGGGTTTTTTCAGGGGGTTGGCGCCCCGAGACGCCCGCCGCCCCAAAACCACCATCAACGCGGCCCGGCCGGAGGCCCCCCCCCCTGCCCGTTTCCGCGGccaat                                                  >  1:117180/1‑147 (MQ=255)
                                                                                                 gggTAGATGCTGAGTGATAAAGTGGCGGGGCAGGGGATTAGCGCGGGGTTTTTTCCGCGGGTTTGCGCTCCGAGGCGCGGGCGGGCCTAATTCATCAATCAGAGCGGCCCGGCCGGAGCCCGCCTCCCGTGGCCGTTTGCGCGGGAATTc                                                 >  1:517610/1‑150 (MQ=255)
                                                                                                  ggTAGATGCTGAGTGATAAAGTGGCGGGGCAGGGAATTAGAGCGGGGGTTATCCAGCGGGTTGCCGCTCCGTGCCGCGGGCGGGCCTAATCCCACCAAAAGGGTGGCCCGGGCCCAGCCCCCCTCCCGTGCCTGGTCGCGCCGCGAATcg                                                >  1:290310/1‑150 (MQ=255)
                                                                                                                            ggTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGCGACTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCCGCTTCGa                      >  1:303170/1‑150 (MQ=255)
                                                                                                                             gTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGAc                     <  2:100300/150‑1 (MQ=255)
                                                                                                                                                 ggTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGgctg      >  1:509218/1‑148 (MQ=255)
                                                                                                                                                 ggTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGgctg      >  2:210576/1‑148 (MQ=255)
                                                                                                                                                  gTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGa  >  2:166758/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGA  >  NC_000913/328659‑328952

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