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%	187.2	51	intergenic (+7/+85)	bdhA/RCAP_rcc01818	3‑hydroxybutyrate dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (23/28):  ref base (0/0):  total (23/28)

GGGGGGAACGACGGTCTTCCTGTGTGCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTGCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC..GCGGGCTCCG.CCCCGGACCCCCGAGGTATTTCAGGCAAGA.TGAAAGAGGGGCGCTCAG  >  NC_014034/1965142‑1965327                                                                                                  |                                                                                              ggggggAACGACGGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTtcgc                                                                                            >  2:856402/1‑101 (MQ=255)  gggggAACGACGGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGcc                                                                                           <  2:16384/101‑1 (MQ=255)   ggggAACGACGGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCt                                                                                          <  2:3245283/101‑1 (MQ=255)     ggAACGACGGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCg                                                                                        <  2:20439/101‑1 (MQ=255)     ggAACGACGGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCg                                                                                        <  1:2956657/101‑1 (MQ=255)      gAACGACGGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCgg                                                                                       <  1:1827122/101‑1 (MQ=255)      gAACGACGGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCgg                                                                                       <  1:3031083/101‑1 (MQ=255)      gAACGACGGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCgg                                                                                       >  1:2649584/1‑101 (MQ=255)       gacgacgGGCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCggg                                                                                      <  2:1337801/100‑1 (MQ=255)       aaCGACGGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCggg                                                                                      >  2:3557073/1‑101 (MQ=255)       aaCGACGGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCggg                                                                                      <  1:1004970/101‑1 (MQ=255)        acggcagTCTTCCTGTGTTCCGAGGCGGCCGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTTGACGGCGCTGTCAGATACTTCGCCTCgggg                                                                                     <  2:1841707/95‑1 (MQ=255)        acgacgGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCgggg                                                                                     <  2:1765125/101‑1 (MQ=255)        acgacgGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCgggg                                                                                     <  2:3804178/101‑1 (MQ=255)        acgacgGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCgggg                                                                                     <  2:3400458/101‑1 (MQ=255)         cgacgGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGc                                                                                    <  2:2397694/101‑1 (MQ=255)         cgacgGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGc                                                                                    >  2:524880/1‑101 (MQ=255)         cgacgGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCCCCTCGGGGc                                                                                    >  1:456536/1‑101 (MQ=255)           acgGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCaa                                                                                  >  2:3013339/1‑101 (MQ=255)            cgGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAg                                                                                 <  1:1780663/101‑1 (MQ=255)             ggTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGc                                                                                <  1:2156780/101‑1 (MQ=255)             ggTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGACGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGc                                                                                <  1:219709/101‑1 (MQ=255)              gTCTTCCTGTGTTCCGAGGCGGTGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGCCTGGACGGCGCTGTGGGAGATTTCGCCTCGGGGCAAGcc                                                                               <  2:759952/101‑1 (MQ=255)              gTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGcc                                                                               >  1:4020905/1‑101 (MQ=255)               tCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGccc                                                                              <  1:2645583/101‑1 (MQ=255)                cTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCg                                                                             >  1:1869557/1‑101 (MQ=255)                 ttCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGa                                                                            >  1:4103750/1‑101 (MQ=255)                 ttCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCAGa                                                                            <  1:3302741/101‑1 (MQ=255)                  tCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAg                                                                           <  2:3466865/101‑1 (MQ=255)                    cTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGc                                                                         <  1:2075682/101‑1 (MQ=255)                     tgtgttCCGAGGCGGCGGCGCAGATCACCGGCACGACGTTTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCg                                                                        <  1:1861211/101‑1 (MQ=255)                     tgtgttCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCg                                                                        <  2:1499291/101‑1 (MQ=255)                      gtgttCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGCGGCGa                                                                       >  2:1577024/1‑101 (MQ=255)                      gtgttCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGa                                                                       >  1:770419/1‑101 (MQ=255)                       tgttCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGAc                                                                      <  2:4263785/101‑1 (MQ=255)                        gttccGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGAcc                                                                     <  2:3570503/98‑1 (MQ=255)                        gttccGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGAcc                                                                     <  2:2822451/98‑1 (MQ=255)                         ttccGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGAccc                                                                    <  1:207583/99‑1 (MQ=255)                         ttccGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCGGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGAccc                                                                    <  1:749270/99‑1 (MQ=255)                            cGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACccg.cc                                                                >  2:3211134/1‑101 (MQ=255)                            cGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACccg.cc                                                                >  1:3624224/1‑101 (MQ=255)                                 cggcggcgCAGATCCCCGGCACGACGTTTTCTGGGGATGGCGGCTGGACGGCGCTGGGAGAGATTTCGCCTGGGGGCAAGCCCGAGGCGCCCCG.CC..GGggg                                                         <  2:3424447/101‑1 (MQ=255)                                  ggcggcgCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGGCCCG.CC..GGgggg                                                        >  2:751653/1‑100 (MQ=255)                                       cgcATATCACCGGCACGACGGTTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGGCCCGCCC..GCgggggcgg                                                    >  2:155574/1‑96 (MQ=255)                                        gcAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC..GCgggggggccc                                                  >  1:364146/1‑94 (MQ=255)                                          aGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTAGCCTCGGGGCAAGCCCGCGGCGACCCG.CC..GGGGGGTCCG.cacc                                               >  2:2136283/1‑101 (MQ=255)                                             tCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC..GGGGGGCCCC.ccccccg                                            >  1:1395224/1‑98 (MQ=37)                                               aCCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC..CCGGGGGCCC.ccccccccc                                          >  1:1939006/1‑96 (MQ=38)                                                ccGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCGCCC..GGGGGCCCCG.cccccggac                                          >  2:1482077/1‑96 (MQ=37)                                                ccGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC..GGGGGGGCCC.cccccccccc                                         >  1:617107/1‑95 (MQ=37)                                                 cGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC..CGGGGGTGCG.ccccccccccc                                        >  2:2511847/1‑94 (MQ=25)                                                                            ggggcgggCGCTGGGAGAGTTTCCCCGGCGGGACCGCCCCGGGGCGCCCCG.CC..GGGGGCTCCGCCCCCGGACCCCCGAGGTATTTCAGGCAAGA.TGAAaga             <  1:2163311/95‑1 (MQ=2)                                                                             tGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC..GGGGGGTCCGCCCCCCGACCCCCGCGGGATTTTCGGCCAGA.TGGAaaag            >  2:2150861/1‑101 (MQ=12)                                                                                     gcTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCG.CCGGGGGGGCCCCC.CCCCCGCGCCCCCAGGAATTTATGGAAAGACAGAAAGAgggggcc      >  2:4147038/1‑98 (MQ=9)                                                                                       tgtgAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC..GGGGGCTCCGCCCCCGGACCCCCGAGGGATTTCAGGCAAGA.TGAAAGAGGGGCGCTCAg  >  2:2198275/1‑101 (MQ=21)                                                                                                  |                                                                                              GGGGGGAACGACGGTCTTCCTGTGTGCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTGCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC..GCGGGCTCCG.CCCCGGACCCCCGAGGTATTTCAGGCAAGA.TGAAAGAGGGGCGCTCAG  >  NC_014034/1965142‑1965327 Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG