Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	642,510	C→G	7.4%	R159R (CGG→CGC)	acoC ←	acetoin cleaving system dihydrolipoyllysine‑residue acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	642,510	0	C	G	7.4%	173.2 / 4.3	54	R159R (CGG→CGC)	acoC	acetoin cleaving system dihydrolipoyllysine‑residue acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (21/29);  new base G (2/2);  total (23/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.25e-01

GCCGGCCAGGTGGGCCTTGGCGATGTCCAGGTGGTCGAGCAGGGCGAGCACGGTTTCGCTCAACTCATCCAGGTCGCCCCGCTGCAGGGCCTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACCCGTGCACCAGCACCAGCGGTGTACCGCCTTCACCCAGCTCGAACCAGCGCAGCAGGCGCCCGCCTACTTCCGCCTTCTGAGGGGCCGGC  >  NC_002947/642375‑642657                                                                                                                                        |                                                                                                                                                     gccgcccAGGTGGGCCTTGGCGATGTCCAGGTGGTCGAGCAGGGCGAGCACGGTTTCGCTCAACTCATCCAGGTCGCCCCGCTGCAGGGCCTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGcagcc‑a                                                                                                                                        <  1:15296/149‑1 (MQ=255)  ccggccAGGTGGGCCTTGGCGATGTCCAGGTGGTCGAGCAGGGCGAGCACGGTTTCGCTCAACTCATCCAGGTCGCCCCGCTGCAGGGCCTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGcagcc‑ag                                                                                                                                       >  3:63408/1‑149 (MQ=255)              gCCTTGGCGATGTCCAGGTGGTCGAGCAGGGCGAGCACGGTTTCGCTCAACTCATCCAGGTCGCCCCGCTGCAGGGCCTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGtt                                                                                                                           >  3:78913/1‑149 (MQ=255)               ccTTGGCGATGTCCAGGTGGTCGAGCAGGGCGAGCACGGTTTCGCTCAACTCATCCAGGTCGCCCCGCTGCAGGGCCTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGGCGCT‑GGCTGCCGAGCGGCGggtcatcg                                                                                                                          <  3:90181/149‑6 (MQ=255)                    gCGATGTCCAGGTGGTCGAGCAGGGCGAGCACGGTTTCGCTCAACTCATCCAGGTCGCCCCGCTGCAGGGCCTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACa                                                                                                                      <  7:295093/148‑1 (MQ=255)                    gCGATGTCCAGGTGGTCGAGCAGGGCGAGCACGGTTTCGCTCAACTCATCCAGGTCGCCCCGCTGCAGGGCCTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACa                                                                                                                      <  7:283648/148‑1 (MQ=255)                     cGATGTCCAGGTGGTCGAGCAGGGCGAGCACGGTTTCGCTCAACTCATCCAGGTCGCCCCGCTGCAGGGCCTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGc                                                                                                                    <  4:63408/149‑1 (MQ=255)                          tCCAGGTGGTCGAGCAGGGCGAGCACGGTTTCGCTCAACTCATCCAGGTCGCCCCGCTGCAGGGCCTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAgtt                                                                                                               <  2:23966/149‑1 (MQ=255)                            cAGGTGGTCGAGCAGGGCGAGCACGGTTTCGCTCAACTCATCCAGGTCGCCCCGCTGCAGGGCCTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAgttgt                                                                                                             >  8:343178/1‑149 (MQ=255)                             aGGTGGTCGAGCAGGGCGAGCACGGTTTCGCTCAACTCATCCAGGTCGCCCCGCTGCAGGGCCTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAgttgtt                                                                                                            <  4:12891/149‑1 (MQ=255)                             aGGTGGTCGAGCAGGGCGAGCACGGTTTCGCTCAACTCATCCAGGTCGCCCCGCTGCAGGGCCTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAgttgtt                                                                                                            <  2:113567/149‑1 (MQ=255)                                tggtCGAGCAGGGCGAGCACGGTTTCGCTCAACTCATCCAGGTCGCCCCGCTGCAGGGCCTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAg                                                                                                         <  5:368403/149‑1 (MQ=255)                                 ggtCGAGCAGGGCGAGCACGGTTTCGCTCAACTCATCCAGGTCGCCCCGCTGCAGGGCCTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAgg                                                                                                        <  5:324219/149‑1 (MQ=255)                                      aGCAGGGCGAGCACGGTTTCGCTCAACTCATCCAGGTCGCCCCGCTGCAGGGCCTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTcacc                                                                                                   <  3:144374/149‑1 (MQ=255)                                        cAGGGCGAGCACGGTTTCGCTCAACTCATCCAGGTCGCCCCGCTGCAGGGCCTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTcaccac                                                                                                 <  3:72039/149‑1 (MQ=255)                                                gCACGGTTTCGCTCAACTCATCCAGGTCGCCCCGCTGCAGGGCCTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGa                                                                                                          <  1:237126/132‑1 (MQ=255)                                                gCACGGTTTCGCTCAACTCATCCAGGTCGCCCCGCTGCAGGGCCTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGa                                                                                                          >  2:237126/1‑132 (MQ=255)                                                      tttCGCTCAACTCATCCAGGTCGCCCCGCTGCAGGGCCTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACCCGTgcacc                                                                                   <  7:343178/149‑1 (MQ=255)                                                        tCGCTCAACTCATCCAGGTCGCCCCGCTGCAGGGCCTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACCCGTgcacca                                                                                  <  8:22090/148‑1 (MQ=255)                                                              aaCTCATCCAGGTCGCCCCGCTGCAGGGCCTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACCCGt                                                                                        >  1:311177/1‑136 (MQ=255)                                                              aaCTCATCCAGGTCGCCCCGCTGCAGGGCCTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACCCGt                                                                                        <  2:311177/136‑1 (MQ=255)                                                                cTCATCCAGGTCGCCCCGCTGCAGGGCCTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGTAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACCCGTGCACCAGCACCAGCg                                                                         <  8:322998/149‑1 (MQ=255)                                                                cTCATCCAGGTCGCCCCGCTGCAGGGCCTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACCCGTGCACCAGCACCAGCg                                                                         >  3:214277/1‑149 (MQ=255)                                                                 tcatcaAGGTCGCCCCGCTGCAGGGCCTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACCCGTGCACCAGCACCAGCg                                                                         >  7:120052/1‑148 (MQ=255)                                                                   aTCCAGGTCGCCCCGCTGCAGGGCCTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACCCGTGCACCAGCACCAGCg                                                                         <  8:120052/146‑1 (MQ=255)                                                                   aTCCAGGTCGCCCCGCTGCAGGGCCTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACCCGTGCACCAGCACCAGCGgtg                                                                      >  3:265427/1‑149 (MQ=255)                                                                         gTCGCCCCGCTGCAGGGCCTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGc                                                                                                                    <  3:139303/97‑1 (MQ=255)                                                                         gTCGCCCCGCTGCAGGGCCTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGc                                                                                                                    >  4:139303/1‑97 (MQ=255)                                                                                cGCTGCAGGGCCTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACCCGTGCACCAGCACCAGCGGTGTACCGCCTTCAcc                                                         <  6:149545/149‑1 (MQ=255)                                                                                     cAGGGCCTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACCCGTGCACCAGCACCAGCGGTGTACCGCCTTCACCCAGCt                                                    <  4:78913/149‑1 (MQ=255)                                                                                     cAGGGCCTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACCCGTGCACCAGCACCAGCGGTGTACCGCCTTCACCCAGCt                                                    >  6:278087/1‑149 (MQ=255)                                                                                         gCCTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACCCGTGCACCAGCACCAGCGGTGTACCGCCTTCACCCAGCTCGaa                                                >  8:44058/1‑149 (MQ=255)                                                                                          ccTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACCCGTGCACCAGCACCAGCGGTGTACCGCCTTCAccc                                                        >  3:97389/1‑140 (MQ=255)                                                                                          ccTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACCCGTGCACCAGCACCAGCGGTGTACCGCCTTCAccc                                                        <  4:97389/140‑1 (MQ=255)                                                                                          ccTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACCCGTGCACCAGCACCAGCGGTGTACCGCCTTCACCCAGCTCGAAc                                               <  4:214277/149‑1 (MQ=255)                                                                                          ccTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACCCGTGCACCAGCACCAGCGGTGTACCGCCTTCACCCAGCTCGAAc                                               <  3:122081/149‑1 (MQ=255)                                                                                          ccTTGGCCGACTCGCCGTGGCCCGCCAGGTCGAGGGCGAAGACCCGCCGCTCGGAAGCC‑AGAGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACCCGTGCACCAGCACCAGCGGAGTACCGCCTTCACCCAGCTCGAAc                                               >  1:141464/1‑149 (MQ=255)                                                                                           cTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACCCGTgcaccag                                                                                 >  1:53781/1‑114 (MQ=255)                                                                                           cTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACCCGTgcaccag                                                                                 <  2:53781/114‑1 (MQ=255)                                                                                                 cGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGGCGCT‑GGCTGCCGAGCGGCggg                                                                                                                                <  1:319627/61‑1 (MQ=18)                                                                                                 cGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGGCGCT‑GGCTGCCGAGCGGCggg                                                                                                                                >  2:319627/1‑61 (MQ=18)                                                                                                               ccGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACCCGTGCACCAGCACCAGCGGTGTACCGCCTTCACCCAGCTCGAACCAGCGCAGCAGGcgcccgcc                           >  3:288635/1‑148 (MQ=255)                                                                                                                cGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACCCGTGCACCAGCACCAGCGGTGTACCGCCTTCACCCAGCTCGAACCAGCGCAGCAGGCGCCCGCCTa                         >  6:114996/1‑149 (MQ=255)                                                                                                                cGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACACGTGCACCAGCACCAGCGGTGTACCGCCTTCACCCAGCTCGAACCAGCGCAGCAGGCGCCCGCCTa                         >  4:365141/1‑149 (MQ=255)                                                                                                                    aGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACCCGTGCACCAGCACCAGCGGTGTACCGCCTTCACCCAGCTCGAACCAGCGCAGCAGGCGCCCGCCTACTTc                     >  4:287379/1‑149 (MQ=255)                                                                                                                          aGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACCCGTGCACCAGCACCAGCGGTGTACCGCCTTCACCCAGCTCGAACCAGCGCAGCAGGCGCCCGCCTACTTCCGCCtt               <  3:287379/149‑1 (MQ=255)                                                                                                                           gggCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACCCGTGCACCAGCACCAGCGGTGTACCGCCTTCACCCAGCTCGAACCAGCGCAGCAGGCGCCCGCCTACTTCCGCCTTc              <  6:319161/149‑1 (MQ=255)                                                                                                                             tcaaccaCCCGGCGCT‑GGCTGCCGAGCGGCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACCCGTGCACCAGCACCAGCGGTGTACCGCCTTCACCCAGCTCGAACCAGCGCAGCAGGCGCCCGCCTACTTCCGCCTTCt             >  5:134791/7‑148 (MQ=255)                                                                                                                               gATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACCCGTGCACCAGCACCAGCGGTGTACCGCCTTCACCCAGCTCGAACCAGCGCAGCAGGCGCCCGCCTACTTCCGCCTTCTGAg          >  3:2767/1‑149 (MQ=255)                                                                                                                               gATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACCCGTGCACCAGCACCAGCGGTGTACCGCCTTCACCCAGCTCGAACCAGCGCAGCAGGCGCCCGCCTACTTCCGCCTTCTGAg          >  2:256271/1‑149 (MQ=255)                                                                                                                                aTGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACCCGTGCACCAGCACCAGCGGTGTACCGCCTTCACCCAGCTCGAACCAGCGCAGCAGGCGCCCGCCTACTTCCGCCTTCTGAgg         <  5:114996/149‑1 (MQ=255)                                                                                                                                 tGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACCCGTGCACCAGCACCAGCGGTGTACCGCCTTCACCCATCTCGAACCAGCGCAGCAGGCGCCCGCCTACTTCCTCCTTCTGAgg         <  2:141464/148‑1 (MQ=255)                                                                                                                                  gACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACCCGTGCACCAGCACCAGCGGTGTACCGCCTTCACCCAGCTCGAACCAGCGCAGCAGGcgcc                               >  6:144342/1‑125 (MQ=255)                                                                                                                                  gACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACCCGTGCACCAGCACCAGCGGTGTACCGCCTTCACCCAGCTCGAACCAGCGCAGCAGGcgcc                               <  5:144342/125‑1 (MQ=255)                                                                                                                                       gccgcTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACCCGTGCACCAGCACCAGCGGTGTACCGCCTTCACCCAGCTCGAACCAGCGCAGCAGGCGCCCGCCTACTTCCGCCTTCTGAGGggccggc  <  4:288635/149‑1 (MQ=255)                                                                                                                                       gccgcTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACCCGTGCACCAGCACCAGCGGTGTACCGCCTTCACCCAGCTCGAACCAGCGCAGCAGGCGCCCGCCTACTTCCGCCTTCTGAGGggccggc  >  1:194506/1‑149 (MQ=255)                                                                                                                                        |                                                                                                                                                     GCCGGCCAGGTGGGCCTTGGCGATGTCCAGGTGGTCGAGCAGGGCGAGCACGGTTTCGCTCAACTCATCCAGGTCGCCCCGCTGCAGGGCCTTGGCCGACTCGCCGTGGCCCGGCAGGTCGAGGGCGATGACCCGCCGCTCGGCAGCC‑AGCGCCGGGTGGTTGAACAGCCAGTTGTTGAGGTCACCACCAAACCCGTGCACCAGCACCAGCGGTGTACCGCCTTCACCCAGCTCGAACCAGCGCAGCAGGCGCCCGCCTACTTCCGCCTTCTGAGGGGCCGGC  >  NC_002947/642375‑642657

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