Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,117,038	C→A	V62L (GTA→TTA)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,117,038	0	C	A	98.2%	161.4 / NA	54	V62L (GTA→TTA)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (25/28);  minor base T (0/1);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGG  >  NC_000913/4116894‑4117181                                                                                                                                                 |                                                                                                                                                cGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGc                                                                                                                                             <  2:408697/149‑1 (MQ=255)   gTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGt                                                                                                                                           <  2:71847/149‑1 (MQ=255)    ttcgacgGCACTGCCAGACTATTGTGTTATAGATTGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTg                                                                                                                                          <  1:250046/147‑1 (MQ=255)    tACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTg                                                                                                                                          <  2:934131/149‑1 (MQ=255)    tACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTg                                                                                                                                          <  2:934069/149‑1 (MQ=255)        acgGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGc                                                                                                                                      >  1:926856/1‑149 (MQ=255)               tGCCCGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCACAGACAATAGTGGTTTAAAGCTGGTACGTAATACCGATAGTTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGt                                                                                                                               >  2:928345/1‑149 (MQ=255)               tGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCTTTTTTTTTAGTGGGTTCACGCGGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCTGCACTTCTAAAAGCGTGGAGCTTTGGGt                                                                                                                               >  2:463406/1‑149 (MQ=255)                 ccAGACACTGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTgg                                                                                                                             <  1:669991/149‑1 (MQ=255)                 ccAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTgg                                                                                                                             >  2:722450/1‑149 (MQ=255)                 ccAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCGAACAGCGTGGAGCTTTGGGTgg                                                                                                                             >  1:630538/1‑149 (MQ=255)                  cAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCCCTTCTAACAGCGTGGAGCTTTGGGTGGc                                                                                                                            >  2:261349/1‑149 (MQ=255)                   aGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACACTAGTGGTTTCACGCTGGTTCGTAGTACCGATAGCTGCAATTTGGTCGGGAATGATATCGGCTTTTGCCAGCACTTCTAACACCATGGAGCTTTGGGTGGcc                                                                                                                           >  2:773164/1‑149 (MQ=255)                   aGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATCTCGGCTTTCGCCAGCACGTCTAACAGCGTGGAGCTTTGGGTGGcc                                                                                                                           >  2:404452/1‑149 (MQ=255)                   aGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTGGcc                                                                                                                           >  2:773167/1‑149 (MQ=255)                       aaTGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGa                                                                                                                       >  1:612280/1‑149 (MQ=255)                                atagatagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTg                                                                                                              <  2:760836/149‑1 (MQ=255)                                 tagatagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAATAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTgg                                                                                                             <  2:239300/149‑1 (MQ=255)                                       ggCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtg                                                                                                       <  1:236173/149‑1 (MQ=255)                                          ttGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTc                                                                                                    >  2:821947/1‑149 (MQ=255)                                           tGCCGGTTTGCTTTTCCCAGCCTATATTGCGTTCACGCTGGTTTGTAATCCCGACAGCTGCAATTTGATCGGAACTGATGTCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCt                                                                                                   <  2:887944/149‑1 (MQ=255)                                             cccgTTTCTTTTTCCCAGCCAATAGTGGTTTCACGCGTGTTCGTAATACCGATGGCCGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTAc                                                                                                 <  2:318024/146‑1 (MQ=255)                                             ccGGTTTCTGTTTCTCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTAc                                                                                                 <  1:589881/149‑1 (MQ=255)                                                gTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTACTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCa                                                                                              >  1:863262/1‑149 (MQ=255)                                                     tttttCCCAGACAGTAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGCTATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTg                                                                                         >  2:655609/1‑149 (MQ=255)                                                     tttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGGAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTg                                                                                         >  1:647948/1‑149 (MQ=255)                                                          cccAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGtttt                                                                                    <  1:821947/149‑1 (MQ=255)                                                            cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                                                  <  1:446137/149‑1 (MQ=255)                                                            cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                                                  >  1:452863/1‑149 (MQ=255)                                                            cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                                                  <  1:608509/149‑1 (MQ=255)                                                                 aaTAGTGGTTTCACGCTGGTTCGTAATACCGATCGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCGACCTGGTTTTGGGTAGa                                                                             >  2:541238/1‑149 (MQ=255)                                                                     gTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTg                                                                         <  2:334535/149‑1 (MQ=255)                                                                         tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCCTTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCa                                                                     <  1:930844/149‑1 (MQ=255)                                                                           tCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGCAGTCCGTTCGTTTTGACGCGGGCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCaaa                                                                   <  1:928345/149‑1 (MQ=255)                                                                             aCGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGGTTTGGGTAGGTTTGCTCAAAtt                                                                 >  2:311482/1‑149 (MQ=255)                                                                               gCTGGTTCGTAATAACGATAGCTGCAATTTGAGCGGATCTGGTATCGGCTTTCGCCAGCACTTCTAACAGCGTGCGGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTcg                                                               <  2:85579/149‑1 (MQ=255)                                                                                    ttCGTAATCCCGATGGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTg                                                          <  2:56461/149‑1 (MQ=255)                                                                                    ttCGCAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTg                                                          <  1:435335/149‑1 (MQ=255)                                                                                              cGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTa                                                >  1:290307/1‑149 (MQ=255)                                                                                                     tGCAATTTGATCGGAACTGGTATCGGCTTTCGCCAGCAATTATCACAGAGTGGAGATTTGACTGCCCCTTATGACCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGGTCAAATTCGCGCTCCGACACGCTAATGatat                                         >  2:436300/1‑149 (MQ=255)                                                                                                      gCAATTTGATCGGAACTGATAGCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATAtt                                        <  2:838678/149‑1 (MQ=255)                                                                                                             gATCGGAACTGATATCGGCTTTCGCCAGCGCTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGATTTGGGTAGATTTGCTCAAAGTGGCGCTGCGACACGCTAAAGATATGGGCATCg                                 >  2:819243/1‑149 (MQ=255)                                                                                                                cGGAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGTCACGCTAATGATATTGGCATCGTGa                              <  1:890673/149‑1 (MQ=255)                                                                                                                 ggAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGAt                             <  2:751615/149‑1 (MQ=255)                                                                                                                 ggAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGGAGATTTGTTCAAATTCGCGCTGCGACACGCTAATGATTTTGGCATAGTggt                             >  2:798109/1‑147 (MQ=255)                                                                                                                  gAACTGATATCGGCTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATc                            <  1:660850/149‑1 (MQ=255)                                                                                                                        tgatCGGCTTTCCCCAGCACTTCTATCAGCGTTGAGTTTTGGGTGGCCCAGTTTTCCACTGGCTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACCCGCTAACGATATTGGCATCGTGACCCATTAc                      <  2:88312/147‑1 (MQ=255)                                                                                                                               cTTTCGCCAGCACTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACcgcg               <  2:630538/149‑1 (MQ=255)                                                                                                                                      cAGCACTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCCTTACGACCGCGCGGGacc        >  1:354025/1‑147 (MQ=255)                                                                                                                                      cAGCACTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGc        >  1:820126/1‑149 (MQ=255)                                                                                                                                          aCTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCAGTAGGGCCGCGCGGGGGCcggc    >  2:331722/1‑148 (MQ=255)                                                                                                                                           cTTCTAACGGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCtggtg   <  2:547311/149‑1 (MQ=255)                                                                                                                                           cTTCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCtggtg   <  1:624050/149‑1 (MQ=255)                                                                                                                                            ttCTAACAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGGGCtggtgg  >  1:481850/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                CGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGG  >  NC_000913/4116894‑4117181

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