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

TTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTGCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC..GCGGG..CT.CCG.CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCGGC  >  NC_014034/1965202‑1965374                                                                                     |                                                                                                tggggtggggggcggcGGGCGGGCGCGGGGAGAGTCGCCGCCCGGGGGCACCCCCGCGCCCCCCCG.CC..GGGGG..CT.CCGCCCCCGGACCCCCGAGGTATTTc                                                                           <  2:1356168/91‑1 (MQ=0)   cTGTGGATGGCGGCTGGACGGCGCTGGGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCGCCC..GGGGGC.CC.CCC.CCCCGGCCCCCCGCGG.ATTTCAg                                                                         >  2:3390175/1‑101 (MQ=16)           ggcggcTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCG.CCCGGGGGG..GC.CCC.CCCCGCCCCCCCCCGGCGTTTCAGGCAAggt                                                                  >  2:4771822/1‑99 (MQ=12)                    aCGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCG.CCGGGGGGG..CC.CCC.CCCCCGACCCCCCCGGGGTTTCAGGCAAGAGGAAAAAggg                                                         >  2:4826728/1‑101 (MQ=16)                              gagagaTTTCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC..GGGGG..CTCCCC.CCCCCGCCCCCCCCGGGGTTTCTGGGAAGATGAGAAGGGGGGGCCCAggcg                                              >  2:3000679/1‑101 (MQ=2)                                   aTTTCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC..GGGGGGCCC.CCC.CCCCCCCCCCCCCGGGGGTTTAAGGCAAAAGAAAAAAGGGGGGCGCAGGGGGccc                                          >  2:1419160/1‑100 (MQ=0)                                            cGGGGCAAGCCCGAGGCGACCCG.CCG.GGGGG..CC.CCC.CCCCCGGCCCCCCAGGGATTTCAGGGCAGAGGAAAAAGGGGGGGCCCGGCGGCCCTCTTcgcccc                                >  1:655753/1‑99 (MQ=12)                                                    gCCCGGGCCGCCCCG.CC..GGGGG..CT.CCGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCg                        <  1:3794863/101‑1 (MQ=25)                                                       gggggccccccG.CC..GGGGG..CT.CCGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGgtg                     <  2:2050063/93‑1 (MQ=25)                                                       gctctccccccG.CC..GGGGG..CT.CCGCCCCCGGTCCCCCGAGGTATTTCAGTCAAGATGAAAGAGGGTCGCTCAGGCGGCCTTTTTCGCCGCCTTGCCCGgtg                     <  2:4937263/93‑1 (MQ=21)                                                       gcacaacccccG.CC..GGGGG..CT.CCGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGgtg                     <  1:1709446/93‑1 (MQ=25)                                                          ccccccccG.CC..GGGGG..CT.CCGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGtct                  <  1:4711328/96‑1 (MQ=25)                                                           cccccccG.CC..GGGGG..CT.CCGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGtctc                 <  2:3226265/97‑1 (MQ=25)                                                           cccccccG.CC..GGGGG..CT.CCGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGtctc                 <  2:2913735/97‑1 (MQ=25)                                                             cccccG.CC..GGGGG..CT.CCGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCgg               <  1:3272171/99‑1 (MQ=255)                                                              ccccG.CC..GGGGG..CT.CCGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCggg              <  1:720428/100‑1 (MQ=255)                                                              ccccG.CC..GGGGG..CT.CCGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCggg              <  1:2170517/100‑1 (MQ=255)                                                              ccccG.CC..GGGGG..CT.CCGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCggg              <  2:154660/100‑1 (MQ=255)                                                               cccG.CC..GGGGG..CT.CCGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGc             <  1:4772145/101‑1 (MQ=255)                                                               cccG.CC..GGGGG..CT.CCGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGc             <  2:1731684/101‑1 (MQ=255)                                                                 cg.cc..ggGGG..CT.CCGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGctt           <  1:4440118/101‑1 (MQ=255)                                                                  g.cc..ggGGG..CT.CCGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGcttc          <  2:4473384/101‑1 (MQ=255)                                                                    cc..ggGGG..CT.CCGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGcttct         <  1:1602031/101‑1 (MQ=255)                                                                        ggggg..GT.CCGCCCCCCGACCCCCCAGGGATTTTAGGCAAGATGAAAGAGGGGGGCTCAGGCGGGCCTCTTCGCCGCCCTGCCCCGGGGCTCGGGgttcttt       >  2:2777972/3‑100 (MQ=255)                                                                         gggg..CT.CCGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGt      <  2:2488534/100‑1 (MQ=255)                                                                         gggg..CT.CCGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGt      >  2:1775634/2‑101 (MQ=255)                                                                         gggg..CT.CCGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGt      >  1:4120829/2‑101 (MQ=255)                                                                         gggg..CT.CCGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGt      <  2:773580/100‑1 (MQ=255)                                                                          ggg..CT.CCGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTc     >  1:4645257/1‑101 (MQ=255)                                                                          ggg..CT.CCGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTc     >  2:4600220/1‑101 (MQ=255)                                                                           gg..CT.CCGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCg    <  1:2791250/101‑1 (MQ=255)                                                                               cT.CCGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCGGc  >  2:1192219/1‑101 (MQ=255)                                                                                     |                                                                                                TTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTGCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC..GCGGG..CT.CCG.CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCGGC  >  NC_014034/1965202‑1965374 Base quality scores:  ATCG  < 3 ≤  ATCG  < 19 ≤  ATCG  < 29 ≤  ATCG  < 34 ≤  ATCG  < 40 ≤  ATCG