Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,756,434	A→G	D246G (GAC→GGC)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,756,434	0	A	G	98.2%	203.0 / ‑4.0	54	D246G (GAC→GGC)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (24/29);  total (24/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

GACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGC  >  NC_000913/1756292‑1756575                                                                                                                                               |                                                                                                                                              gACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGgtg                                                                                                                                       >  1:903621/1‑151 (MQ=255) gACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGgtg                                                                                                                                       >  2:841850/1‑151 (MQ=255) gACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGgtg                                                                                                                                       >  1:595794/1‑151 (MQ=255)  aCGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGgtgt                                                                                                                                      >  1:411182/1‑151 (MQ=255)    gTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAg                                                                                                                                    <  2:746432/151‑1 (MQ=255)      tATCGAAATCCGTTAGCACCTGAAAGCGCCCGGGGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGa                                                                                                                                   <  2:149033/150‑1 (MQ=255)       aTCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGaaa                                                                                                                                 <  2:933255/151‑1 (MQ=255)         cGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAAt                                                                                                                                <  2:502104/150‑1 (MQ=255)                cGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGt                                                                                                                          >  2:52815/1‑149 (MQ=255)                      cACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGaagaag                                                                                                                  <  1:194630/151‑1 (MQ=255)                         cTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGtt                                                                                                                >  1:767248/1‑150 (MQ=255)                           gAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTa                                                                                                               <  2:434464/149‑1 (MQ=255)                            aaaGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCt                                                                                                            <  1:209826/151‑1 (MQ=255)                             aaGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCtt                                                                                                           >  1:825061/1‑151 (MQ=255)                                   aCGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGccc                                                                                                      >  2:488502/1‑150 (MQ=255)                                     ggcggcGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCaga                                                                                                   <  1:921701/151‑1 (MQ=255)                                       cggcgAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCagaag                                                                                                 <  1:868949/151‑1 (MQ=255)                                        ggcgAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCaga                                                                                                   >  2:114692/1‑148 (MQ=255)                                           gAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAgatg                                                                                              <  2:528119/150‑1 (MQ=255)                                              aaCATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAgatgatga                                                                                          <  1:29119/151‑1 (MQ=255)                                              aaCATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAgatgatg                                                                                           >  2:773137/1‑150 (MQ=255)                                                cATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATc                                                                                        <  2:1022996/151‑1 (MQ=255)                                                 aTCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCg                                                                                       >  1:671581/1‑151 (MQ=255)                                                    cacaTCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGaaa                                                                                    <  2:639146/151‑1 (MQ=255)                                                         catcTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGt                                                                               <  1:143711/151‑1 (MQ=255)                                                              ccAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATcc                                                                           <  1:1020852/150‑1 (MQ=255)                                                                aaaaTCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTg                                                                        <  2:38967/151‑1 (MQ=255)                                                                  aaTCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCa                                                                      <  2:384827/151‑1 (MQ=255)                                                                   aTCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCAc                                                                     <  2:595891/151‑1 (MQ=255)                                                                    tCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACg                                                                    >  2:945989/1‑151 (MQ=255)                                                                    tCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACg                                                                    >  2:1108079/1‑151 (MQ=255)                                                                        aaaCCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTaa                                                                 >  1:794487/1‑150 (MQ=255)                                                                         aaCCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTaa                                                                 >  1:663061/1‑149 (MQ=255)                                                                           ccAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTc                                                             <  2:20767/151‑1 (MQ=255)                                                                             aGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGt                                                           >  2:783645/1‑151 (MQ=255)                                                                               gAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGtt                                                          >  1:921037/1‑150 (MQ=255)                                                                               gAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTa                                                         <  2:937953/151‑1 (MQ=255)                                                                                  ggCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCa                                                      <  1:721788/151‑1 (MQ=255)                                                                                            acTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGAccc                                             >  2:39611/1‑150 (MQ=255)                                                                                               tCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGAccc                                             >  2:831882/1‑147 (MQ=255)                                                                                                cgacgaAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGAt                                         <  1:49068/150‑1 (MQ=255)                                                                                                  acgaAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCt                                      <  1:202560/151‑1 (MQ=255)                                                                                                     aaaTCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTgg                                    >  1:951803/1‑150 (MQ=255)                                                                                                             gAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGa                           <  2:554419/151‑1 (MQ=255)                                                                                                              aaGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGAt                          <  2:399924/151‑1 (MQ=255)                                                                                                               aGCCTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATc                         <  1:1030282/151‑1 (MQ=255)                                                                                                                 ccTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCaa                       <  1:893143/151‑1 (MQ=255)                                                                                                                 ccTCTGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCa                        >  2:254058/1‑150 (MQ=255)                                                                                                                    ctGACGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCaaaa                     >  1:680400/1‑150 (MQ=255)                                                                                                                        cGGCATCATGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAAccc                 >  2:601703/1‑150 (MQ=255)                                                                                                                               atGGTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGa         <  1:965428/151‑1 (MQ=255)                                                                                                                                 ggTTGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGAc        <  2:648902/150‑1 (MQ=255)                                                                                                                                    tGCGCGTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTcg     <  2:613462/150‑1 (MQ=255)                                                                                                                                        cgTGGCGGCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCgcgc  >  2:382843/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                              GACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGC  >  NC_000913/1756292‑1756575

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