Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913328,8160TG52.0% 30.1 / 38.0 25K173N (AAA→AAC) betIcholine‑inducible betIBA‑betT divergent operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (7/5);  new base G (13/0);  total (20/5)
Fisher's exact test for biased strand distribution p-value = 1.49e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.51e-04

GTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCA‑GGGAAT‑TAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACG  >  NC_000913/328687‑328965
                                                                                                                                    |                                                                                                                                                     
gTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCC‑GGGGAT‑TAGCGC‑GGGGGTTTTCCAGCGGGTTGCCGCTCCCAGCcccgcg                                                                                                                                   >  1:295890/1‑150 (MQ=255)
   taACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCA‑GGGAAT‑TAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGc                                                                                                                               <  1:26259/151‑1 (MQ=255)
        accaTGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCA‑GGGGAA‑TTGCGC‑GGGGGTTTTTCAGCGGTGTTGCGCCCCGCGCCGCGCgccgccccca                                                                                                                          >  1:417656/1‑147 (MQ=255)
        accaTGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCA‑GGGAAT‑TAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATaa                                                                                                                          <  1:34920/151‑1 (MQ=255)
         ccaTGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCAGGGGAATATAGCG‑‑GGGTTTTTTCCACCGGTTTTCC‑CCCCCAGCCGCGCgccgcccaaaa                                                                                                                         >  2:366395/1‑146 (MQ=255)
           aTGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCA‑GGGAAT‑TAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATcc                                                                                                                       <  2:24938/151‑1 (MQ=255)
            tGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGGC‑GGGGAT‑TAGCGG‑GGGGTTTTTTTCCCGGGTTGCCGCCCCCAGCCccgcgcgggcccaaatc                                                                                                                       >  1:183770/1‑139 (MQ=255)
               atataAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCG‑GGGGAT‑TAGTGG‑GGGGTTTTTTCAACGGGTTGTCGCCCCGAGCCGCGCGCAGCCCTAATccacca                                                                                                                   >  1:245014/1‑151 (MQ=255)
               atataAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCA‑GGGGAT‑TAGCGC‑GGGGTTTTTCCAGCGGGTTGCCGCTCAGAGCCGCGCGCCGCCCTAATCCatca                                                                                                                   >  2:156961/1‑151 (MQ=255)
               atataAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGGA‑GGGGAA‑TTGGGC‑GGGGGTTTTTCCCCGGGTTGCCGCTCAGAGGCGCGCGCCGCCCAAAACCatc                                                                                                                    >  2:51444/1‑150 (MQ=255)
                       ctgctgTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCA‑GGGAAT‑TAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAgcgc                                                                                                           <  1:133340/151‑1 (MQ=255)
                            gTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCA‑GGGGAT‑TAGGGC‑GGGGTTTTTTCAGCGGGTTTTCGCTCCGAGCCGCGCGGAGCCCTAATACATCAATAAACGCGGcc                                                                                                       >  2:139511/1‑150 (MQ=255)
                            gTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCA‑GGGGAA‑TAGCGCGGGGGTTTTTCCA‑CGGGTTGTCGCCCCGAGCCCCGCGCAGGCCCAATCCATCAATCAGCgccgcc                                                                                                       >  1:452455/1‑150 (MQ=255)
                            gTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCA‑GGGGAA‑TAGCGG‑GGGGTTTTTTCAGCGGGTTGCCGCCCCGAGCCCCCCGCAGCCCCAATACATCCAACAGCGCGGcc                                                                                                       >  1:578822/1‑150 (MQ=255)
                                       cGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCA‑GGGGAA‑TAGCGG‑GGGGGTTTTTCAGCGGGTTGCCGCCCCGAGCCCCGCGCGGCCCCAAACCATCAAACAACGcggcccggccccggcc                                                                                            >  1:274856/1‑144 (MQ=255)
                                              tCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCC‑GGGGAT‑TAACGC‑GGGGTTTTTTCAACGGGTTTTCGCCCCCAGCCGCGCGCGGCCCTAATTCATCAATCAGCGGGGCCCGGCCGGCGcccccctc                                                                                      >  2:148404/1‑144 (MQ=255)
                                                                   gTGGGTAGATGCTGAGTGATAAAGTGGCGGGGGA‑GGGGAA‑TAGCGG‑GGGGTTTTTTCAGGGGGTTGGCGCCCCGAGACGCCCGCCGCCCCAAAACCACCATCAACGCGGCCCGGCCGGAGGCCCCCCCCCCTGCCCGTTTCCGCGGccaat                                                               >  1:117180/1‑147 (MQ=255)
                                                                     gggTAGATGCTGAGTGATAAAGTGGCGGGGCA‑GGGGAT‑TAGCGC‑GGGGTTTTTTCCGCGGGTTTGCGCTCCGAGGCGCGGGCGGGCCTAATTCATCAATCAGAGCGGCCCGGCCGGAGCCCGCCTCCCGTGGCCGTTTGCGCGGGAATTc                                                              >  1:517610/1‑150 (MQ=255)
                                                                      ggTAGATGCTGAGTGATAAAGTGGCGGGGCA‑GGGAAT‑TAGAGC‑GGGGGTTATCCAGCGGGTTGCCGCTCCGTGCCGCGGGCGGGCCTAATCCCACCAAAAGGGTGGCCCGGGCCCAGCCCCCCTCCCGTGCCTGGTCGCGCCGCGAATcg                                                             >  1:290310/1‑150 (MQ=255)
                                                                                                ggTCA‑GGGAAT‑TAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGCGACTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCCGCTTCGa                                   >  1:303170/1‑150 (MQ=255)
                                                                                                 gTCA‑GGGAAT‑TAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGAc                                  <  2:100300/150‑1 (MQ=255)
                                                                                                                     ggTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGgctg                   >  2:210576/1‑148 (MQ=255)
                                                                                                                     ggTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGgctg                   >  1:509218/1‑148 (MQ=255)
                                                                                                                      gTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGa               >  2:166758/1‑151 (MQ=255)
                                                                                                                                   gTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACg  >  2:605092/1‑151 (MQ=255)
                                                                                                                                    |                                                                                                                                                     
GTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCA‑GGGAAT‑TAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACG  >  NC_000913/328687‑328965

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: ‑