Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	1,965,271	C→G	100%	intergenic (+39/+53)	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,271	0	C→G	100.0%	40.1	14	intergenic (+39/+53)	bdhA/RCAP_rcc01818	3‑hydroxybutyrate dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (5/9):  ref base (0/0):  total (5/9)

CGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTGCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC..GCGGGCTCCG.CCCCGGACCCC.CGAGGTATTTCAGGCAAGA.TGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTT  >  NC_014034/1965174‑1965365                                                                                                     |                                                                                                  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..GGGGGCTCCGCCCCCGGACCCC.CGAGGTATTTCAGGCAAGA.TGAAaga                                                   <  1:2163311/95‑1 (MQ=2)                                             tGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC..GGGGGGTCCGCCCCCCGACCCC.CGCGGGATTTTCGGCCAGA.TGGAaaag                                                  >  2:2150861/1‑101 (MQ=12)                                                     gcTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCG.CCGGGGGGGCCCCC.CCCCCGCGCCC.CCAGGAATTTATGGAAAGACAGAAAGAgggggcc                                            >  2:4147038/1‑98 (MQ=9)                                                       tgtgAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC..GGGGGCTCCGCCCCCGGACCCC.CGAGGGATTTCAGGCAAGA.TGAAAGAGGGGCGCTCAg                                        >  2:2198275/1‑101 (MQ=21)                                                                         ggggCAAGCCCGAGGCGACCCG.CCGGGGGGGCCCCC.CCCCCGCCCCCGCGGGGT.TTTCGGGAAAAA.TGAAGAAGGGGGGCCCGGGGGGCcctcttgccccc                       >  2:2629012/1‑95 (MQ=2)                                                                                 gcggggcccccccG.CC..GGGGGCTCCGCCCCCGGACCCC.CGAGGTATTTCAGGCAAGA.TGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCgg              <  1:3619264/91‑1 (MQ=25)                                                                                     cgcccccccG.CC..GGGGGCTCCGCCCCCGGACCCC.CGAGGTATTTCAGGCAAGA.TGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGtc          <  1:785687/95‑1 (MQ=25)                                                                                       cccccccG.CC..GGGGGCTCCGCCCCCGGACCCC.CGAGGTATTTCAGGCAAGA.TGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGtctc        <  1:699280/97‑1 (MQ=25)                                                                                        ccccccG.CC..GGGGGCTCCGCCCCCGGACCCC.CGAGGTATTTCAGGCAAGA.TGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCg       <  2:1395117/98‑1 (MQ=37)                                                                                         cccccG.CC..GGGGGCTCCGCCCCCGGACCCC.CGAGGTATTTCAGGCAAGA.TGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCgg      <  1:3211363/99‑1 (MQ=255)                                                                                           cccG.CC..GGGGGCTCCGCCCCCGGACCCC.CGAGGTATTTCAGGCAAGA.TGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGc    <  2:3981054/101‑1 (MQ=255)                                                                                            ccg.cc..gGGGGCTCCGCCCCCGGACCCC.CGAGGTATTTCAGGCAAGA.TGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGct   <  1:155595/101‑1 (MQ=255)                                                                                             cg.cc..ggGGGCTCCGCCCCCGGACCCC.CGAGGTATTTCAGGCAAGA.TGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGctt  <  2:3623971/101‑1 (MQ=255)                                                                                                     |                                                                                                  CGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTGCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC..GCGGGCTCCG.CCCCGGACCCC.CGAGGTATTTCAGGCAAGA.TGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTT  >  NC_014034/1965174‑1965365 Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 29 ≤  ATCG  < 34 ≤  ATCG  < 40 ≤  ATCG