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→AACbetIcholine‑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

GTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACG  >  NC_000913/328687‑328965
                                                                                                                                    |                                                                                                                                                     
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                   >  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)
                                                                                                                                    |                                                                                                                                                     
GTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACG  >  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: