Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913328,8030TG50.0% 30.5 / 22.8 36T178P (ACC→CCC) betIcholine‑inducible betIBA‑betT divergent operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (4/14);  new base G (18/0);  total (22/14)
Fisher's exact test for biased strand distribution p-value = 1.61e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTG  >  NC_000913/328657‑328951
                                                                                                                                                   |                                                                                                                                                     
gTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTa                                                                                                                                                   <  2:465822/151‑1 (MQ=255)
     gAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAg                                                                                                                                              <  1:543177/151‑1 (MQ=255)
        ggTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGGAGGGGGTTAGAGCGGGGTTTTTTCAGCgg                                                                                                                                           >  2:658913/1‑151 (MQ=255)
             gACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGGAATAGCGCGGGGTTTTTTCAGCGGgtt                                                                                                                                        >  2:33470/1‑149 (MQ=255)
             gACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGc                                                                                                                                      <  2:341664/151‑1 (MQ=255)
                  cTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCAGGGGAATAAGGCGGGGTTTTTTCAACGGGTTTcc                                                                                                                                     >  2:316390/1‑147 (MQ=255)
                    ggtggCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCa                                                                                                                                <  2:233143/150‑1 (MQ=255)
                      tggCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGGATTAGGGCGGGGTTTTTTCAGCGGTTTg                                                                                                                                       >  2:130845/1‑141 (MQ=255)
                        gCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCAGGGGAGTCGCGCGGGGTTTTTCGCGCCGGGTGCCGCCCag                                                                                                                               >  2:590653/1‑147 (MQ=255)
                                     caccaTGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGCGCAGGGGAGTTGCGCGGGGGTTTTTCAGCGGTTTgttgcctacgcccgc                                                                                                                         >  1:254501/1‑126 (MQ=255)
                                         aTGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTGAGGGGAATAGGGGGGGGGTTTTTCACCGGGTTGCCGCCCCGAGGCCCGCGCCGCCCTaa                                                                                                             >  2:411688/1‑148 (MQ=255)
                                         aTGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCCGGGGCATAGCGGGGGGGTTATCCAGCCGTGTGCCCCGCCCCGCCGCGCGCAGCCCTAATcc                                                                                                          >  1:129848/1‑151 (MQ=255)
                                          tGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCGGGGGGGTAACGCGGGGGTTTTTCACACGGTTTCCGCCCCGAGCCCCGCgccgcccaaaatc                                                                                                          >  2:511118/1‑143 (MQ=255)
                                           gTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCa                                                                                                         <  2:135475/150‑1 (MQ=255)
                                               ataAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGGAGGGGAATAAAGCGGGGGTTTTTCACCGGGTTTGCGCCCCCAGCCGCGCgccgccccaa                                                                                                             >  2:651940/1‑138 (MQ=255)
                                                                    tCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCAGGGGAATAACGCGGGGTTTTTTCTGCGGTTTTGCCCTCCGAGCCGCGCGCCGCCCCAATCCATCAAACAACGCGGCCCGGCCgcacc                                                                                >  1:65583/1‑146 (MQ=255)
                                                                     cGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGccc                                                                              <  2:631820/151‑1 (MQ=255)
                                                                       ggACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTTATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGc                                                                            <  1:340859/151‑1 (MQ=255)
                                                                        gACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCCGGGGGGTAACGCGGGGGTTTTTTAACCGGTTTTCCCTCCGAGCCGCGCGCCGCCCCAATCCATCAATAAACGCGGGCAGGCCggcaccccc                                                                            >  2:699214/1‑142 (MQ=255)
                                                                                  ctcctcGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCAGGGGAATAGAGCGGGGTTTTTTCAACGGGTTGTCGCCCAGAGCCGCGCGCAGCCCCAAAACATCAAACAACGCCGGCCGGCCGCAGCCCCCCTCCCCTGcc                                                                  >  2:592619/1‑150 (MQ=255)
                                                                                        ggTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTGAGGGGATTTGAGCGGGGGTTTTTCAGCCGTTTGGCGCCCAGAGCCGCGCCCCGCCCTCAAACATCCCTACACGCGGCCCGGCCGGCGCCCGCGTCCCCTGCCTGTTct                                                            >  1:40403/1‑149 (MQ=255)
                                                                                        ggTTAATCGGGGGGTAGAGGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGtt                                                              <  1:171875/148‑1 (MQ=255)
                                                                                           tAATCGGGGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGTTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTcgcgc                                                         <  1:56346/150‑1 (MQ=255)
                                                                                              tCGGAGGGTAGATGCTGAGTGATAAAGTGGCGGGGGATGGGAATAGGGGGGGGGTTTTCTAGCGGGTTGTCGCGCCGAGCCGCGCGCCGCCCCAATCCCTCCATCAACGGGGGCAGGGCCCAGCCCGCGTCCCGCGCCCGTTGGCCCgg                                                       >  1:48422/1‑149 (MQ=255)
                                                                                                 gTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAAtt                                                 <  2:410322/151‑1 (MQ=255)
                                                                                                      tAGATGCTGAGTGATAAAGTGGCGGGGCCGGGGAATAACGTGGGGGTTTTTCCGCGGGTTGGCGCTCCGAGACCCCCGCCGCCACAAACCATCAACAAACGCGGCCCGGCCGGTACCCGCCTCCTGTGCCTGTTCTCGCGGcgactc                                                >  1:457276/1‑142 (MQ=255)
                                                                                                       aGATGCTGAGTGATAAAGTGGCGGGTCAGGGGAATAACGCGGGGTTTTTTCCTCGGGTTTTCGCCCCGAACCGCGCGCCGCCCAAAACCCTCAATCAACGCGGCCCGGCCGTAGCCCCCCTTCTTTGCCTGTTTGCGCGGCCAATTcc                                               >  2:342603/1‑147 (MQ=255)
                                                                                                         aTGCTGAGTGATAAAGTGGCGGGGCCGGGAAATAGCGCGGGGGTTTTTCAGCGGTTTTCCCCCCACAGCCGCGCGCCGCCCCAAACCATCCATCAACGCGGCCCGGCCGGAGCCCCCCTCCTCTGCCTGTTTGCGCCGCCATTCGcgac                                           >  2:495473/1‑149 (MQ=255)
                                                                                                          tGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTGGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGcgacga                                         <  1:140240/150‑1 (MQ=255)
                                                                                                                        gTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGa                          <  2:351936/151‑1 (MQ=255)
                                                                                                                          ggCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGAtt                        <  1:447751/151‑1 (MQ=255)
                                                                                                                             gggTCGGGGAATTAGCGCGGTTTTTATCCAGCGTTTTGCCCCTCAGAGCCGCGCCCAGCCATAATCCATCAAGCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTc                       <  1:27659/149‑1 (MQ=255)
                                                                                                                              ggTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGAc                    >  1:545237/1‑151 (MQ=255)
                                                                                                                                    ggAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCa                >  2:510630/1‑149 (MQ=255)
                                                                                                                                             gcgGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGgctg     >  1:507395/1‑151 (MQ=255)
                                                                                                                                                  gTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGctgctg  >  1:230905/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTG  >  NC_000913/328657‑328951

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