Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	328,803	0	T	G	42.5%	65.5 / 30.4	47	T178P (ACC→CCC)	betI	choline‑inducible betIBA‑betT divergent operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (3/24);  new base G (20/0);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 1.81e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTC‑AGGGAATTAGCGCGGGTTTTATCCAGCGG‑TTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGAC  >  NC_000913/328661‑328953                                                                                                                                                |                                                                                                                                                        cGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTC‑AGGGGAATAGAGCGGGgttttt                                                                                                                                                    >  2:878467/1‑146 (MQ=255)  gAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGC‑AGGGGAGTAACGCGGGgttttt                                                                                                                                                    >  2:940423/1‑145 (MQ=255)      gTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTC‑AGGGAATTAGCGCGGGTTTTATCCAGCGG‑t                                                                                                                                           <  1:107812/151‑1 (MQ=255)        gCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTC‑AGGGAATTAGCGCGGGTTTTATCCAGCGG‑tt                                                                                                                                          <  2:586486/150‑1 (MQ=255)            ccGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTC‑AGGGAATTAGCGCGGGTTTTATCCAGCGG‑TTTGCCg                                                                                                                                     <  2:542048/151‑1 (MQ=255)             cGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTC‑AGGGAATTAGCGCGGGTTTTATCCAGCGG‑TTTGCCGc                                                                                                                                    <  2:302848/151‑1 (MQ=255)              gCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTC‑AGGGAATTAGCGCGGGTTTTATCCAGCGG‑TTTGCCGCt                                                                                                                                   <  2:837378/151‑1 (MQ=255)                 ggtggCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTC‑CGGGGATTTGAGAGGGGTTTTTTCAGGGG‑TTTTcc                                                                                                                                      >  2:383079/1‑145 (MQ=255)                 ggtggCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTC‑AGGGAATTAGCGCGGGTTTTATCCAGCGG‑TTTGCCGCTc                                                                                                                                  <  2:742790/149‑1 (MQ=255)                 ggtggCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTC‑AGGGAATTAGCGCGGGTTTTATCCAGCGG‑TTTGCCGCTCa                                                                                                                                 <  2:332896/150‑1 (MQ=255)                   tggCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGC‑AGGGGATTAGCGCGGGGTTTTTTCAACGG‑GTTGCCGCTCa                                                                                                                                 >  2:582965/1‑148 (MQ=255)                       ggAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTC‑AGGGAATTAGCGCGGGTTTTATCCAGCGG‑TTTGCCGCTCAGAGCcgc                                                                                                                          <  2:556909/151‑1 (MQ=255)                             ataACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGC‑TGGGGATTAGCGCGGGGTTTTTTCAGCGG‑GTTGTCGCCCACAGAcgcggcc                                                                                                                      >  1:784949/1‑146 (MQ=255)                                           tataAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGT‑GGGGGAGTAGCGCGGGGGTTTTCAAGCGGCTTTGCCGC‑CCCAGGCGCGCGCCGCCACAAAccccc                                                                                                        >  1:533428/1‑146 (MQ=255)                                              aaaGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTC‑AGGGAATTAGCGCGGGTTTTATCCAGCGG‑TTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCatcaatca                                                                                                   <  2:170779/151‑1 (MQ=255)                                                aGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTC‑AGGGAATTAGCGCGGGTTTTATCCAGCGG‑TTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAgc                                                                                                 <  1:233420/151‑1 (MQ=255)                                                  ctgctgTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTG‑AGGGAAATAAAGCGGGGGTTTTTCAGCGG‑TTTGGCGCTCCGAGGCGCCCGCCGCCCCAATCCaacaatcacc                                                                                                 >  1:553304/1‑147 (MQ=255)                                                  ctgatgTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGG‑AGGGGGTTAGCGGGGGGTTTTTTCAGCGG‑GTTGCCGCTCCGAGCCGCGCGCCGGCATAATACctcaatcaac                                                                                                 >  1:591543/1‑147 (MQ=255)                                                     ctgTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGC‑CGGGGATTAGCGCGGGGTTTTTTCAACGG‑GTTTGCGCCCCGAACCGCGCGCCGGCATAAACCATCAAACAacgccg                                                                                             >  2:889491/1‑148 (MQ=255)                                                       gTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCGGGGGAATAAGGAGGGGTTTTTACAAGGGT‑TTGGGCGC‑CAGAGCCGGGCGCTGCCCTAAGCCAAAAATAAgcg                                                                                                >  2:392842/1‑145 (MQ=255)                                                       gTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTC‑AGGGAATTAGCGCGGGTTTTATCCAGCGG‑TTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCa                                                                                          <  2:160079/151‑1 (MQ=255)                                                         tCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTC‑AGGGGATTAGCGCGGGGTTTTTTCCGCGG‑GTTTCCGCTCCCAACCGCGCGCCGACACAAAACATCAATAAAAGCGGcc                                                                                           >  1:378465/1‑148 (MQ=255)                                                             ccATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGC‑AGGGGATTAGCGCGGGGTTTATTCAGCGG‑GTTGCCGCCCAGAGCCGCGCGCAGCCATAATCCATCAAACAACGcggcccggcc                                                                                      >  2:833423/1‑149 (MQ=255)                                                             ccATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGC‑AGGGGAATAACGCGGGGTTTTTTCACCGG‑GTTGCCGCCCCGAGCCGCGCGCCGCCCCAATCCCTCCATCAACGcggcccggcc                                                                                      >  2:710742/1‑149 (MQ=255)                                                                ttCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGC‑CGGGGAATAAGGCGGGGTTTTTTCCGCGG‑GTTGGCCCCCCGAGACGCGCGCCGCCCAAATTCCTCCAACACAGcggcccggcc                                                                                      >  1:355546/1‑146 (MQ=255)                                                                  cGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTC‑AGGGAATTAGCGCGGGTTTTATCCAGCGG‑TTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGcc                                                                                <  2:173407/150‑1 (MQ=255)                                                                  cGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTC‑AGGGAATTAGCGCGGGGTTTTTCCAGCGG‑GTTGCCGCTCAGAGCCGCGCGCCGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGccc                                                                               >  2:796730/1‑151 (MQ=255)                                                                   gggACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTC‑AGGGGAATAGGGCGGGGTTTTTCCAACGG‑GTTTCCGCTCCGAGCCGCGGGCCGCCCTAAACCATCAATCAACCcggcccggccccggccc                                                                               >  2:593503/1‑143 (MQ=255)                                                                   gggACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTC‑AGGGAATTAGCGCGGGTTTTATCCAGCGG‑TTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGccc                                                                               <  1:662086/150‑1 (MQ=255)                                                                    ggACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTC‑AGGGAATTAGCGCGGGTTTTATCCAGCGG‑TTGGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGcc                                                                                <  1:435284/148‑1 (MQ=255)                                                                        aTCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTC‑AGGGAATTAGCGCGGGTTTTATCCAGCGG‑TTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTcc                                                                         <  2:124570/151‑1 (MQ=255)                                                                            cGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTC‑AGGGAATTAGCGCGGGTTTTATCCAGCGG‑TTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCtgtg                                                                     <  1:765361/151‑1 (MQ=255)                                                                                   tcGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGG‑AGGGAAATAGGGCGGGGGTTTTTCAGCGG‑GTTGGCGCCCCGGGCCGCGCGCCGCCATAATACATCAAAAACCGCGGCCGGGGCGCCGCCCGCCTCCTCTGGCTgct                                                               >  1:610108/1‑148 (MQ=255)                                                                                   tcGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGC‑AGGGGATTAACGCGGGGTTTTTTCAACGG‑GTTGCCGCCCAGAGCCGCGCGCCGCCATAATACCTCAATCAGCGCGGCCCGGGCCGAGCCCCCCTCCTGTGCCTGtt                                                               >  2:542031/1‑150 (MQ=255)                                                                                    cGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTC‑AGGGAATTAGCGCGGGTTTTATCCAGCGG‑TTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGtt                                                               <  2:468086/149‑1 (MQ=255)                                                                                       ttAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGC‑CGGGGAGTAGGGCGGGGGTTTTCCCGCGG‑GTTGGCGCTCCCAGGCGCGCGGAGCCAAAAACACTCAATCATCCCGGGCACGCCGGACCCCCCCTCCTCTGCCCGTTctccc                                                          >  1:431984/1‑147 (MQ=255)                                                                                       ttAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGC‑AGGGAAATAGGGCGGGGTTTTTTCACCGG‑GTTGGCGCCCAGAGGCGCGCGCAGCCCTAAACCATCAAACAACGCGGCCAGGCCGCCGCCCGCCTCCTCTGcc                                                                   >  1:635203/1‑142 (MQ=255)                                                                                                gggTAGATGCTGAGTGATAAAGTGGCGGGTC‑AGGGAATTAGCGCGGGTTTTATCCAGCGG‑TTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTc                                                  <  2:731448/150‑1 (MQ=255)                                                                                                    aGATGCTGAGTGATAAAGTGGCGGGGC‑AGGGGATTAGCGCGGGGTTTTTTCCGCGG‑GTTTCCGCCCCCAGCCGCGCGCCGCCCTAAACCCTCAAACAACGCGGCCCGGGCCGAGCCCCCCTCCCGTGCCCGTTTTCGCgggcaat                                                   >  2:868117/1‑140 (MQ=255)                                                                                                        gCTGAGTGATAAAGTGGCGGGTC‑AGGGGATTATCGCGGGGTTTATTCACCGG‑TTTGTCGCCCAGAGCCGCGCGCCGCCATAATCCATCAAACAGCGCGGCCCCGCCGTAGCCCGCGTCCTGTGCCTGTTCGCGCGGCGATTc                                                  >  1:900757/1‑142 (MQ=255)                                                                                                               gATAAAGTGGCGGGTC‑AGGGAATTAGCGCGGGTTTTATCCAGCGG‑TTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCt                                   <  2:946804/150‑1 (MQ=255)                                                                                                                   aaGTGGCGGGTC‑AGGGAATTAGCGCGGGTTTTATCCAGCGG‑TTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCa                              <  2:456565/151‑1 (MQ=255)                                                                                                                      tGGCGGGTC‑AGGGAATTAGCGCGGGGTTTATCCAGCGG‑TTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCAcc                               >  1:205846/1‑147 (MQ=255)                                                                                                                        gtgggTC‑AGGGAATTAGCGCGGTTTTTACCCAGCGG‑TTTGCCGCCAAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTc                         <  1:164706/149‑1 (MQ=255)                                                                                                                              c‑AGGGAATTAGCGCGGGTTTTATCCAGCGG‑TTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAg                    <  2:794575/150‑1 (MQ=255)                                                                                                                              c‑AGGGAATTAGCGCGGGTTTTATCCAGCGG‑TTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGc                   >  2:217900/1‑151 (MQ=255)                                                                                                                                 gggAATTAGCGCGGGTTTTATCCAGCGG‑TTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCaa                 <  2:516401/151‑1 (MQ=255)                                                                                                                                       tAGCGCGGGTTTTATCCAGCGG‑TTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAgcggc            <  1:940313/150‑1 (MQ=255)                                                                                                                                            cgGGTTTTATCCAGCGG‑TTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGgctgc      >  1:505485/1‑151 (MQ=255)                                                                                                                                                ttttATCCAGNGG‑TTNNNNNCNCNNNNNCGCNCGCANNNANNATCCATCAATNAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGAc  >  1:100414/1‑151 (MQ=255)                                                                                                                                                ttttATCCAGCGG‑TTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGAc  >  2:219004/1‑151 (MQ=255)                                                                                                                                                |                                                                                                                                                        CGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTC‑AGGGAATTAGCGCGGGTTTTATCCAGCGG‑TTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGAC  >  NC_000913/328661‑328953

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