Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	348,710	G→A	8.0%	P199P (CCG→CCA)	PP_0288 →	oxalate/formate antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	348,710	0	G	A	8.0%	141.3 / 2.3	50	P199P (CCG→CCA)	PP_0288	oxalate/formate antiporter
Reads supporting (aligned to +/- strand):  ref base G (29/17);  new base A (2/2);  total (31/19)
Fisher's exact test for biased strand distribution p-value = 6.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.49e-01

GCTCGGGCGTGATTGGCGGTATCGGTCTGGGCCTGGGCTACATCTCACCCGTGTCGACCTTGATCAAGTGGTTCCCGGACAAACGCGGCATGGCCACCGGTATGGCGATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTATGGCAGAGCTTTGTCGCCATGGCGGCGATCTACTTCGTGTTCATGACCGCTGGCGCCCTGGCGTATCGCGTACCACCGACCGGCTGGAAGCC  >  NC_002947/348565‑348853                                                                                                                                                  |                                                                                                                                                 gctCGGGCGTGATTGGCGGTATCGGTCTGGGCCTGGGCTACATCTCACCCGTGTCGACCTTGATCAAGTGGTTCCCGGACAAACGCGGCATGGCCACCGGTATGGCGATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTg                                                                                                                                               <  4:164055/149‑1 (MQ=255)      ggCGTGATTGGCGGTATCGGTCTGGGCCTGGGCTACATCTCACCCGTGTCGACCTTGATCAAGTGGTTCCCGGACAAACGCGGCATGGCCACCGGTATGGCGATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑c                                                                                                                                         >  4:241377/1‑149 (MQ=255)       gCGTGATTGGCGGTATCGGTCTGGGCCTGGGCTACATCTCACCCGTGTCGACCTTGATCAAGTGGTTCCCGGACAAACGCGGCATGGCCACCGGTATGGCGATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑cc                                                                                                                                        >  2:288635/1‑149 (MQ=255)                    tATCGGTCTGGGCCTGGGCTACATCTCACCCGTGTCGACCTTGATCAAGTGGTTCCCGGACAAACGCGGCATGGCCACCGGTATGGCGATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATgggg                                                                                                                            >  3:179944/1‑148 (MQ=255)                                      cTACATCTCACCCGTGTCGACCTTGATCAAGTGGTTCCCGGACAAACGCGGCATGGCCACCGGTATGGCGATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCAc                                                                                                          >  3:303081/1‑148 (MQ=255)                                         cATCTCACCCGTGTCGACCTTGATCAAGTGGTTCCCGGACAAACGCGGCATGGCCACCGGTATGGCGATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAgg                                                                                                      <  1:288635/149‑1 (MQ=255)                                                 ccGTGTCGACCTTGATCAAGTGGTTCCCGGACAAACGCGGCATGGCCACCGGTATGGCGATCATGGGCTTCGGCGGCGGGGCCATTGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTg                                                                                                   <  5:178448/144‑1 (MQ=255)                                                 ccGTGTCGACCTTGATCAAGTGGTTCCCGGACAAACGCGGCATGGCCACCGGTATGGCGATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTg                                                                                                   >  6:178448/1‑144 (MQ=255)                                                 ccGTGTCGACCTTGATCAAGTGGTTCCCGGACAAACGCGGCATGGCCACCGGTATGGCGATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTa                                                                                              >  6:305670/1‑149 (MQ=255)                                                              gATCAAGTGGTTCCCGGACAAACGCGGCATGGCCACCGGTATGGCGATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTATGGCAGAGCttt                                                                                  <  4:179944/148‑1 (MQ=255)                                                                tCAAGTGGTTCCCGGACAAACGCGGCATGGCCACCGGTATGGCGATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTATGGCAGAGCTTTGTc                                                                               >  5:69255/1‑149 (MQ=255)                                                                  aaGTGGTTCCCGGACAAACGCGGCATGGCCACCGGTATGGCGATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTATGGCAGAGCTTTGTCGc                                                                             >  8:228239/1‑149 (MQ=255)                                                                    gTGGTTCCCGGACAAACGCGGCATGGCCACCGGTATGGCGATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAgg                                                                                                      >  7:76422/1‑122 (MQ=255)                                                                    gTGGTTCCCGGACAAACGCGGCATGGCCACCGGTATGGCGATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAgg                                                                                                      <  8:76422/122‑1 (MQ=255)                                                                         tCCCGGACAAACGCGGCATGGCCACCGGTATGGCGATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTATGGCAGAGCTTTGTCGCCATggcg                                                                      >  3:216378/1‑149 (MQ=255)                                                                         tCCCGGACAAACGCGGCATGGCCACCGGTATGGCGATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTATGGCAGAGCTTTGTCGCCATggcg                                                                      >  2:164396/1‑149 (MQ=255)                                                                           ccGGACAAACGCGGCATGGCCACCGGTATGGCGATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTATGGCAGAGCTTTGTCGCCATggcggc                                                                    >  6:166893/1‑149 (MQ=255)                                                                             ggACAAACGCGGCATGGCCACCGGTATGGCGATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTATGGCAGAGCTTTGTCGCCATGGCGGCGa                                                                  >  5:50990/1‑149 (MQ=255)                                                                                  aaCGCGGCATGGCCACCGGTATGGCGATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTATGGCAGAGCTTTGTCGCCATGGCGGCGATCTAc                                                             >  8:123091/1‑149 (MQ=255)                                                                                   aCGCGGCATGGCCACCGGTATGGCGATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTATGGCAGAGCTTTGTCGCCATGGCGGCGATCTACt                                                            >  8:184905/1‑149 (MQ=255)                                                                                       ggCATGGCCACCGGTATGGCGATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTAt                                                                                             >  6:387546/1‑112 (MQ=255)                                                                                       ggCATGGCCACCGGTATGGCGATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTAt                                                                                             <  5:387546/112‑1 (MQ=255)                                                                                        gCATGGCCACCGGTATGGCGATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATggg                                                                                                                             <  5:81794/79‑1 (MQ=255)                                                                                        gCATGGCCACCGGTATGGCGATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATggg                                                                                                                             >  6:81794/1‑79 (MQ=255)                                                                                             gCCACCGGTATGGCGATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTATGGCAGAGCTTTGTCGCCATGGCGGCGATCTACTTCGTGTTCa                                                   <  7:184905/148‑1 (MQ=255)                                                                                             gCCACCGGTATGGCGATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTATGGCAGAGCTTTGTCGCCATGGCGGCGATCTACTTCGTGTTCAt                                                  >  2:190099/1‑149 (MQ=255)                                                                                             gCCACCGGTATGGCGATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCACTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTATGGCAGAGCTTTGTCGCCATGGCGGCGATCTa                                                              >  1:75973/1‑137 (MQ=255)                                                                                             gCCACCGGTATGGCGATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCACTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTATGGCAGAGCTTTGTCGCCATGGCGGCGATCTa                                                              >  5:11370/1‑137 (MQ=255)                                                                                             gCCACCGGTATGGCGATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCACTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTATGGCAGAGCTTTGTCGCCATGGCGGCGATCTa                                                              <  6:11370/137‑1 (MQ=255)                                                                                             gCCACCGGTATGGCGATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCACTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTATGGCAGAGCTTTGTCGCCATGGCGGCGATCTa                                                              <  2:75973/137‑1 (MQ=255)                                                                                              ccACCGGTATGGCGATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTATGGCAGAGCTTTGTCGCCATGGCGGCGATCTACTTCGTGTTCATg                                                 >  8:264233/1‑149 (MQ=255)                                                                                               cACCGGTATGGCGATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTATGGCAGAGCTTTGTCGCCATGGCGGCGATCTACTTCGTGTTCATGa                                                <  3:241377/149‑1 (MQ=255)                                                                                                    gTATGGCGATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTATGGCAGAGCTTTGTCGCCATGGCGGCGATCTACTTCGTGTTCATGACCGCt                                           >  3:82607/1‑149 (MQ=255)                                                                                                       tGGCGATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTATGGCAGAGCTTTGTCGCCATGGCGGCGATCTACTTCGTGTTCATGACCGCTGgc                                        <  6:69255/149‑1 (MQ=255)                                                                                                         gCGATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTATGGCAGAGCTTTGTCGCCATGGCGGCGATCTACTTCGTGTTCATGACCGCTGgcgc                                      >  5:294886/1‑149 (MQ=255)                                                                                                           gATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTATGGCAGAGCTTTGTCGCCATGGCGGCGATCTACTTCGTGTTCATGACCGCTGGCGccc                                    >  1:164997/1‑149 (MQ=255)                                                                                                             tCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGCCGTATGGCAGAGCTTTGTCGCCATGGCGGCGATCTACTTCGTGTTCATGACCGCTGGCGCCCTg                                  <  6:50990/149‑1 (MQ=255)                                                                                                                 gggCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTATGGCAGAGCTTTGTCGCCATGGCGGCGATCTACTTCGTGTTCATGACCGCTGGCGCCCTGGCg                               >  6:366369/1‑148 (MQ=255)                                                                                                                  ggCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTATGGCAGAGCTTTGTCGCCATGGCGGCGATCTACTTCGTGTTCATGACCGCTGGCGCCCTGGCGt                              <  7:264233/148‑1 (MQ=255)                                                                                                                   gCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTATGGCAGAGCTTTGTCGCCATGGCGGCGATCTACTTCGTGTTCATGACCGCTGGCGCCCTGGCGTAt                            >  2:67262/1‑149 (MQ=255)                                                                                                                   gCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTATGGCAGAGCTTTGTCGCCATGGCGGCGATCTACTTCGTGTTCATGACCGCTGGCGCCCTGGCGTAt                            >  3:150461/1‑149 (MQ=255)                                                                                                                    cTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTATGGCAGAGCTTTGTCGCCATGGCGGCGATCTACTTCGTGTTCATGACCGCTGGCGCCCTGGCGTATc                           <  5:305670/149‑1 (MQ=255)                                                                                                                           ggcggGGCCATGGTCGGGGCGCCGCTGGCCATCCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTATGGCAGAGCTTTGTCGCCATGGCGGCGATCTACTTCGTGTTCATGACCGCTGGCGCCCTGGCGTATCGCGTac                     >  1:40299/1‑149 (MQ=255)                                                                                                                           ggcggGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTATGGCAGAGCTTTGTCGCCATGGCGGCGATCTACTTCGTGTTCATGACCGCTGGCGCCCTGGCGTATCGCGTac                     >  2:269866/1‑148 (MQ=255)                                                                                                                                 gCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTATGGCAGAGCTTTGTCGCCATGGCGGCGATCTACTTCGTGTTCATGACCGCTGGCGCCCTGGCGTATCGCGTACCaccgac              <  6:294886/149‑1 (MQ=255)                                                                                                                                 gCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTATGGCAGAGCTTTGTCGCCATGGCGGCGATCTACTTCGTGTTCATGACCGCTGGCGCCCTGGCGTATCGCGTACCaccgac              <  1:164396/149‑1 (MQ=255)                                                                                                                                  ccATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTATGGCAGAGCTTTGTCGCCATGGCGGCGATCTACTTCGTGTTCATGACCGCTGGCGCCCTGGCGTATCGCGTACCaccgacc             >  5:272229/1‑149 (MQ=255)                                                                                                                                       gTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTATGGCAGAGCTTTGTCGCCATGGCGGCGATCTACTTCGTGTTCATGACCGCTGGCGCCCTGGCGTATCGCGTACCACCGACCGGCTg        >  8:116818/1‑149 (MQ=255)                                                                                                                                          ggggCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTATGGCAGAGCTTTGTCGCCATGGCGGCGATCTACTTCGTGTTCATGACCGCTGGCGCCCTGGCGTATCGCGTACCACCGACCGGCTGGaa     <  6:272229/149‑1 (MQ=255)                                                                                                                                             gcgcCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTATGGCAGAGCTTTGTCGCCATGGCGGCGATCTACTTCGTGTTCATGACCGCTGGCGCCCTGGCGTATCGCGTACCACCGACCGGCTGGAAGcc  <  7:116818/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                 GCTCGGGCGTGATTGGCGGTATCGGTCTGGGCCTGGGCTACATCTCACCCGTGTCGACCTTGATCAAGTGGTTCCCGGACAAACGCGGCATGGCCACCGGTATGGCGATCATGGGCTTCGGCGGCGGGGCCATGGTCGGGGCGCCGCTGGCCA‑CCGCGCTTATGGGGCATTTTGGCAATGAGCACGAGGTTGGCGTATGGCAGAGCTTTGTCGCCATGGCGGCGATCTACTTCGTGTTCATGACCGCTGGCGCCCTGGCGTATCGCGTACCACCGACCGGCTGGAAGCC  >  NC_002947/348565‑348853

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