Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	1,965,239	G→T	100%	intergenic (+7/+85)	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,239	0	G→T	100.0%	255.8	69	intergenic (+7/+85)	bdhA/RCAP_rcc01818	3‑hydroxybutyrate dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (32/37):  ref base (0/0):  total (32/37)

TGGGGGGAACGACGGTCTTCCTGTGTGCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTGCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC..GCGGG..CT.CCG.CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCT  >  NC_014034/1965141‑1965334                                                                                                   |                                                                                                       tGGGGGGAACGACGGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGGGGCTGGACGGCGCTGTGAGAGATTtcg                                                                                                      >  1:105138/1‑101 (MQ=255) tGGGGGGAACGACGGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGGGAGATTtcg                                                                                                      >  1:4533767/1‑101 (MQ=255) tGGGGGGAACGACGGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTtcg                                                                                                      <  1:1808053/101‑1 (MQ=255) tGGGGGGAACGACGGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGGGAGAGATTtcg                                                                                                      >  1:4759985/1‑101 (MQ=255)  ggggggAACGACGGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTtcgc                                                                                                     <  1:739263/101‑1 (MQ=255)  ggggggAACGACGGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTtcgc                                                                                                     <  1:4670284/101‑1 (MQ=255)   gggggAACGACGGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGGTGGGCGGGGCTGTGAGAGAGTTCGcc                                                                                                    >  2:2316087/1‑101 (MQ=255)    ggggAACGACGGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCt                                                                                                   >  1:781643/1‑101 (MQ=255)     gggAACGACGGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTc                                                                                                  <  2:277430/101‑1 (MQ=255)     gggAACGACGGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTc                                                                                                  >  2:2856699/1‑101 (MQ=255)      ggAACGACGGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCg                                                                                                 <  1:1497850/101‑1 (MQ=255)      ggAACGACGGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCg                                                                                                 <  2:3714296/101‑1 (MQ=255)      ggAACGACGGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCg                                                                                                 <  1:1110047/101‑1 (MQ=255)        aaCGACGGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCggg                                                                                               <  1:1617637/101‑1 (MQ=255)         acgacgGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCgggg                                                                                              <  1:5039949/101‑1 (MQ=255)         acgacgGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCgggg                                                                                              <  1:4963743/101‑1 (MQ=255)           gacgGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCa                                                                                            <  2:1913274/101‑1 (MQ=255)           gacgGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCa                                                                                            <  2:1254601/101‑1 (MQ=255)             cgGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAg                                                                                          <  1:491609/101‑1 (MQ=255)             cgGTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAg                                                                                          >  1:4269743/1‑101 (MQ=255)              ggTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGc                                                                                         >  2:4036114/1‑101 (MQ=255)              ggTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGc                                                                                         <  1:329093/101‑1 (MQ=255)              ggTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGc                                                                                         >  1:1258133/1‑101 (MQ=255)               gTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTTGACGGCGCTTTGAGAGATTTCGCCTCGGGGCAAGcc                                                                                        >  2:3933926/1‑101 (MQ=255)               gTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGcc                                                                                        >  1:3800793/1‑101 (MQ=255)               gTCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGcc                                                                                        <  1:719275/101‑1 (MQ=255)                tCTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGTCGGCTGGACGGCGCTGGGAGAGATTTCGCCTCGGGGCAAGccc                                                                                       <  2:415195/101‑1 (MQ=255)                 cTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCg                                                                                      <  2:2483767/101‑1 (MQ=255)                 cTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCg                                                                                      >  2:1896804/1‑101 (MQ=255)                 cTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCg                                                                                      <  1:108100/101‑1 (MQ=255)                 cTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCg                                                                                      >  1:580089/1‑101 (MQ=255)                 cTTCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCg                                                                                      >  1:2715301/1‑101 (MQ=255)                  ttCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGa                                                                                     >  1:2729843/1‑101 (MQ=255)                  ttCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGa                                                                                     <  2:3266270/101‑1 (MQ=255)                  ttCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGa                                                                                     <  1:3231068/101‑1 (MQ=255)                  ttCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGa                                                                                     >  1:1220020/1‑101 (MQ=255)                   tCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAg                                                                                    >  1:1133471/1‑101 (MQ=255)                   tCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAg                                                                                    <  2:429954/101‑1 (MQ=255)                   tCCTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAg                                                                                    <  1:4039986/101‑1 (MQ=255)                    ccTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAgg                                                                                   <  1:3459300/101‑1 (MQ=255)                    ccTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAgg                                                                                   >  1:1291662/1‑101 (MQ=255)                    ccTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAgg                                                                                   <  1:70058/101‑1 (MQ=255)                     cTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAggg                                                                                  >  2:1521352/1‑100 (MQ=255)                     cTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAggg                                                                                  >  1:1504611/1‑100 (MQ=255)                     cTGTGTTCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGc                                                                                  <  1:3813334/101‑1 (MQ=255)                      tgtgttCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCg                                                                                 <  2:3372623/101‑1 (MQ=255)                      ngtgttCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCg                                                                                 <  1:2905776/100‑1 (MQ=255)                       gtgttCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGa                                                                                <  2:4836817/101‑1 (MQ=255)                       gtgttCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGa                                                                                >  2:1671730/1‑101 (MQ=255)                       gtgttCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGa                                                                                >  1:2102485/1‑101 (MQ=255)                        tgttCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGAc                                                                               <  1:1772465/101‑1 (MQ=255)                         gttccGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGAcc                                                                              <  2:4533447/98‑1 (MQ=255)                         gttccGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGAcc                                                                              <  2:4527292/98‑1 (MQ=255)                         gttccGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGAcc                                                                              <  2:1216687/98‑1 (MQ=255)                         gatctgAGGCGGCGGCGCAGATCACCGGCACGATGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGAcc                                                                              <  1:4331744/96‑1 (MQ=255)                          ttccGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGAccc                                                                             <  2:806916/99‑1 (MQ=255)                          ttccGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGAccc                                                                             <  2:3180729/99‑1 (MQ=255)                          ttccGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGAccc                                                                             <  2:387420/99‑1 (MQ=255)                          ttccGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGAccc                                                                             <  1:4849598/99‑1 (MQ=255)                           tccGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACccg                                                                            <  2:1708751/100‑1 (MQ=255)                           tccGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACccg                                                                            >  1:1650947/2‑101 (MQ=255)                           tccGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACccg                                                                            <  1:3889394/100‑1 (MQ=255)                            ccGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACccg.c                                                                          <  2:2615570/101‑1 (MQ=255)                              gAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGAccagatc                                                                         >  2:4331427/1‑96 (MQ=255)                              gAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACccg.ccc                                                                        >  1:94836/1‑100 (MQ=255)                              gAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACccg.cc..g                                                                      >  2:2440653/1‑101 (MQ=255)                               aGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACccg.cc..gg                                                                     >  2:4940202/1‑100 (MQ=255)                               aGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCCAGGCGACccggccg                                                                        >  2:998128/1‑97 (MQ=255)                                ggcggcggcGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACccg.cccgc                                                                      >  1:2072929/1‑98 (MQ=255)                                   ggcggcgCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACccg.cccggggg                                                                   >  2:4317165/1‑95 (MQ=255)                                   ggcggcgCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC..GCgggg                                                                 >  1:1154506/1‑100 (MQ=255)                                      ggcgCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC..GGgggggcc                                                              >  1:2286917/1‑97 (MQ=255)                                                              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)                                                                                                   |                                                                                                       TGGGGGGAACGACGGTCTTCCTGTGTGCCGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTGCGCCTCGGGGCAAGCCCGAGGCGACCCG.CC..GCGGG..CT.CCG.CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCT  >  NC_014034/1965141‑1965334 Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG