Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,234,335	T→A	7.8%	E317D (GAA→GAT)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,234,335	0	T	A	7.8%	143.6 / 2.7	51	E317D (GAA→GAT)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand):  ref base T (25/22);  new base A (2/2);  total (27/24)
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.41e-01

AACAACGCACCACGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCGTGGGACGACTTCACCGCCTGGTCGCCCGTGTTCATCAACAGGA  >  NC_002947/1234198‑1234483                                                                                                                                          |                                                                                                                                                     aacaacGCACCACGGGCATACGGGTCCCAGTACGGGGCACACAGGCCGGTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAg                                                                                                                                           >  5:202069/1‑149 (MQ=255)    aacGCACCACGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGAACTGAAGATCGTCaa                                                                                                                                               <  2:35908/142‑2 (MQ=255)    aacGCACCACGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGAACTGAAGATCGTCaa                                                                                                                                               >  1:35908/1‑141 (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                                                                                                                            >  3:152482/1‑144 (MQ=255)                     cGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTaga                                                                                                                            <  4:152482/144‑1 (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)                                                ggTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTa                                                                                                                              <  2:291717/115‑1 (MQ=255)                                                ggTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTa                                                                                                                              >  1:291717/1‑115 (MQ=255)                                                 gttaaGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTa                                                                                                                              <  6:265113/111‑1 (MQ=255)                                                 gTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATTTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATa                                                                                           <  2:318562/149‑1 (MQ=255)                                                 gTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTa                                                                                                                              >  5:265113/1‑114 (MQ=255)                                                        cTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATa                                                                                           >  3:164731/1‑142 (MQ=255)                                                        cTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATa                                                                                           <  4:164731/142‑1 (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                                          <  7:143784/122‑1 (MQ=255)                                                                                                                         tgcgtGACGAACTGAAGATCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCGTgg                                          >  8:143784/5‑126 (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)                                                                                                                                          ttCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCGTGGGACGACTTCACCGCCTGGTCGCCCGTGTTCATCAACAGGa  <  3:58787/149‑1 (MQ=255)                                                                                                                                          |                                                                                                                                                     AACAACGCACCACGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTACACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCGTGGGACGACTTCACCGCCTGGTCGCCCGTGTTCATCAACAGGA  >  NC_002947/1234198‑1234483

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