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%	91.8	27	intergenic (+47/+45)	bdhA/RCAP_rcc01818	3‑hydroxybutyrate dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (13/14):  ref base (0/0):  total (13/14)

AGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTGC.GCCTCGGGGCAAGCCCGAGGCGACCCG.CC..GCGGGCT.CCG.C.CCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCGGC  >  NC_014034/1965183‑1965374                                                                                                       |                                                                                                 aGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTC.GCCTCGGGGCAAGCCCGAGGCGACCCG.CC..GGGGGGG.GCG.c.ccc                                                                                             >  1:2705529/1‑101 (MQ=38)       ccGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTC.GCCTCGGGGCAAGCCCGAGGCGACCCG.CC..CCGGGGG.GCG.c.cccccgcgc                                                                                       >  1:1178794/1‑95 (MQ=37)                cgaTTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTC.GCCTCGGGGCAAGCCCGAGGCGACCCGCCC..GGGGGTT.CCC.C.CCCCGGccccccggggc                                                                               >  2:1218062/1‑95 (MQ=25)                     tCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTC.GCCTCGGGGCAAGCCCGAGGCGACCCG.CC..GGGGGCT.CCGCC.CCCCGGCCCCCCAGGGATTTCa                                                                          >  2:2595796/1‑101 (MQ=21)                             ggggggcGGGGCGGCGCTGGGAG.TATTGCGGGCGCGGGCCAGGCCCGAGGGCCCCCG.CC..GGGGGCT.CCGCC.CCCGGACCCCCGAGGTATTTCAGGCAAGAt                                                                  <  1:559281/97‑1 (MQ=2)                                      gccGGCGCTGGGAGATTTCGC.CCCGGGGGGCACGCCCGAGCCCCCCCG.CC..GGGGGCT.CCGCC.CCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAggg                                                         <  2:459841/99‑1 (MQ=9)                                            gttgtgAGAGTTTCC.GCCGCGGCGCACGCCCGGGCCCCCCCG.CC..GGGGGCT.CCGCC.CCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTc                                                   <  2:2705323/99‑1 (MQ=16)                                                  agagaTTTC.GCCTCGGGGCAAGCCCGAGGCGACCCG.CCGGGGGGGCC.CCC.C.CCCGGCCCCCCCGGGGGTTTTGGGCAGAAGAAAAGAGGGGGGC.CAggcgg                                             >  2:2464724/1‑101 (MQ=2)                                                            gcCTCGGGGCAAGCCCGAGGCGACCCG.CC..GGGGGCTCCCC.C.CCCCCGCCCCCCCGGGGTTTTAGGCAAGAAGAAAGAGGGGCGGCCAGGCGGCCtgcttct                                    >  2:629683/1‑100 (MQ=12)                                                              ggggggccgcgCCCGGGGGGCCCCG.CC..GGGGGCT.CCGCC.CCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTcgcc                                  <  1:2796953/91‑1 (MQ=21)                                                                       ggCCCGGGCCGCCCCG.CC..GGGGGCT.CCGCC.CCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGccc                         <  2:4075643/100‑1 (MQ=25)                                                                           gggccccccccG.CC..GGGGGCT.CCGCC.CCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGgtg                     <  1:2247624/93‑1 (MQ=25)                                                                             gccccccccG.CC..GGGGGCT.CCGCC.CCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGtc                   <  1:1711898/95‑1 (MQ=25)                                                                             agcagtcacg.cc..gaGGGCT.CCG.CGCCCGGCCCCCCGAGGTATTTCAGGCAAGACGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTACCCGgtgcc                   <  2:4091640/93‑3 (MQ=255)                                                                               cccccccG.CC..GGGGGCT.CCGCC.CCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGtctc                 <  1:1218156/97‑1 (MQ=25)                                                                                 cccccG.CC..GGGGGCT.CCGCC.CCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCgg               <  2:433460/99‑1 (MQ=255)                                                                                    ccg.cc..gGGGGCT.CCGCC.CCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGct            <  1:373268/101‑1 (MQ=255)                                                                                            gggggCT.CCGCC.CCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTg       <  2:830715/99‑1 (MQ=255)                                                                                             ggggCT.CCGCC.CCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGt      >  2:4181341/2‑101 (MQ=255)                                                                                             ggggCT.CCGCC.CCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGt      >  1:1916645/2‑101 (MQ=255)                                                                                              gggCT.CCGCC.CCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTc     <  1:84940/101‑1 (MQ=255)                                                                                              gggCT.CCGCC.CCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTc     >  2:3786807/1‑101 (MQ=255)                                                                                              gggCT.CCGCC.CCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTc     >  1:2069751/1‑101 (MQ=255)                                                                                              gggCT.CCGCC.CCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTc     >  1:1089680/1‑101 (MQ=255)                                                                                               ggCT.CCGCC.CCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCg    >  1:754268/1‑101 (MQ=255)                                                                                               ggCT.CCGCC.CCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCg    >  1:538148/1‑101 (MQ=255)                                                                                               ggCT.CCGCC.CCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCg    <  1:3290677/101‑1 (MQ=255)                                                                                                gCT.CCGCC.CCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCgg   >  2:2252461/1‑101 (MQ=255)                                                                                                gCT.CCGCC.CCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCgg   >  2:3048821/1‑101 (MQ=255)                                                                                                gCT.CCGCC.CCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCgg   >  2:493105/1‑101 (MQ=255)                                                                                                gCT.CCGCC.CCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCgg   >  1:1597587/1‑101 (MQ=255)                                                                                                 cT.CCGCC.CCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCGGc  >  2:1558585/1‑101 (MQ=255)                                                                                                 cT.CCGCC.CCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCGGc  <  2:3480982/101‑1 (MQ=255)                                                                                                       |                                                                                                 AGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTGC.GCCTCGGGGCAAGCCCGAGGCGACCCG.CC..GCGGGCT.CCG.C.CCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCGGC  >  NC_014034/1965183‑1965374 Base quality scores:  ATCG  < 3 ≤  ATCG  < 24 ≤  ATCG  < 31 ≤  ATCG  < 35 ≤  ATCG  < 41 ≤  ATCG