Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,233,767	A→C	11.2%	intergenic (‑83/+19)	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,767	0	A	C	11.2%	151.7 / 7.7	53	intergenic (‑83/+19)	glpR/glpK	DNA‑binding transcriptional repressor GlpR/glycerol kinase
Reads supporting (aligned to +/- strand):  ref base A (24/23);  new base C (2/4);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

ACATAGCCGCGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGACAA‑ATTCTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTCCCAGTCACGGGTACGCTCGACCGCCTTCTTCCAGCCGGCATACAGCTTTTCTTTCTGCGTTTCGTCCAGTTGCGGGCTGAACTCACGCTCGATGATCGCCTTGTCGC  >  NC_002947/1233626‑1233910                                                                                                                                                |                                                                                                                                                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                                                                                                                         <  2:301239/149‑1 (MQ=255)                  tACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGACAA‑ATTCTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGc                                                                                                                         <  6:229270/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:218495/148‑1 (MQ=255)                                                         cATGGGCGGTCCTTTGGGGCTGCCGGACAA‑ATTCTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTCCCAGTCACGGGTACGCTCGACCGCct                                                                                   <  5:311443/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                                                                                                                      >  8:94169/1‑88 (MQ=255)                                                                                  gACAA‑ATTCTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTc                                                                                                                      <  7:94169/88‑1 (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)                                                                                                                                                |                                                                                                                                                ACATAGCCGCGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGACAA‑ATTCTGCAATGATGCCGTAGGAAGGGGATG‑ACGGCCAGTAGACGATTGGGGCCGCATGGCGGCCCCAACAGCTATCATGCCTCGTGATCTTCCCAGTCACGGGTACGCTCGACCGCCTTCTTCCAGCCGGCATACAGCTTTTCTTTCTGCGTTTCGTCCAGTTGCGGGCTGAACTCACGCTCGATGATCGCCTTGTCGC  >  NC_002947/1233626‑1233910

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