Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,743,657	0	A	C	44.0%	77.5 / 15.2	50	T309P (ACC→CCC)	yiaK	2,3‑diketo‑L‑gulonate reductase, NADH‑dependent
Reads supporting (aligned to +/- strand):  ref base A (24/4);  new base C (0/22);  total (24/26)
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

TTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGGCGC  >  NC_000913/3743509‑3743806                                                                                                                                                     |                                                                                                                                                      tttacaAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTAcc                                                                                                                                                     >  2:790233/4‑151 (MQ=255) ttCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTtactac                                                                                                                                                      <  1:398448/150‑1 (MQ=255) ttCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTtactac                                                                                                                                                      <  2:413991/150‑1 (MQ=255) ttCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTAcc                                                                                                                                                     >  1:356508/1‑151 (MQ=255) ttCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTAcc                                                                                                                                                     >  1:435043/1‑151 (MQ=255)     caAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCctgc                                                                                                                                                 >  1:186641/1‑151 (MQ=255)     caAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCctgc                                                                                                                                                 >  1:322503/1‑151 (MQ=255)           tttATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGcc                                                                                                                                            <  1:163695/150‑1 (MQ=255)           tttATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCg                                                                                                                                           >  1:492858/1‑151 (MQ=255)            ttATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCg                                                                                                                                           >  2:459638/1‑150 (MQ=255)              aTTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGaa                                                                                                                                         <  2:463507/150‑1 (MQ=255)              aTTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGaa                                                                                                                                         >  1:624168/1‑150 (MQ=255)                 gCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAAcc                                                                                                                                     >  2:807892/1‑151 (MQ=255)                    aTTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAAccgcc                                                                                                                                  <  2:426492/151‑1 (MQ=255)                            ggACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTATGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGc                                                                                                                          >  2:679165/1‑151 (MQ=255)                             gACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGc                                                                                                                          >  2:958167/1‑150 (MQ=255)                                    ttATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGt                                                                                                                  >  2:158378/1‑151 (MQ=255)                                    tatccGACGCTCCCCCCCGCGATGCCAAGCTCCAACCCATCATGGTTTACGTTACTAGTGCGGGGCGTGCTGACAAAAATCAGGCCATTCGCTTCCCCGCCCATGAATTTCCTCCCCTGCGGGCAAAAAACCGCCGTAACGGCATCACt                                                                                                                    <  1:705375/145‑1 (MQ=255)                                     tATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGtt                                                                                                                 >  2:870914/1‑151 (MQ=255)                                        cgccGGCCCCCCCCGCGATGCCAAGCTGCAACGCATCATGTATTCCTTTACTGGTGCCGGGCGTGCTGACAAAAATCGGCCCTTTCGCTTCCCCGGCCATGATTTTACTCCCCTGCGGGCGAAAAACCGCCGTAACGGCATCACTGTtga                                                                                                               <  1:749891/147‑1 (MQ=255)                                        cGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTtgat                                                                                                              >  2:268035/1‑151 (MQ=255)                                         gacgtcccccccCGCGAGGCCAAGCTGCAACGCATCATGGTTTACTTTACTGGTGCGGAGCGTGCTGACAAAATTCAGGCCTTTCGCTTCCCCGGCCATAAATTTACTCCCCTGCGGCCCAAAAACCGCCGTAACGGCATCACTGt                                                                                                                  <  1:933012/140‑1 (MQ=255)                                            ggcccccccCGCGATGCCAAGCTGCACCGCATCATGGATTACTTTACTAGTGCCGGGCGTGCTGACGAAAATCAGGCCATCCCCTTCCCCGGCAATGAATTTACTCCCCGGCTGCCCAAAAACCGCCGTAACGGCATCACTGTtgatga                                                                                                            <  1:169493/146‑1 (MQ=255)                                               ccccccCGCGATGCCAAGCTGCACCGAATAATGGATTACGTTACTAGGGCCGGGCGTGCTGACAAAAATCAGGCCATTCGCTTCCCCGGCAATGAATTTACTCCCCTGCTGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACa                                                                                                          <  1:220177/148‑1 (MQ=255)                                                ccccccGCGATGCAAAGCTGAACCGCATCTTGGATTCCGTTACTAGTGCCGGGCGTGCTGACAAAAACCAGGCCATCCGCTTCCCCGGCCATGATTTTCCTCCCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCgt                                                                                                      <  1:86996/148‑1 (MQ=255)                                                        gATCCCAAGCTGCAACCCTCCATGGTTTACGTTACTTGTGCCGAGCTTGCTGACAAAAACCAGGCCATTCGCTTCCCCGGCCATGAATTTCCTCCCCTCCGGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGc                                                                                                <  1:817309/149‑1 (MQ=255)                                                                 cTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAgg                                                                                     >  2:540981/1‑151 (MQ=255)                                                                  tGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGc                                                                                    >  2:951992/1‑151 (MQ=255)                                                                     accGCATCATGGATTACTTTACTAGTGCCGGGGGTGCTGACAAAAACCAGGCCTTCCGCTTCCCCGCCAATAATTTTACTCCCCTGCTGGCGAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGtt                                                                                 <  2:569744/149‑1 (MQ=255)                                                                     aacacctcatGGATTACTTTACTGGTGCCGGGGGTGCTGACAAAAATCGGCCCCTCCGTTTCCCCGGCCATGTTTTTCCTCCCCTGTTGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAgg                                                                                     <  2:926113/141‑1 (MQ=255)                                                                     aaCGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGtt                                                                                 >  2:85059/1‑151 (MQ=255)                                                                         catcatGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGCTGATGACAGCGTGTGGGCCAACATCCAGGCGTTAt                                                                               >  1:794134/1‑149 (MQ=255)                                                                                      gTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCAt                                                                >  2:280162/1‑151 (MQ=255)                                                                                       ttccTAGTGCCGGGCGTGCTGACAAAAACCAGCCCATTCCCTTCCCCGGCAATAAATTTCCTCCCCTGCTGCCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCa                                                                 <  2:204969/146‑1 (MQ=255)                                                                                          cTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGAtt                                                            >  2:483563/1‑151 (MQ=255)                                                                                             gtccggggCGTGCTGAAAAAAAGCAGCCCTTCCGCTTCCCCGGCATTGATTTTACTCCCCTGTGGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGAtt                                                            <  2:798314/141‑1 (MQ=255)                                                                                                cggggCGTGCGGACAAAAACCGGGCCTTTCGCTTCCCCGGCCATAATTTTACTCCCCTGCTGGCGAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTg                                                        <  1:536727/145‑1 (MQ=255)                                                                                                cggAGCGTGCTAACGAAAACCAGCCCATTCGTTTCCCCGGCAATAATTTTCCGCCCCTCCTGCCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATG‑GGAGATAAGTCATGATTTTTgg                                                       <  1:82154/147‑1 (MQ=255)                                                                                                ccgggggTGCTGAAAAAAATCAGCCCATCCGCTTCCCCGCCCAGAAATTTACTACCCGGCGGGCAAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTgg                                                       <  1:952360/144‑1 (MQ=255)                                                                                                ccGGGCGTGCTAACAAAAACCGGCCCATCCCCTTCCCCGCCCATAATTTTACTCCCCTGCGGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGAt                                                             <  1:897148/144‑1 (MQ=255)                                                                                                    gCGTGCGGACAAAAATCAGGCCTTTCCCTTCCCCGCCCATGAATTTACTCCCCTGCTGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACata                                                  <  1:967462/151‑1 (MQ=255)                                                                                                          tGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCgcgc                                            >  1:35132/1‑151 (MQ=255)                                                                                                               aaaaCCAGGCCTTTCCCTTCCCCGGCCATGATTTTACTCCCCGGCGGGCAAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCgcgc                                            <  2:623307/146‑1 (MQ=255)                                                                                                               aaaaCCAGCCCTTCCTCTTCCCCGGCAATAATTTTACTCCCCGGCGGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACa                                                    <  2:842202/138‑1 (MQ=255)                                                                                                                aaaccggcccATTCGTTTCCCCGCCCATAATTTTCCTCCCCTGCTGCCAAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCa                                      <  1:908355/143‑1 (MQ=255)                                                                                                                      ggCCTTTCCCTTCCCCGGCAATAATTTTACTCCCCGGCGGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAAt                                    <  2:186626/147‑1 (MQ=255)                                                                                                                          aTTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACAGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCg                            >  1:223494/1‑151 (MQ=255)                                                                                                                             cgtttccccGCCCATAATTTTCCTCCCCTGCGGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGttt                         <  2:954849/145‑1 (MQ=255)                                                                                                                             cGCTTCCCCGGCCATAAATTTACTCCCCTGCCGGCAAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGttt                         <  2:937638/151‑1 (MQ=255)                                                                                                                                   cccGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCccgc                   >  1:157918/1‑151 (MQ=255)                                                                                                                                        ccATAATTTTACTCCCCTGCTGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCa              <  1:483351/151‑1 (MQ=255)                                                                                                                                             aaTTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGa          >  1:318135/1‑150 (MQ=255)                                                                                                                                                    tacCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGgcgc  >  2:171166/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                      TTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGGCGC  >  NC_000913/3743509‑3743806

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