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%	94.9	35	intergenic (+48/+124)	dcd/guaB	deoxycytidine triphosphate deaminase/inosine‑5'‑monophosphate dehydrogenase
Reads supporting (aligned to +/- strand):  new base (3/32):  ref base (0/0):  total (3/32)

CGAACTATCAGGGCCAGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGGC.GCGCGGATGGATTTCGCCGCGG.CTTTG.CCTGCGGCGACCCGCC.GGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAGGGCCGGGG  >  NC_014034/2070533‑2070721                                                                                                         |                                                                                         cGAACTATCAGGGCCAGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGGC.GCGCGGATGGATTTCGCCGCGG.CTTTG.CCTGCGGCGACccgcc                                                                                           >  1:3836589/1‑101 (MQ=255)     cTATCAGGGCCAGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGGC.GCGCGGATGGATTTCGCCGCGG.CTTTG.CCTGCGGCGACCCGCC.cggg                                                                                      >  2:1519311/1‑101 (MQ=255)      tATCAGGGCCAGGGGCTGAAGCTCGCCAAGCATTTCCGCGCTGGGTGAGGGGC.GCGCGGATGGATTTCGCCGCGG.CTTTG.CCTGCGGTGACccgccctggg                                                                                      >  2:1902116/1‑96 (MQ=255)      tATCAGGGCCAGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGGC.GCGCGGATGGATTTCGCCGCGG.CTTTG.CCTGCGGCGACccgccccgcg                                                                                      >  1:4764753/1‑96 (MQ=255)       aTCAGGGCCAGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGGC.GCGCGGATGGATTTCGCCGCGG.CTTTG.CCTGCGGCGACccgccccgccc                                                                                     >  1:3976412/1‑95 (MQ=255)       aTCAGGGCCAGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGGC.GCGCGGATGGATTTCGCCGCGG.CTTTG.CCTGCGGCGACccgcccccgcg                                                                                     >  1:3950804/1‑95 (MQ=255)          aGGGCCAGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGGC.GCGCGGATGGATTTCGCCGCGG.CTTTG.CCTGCGGCGACCCGCC.ggggggggg                                                                                 >  2:2469311/1‑97 (MQ=255)          aGGGCCAGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGGC.GCGCGGATGGATTTCGCCGCGG.CTTTG.CCTGCGGCGACCCGCC.cgggggggg                                                                                 >  1:5464348/1‑97 (MQ=255)             gCCAGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGGC.GCGCGGATGGATTTCGCCGCGG.CTTTG.CCTGCGGCGACCCGCC.cggggggggggg                                                                              >  2:2369573/1‑94 (MQ=255)               cAGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGGC.GCGCGGATGGATTTCGCCGCGG.CTTTG.CCTGCGGCGACCCGCCGGGGGGCGCTGcca                                                                             >  2:4066628/1‑100 (MQ=255)                aGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGGC.GCGCGGATGGATTTCGCCGCGG.CTTTG.CCTGCGGCGACCCGCC.cgggggggggggggg                                                                           >  1:5222116/1‑91 (MQ=255)                  gggCTGAGGCGCGCCAGGCGTTCCGGGGCGGGGGGGGGGCG.GCGGGGGGGGTTCCCCCCGCGG.TTTTG.CCGGCGGCGCCCCGCCGGGGGGCGCTGcccccc                                                                          <  1:905966/101‑1 (MQ=255)                   ggcggacccTCGCAAAGCTTTCCCGCGGGGGGGGGGGGGC.GCGGGGATGGTTTCCCCCGCGG.CTTTG.CCGGCGGCGCCCCGCCGGGGGGCGCTGCCCCCCg                                                                         <  1:5647602/93‑1 (MQ=255)                   ggCTGAGGCGCGCCAAGCTTTCCGGGGGGGGGGGGGGGGG.GCGGGGATGG.TTTCCCCGGGGTTTTTG.CCGGCGGCGCCCCGCCGGGGGGCGCTGCCCCCCg                                                                         <  1:2736662/101‑1 (MQ=255)                    ggtGAAGCCCGCAAAGTTTTCCCGCGGGGGGGGGGGGGC.GGGGGGATGTTTTCCGCGGCGG.TTTTC.CCGGGGGCGCCCCGCCGGGGGGCGCTGCCCCCCgg                                                                        <  2:328457/99‑1 (MQ=255)                    gCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGGC.GCGCGGATGGATTTCGCCGCGG.CTTTG.CCTGCGGCGACCCGCCGGGGGGCGCTGcagatcgg                                                                        >  1:5580273/1‑94 (MQ=38)                            ccGCCAAGCTTGCCCGGGGGGGGGGGGGGGG.GCGCGGTTGGTTCCCGCCGCGT.CTTGC.CCTGCGGCGCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCg                                                                <  1:3987742/100‑1 (MQ=12)                            cagcaaaGCA.TTCCGGGGGGGGGGGGGGGCGGCGGGGGTGGTTGCCGCCGCGG.TTTTC.CCGGCGGCCCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCg                                                                <  1:933697/96‑1 (MQ=2)                               ccAAGCTTTCCGGGGGGGGGGGGGGGGG.GCGCGGGTGTTTTCCGCCGCGG.CTTTG.CCTGGGGCCCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGa                                                             <  1:4179958/101‑1 (MQ=12)                                caggcttgcccGCGGGGGGTGGGGGGC.GCGGGGATGTTTTGCGCCGCGC.TTTTG.CCTGCGGCCCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGat                                                            <  2:4698091/92‑1 (MQ=16)                                 agccttccccGCGGGGGGGGGGGGGC.GCGCGGATGGTTTCCCCCGCGT.CTTTG.CGTGGGCCCCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGata                                                           <  2:1506867/93‑1 (MQ=12)                                  aGCATTTCCGCGCGGGGTGAGGGGC.GCGCGGATGGATTTCGCCGCGG.CTTTG.CCTGCGGCGACCCGCCGGGGGGGGCTGCCCCCCGGACCCCCCCGGatat                                                          >  2:3245883/1‑101 (MQ=39)                                          cgcgggGGGTGGGTGGC.GCGCGTTGGGTTTTCGCCGCGT.CTTTG.CCTGCGGCGCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCa                                                  <  2:4764472/101‑1 (MQ=21)                                             ggggggTGGGGGGG.GCGGGGATGGTTGCCGCCGCGG.CTTTC.CCGGGGGCGCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGa                                               <  2:4756300/99‑1 (MQ=16)                                             ggggggGGGGGGGC.GGGCGGTTGGGTCCCGCCGCTG.CTTCC.CCTGCGGCGCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGa                                               <  2:1114872/99‑1 (MQ=12)                                                ggggggggggC.GGGGGGGTGGTTTCCCCCGCGG.CTTGC.CCGGCGGCCCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGa                                            <  1:4627803/95‑1 (MQ=9)                                                    ggggggC.GGGGGGGTGGTTTCCCCCGCGG.TTTTC.CCGGGGGCGCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAggg                                        <  1:4200842/99‑1 (MQ=12)                                                    ggggggC.GCGGGGGTGGTTCCCCCCGCGC.TTGTG.CCGGCGGCCCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAggg                                        <  2:3258629/99‑1 (MQ=9)                                                      ggggC.GCGCGGAGGAGTTACCCGGCGG.CTTCG.CCGGCGGCGCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGa                                      <  2:3836347/101‑1 (MQ=18)                                                            gcgggGGTGGCTTCCCCCGCGG.CTTGC.CCTGCGGCCCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCgg                                 <  2:2621221/101‑1 (MQ=17)                                                                 gtgttttcccgccgcGC.TTTCG.CGTGCGGCGCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGttt                            <  2:4601915/92‑1 (MQ=18)                                                                  gtggttTCCGCCGCGG.CTTTC.CCGGCGGCCCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTg                           <  2:5464058/100‑1 (MQ=18)                                                                     tttcccgccgcGG.CTGCC.CCGGCGGCGCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCg                        <  2:3950560/96‑1 (MQ=18)                                                                            ccgcGG.TTTCG.CCGGCGCCCCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTAcaca                 <  1:1519404/101‑1 (MQ=21)                                                                               ccg.CTTTGCCCGGCGGCCCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGa               <  1:4574516/99‑1 (MQ=21)                                                                                   tgtcg.CCTGCGGCCCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGaaaaa           <  2:5145667/97‑1 (MQ=25)                                                                                     ccg.CCGGCGGCCCCCCGCCGGGGGGCGCTGCCCTCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAgg         <  1:4931681/99‑1 (MQ=25)                                                                                         ccgcgggccccccGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAGGGcc      <  1:3469474/91‑1 (MQ=25)                                                                                          gggcggccCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAGGGCCg     <  1:5449030/99‑1 (MQ=25)                                                                                           ggcgcccccccGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAGGGCCgg    <  2:5221829/93‑1 (MQ=25)                                                                                            ggggccccccGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAGGGCCggg   <  1:5231944/94‑1 (MQ=25)                                                                                            gcgccgCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAGGGCCggg   <  2:4326817/101‑1 (MQ=25)                                                                                             ggcccccccGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAGGGCCgggg  <  1:3246083/95‑1 (MQ=25)                                                                                             ggcccccccGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAGGGCCgggg  <  1:1458026/95‑1 (MQ=25)                                                                                             cccccccccGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAGGGCCgggg  <  1:4173853/95‑1 (MQ=25)                                                                                                         |                                                                                         CGAACTATCAGGGCCAGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGGC.GCGCGGATGGATTTCGCCGCGG.CTTTG.CCTGCGGCGACCCGCC.GGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAGGGCCGGGG  >  NC_014034/2070533‑2070721 Base quality scores:  ATCG  < 3 ≤  ATCG  < 23 ≤  ATCG  < 32 ≤  ATCG  < 35 ≤  ATCG  < 41 ≤  ATCG