Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,457,578	G→A	7.7%	S371S (TCC→TCT)	bhbP ←	D‑beta‑hydroxybutyrate permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,457,578	0	G	A	7.7%	150.4 / 2.1	52	S371S (TCC→TCT)	bhbP	D‑beta‑hydroxybutyrate permease
Reads supporting (aligned to +/- strand):  ref base G (23/25);  new base A (2/2);  total (25/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.30e-01

CCGCGCCATTGTGCGGCAGGGTGTCCATGCCACCGCTGGCCATGGCGGCTACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAACCCGGCAGCGAGGCGATTACCGCGCCAAAGCCGTATTCCGAAGCGGTGTTCATG  >  NC_002947/3457430‑3457725                                                                                                                                                      |                                                                                                                                                    ccGCGCCATTGTGCGGCAGGGTGTCCATGCCACCGCTGGCCATGGCGGCTACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCg                                                                                                                                                     <  4:305006/149‑1 (MQ=255)  cgcgCCATTGTGCGGCAGGGTGTCCATGCCACCGCTGGCCATGGCGGCTACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCg                                                                                                                                                     >  6:46151/1‑148 (MQ=255)           gtgCGGCAGGGTGTCCATGCCACCGCTGGCCATGGCGGCTACCCGGTGTAGCACTTCCAGCGGGATATCTGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTg                                                                                                                                            <  2:170148/148‑1 (MQ=255)              cGGCAGGGTGTCCATGCCACCGCTGGCCATGGCGGCTACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGc                                                                                                                                        <  8:147751/149‑1 (MQ=255)                gCAGGGTGTCCATGCCACCGCTGGCCATGGCGGCTACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGcc                                                                                                                                       >  5:252678/1‑148 (MQ=255)                            tGCCACCGCTGGCCATGGCGGCTACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTc                                                                                                                          >  3:45825/1‑149 (MQ=255)                               cACCGCTGGCCATGGCGGCTACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACg                                                                                                                       >  4:27623/1‑149 (MQ=255)                                   gCTGGCCATGGCGGCTACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGa                                                                                                                   <  5:251819/149‑1 (MQ=255)                                    cTGGCCATGGCGGCTACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGATTGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGa                                                                                                                   <  1:92314/149‑1 (MQ=255)                                    cTGGCCATGGCGGCTACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGAt                                                                                                                  >  2:193442/1‑149 (MQ=255)                                          aTGGCGGCTACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTc                                                                                                            <  1:193442/149‑1 (MQ=255)                                             gcggcTACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGtt                                                                                                         <  3:27623/149‑1 (MQ=255)                                                    cccGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAg                                                                                                  <  3:319912/149‑1 (MQ=255)                                                        gtgtAGCACTTCCAGCGGGATATCGGCCGCGAGGGCGGCCCCGATGAACTTGTCGCACATGGCTGCCAAGGCGA‑AGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCg                                                                                                                                      >  5:228431/1‑109 (MQ=255)                                                             gCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTgggg                                                                                         <  6:155567/149‑1 (MQ=255)                                                             gCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTggg                                                                                          >  8:186886/1‑148 (MQ=255)                                                             gCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTggg                                                                                          >  8:86143/1‑148 (MQ=255)                                                             gCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTggg                                                                                          <  2:9560/148‑1 (MQ=255)                                                              cACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTgggg                                                                                         <  4:45825/148‑1 (MQ=255)                                                              cACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGa                                                                                        >  4:4664/1‑149 (MQ=255)                                                               aCTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCATGGTGATGGCTTCGTTCACCAAGGGGTTGGGGAt                                                                                       >  1:346913/1‑149 (MQ=255)                                                                         ggATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAgg                                                                             <  3:4664/149‑1 (MQ=255)                                                                          gATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAggg                                                                            <  7:86143/149‑1 (MQ=255)                                                                             atCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCAt                                                                         >  2:26035/1‑149 (MQ=255)                                                                                 gCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCAGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGc                                                                     >  6:18266/1‑149 (MQ=255)                                                                                  ccGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGcc                                                                    <  5:46151/149‑1 (MQ=255)                                                                                         ggcggcggcGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCAc                                                               >  1:318073/1‑147 (MQ=255)                                                                                             gcggcgATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGaa                                                          >  3:128245/1‑148 (MQ=255)                                                                                                  gATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCAGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCa                                                                                                       <  6:24309/98‑1 (MQ=255)                                                                                                  gATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCAGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCa                                                                                                       >  5:24309/1‑98 (MQ=255)                                                                                                      aaCGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAACCCGGCa                                                 <  2:318073/148‑1 (MQ=255)                                                                                                         gtgtCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAACCCGGCAGc                                               >  6:318564/1‑147 (MQ=255)                                                                                                          tgtCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAACCCGGCAGCGAg                                            <  1:26035/149‑1 (MQ=255)                                                                                                                ccATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGa                                                           >  2:296902/1‑128 (MQ=255)                                                                                                                ccATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAACCCGGCAGCGAgg                                           <  4:196169/144‑1 (MQ=255)                                                                                                                ccATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAACCCGGCAGCGAgg                                           >  3:196169/1‑144 (MQ=255)                                                                                                                  aTGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAg                                                            <  2:163059/125‑1 (MQ=255)                                                                                                                  aTGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAg                                                            >  1:163059/1‑125 (MQ=255)                                                                                                                    ggCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAACCCGGCAGCGAgg                                           <  8:300647/140‑1 (MQ=255)                                                                                                                    ggCTGCCAAGGCGA‑GGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAACCCGGCAGCGAgg                                           >  7:300647/1‑140 (MQ=255)                                                                                                                         ccAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCATGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAACCCGGCAGCGAGGCGATTACCgcgc                               <  2:346913/147‑1 (MQ=255)                                                                                                                         ccAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAACCCGGCAGCGAGGCGATTACCGCGcc                              >  5:177715/1‑148 (MQ=255)                                                                                                                          cAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAACCCGGCAGCGAGGCGATTACCGCGCCa                             >  4:151973/1‑148 (MQ=255)                                                                                                                             ggCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAACCCGGCAGCGAg                                            <  8:18088/130‑1 (MQ=255)                                                                                                                             ggCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAACCCGGCAGCGAg                                            >  7:18088/1‑130 (MQ=255)                                                                                                                              gCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAACCCGGCAGCGAGGCGATTACCGCGCCAAAGcc                        <  6:177715/149‑1 (MQ=255)                                                                                                                               cGA‑TGCTCATGCCGCCCGAAGCAGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAACCCGGCAGCGAGGCGATTACCGCGCCAAAGCCg                       <  5:18266/149‑1 (MQ=255)                                                                                                                                               ccGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAACCCGGCAGCGAGGCGATTACCGCGCCAAAGCCGTATTCCGAAGCGgtg        >  1:105681/1‑149 (MQ=255)                                                                                                                                               ccGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAACCCGGCAGCGAGGCGATTACCGCGCCAAAGCCGTATTCCGAAGCGgtg        <  2:105681/149‑1 (MQ=255)                                                                                                                                                  aaGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAACCCGGCAGCGAGGCGATTACCGCGCCAAAGCCGTATTCCGAAGCGGTGtt      >  1:6245/1‑148 (MQ=255)                                                                                                                                                  aaGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAACCCGGCAGCGAGGCGATTACCGCGCCAAAGCCGTATTCCGAAGCGGTGtt      >  4:218866/1‑148 (MQ=255)                                                                                                                                                   aGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAACCCGGCAGCGAGGCGATTACCGCGCCAAAGCCGTATTCCGAAGCGGTGTTCa    <  7:186886/149‑1 (MQ=255)                                                                                                                                                     cGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAACCCGGCAGCGAGGCGATTACCGCGCCAAAGCCGTATTCCGAAGCGGTGTTCATg  <  6:223178/149‑1 (MQ=255)                                                                                                                                                     cGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAACCCGGCAGCGAGGCGATTACCGCGCCAAAGCCGTATTCCGAAGCGGTGTTCATg  <  4:128245/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                    CCGCGCCATTGTGCGGCAGGGTGTCCATGCCACCGCTGGCCATGGCGGCTACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAACCCGGCAGCGAGGCGATTACCGCGCCAAAGCCGTATTCCGAAGCGGTGTTCATG  >  NC_002947/3457430‑3457725

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