Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	328,791	0	A	G	44.0%	58.8 / 19.3	50	S182P (TCC→CCC)	betI	choline‑inducible betIBA‑betT divergent operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (4/24);  new base G (22/0);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 1.72e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAATTAG‑CGC‑GGGTTTTATCCAGC‑GGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGC  >  NC_000913/328645‑328940                                                                                                                                                    |                                                                                                                                                         gCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAATTa                                                                                                                                                      <  1:228935/151‑1 (MQ=255)     ggTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGGCCGGGGAGTAg‑cg                                                                                                                                                  >  1:297899/1‑150 (MQ=255)               ctcGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGGCCGGGGAATAG‑CGC‑GGGgtt                                                                                                                                          >  2:773548/1‑147 (MQ=255)                 cGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGGCCGGGGAATAG‑TGC‑GGGgttt                                                                                                                                         >  2:254451/1‑146 (MQ=255)                  gAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGGCAGGGGGTTTG‑GGC‑gggggtttttca                                                                                                                                    >  2:71037/1‑141 (MQ=255)                      gTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGGGCGGGGATTAGAGGG‑GGGTTTTTTCCacgg                                                                                                                                 >  2:20958/1‑147 (MQ=255)                       tGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGGAGTAG‑GGC‑gggggtttttc                                                                                                                                     >  1:70613/1‑136 (MQ=255)                       tGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTAAGTGATAAAGTGGC‑GGGTCGGGGGATGAG‑GGC‑gggggtttttccccg                                                                                                                                 >  1:215721/1‑136 (MQ=255)                         cGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAATTAG‑CGC‑GGGTTTTATCCAGC‑GGTTTg                                                                                                                           <  1:532190/151‑1 (MQ=255)                            ccGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAATTAG‑CGC‑GGGTTTTATCCAGC‑GGTTTGcc                                                                                                                         <  1:903873/150‑1 (MQ=255)                             cGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCCGGGGAATTGG‑GGG‑GGGTTTTTTctgcgggtt                                                                                                                              >  1:256362/1‑138 (MQ=255)                             cGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAATTAG‑CGC‑GGGTTTTATCCAGC‑GGTTTGCCGc                                                                                                                       <  2:649903/151‑1 (MQ=255)                                   tggCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTTAGGGGAGTAG‑GGC‑GGGGGTTTTTCAGC‑GGGTtgctgcgccc                                                                                                                   >  1:950922/1‑145 (MQ=255)                                     gCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAATTAG‑CGC‑GGGTTTTATCCAGC‑GGTTTGCCGCTCAGAGcc                                                                                                               <  1:90618/151‑1 (MQ=255)                                          ggTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAATTAG‑CGC‑GGGTTTTATCCAGC‑GGTTTGCCGCTCAGAGCcgcgcg                                                                                                          <  1:70657/151‑1 (MQ=255)                                              taACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAATTAG‑CGC‑GGGTTTTATCCAGC‑GGTTTGCCGCTCAGAGCCGCGCGCAg                                                                                                       <  1:18699/150‑1 (MQ=255)                                                  caccaTGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGGCAGGGGAGTAG‑CGC‑GGGGGTTTTTTCGC‑CGGTTGTCGCTCCGAGCcgcccggcgcc                                                                                                     >  1:831251/1‑143 (MQ=255)                                                     caTGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGGCAGGGGAATAGAGGC‑GGGTTTT‑TTCAAC‑GGGTTGCCGCTCCGAAGCCCGCGCGGCCCTaaa                                                                                                >  2:679917/1‑149 (MQ=255)                                                        gTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAATTAG‑CGC‑GGGTTTTATCCAGC‑GGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCat                                                                                            <  2:788420/151‑1 (MQ=255)                                                        gTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGGCAGGGGAAGAT‑GGC‑GGGGGTTTTTTAAC‑GGGTTGGCGCGCACAGCCCCGCGCCGCCaaaa                                                                                                 >  2:441429/1‑146 (MQ=255)                                                            ataAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGGATTAG‑CGC‑GGGGGTTTTTCAAC‑GGGTTGCCGCCCCGAGCCGCGCGCAGCCCTAATCCctcaa                                                                                         >  2:351423/1‑150 (MQ=255)                                                             taAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGGATTAG‑CGC‑GGGGTTTTTTCAGC‑GGGTTTGCGCCCCGAGGCGCGCGCCGCCATAATCCatcaa                                                                                         >  2:555310/1‑149 (MQ=255)                                                               aaGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAATTAG‑CGC‑GGGTTTTATCCAGC‑GGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCAACAATCAg                                                                                     <  2:408125/151‑1 (MQ=255)                                                                 gctgctGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGGCAGGGGGATAG‑AGG‑GGGGTTTTTTCACC‑GGGTTTCCGCCCCAAGCCCCCCgccgccccaa                                                                                                 >  1:1111858/1‑133 (MQ=255)                                                                    gctgTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAATTAG‑CGC‑GGGTTTTATCCAGC‑GGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAgcgc                                                                                  <  2:569666/149‑1 (MQ=255)                                                                    gctgTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAATTAG‑CGC‑GGGTTTTATCCAGC‑GGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGcgcg                                                                                 <  2:886789/150‑1 (MQ=255)                                                                          cTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGGCAGGGGAGTAA‑CGC‑GGGGGTTATCCAAC‑GGGTTGGCGCTCCTAGCCGCCCGCCGCCACAATCAATCCAACAccgc                                                                                  >  1:728487/1‑143 (MQ=255)                                                                            gCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCCGGGGAATAA‑CGCGGGGGTTTTTCAACC‑GGTTT‑CCGCTCAGAGCCGCGCGCCGCCCTAATCCATCAATCAACGCgg                                                                                >  2:867384/1‑143 (MQ=255)                                                                            gCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAATTAG‑CGC‑GGGTTTTATCCAGC‑GGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAgg                                                                           <  2:864455/148‑1 (MQ=255)                                                                                   gggACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGGCAAGGGATTAG‑CGC‑GGGGTTTATTCAAC‑CGCTTTGCGCCCCAAGACGCGCGCCGCCCCAAACCAAAAAACAGCgggggccagccgcaacccc                                                                 >  2:947273/1‑136 (MQ=255)                                                                                    ggACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAATTAG‑CGC‑GGGTTTTATCCAGC‑GGTTTGCCGCGCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGc                                                                <  1:805554/151‑1 (MQ=255)                                                                                      aCATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGCGGC‑GGGGCAGGGGATTAG‑CGC‑GGGGTTTATTCAGC‑GGGTTTCCGCTCAGAGCCGCGCGCCGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGcct                                                              >  2:375588/1‑149 (MQ=255)                                                                                          cgcGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAATTAG‑CGC‑GGGTTTTATCCAGC‑GGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCt                                                           <  1:589167/148‑1 (MQ=255)                                                                                                 cctcGTTTAATCGTTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGATTTAG‑CGC‑GGGTTTTATCCAGC‑GGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCt                                                     <  1:946265/149‑1 (MQ=255)                                                                                                 cctcGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAATTAG‑CGC‑GGGTTTTATCCAGC‑GGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGt                                                   <  2:597456/151‑1 (MQ=255)                                                                                                     ggTTAATCGGTGGGTAGATGCTAAGTGATAAAGGGGC‑GGTTCAGGGAATTAG‑CGC‑GGGTTTTATCCAGC‑GTTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTcg                                                <  1:1090911/150‑1 (MQ=255)                                                                                                        tAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGGGCGGGGATTAG‑CGC‑GGGGTTTTTTCAGC‑GGGTTGTCCCTCAGAGCCCCGCGCCGCCATAATCCCCTAAACAACGCGGCCCCGCCCCAGCCCGCCTCCTCTGcccg                                                    >  2:446728/1‑141 (MQ=255)                                                                                                         aaTCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGGGTTAG‑CGC‑GGGGTTTTTTCAGCGGGTTTGCC‑CCCCGAGCCGCGCGCCGCCCTAAACCATCAAACAACGCGggccggccccaagcccccttc                                                            >  2:354593/1‑121 (MQ=255)                                                                                                         aaTCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAATTAG‑CGC‑GGGTTTTATCCAGC‑GGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCgg                                           <  2:1100731/151‑1 (MQ=255)                                                                                                               tGGGTAGATGCTGAGTGATAAAGTGGC‑GGGGCAGGGGAATAA‑TGC‑GGGGTTTTTTCAGC‑GGGTTGTCGCCCACAGGCGCGCGCCGCCACAATCCCTCAAACACCGCGGCCAGGCCGGAGCCCCCCTCCCTTGCCCGTTTCCGCGGcgaatc                                     >  1:151670/1‑146 (MQ=255)                                                                                                                   tAGATGCTGAGTGATAAAGTGGC‑GGGGCCGGGGAATAA‑GGC‑GGGGGTTTTTCCGC‑GGGTTGGCGCTCAGAGCCGCGCGGACGCATAAACCCACAATCAGCGCGGGCCGGCCCGCGCCCGCCTCCCTTGCCTGTTCGCGCgggcaat                                      >  1:956539/1‑141 (MQ=255)                                                                                                                     gATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAATTAG‑CGC‑GGGTTTTATCCAGC‑GGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGcgacg                               <  2:118741/151‑1 (MQ=255)                                                                                                                           gAGTGATAAAGTGGC‑GGGTCAGGGAATTAG‑CGC‑GGGTTTTATCCAGC‑GGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTc                         <  2:70613/151‑1 (MQ=255)                                                                                                                            aGTGATAAAGTGGC‑GGGTCAGGGAATTAG‑CGC‑GGGTTTTATCCAGC‑GGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCg                        <  2:215725/151‑1 (MQ=255)                                                                                                                                   aaGTGGC‑GGGGCAGGGAAATAG‑AGC‑GGGGTTTTTTCAGC‑GGGTTTCCGCCCAGAGCCGCGCGCAGCCCTAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCCGTGCCTGTTCGCGCGGCCATTCGCGGCGAAACTCCCTcccc                  >  1:832716/1‑150 (MQ=255)                                                                                                                                     gTGGC‑GGGTCAGGGAATTAG‑CGC‑GGGTTTTATCCAGC‑GGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGa               <  1:254447/151‑1 (MQ=255)                                                                                                                                       ggC‑GGGTCAGGGAATTAG‑CGC‑GGGTTTTATCCAGC‑GGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGAt              <  2:102536/150‑1 (MQ=255)                                                                                                                                            ggTCAGGGAATTAG‑CGC‑GGGTTTTATCCAGC‑GGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGAc         >  1:442588/1‑151 (MQ=255)                                                                                                                                                 gggAATTAG‑CGC‑GGGTTTTATCCAGC‑GGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCaa    <  1:1028645/151‑1 (MQ=255)                                                                                                                                                   gAATTAG‑CGC‑GGGTTTTATCCAGC‑GGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAgc  >  1:589950/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                         GCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAATTAG‑CGC‑GGGTTTTATCCAGC‑GGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGC  >  NC_000913/328645‑328940

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