Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	2,070,633	+G	100%	intergenic (+48/+124)	dcd → / ← guaB	deoxycytidine triphosphate deaminase/inosine‑5'‑monophosphate dehydrogenase

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_014034	2,070,633	1	.→G	100.0%	75.0	28	intergenic (+48/+124)	dcd/guaB	deoxycytidine triphosphate deaminase/inosine‑5'‑monophosphate dehydrogenase
Reads supporting (aligned to +/- strand):  new base (1/27):  ref base (0/0):  total (1/27)

GAACTATCAGGGCCAGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGG.CGCGCGGATGGATTTCGCCGCGGCTTT.GCCTGC.GGCGACCCGCC.GGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAGGGCCGGGGTTTCC  >  NC_014034/2070534‑2070726                                                                                                        |                                                                                              gAAATATCAGGGCCAGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGG.CGCGCGGATGGATTTCGCCGCGGCTTT.GCCTGC.GGCGACccgccc                                                                                               >  2:614375/1‑100 (MQ=255)     tATCAGGGCCAGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGG.CGCGCGGATGGATTTCGCCGCGGCTTT.GCCTGC.GGCGACccgccccccc                                                                                           >  1:17544/1‑96 (MQ=255)     tATCAGGGCCAGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGG.CGCGCGGATGGATTTCGCCGCGGCTTT.GCCTGC.GGCGACCCGCC.ggggg                                                                                          >  2:275596/1‑101 (MQ=255)      aTCAGGGCCAGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGG.CGCGCGGATGGATTTCGCCGCGGCTTT.GCCTGC.GGCGACccgcccccgcc                                                                                          >  1:831956/1‑95 (MQ=255)          gggCCAGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGG.CGCGCGGATGGATTTCGCCGCGGCTTT.GCCTGC.GGCGACCCGCC.CCGGgcgggc                                                                                     >  1:827624/1‑98 (MQ=255)             ccgggggCGGAGGCTCGCCAAGCTTTCCCGCGGGGGGGGGGGGG.CGCGGGGATGGTTTCCGCCGCGGCTTT.GCCGGC.GGCGCCCCGCCGGGGGGCGCTGcc                                                                                   <  2:1056688/98‑1 (MQ=255)              cAGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGG.CGCGCGGATGGATTTCGCCGCGGCTTT.GCCTGC.GGCGACCCGCCGGGGGGCGCTGccc                                                                                  >  1:1070549/1‑101 (MQ=255)                  gccTGAAGCCCGCCAAGCTTTCCGGGGGGGGGGGGGGGG.CGCGGGGTGGGTTTCCGCCGCGGTTTG.CCCGGC.GGCGCCCCGCCGGGGGGCGCTGCCCCCCg                                                                              <  2:1028199/99‑1 (MQ=255)                   gCTGACGCCCGCCACGCTTTCCCGGGGGGGGGGGGGGG.CGGGGGGGGGGTTTCCGCCGGGTTTTT.GCCTGC.GGCGCCCCGCCGGGGGGCGCTGCCCCCCgg                                                                             <  2:4031710/101‑1 (MQ=255)                      gacgccccccAGGCTTTCCCGGGGGGGGGGGGGGG.CGGGGGGAGGGTTTTCGCGGGGGTTTT.GCCGGC.GGCGCCCCGCCGGGGGGCGCTGCCCCCCGGAcc                                                                          <  1:590753/93‑1 (MQ=255)                          cccGCCAAGCTTTCCGGGGGGGGGGGGGGGG.CGGGCGGTTGCTTTCCGCCG.GGTTTT.TCCCGCGGGCCCCCCGCCGGGGGGCGCTGCCCCCCGGAcccccc                                                                      <  2:2136990/99‑1 (MQ=2)                           ccGCCAAGCTTTGCCGGGGGGGGGGGGGGG.CGGGGGGATGGTTTCCGCGGCTGTTTT.CCCTGC.GCCACCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCg                                                                     <  2:3941462/100‑1 (MQ=9)                             cccAAGCTTTTCCGGGGGGGGGGGGGGG.GGCGGGGAGGTTTTCCGCCGCGGTTTC.CCCGGC.GGCGCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCggg                                                                   <  1:3070898/100‑1 (MQ=12)                               cAAGCCTTTCGGGGGGGGGGGGGGGG.CGCGCGGGGGGTTTCCCCCGCGGTTTT.CCCTGC.GGCGCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGat                                                                 <  1:3403679/101‑1 (MQ=16)                                caGCCTTTCCGGGGGGGGGGGGGGGCCCCGGGGATGGTTTCCGCCGGGCTTTT.GCCGGC.GGCGCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGat                                                                 <  2:17531/100‑1 (MQ=12)                                    gttcccggggggggTGGGGGG.CGGGGGGTTGGTTCCCCCCGGGGTTTT.GCCCGC.GGCGCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATAtttt                                                            <  2:827543/91‑1 (MQ=9)                                     tttCCGCGCGGGGTGAGGGG.CGCGCGGATGGATTTCGCCGCGGCTTT.GCCTGC.GGCGACCCGCCGGGGGGGGCGGCCCCCCCGCCCCCCCCGGGCAttttt                                                           >  2:2042368/1‑100 (MQ=25)                                         gggggggggGGGGGGG.GGCGCGGGTTGTTTCCGCCGGGTTTTT.GCCGGC.GGCCCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCa                                                       <  1:122586/96‑1 (MQ=9)                                           ggcgGGGGGGGGGG.CGCGCGTATGGTTCCCGCCGCGTCTTG.GCCGGC.GCCCCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAg                                                     <  2:112545/100‑1 (MQ=12)                                               gggggggggg.GGCGCGGTTGGTTTCCGCCGCGGCTTG.CCCGGC.GGCCCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGa                                                 <  1:2963998/95‑1 (MQ=12)                                                         ggggcggggtgggtTTCCCCGGCGCTTTGGCCGGC.GGCGCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAg                                         <  1:2538679/92‑1 (MQ=12)                                                            cgcgGGTGGTTTCCGCCGCGGTTTC.GCGGGC.GCCGCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGgt                                     <  2:439803/101‑1 (MQ=17)                                                                 ttggttTCCGCCGCGGCTTG.CCCGGC.GGCCCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTg                                <  2:513793/100‑1 (MQ=18)                                                                  ttgttctCGCCGGGGCTTC.GCCGGC.GCCCCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTCGc                               <  1:804434/95‑1 (MQ=17)                                                                  tggttccccccgcGTCTTC.GCCTGC.GGCCCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGc                               <  1:804435/92‑1 (MQ=18)                                                                    gttccccccgcGTCTTC.GCCCGC.GGCGCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCg                             <  1:1263605/94‑1 (MQ=21)                                                                              cttcTTT.GCCTGC.GGCCCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGaa                   <  2:3737758/98‑1 (MQ=25)                                                                                       ccggc.ggccCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAGGGcc           <  1:121700/98‑1 (MQ=25)                                                                                       cccgc.ggccCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAGGGcc           <  2:1728311/98‑1 (MQ=25)                                                                                         cgcgccccccccGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAGGGCCg          <  1:2995924/92‑1 (MQ=25)                                                                                           c.ggcgCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAGGGCCgggg       <  1:2240803/101‑1 (MQ=255)                                                                                               ccccccGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAGGGCCGGGGttt    <  2:1132203/98‑1 (MQ=255)                                                                                               ccccccGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAGGGCCGGGGttt    <  2:3963851/98‑1 (MQ=255)                                                                                                 ccccGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAGGGCCGGGGTTTcc  <  1:1627902/100‑1 (MQ=255)                                                                                                        |                                                                                              GAACTATCAGGGCCAGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGG.CGCGCGGATGGATTTCGCCGCGGCTTT.GCCTGC.GGCGACCCGCC.GGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAGGGCCGGGGTTTCC  >  NC_014034/2070534‑2070726 Base quality scores:  ATCG  < 3 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 35 ≤  ATCG  < 41 ≤  ATCG