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%	100.7	31	intergenic (+47/+45)	bdhA/RCAP_rcc01818	3‑hydroxybutyrate dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (9/21):  ref base (1/0):  total (10/21)

GATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTGCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC...GCGGG..CT..CCG..CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCGGC  >  NC_014034/1965184‑1965374                                                                                                         |                                                                                                 gATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCGCCC...GGGGG..CG..CCC..cccc                                                                                             >  2:419049/1‑101 (MQ=38)  aTCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC...GGGGG..GG..CCT..cccccg                                                                                           >  2:1708443/1‑99 (MQ=38)   tCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC...CCGGGC.CT..CCC..cccccc                                                                                           >  2:4682411/1‑99 (MQ=38)      ccGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC...GGGGG..GG..CCG..cgccccccgc                                                                                       >  2:3308862/1‑95 (MQ=37)       cGGAACGACGATTTCTGTGGATGCCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCG.CG...GGGGG..CG..CAG..ccccccggacc                                                                                      >  2:939848/1‑94 (MQ=21)           acgacgATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC...GGGGG..GG..CCGC.CCCCCGAccccccg                                                                                   >  2:3781504/1‑99 (MQ=25)               cgaTTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC...GGGGG..GTGCCCC..CCCCGGCccccccggga                                                                                >  2:5703908/1‑96 (MQ=21)                   ttgtgggggtGGCGGGTGGGCGGGGCGGGGAGAGATTCCCCCGCGGCGCAACCCGGGGCCCCCCCG.CC...GGGGG..CT..CCGC.CCCCGGACCCCCGAGGTATTTc                                                                           <  2:4655465/92‑1 (MQ=2)                                        ggCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC...GGGGG..CT..CCGC.CCCCGGACCCCCGAGGGATTTCAGGCAAGATGAAAGAGGGgcg                                                      >  2:2864647/1‑101 (MQ=21)                                        ggCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC...CGGGG..GT..CCGCCCCCCGGACCCCC.AGGGATTTCCGGCAAGATGAAAGAgggggg                                                      >  2:5706962/1‑99 (MQ=18)                                            cTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCG.CCCGGGGGGG..CG..CCG..CCCC.CCCCCCCCAGCCATTTCGGCCAGGAAGAAAGAggggggggc                                                   >  1:4745001/1‑96 (MQ=9)                                                gagagaTTTCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC...CGGGG..CT..CCGC.CCCCCCACCCCCGAGGCATTTCAGGCAAGATGAAAGCAGGGCGCTcgcgcg                                              >  2:3861534/1‑96 (MQ=18)                                                gaaactctccgcCGCGGGACAACCCCGGGCCGCCCCG.CC...GGGGG..CT..CCGC.CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAggcg                                              <  1:3086230/92‑1 (MQ=17)                                                  gagaTTTCGCCTCGGGGCAAGCCCGAGGCGACCCG.CCCG.GGGGG..CC..CCC..CCCCGCGCCCCCCCGGGGTTTCCGGCAGAAGGAAAAAGGGGGGCCCGggcgc                                             >  2:4722940/1‑100 (MQ=2)                                                        ccgcCGCGGGCCACGCGCGGGCCGCCCCG.CC...GGGGG..CT..CCGC.CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCcttctt                                      <  2:2382752/100‑1 (MQ=18)                                                         cccggcGGGCCAGGCCCGGGCCCCCCCG.CC...GGGGG..CT..CCGC.CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCttcttc                                     <  2:3175623/96‑1 (MQ=18)                                                              gggggCACGCCCGGGCCGCCCCG.CC...GGGGG..CT..CCGC.CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTcgccgc                                <  2:655900/100‑1 (MQ=21)                                                               ggggCAAGCCCGAGGCGACCCG.CC...GGGGGGTCC..CCC..CCCCGGACCCCC.CGGGATTTCAGGCAAGATGGAAAAGGGGCGGCCAGGGGGCCTTCTTCgccgcc                               >  2:655000/1‑101 (MQ=16)                                                                           gccccccccG.CC...GGGGG..CT..CCGC.CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGtc                   <  1:714661/95‑1 (MQ=25)                                                                               gccccG.CC...GGGGG..CT..CCGC.CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCgg               <  1:1516874/99‑1 (MQ=255)                                                                               cccccG.CC...GGGGG..CT..CCGC.CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCgg               <  1:4470079/99‑1 (MQ=255)                                                                               cccccG.CC...GGGGG..CT..CCGC.CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCgg               <  1:419093/99‑1 (MQ=255)                                                                               cccccG.CC...GGGGG..CT..CCGC.CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCgg               <  1:2136738/99‑1 (MQ=255)                                                                                ccccG.CC...GGGGG..CT..CCGC.CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCggg              <  1:4105728/100‑1 (MQ=255)                                                                                  ccg.cc...gGGGG..CT..CCGC.CCCCGGACCCCCGAGGGATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGct            >  2:2748579/1‑101 (MQ=37)                                                                                   cg.cc...ggGGG..CTC.CCC..CCCCGGACCCCCTAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGctt           >  1:1871959/1‑101 (MQ=255)                                                                                    g.cc...ggGGG..CT..CCGC.CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAACAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGcttc          <  1:1378221/101‑1 (MQ=37)                                                                                          cggggg..CT..CCGC.CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGcttctt        <  1:5658511/98‑1 (MQ=255)                                                                                           ggggg..CT..CCGC.CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTg       >  2:5725639/3‑101 (MQ=255)                                                                                           ggggg..CT..CCGC.CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTg       <  2:5098241/99‑1 (MQ=255)                                                                                           ggggg..CT..CCGC.CCCCGGACCCCCGAGGGATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTg       >  1:460325/3‑101 (MQ=255)                                                                                           ggggg..CT..CCGC.CCCCGGACCCCCGAGGGATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGGGTCTCCGGCTTCTTg       >  2:631609/3‑101 (MQ=37)                                                                                            gggg..CT..CCGC.CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGt      <  1:1368939/100‑1 (MQ=255)                                                                                            gggg..CT..CCGC.CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGt      <  1:3891713/100‑1 (MQ=255)                                                                                            gggg..CT..CCGC.CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGt      <  1:677869/100‑1 (MQ=255)                                                                                            gggg..CT..CCGC.CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGt      <  2:1139300/100‑1 (MQ=255)                                                                                             ggg..CT..CCGC.CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTc     <  1:4014981/101‑1 (MQ=255)                                                                                             ggg..CT..CCGC.CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTc     >  2:3949289/1‑101 (MQ=255)                                                                                              gg..CT..CCGC.CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCg    >  2:2896167/1‑101 (MQ=255)                                                                                               g..CT..CCGC.CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCgg   <  1:390783/101‑1 (MQ=255)                                                                                                gccT..CCGC.CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCg    <  2:5738478/99‑1 (MQ=255)                                                                                                  cT..CCGC.CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCGTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCGGc  >  1:4887037/1‑101 (MQ=255)                                                                                                         |                                                                                                 GATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTGCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC...GCGGG..CT..CCG..CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCGGC  >  NC_014034/1965184‑1965374 Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG