Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	1,965,279	+C	100%	intergenic (+47/+45)	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,279	1	.→C	100.0%	62.4	21	intergenic (+47/+45)	bdhA/RCAP_rcc01818	3‑hydroxybutyrate dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (3/17):  ref base (1/0):  total (4/17)

AGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTGCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC..GCGGGCTC.CG.CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCGGC  >  NC_014034/1965183‑1965374                                                                                                      |                                                                                                aGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGAGCCG.CC..GCGGGGGC.CG.cccc                                                                                             >  2:718153/1‑101 (MQ=39)   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)                                                                                          cgggggCTC.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGcttctt        <  2:4379091/98‑1 (MQ=255)                                                                                          cgggggCTC.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGcttctt        <  1:518238/98‑1 (MQ=255)                                                                                             gggCTC.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTc     >  1:3690274/1‑101 (MQ=255)                                                                                               gCTC.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCgg   >  2:1107332/1‑101 (MQ=255)                                                                                                cTC.CGCCCCCGGACCCCCGAGGCATTTCAGGCAAGATGAAAAAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCGGc  >  1:1352948/1‑101 (MQ=38)                                                                                                      |                                                                                                AGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTGCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC..GCGGGCTC.CG.CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCGGC  >  NC_014034/1965183‑1965374 Base quality scores:  ATCG  < 3 ≤  ATCG  < 23 ≤  ATCG  < 31 ≤  ATCG  < 35 ≤  ATCG  < 41 ≤  ATCG