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

AGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTGCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC..GCGGGCTCCG.CCCCGGACCCC.CGAGGTATTTCAGGCAAGA.TGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTC  >  NC_014034/1965183‑1965371                                                                                                     |                                                                                               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)                                                                                           gggggCTCCCCCCCCGGACCCC.CCGGGGATTTCAGGCAAGA.TGAAAGAGGGGCGCTCAGGGGGGCTTCTTCGCCCCCCTGCCCCGGGGCTCCGGcttcttt    >  2:3622423/3‑100 (MQ=255)                                                                                            ggggCTCCGCCCCCGGACCCC.CGAGGTATTTCAGGCAAGA.TGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGt   >  1:4333674/2‑101 (MQ=255)                                                                                            ggggCTCCGCCCCCGGACCCC.CGAGGTATTTCAGGCAAGA.TGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGt   <  2:2649412/100‑1 (MQ=255)                                                                                            ggggCTCCGCCCCCGGACCCC.CGAGGTATTTCAGGCAAGA.TGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGt   <  1:1138590/100‑1 (MQ=255)                                                                                             gggCTCCGCCCCCGGACCCC.CGAGGTATTTCAGGCAAGA.TGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTc  >  1:4236644/1‑101 (MQ=255)                                                                                                     |                                                                                               AGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTGCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC..GCGGGCTCCG.CCCCGGACCCC.CGAGGTATTTCAGGCAAGA.TGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTC  >  NC_014034/1965183‑1965371 Base quality scores:  ATCG  < 3 ≤  ATCG  < 21 ≤  ATCG  < 31 ≤  ATCG  < 34 ≤  ATCG  < 40 ≤  ATCG