Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	1,965,239	G→T	100%	intergenic (+7/+85)	bdhA → / ← RCAP_rcc01818	3‑hydroxybutyrate dehydrogenase/hypothetical protein

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_014034	1,965,239	0	G→T	100.0%	211.7	55	intergenic (+7/+85)	bdhA/RCAP_rcc01818	3‑hydroxybutyrate dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (18/37):  ref base (0/0):  total (18/37)

TGGGGGGAACGACGGTCTTCCTGTGTGCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTGCGCCTCGGGGCAAGCCCGAGGCGACCCGCC..GCGGGCTCCG.CC.CCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGC  >  NC_014034/1965141‑1965332                                                                                                   |                                                                                                  tGGGGGGAACGACGGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTtcg                                                                                                 <  1:189748/101‑1 (MQ=255)    ggggAACGACGGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCt                                                                                              <  2:1826687/101‑1 (MQ=255)    ggggAACGACGGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCt                                                                                              <  2:1930484/101‑1 (MQ=255)    ggggAACGACGGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCt                                                                                              <  2:2369724/101‑1 (MQ=255)     gggAACGACGGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGGTTTCGGCTc                                                                                             >  2:1838112/1‑101 (MQ=255)     gggAACGACGGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTc                                                                                             <  1:3293075/101‑1 (MQ=255)     gggAACGACGGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTc                                                                                             <  1:2816309/101‑1 (MQ=255)      ggAACGACGGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGGGAGATTTCGCCTCg                                                                                            <  2:4216719/101‑1 (MQ=255)      ggAACGACGGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCg                                                                                            <  2:1140185/101‑1 (MQ=255)         acgacgGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCgggg                                                                                         <  1:1644215/101‑1 (MQ=255)          cgacgGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGc                                                                                        >  2:4336776/1‑101 (MQ=255)          cgacgGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGc                                                                                        <  2:3514520/101‑1 (MQ=255)           gacgGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCa                                                                                       <  2:228164/101‑1 (MQ=255)           gacgGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCa                                                                                       <  1:1999823/101‑1 (MQ=255)           gacgGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCa                                                                                       <  1:2221441/101‑1 (MQ=255)           gacgGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCa                                                                                       >  2:380306/1‑101 (MQ=255)            acgGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCaa                                                                                      >  2:3624738/1‑101 (MQ=255)            acgGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCaa                                                                                      <  2:2324147/101‑1 (MQ=255)             cgGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAg                                                                                     <  2:3170983/101‑1 (MQ=255)             cgGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAg                                                                                     >  2:3689686/1‑101 (MQ=255)               gTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGcc                                                                                   >  2:2462156/1‑101 (MQ=255)               gTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGcc                                                                                   <  1:86186/101‑1 (MQ=255)                tCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGccc                                                                                  <  2:2196700/101‑1 (MQ=255)                tCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGccc                                                                                  <  1:824126/101‑1 (MQ=255)                 cTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCg                                                                                 <  2:4596469/101‑1 (MQ=255)                  ttCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGa                                                                                <  2:3928835/101‑1 (MQ=255)                  ttCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGa                                                                                <  1:1255363/101‑1 (MQ=255)                  ttCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGa                                                                                >  1:1812584/1‑101 (MQ=255)                   tCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAg                                                                               <  1:1383696/101‑1 (MQ=255)                   tCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAg                                                                               <  2:3441377/101‑1 (MQ=255)                    ccTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAgg                                                                              <  2:513996/101‑1 (MQ=255)                    ccTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAgg                                                                              <  2:2891047/101‑1 (MQ=255)                    ccTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAgg                                                                              <  1:3348130/101‑1 (MQ=255)                    ccTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAgg                                                                              <  2:3476186/101‑1 (MQ=255)                     cTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGc                                                                             <  1:566871/101‑1 (MQ=255)                     cTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGc                                                                             <  2:3082505/101‑1 (MQ=255)                      tgtgttCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCg                                                                            >  1:3673338/1‑101 (MQ=255)                       gtgttCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCACGCCAGAGGCgc                                                                           >  2:1091875/1‑100 (MQ=255)                       gtgttCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGa                                                                           <  2:3948864/101‑1 (MQ=255)                         gttccGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGAcc                                                                         <  1:4507892/98‑1 (MQ=255)                         gttccGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGAcc                                                                         <  2:3535799/98‑1 (MQ=255)                         atctcGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGAcc                                                                         <  2:823886/97‑1 (MQ=255)                          ttccGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGAccc                                                                        <  2:1705588/99‑1 (MQ=255)                          ttccGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGAccc                                                                        <  1:4596251/99‑1 (MQ=255)                          ttccGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGAccc                                                                        <  1:3152875/99‑1 (MQ=255)                          ttccGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGAccc                                                                        <  1:3058437/99‑1 (MQ=255)                            ccGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACccgc                                                                      <  1:2443211/101‑1 (MQ=255)                             cGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGGGACccggc                                                                     >  2:2642539/1‑99 (MQ=255)                             cGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATATCGCCTCGGGGCAAGCCCGAGGCGAccagat                                                                     >  1:823962/1‑97 (MQ=255)                                ggcggcggcGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACccgcc..ggg                                                                >  2:580021/1‑99 (MQ=255)                                      ggcgCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTTAGAGATTTCGCCTCGGGGCAAGGCCGAGGCGAGccgcccggggggac                                                            >  2:4298098/1‑92 (MQ=255)                                      ggcgCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGGGCACCGCC...CGGGcgggc                                                         >  2:4244043/1‑97 (MQ=255)                                      ggcgCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACccgcccgggggggc                                                            >  1:364286/1‑92 (MQ=39)                                        cgcAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCGCC..GGgggggcgcc                                                        >  2:365681/1‑95 (MQ=38)                                              tCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCGCC..GGGGGGGCCC.cc.ccccc                                                >  1:1399186/1‑98 (MQ=37)                                                  cGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAAGCGACCCGCC..GCGGGGGCCC.cc.cccccgccc                                            >  2:4261370/1‑94 (MQ=255)                                                        gtcgttttGTGTGGGGGGCGGCGGGCCGGCGCTGTGAGAGTTTCCGCGGCGGGCCCACCCGGGGCCGACCCGCC..GGGGGCTCCGCCC.CCGGACCCCCGAgg                                       <  2:3753602/96‑1 (MQ=9)                                                                  gtgGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCGCCGGGGGGGCCCCC.CC.CCGGCCCCCCCGGGGAttttggg                              >  1:2204294/1‑97 (MQ=16)                                                                          cggcTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCGCC..GGGGGCTCCG.CCACCGGACCCCCGAGGTATTTCAGGCAAGATGaa                     >  2:2817148/1‑101 (MQ=25)                                                                                          tgAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCGCC..GGGGGGTCCGCCC.CCCGACCCCCCAGGGATTTTCGGGAAGATGAAAAAgggggggccaggg     >  2:2971960/1‑92 (MQ=12)                                                                                             gagaTTTCGCCTCGGGGCAAGCCCGAGGCGACCCGCC..GGGGGCTCCGCCC.CCGGACCCCCGAGGGATTTCAGGCAAGATGAAAGAGGGGGGCGCAggcggc  >  2:366990/1‑101 (MQ=21)                                                                                                   |                                                                                                  TGGGGGGAACGACGGTCTTCCTGTGTGCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTGCGCCTCGGGGCAAGCCCGAGGCGACCCGCC..GCGGGCTCCG.CC.CCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGC  >  NC_014034/1965141‑1965332 Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 28 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG