Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,234,331	T→G	7.4%	K319Q (AAG→CAG)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,234,331	0	T	G	7.4%	161.7 / 2.8	53	K319Q (AAG→CAG)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand):  ref base T (26/23);  new base G (2/2);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

CCACGTGTCAGGCCGAACAACGCACCACGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCGTGGGACGACTTCACCGCCTGGTCGCCCGTGT  >  NC_002947/1234183‑1234471                                                                                                                                                     |                                                                                                                                             ccACGTGTCAGGCCGAACAACGCACCACGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCtt                                                                                                                                              >  4:71576/1‑149 (MQ=255)               gAACAACGCACCACGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAg                                                                                                                                           >  5:279293/1‑138 (MQ=255)               gAACAACGCACCACGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAg                                                                                                                                           <  6:279293/138‑1 (MQ=255)                aacaacGCACCACGGGCATACGGGTCCCAGTACGGGGCACACAGGCCGGTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAg                                                                                                                               >  5:202069/1‑149 (MQ=255)                   aacGCACCACGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGAACTGAAGATCGTCaa                                                                                                                                   >  1:35908/1‑141 (MQ=255)                   aacGCACCACGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGAACTGAAGATCGTCaa                                                                                                                                   <  2:35908/142‑2 (MQ=255)                    acGCACCACGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCAt                                                                                                                           >  5:329303/1‑149 (MQ=255)                        accacGGGCACACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGGCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCa                                                                                                                       <  7:337880/149‑1 (MQ=255)                          cacGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACg                                                                                                                     >  6:21963/1‑149 (MQ=255)                            cGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGt                                                                                                                   <  2:296774/149‑1 (MQ=255)                               gCATACGGGTCCCAGTACGGGGCACCCAGGCCGTTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTAGGTATCCAGGGAGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTagc                                                                                                                >  3:300405/1‑148 (MQ=255)                                  tACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGcc                                                                                                             >  1:173567/1‑149 (MQ=255)                                    cGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTaga                                                                                                                <  4:152482/144‑1 (MQ=255)                                    cGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTaga                                                                                                                >  3:152482/1‑144 (MQ=255)                                     gggTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGccac                                                                                                           >  6:106563/1‑148 (MQ=255)                                       gTCCCAGTACGGGTCACCCAGGCCTGTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCAt                                                                                                        <  4:300405/149‑1 (MQ=255)                                          ccAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTg                                                                                                      >  7:157321/1‑148 (MQ=255)                                           cAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGaa                                                                                                    >  3:270663/1‑149 (MQ=255)                                                    gCACCCAGGCCGGTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATTTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGcc                                                                                           >  1:318562/1‑149 (MQ=255)                                                     cACCCAGGCCGGTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCt                                                                                          <  3:71576/149‑1 (MQ=255)                                                       cccAGGCCGGTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTc                                                                                        <  2:173567/149‑1 (MQ=255)                                                         cAGGCCGGTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCa                                                                                      <  5:21963/149‑1 (MQ=255)                                                          aGGCCGGTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAg                                                                                     >  1:116371/1‑149 (MQ=255)                                                           ggCCGGTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGt                                                                                    <  6:329303/149‑1 (MQ=255)                                                              cGGTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCtt                                                                                                                                              <  4:292380/88‑1 (MQ=255)                                                              cGGTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCtt                                                                                                                                              >  3:292380/1‑88 (MQ=255)                                                               ggTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTa                                                                                                                  >  1:291717/1‑115 (MQ=255)                                                               ggTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTa                                                                                                                  <  2:291717/115‑1 (MQ=255)                                                                gttaaGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTa                                                                                                                  <  6:265113/111‑1 (MQ=255)                                                                gTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATTTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATa                                                                               <  2:318562/149‑1 (MQ=255)                                                                gTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTa                                                                                                                  >  5:265113/1‑114 (MQ=255)                                                                       cTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATa                                                                               <  4:164731/142‑1 (MQ=255)                                                                       cTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATa                                                                               >  3:164731/1‑142 (MQ=255)                                                                            aCCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTcac                                                                    >  3:34414/1‑148 (MQ=255)                                                                                 gTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCCCGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCAcgcg                                                              >  7:272655/1‑149 (MQ=255)                                                                                    cacGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGc                                                           <  2:116371/149‑1 (MQ=255)                                                                                             gCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGAGCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCaa                                                  >  4:58787/1‑149 (MQ=255)                                                                                                 tCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCATGCGGGCCGCAGGCAATg                                                <  4:270663/147‑1 (MQ=255)                                                                                                 tCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATg                                                <  4:34414/147‑1 (MQ=255)                                                                                                       tacTTGCTGGCGATATACTCGGTATAGAGGGCGTAGTTGACGATCTTCAGTTCGTCAAGGAGCCATTGCACGGTAGATTAACCATTGAATAACGCGCCTTCCTTTGCATAGGCCACCTCTCCACGCGGGACTCAGGCAATTGTCGTGa                                         <  4:217601/146‑1 (MQ=255)                                                                                                       aCCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCTATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAg                                        >  6:45402/1‑149 (MQ=255)                                                                                                       aCCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAg                                        <  8:157321/149‑1 (MQ=255)                                                                                                                gCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCGTg                               <  5:106563/149‑1 (MQ=255)                                                                                                                  gAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCGTgg                              >  6:145328/1‑148 (MQ=255)                                                                                                                   aaaTACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCGTGGga                            >  4:103903/1‑149 (MQ=255)                                                                                                                            gTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAg                                                                              >  2:298321/1‑90 (MQ=255)                                                                                                                                 cAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCgg                                                             >  5:31705/1‑102 (MQ=255)                                                                                                                                 cAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCgg                                                             <  6:31705/102‑1 (MQ=255)                                                                                                                                    ggCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCGTGGGACGACTTCACCGCCTGGt           <  5:145328/149‑1 (MQ=255)                                                                                                                                        tgcgtGACGAACTGAAGATCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCGTgg                              >  8:143784/5‑126 (MQ=255)                                                                                                                                        tgcgtGACGAACTGAAGATCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCGTgg                              <  7:143784/122‑1 (MQ=255)                                                                                                                                          gtTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCGTGGGACGACTTCACCGCCTGGTCGCCCg     >  7:315755/1‑149 (MQ=255)                                                                                                                                             gACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCGTGGGACGACTTCACCGCCTGGTCGCCCgtgt  >  1:134062/1‑149 (MQ=255)                                                                                                                                              aCGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCGTGGGACGACTTCAcc                  >  7:38176/1‑132 (MQ=255)                                                                                                                                              aCGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCGTGGGACGACTTCAcc                  <  8:38176/132‑1 (MQ=255)                                                                                                                                                gATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGt                                           <  3:338400/105‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                             CCACGTGTCAGGCCGAACAACGCACCACGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCGTGGGACGACTTCACCGCCTGGTCGCCCGTGT  >  NC_002947/1234183‑1234471

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