Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,233,768	T→A	13.4%	intergenic (‑84/+18)	glpR ← / ← glpK	DNA‑binding transcriptional repressor GlpR/glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,233,768	0	T	A	13.4%	126.5 / 9.0	51	intergenic (‑84/+18)	glpR/glpK	DNA‑binding transcriptional repressor GlpR/glycerol kinase
Reads supporting (aligned to +/- strand):  ref base T (22/22);  new base A (3/4);  total (25/26)
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

ACATAGCCGCGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGACAA‑ATTCTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTCCCAGTCACGGGTACGCTCGACCGCCTTCTTCCAGCCGGCATACAGCTTTTCTTTCTGCGTTTCGTCCAGTTGCGGGCTGAACTCACGCTCGATGATCGCCTTGTCGCGCAGCT  >  NC_002947/1233626‑1233916                                                                                                                                                 |                                                                                                                                                     acacAGCCGCCTTCGCGAACCTCTACTCGCATTTGTCCTTCCCGGGGGCCCCGACTCATGCCCGGTCCTTTGGGGCTGCCGGACAA‑ATTCTGCAATGATCCCGTAGGAAGGGGATC‑ACGCCCAGTAGACGATTCGGGCCGCATggcg                                                                                                                                                  <  8:209722/147‑1 (MQ=255)  cATAGCCGCGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGACAA‑ATTCTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCCATGCGGcc                                                                                                                                               <  6:2417/149‑1 (MQ=255)     aGCCGGGTTCGCGTACCAGTTCGCGGATTGGTTGTGGGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCGGCCGGACAA‑ATTCTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGTTTGGGGCCGCCATGCGGCCCCa                                                                                                                                            <  7:216212/149‑1 (MQ=255)      gCCCCGTCCGCGTACCAGGTCGAGGATTTGTTGTTGGCGGGGCGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGACAA‑ATTCTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCCATGCGGCCCCaa                                                                                                                                           <  3:142152/149‑1 (MQ=255)       ccGCGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGACAA‑ATTCTGCAATGATGCCGTAGGAAGGGGATGAACGGCCAGTAGACGATTGGGGCCGCCATGCGGCCCCaa                                                                                                                                           <  1:61349/149‑1 (MQ=255)       ccGCGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGACAA‑ATTCTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCaa                                                                                                                                           >  6:200781/1‑148 (MQ=255)        cgcgTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGACAA‑ATTCTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACa                                                                                                                                         <  5:200781/149‑1 (MQ=255)           gTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGACAA‑ATTCTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGc                                                                                                                                       >  6:297677/1‑148 (MQ=255)               gcgTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGACAA‑ATTCTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCa                                                                                                                                  >  1:301239/1‑149 (MQ=255)                  tACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGACAA‑ATTCTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGc                                                                                                                               <  6:229270/149‑1 (MQ=255)                  tACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGACAA‑ATTCTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGc                                                                                                                               <  2:301239/149‑1 (MQ=255)                       gTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGACAA‑ATTCTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCCATGCGGCCCCaatc                                                                                                                                         >  3:212281/1‑132 (MQ=255)                         tCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGACAA‑ATTCTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGa                                                                                                                        >  8:251928/1‑149 (MQ=255)                            aGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGACAA‑ATTCTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCt                                                                                                                     <  6:58611/149‑1 (MQ=255)                                       ttgGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGACAA‑ATTCTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTCCCAGTCACg                                                                                                          >  6:215518/1‑149 (MQ=255)                                                   cAGATTCATGGGCGGTCCTTTGGGGCTGCCGGACAA‑ATTCTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTCCCAGTCACGGGTACGCTCGAc                                                                                              <  7:251928/149‑1 (MQ=255)                                                      aTTCATGGGCGGTCCTTTGGGGCTGCCGGACAA‑ATTCTGCAATGATGTCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTAt                                                                                                                                    >  1:45778/1‑108 (MQ=255)                                                      aTTCATGGGCGGTCCTTTGGGGCTGCCGGACAA‑ATTCTGCAATGATGTCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTAt                                                                                                                                    <  2:45778/108‑1 (MQ=255)                                                      aTTCATGGGCGGTCCTTTGGGGCTGCCGGACAA‑ATTCTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTCCCAGTCACGGGTACGCTCGACCGc                                                                                           >  3:129583/1‑149 (MQ=255)                                                         cATGGGCGGTCCTTTGGGGCTGCCGGACAA‑ATTCTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTCCCAGTCACGGGTACGCTCGACCGCct                                                                                         <  5:311443/148‑1 (MQ=255)                                                         cATGGGCGGTCCTTTGGGGCTGCCGGACAA‑ATTCTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTCCCAGTCACGGGTACGCTCGACCGCct                                                                                         <  5:218495/148‑1 (MQ=255)                                                            gggCGGTCCTTTGGGGCTGCCGGACAATATTCTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCAt                                                                                                                                                      >  1:14266/1‑85 (MQ=255)                                                            gggCGGTCCTTTGGGGCTGCCGGACAA‑ATTCTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCAt                                                                                                                                                      <  2:14266/84‑1 (MQ=255)                                                               cGGTCCTTTGGGGCTGCCGGACAA‑ATTCTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTCCCAGTCACGGGTACGCTCGACTGCCttcttc                                                                                    >  8:180891/1‑147 (MQ=255)                                                                  tCCTTTGGGGCTGCCGGACAA‑ATTCTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTccc                                                                                                                 <  7:291697/115‑1 (MQ=255)                                                                  tCCTTTGGGGCTGCCGGACAA‑ATTCTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTccc                                                                                                                 >  8:291697/1‑115 (MQ=255)                                                                  tCCTTTGGGGCTGCCGGACAA‑ATTCTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTCCCAGTCACGGGTACGCTCGACCGCCTTCTTCCAGcc                                                                               >  2:307635/1‑149 (MQ=255)                                                                     tttGGGGCTGCCGGACAA‑ATTCTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTCCCAGTCACGGGTACGCTCGACCGCCTTCTTCCAGCCgg                                                                             >  7:116416/1‑148 (MQ=255)                                                                                cGGACAA‑ATTCTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTCCCAGTCACGGGTACGCTCGACCGCCTTCTTCCAGCCGGCATACAGCtttt                                                                 <  1:307635/149‑1 (MQ=255)                                                                                 ggACAA‑ATTCTGCAATGATGCCGTAGGAATGGGATT‑ACGGACAGAAGACGATTGGGGCCGCAAGGCGGCCCCAACAGCTATCATGCCTAGTGATCTTCCCAGACACGTGTACGATCGAACTCATTCTTCCACCCGGAATACAGAtttac                                                                >  7:60632/1‑147 (MQ=255)                                                                                 ggACAA‑ATTCTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTCCCAGTCACGGGTACGCTCGACCGCCTTCTTCCAGCCGGCATACAGCTtttc                                                                >  2:174489/1‑149 (MQ=255)                                                                                  gACAA‑ATTCTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTc                                                                                                                            <  7:94169/88‑1 (MQ=255)                                                                                  gACAA‑ATTCTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTc                                                                                                                            >  8:94169/1‑88 (MQ=255)                                                                                         ttCTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTCCCAGTCACGGGTACGCTCGACCGCCTTCTTCCAGCCGGCATACAGCTTTTCTTTCTGc                                                         <  5:297677/149‑1 (MQ=255)                                                                                           cTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTCCCAGTCACGGGTACGCTCGACCGCCTTCTTCCAGCCGGCATACAGCTTTTCTTTCTGCGt                                                       >  2:228033/1‑149 (MQ=255)                                                                                            tGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTCCCAGTCACGGGTACGCTCGACCGCCTTCTTCCAGCCGGCATACAGCTTTTCTTTCTGCGtt                                                      <  1:174489/149‑1 (MQ=255)                                                                                               aaTGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTCCCAGTCACGGGTACGCTCGACCGCCTTCTTCCAGCCGGCATACAGCttt                                                                  >  6:315733/1‑134 (MQ=255)                                                                                               aaTGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTCCCAGTCACGGGTACGCTCGACCGCCTTCTTCCAGCCGGCATACAGCttt                                                                  <  5:315733/134‑1 (MQ=255)                                                                                                       cGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTCCCAg                                                                                                               <  7:126269/81‑1 (MQ=255)                                                                                                       cGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTCCCAg                                                                                                               >  8:126269/1‑81 (MQ=255)                                                                                                       cGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTCCCAGTCACGGGTACGCTCGACCGCCTTCTTCCAGCCGGCATACAGCTTTTCTTTCTGCGTTTCGTCCAGTTg                                           >  4:175130/1‑149 (MQ=255)                                                                                                             aaGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTCCCAGTCACGGGTACGCTCGACCGCCTTCTTCCAGCCGGCATACAGCTTTTCTTTCTGCGTTTCGTCCAGTTGCGGGCt                                     >  7:315925/1‑149 (MQ=255)                                                                                                              aGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTCCCAGTCACGGGTACGCTCGACCGCCTTCTTCCAGCCGGCATACAGCTTTTCTTTCTGCGTTTCGTCCAGTTGCGGGCTg                                    <  8:116416/149‑1 (MQ=255)                                                                                                                 ggATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTCCCAGTCACGGGTACGCTCGACCGCCTTCTTCCAGCCGGCATACAGCTTTTCTTTCTGCGTTTCGTCCAGTTGCGGGCTGAAc                                 >  5:253821/1‑149 (MQ=255)                                                                                                                   aTG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTCCCAGTCACGGGTACGCTCGACCGCCTTCTTCCAGCCGGCATACAGCTTTTCTTTCTGCGTTTCGTCCAGTTGCGGGCTGAACTc                               <  3:175130/149‑1 (MQ=255)                                                                                                                     g‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTCCCAGTCACGGGTACGCTCGACCGCCTTCTTCCTGCCGGCATACAGCTTTTCTTTCTGCGTTTCGTCCAGTTGCGGGCTGAACTCAc                             >  1:280991/1‑149 (MQ=255)                                                                                                                        cGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTCCCAGTCACGGGTACGCTCGACCGCCTTCTTCCAGCCGGCATACAGCTTTTCTTTCTGCGTTTCGTCCAGTTGCGGGCTGAACTCACGc                           <  4:129583/149‑1 (MQ=255)                                                                                                                          gCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTCCCAGTCACGGGTACGCTCGACCGCCTTCTTCCAGCCGGCATACAGCTTTTCTTTCTGCGTTTCGTCCAGTTGCGGGCTGAACTCACGCTc                         <  1:228033/149‑1 (MQ=255)                                                                                                                             aGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTCCCAGTCACGGGTACGCTCGACCGCCTTCTTCCAGCCGGCATACAGCTTTTCTTTCTGCGTTTCGTCCAGTTGCGGGCTGAACTCACGCTCgat                      <  6:253821/149‑1 (MQ=255)                                                                                                                               tAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTCCCAGTCACGGGTACGCTCGACCGCCTTCTTCCAGCCGGCATACAGCTTTTCTTTCTGCGTTTCGTCCAGTTGCGGGCTGAACTCACGCTCgatga                    >  2:132397/1‑149 (MQ=255)                                                                                                                                        ggggCCGCCATGCGGCCCCAACAGCTATCATGCCTCGTGATCTTCCCAGTCACGGGTACGCTCGACCGCCTTCTTCCAGCCGGCATACAGCTTTTCTTTCTGCGTTTCGTCCAGTTGCGGGCTGAACTCACGCTCGATGATCGCCTTGt           >  2:29249/1‑149 (MQ=255)                                                                                                                                         gggCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTCCCAGTCACGGGTACGCTCGACCGCCTTCTTCCAGCCGGCATACAGCTTTTCTTTCTGCGTTTCGTCCAGTTGCGGGCTGAACTCACGCTCGATGATCGCCTTGTc          <  8:315925/149‑1 (MQ=255)                                                                                                                                           gCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTCCCAGTCACGGGTACGCTCGACCGCCTTCTTCCAGCCGGCATACAGCTTTTCTTTCTGCGTTTCGTCCAGTTGCGGGCTGAACTCACGCTCGATGATCGCCTTGTcgc        >  5:354567/1‑149 (MQ=255)                                                                                                                                                 tGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTCCCAGTCACGGGTACGCTCGACCGCCTTCTTCCAGCCGGCATACAGCTTTTCTTTCTGCGTTTCGTCCAGTTGCGGGCTGAACTCACGCTCGCTGATCGCCTTGTCGCGCAGCt  <  6:354567/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                     ACATAGCCGCGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGACAA‑ATTCTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTCCCAGTCACGGGTACGCTCGACCGCCTTCTTCCAGCCGGCATACAGCTTTTCTTTCTGCGTTTCGTCCAGTTGCGGGCTGAACTCACGCTCGATGATCGCCTTGTCGCGCAGCT  >  NC_002947/1233626‑1233916

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