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%	45.0	17	intergenic (+39/+53)	bdhA/RCAP_rcc01818	3‑hydroxybutyrate dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (3/14):  ref base (0/0):  total (3/14)

CAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTGCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC..GCGGGCTC.CG.CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTC  >  NC_014034/1965182‑1965366                                                                                             |                                                                                                  cAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC..CCggctggccgtt                                                                                        >  1:3827920/1‑92 (MQ=38)  aGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGAGCCG.CC..GCGGGGGC.CG.cccc                                                                                     >  2:718153/1‑101 (MQ=39)    aTCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCGCCC..GCggggggggccccc                                                                                      >  1:1027284/1‑91 (MQ=38)    aTCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC..GGGGGGGG.CG.cccccc                                                                                   >  2:2817090/1‑99 (MQ=37)            cACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC..GGGGGGGC.CC.cccccccgcccccc                                                                           >  2:4430244/1‑91 (MQ=25)                acgaTTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCG.CCG.GGGGGCGC.CC.CCCCCGCccccccccgg                                                                        >  2:4631117/1‑96 (MQ=21)                    tgtCTGTGGGGGGGGGCGGGCCGGCGCTGTGAGATTTTCCGCCGCGGGGCAACCCCGGGGCGCCCCG.CC..GGGGGCTC.CGCCCCCCGACCCCCGAGGTAttt                                                                    <  1:4147859/99‑1 (MQ=12)                       gtggggggggcggcGGGGCGGCGCGGTGAGAGTTTTCCCCTGGGGGCACCCCGGGGGCCCCCCG.CC..GGGGGCTC.CGCCCCCGGACCCCCGAGGTATTTCAg                                                                 <  2:530586/93‑1 (MQ=9)                                 cggcTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC..GGGGGGTC.CGCCCCCGGACCCCCGAGGGATTTCAGGCAAGATGaa                                                       >  2:2331316/1‑101 (MQ=21)                                          gggcgtgggagagtTTGCGCCCGGGGGCACGCCCGGGGCCCCCCG.CC..GGGGGCTC.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGgcg                                              <  1:673868/93‑1 (MQ=16)                                          ggCGGTGTGGGAGTTTGCGCCGGGGGCCACCCCCGGGGCGCCCCG.CC..GGGGGCTC.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGgcg                                              <  1:2150355/101‑1 (MQ=12)                                                gtgAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCGCCC..GGGGGCCCGCC.CCCCCGCCCCCCCCGGGGTTTCGGGCAAAGGGAAAGAGGGgggccccc                                         >  1:2093022/1‑97 (MQ=2)                                                       ttgcccccggggggCACCCCCGGGGCCCCCCG.CC..GGGGGCTC.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCctt                                 <  1:1177988/91‑1 (MQ=18)                                                         ttcgcCTCGGGGCAAGCCCGAGGCGACCCG.CCGGGGGGGCCC.CC.CCCCGGCCCCCCCGGGAATTT.AGGCAAAATGAAAAAGGGGGGCCCAGGCGGCcctct                               >  1:1473568/3‑101 (MQ=12)                                                          cggcgccGGGGCACCCCCGCGCCCCCCCG.CC..GGGGGCTC.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCcttctt                              <  2:1781449/95‑1 (MQ=18)                                                                 ggggCAAGCCCGAGGCGACCCG.CCCGGGGGGGCC.CC.CCCCCGCCCCCCCGGGGATTTTAGGAAAAAGGAAAAAAGGGGGCCCCGGGGGCCTGCTTcccccc                        >  2:1999341/1‑99 (MQ=2)                                                                  gggcgcgCCCGGGGCCCCCCG.CC..GGGGGCTC.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCt                      <  2:4456393/95‑1 (MQ=21)                                                                         ccggggcccccccG.CC..GGGGGCTC.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCgg               <  2:1473474/91‑1 (MQ=25)                                                                         cccacggcaccccG.CC..GGGGGCTC.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCgg               <  2:425763/91‑1 (MQ=25)                                                                                cgccccG.CC..GGGGGCTC.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCg        <  1:4202949/98‑1 (MQ=37)                                                                                ccccccG.CC..GGGGGCTC.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCg        <  1:2537164/98‑1 (MQ=37)                                                                                  ccccG.CC..GGGGGCTC.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCggg      <  1:2817255/100‑1 (MQ=255)                                                                                   cccG.CC..GGGGGCTC.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGc     <  1:1339214/101‑1 (MQ=255)                                                                                    ccg.cc..gGGGGCTC.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGct    <  2:4305387/101‑1 (MQ=255)                                                                                      g.cc..ggGGGCTC.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGcttc  <  1:1869512/101‑1 (MQ=255)                                                                                             |                                                                                                  CAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTGCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC..GCGGGCTC.CG.CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTC  >  NC_014034/1965182‑1965366 Base quality scores:  ATCG  < 3 ≤  ATCG  < 22 ≤  ATCG  < 31 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG