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

CGAACTATCAGGGCCAGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGGCGCGCGGATGGATTT.CGCCGCGGCTTTGCCTGCGGCGACCCGCC.GGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAGGGCCGGGGTTTCC  >  NC_014034/2070533‑2070726                                                                                                       |                                                                                              cGAACTATCAGGGCCAGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGGCGCGCGGATGGATTT.CGCCGCGGCTTTGCCTGCGGCGACccgcc                                                                                                >  2:2677627/1‑101 (MQ=255) cGAACTATCAGGGCCAGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGGCGCGCGGATGGATTT.CGCCGCGGCTTTGCCTGCGGCGACccgcc                                                                                                >  2:1807962/1‑101 (MQ=255) cGAACTATCAGGGCCAGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGGCGCGCGGATGGATTT.CGCCGCGGCTTTGCCTGCGGCGACccgcc                                                                                                >  2:1787483/1‑101 (MQ=255)  gAACTATCAGGGCCAGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGGCGCGCGGATGGATTT.CGCCGCGGCTTTGCCTGCGGCGACccgcc.g                                                                                              >  2:2058854/1‑101 (MQ=255)    aCTATCAGGGCCAGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGGCGCGCGGATGGATTT.CGCCGCGGCTTTGCCTGCGGCGACccgccccc                                                                                             >  1:4138933/1‑98 (MQ=255)    aCTATCAGGGCCAGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGGCGCGCGGATGGATTT.CGCCGCGGCTTTGCCTGCGGCGACCCGCC.cgg                                                                                            >  1:306271/1‑101 (MQ=255)      tATCAGGGCCAGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGGCGCGCGGATGGATTT.CGCCGCGGCTTTGCCTGCGGCGACCCGCC.ggggg                                                                                          >  2:1528964/1‑101 (MQ=255)      tATCAGGGCCAGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGGCGCGCGGATGGATTT.CGCCGCGGCTTTGCCTGCGGCGACCCGCC.cgggg                                                                                          >  2:3459485/1‑101 (MQ=255)      tATCAGGGCCAGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGGCGCGCGGATGGATTT.CGCCGCGGCTTTGCCTGCGGCGACCCGCC.cgggg                                                                                          >  1:2843467/1‑101 (MQ=255)       aTCAGGGCCAGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGGCGCGCGGATGGATTT.CGCCGCGGCTTTGCCTGCGGCGACccgccctgcgg                                                                                          >  1:4458899/1‑95 (MQ=255)        tCAGGGCCAGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGGCGCGCGGATGGATTT.CGCCGCGGCTTTGCCTGCGGCGACCCGCC.cgggggg                                                                                        >  1:4489923/1‑99 (MQ=255)        tCAGGGCCAGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGGCGCGCGGATGGATTT.CGCCGCGGCTTTGCCTGCGGCGACCCGCC.CGGGgcg                                                                                        >  2:3146186/1‑101 (MQ=255)          aGGGCCAGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGGCGCGCGGATGGATTT.CGCCGCGGCTTTGCCTGCGGCGACccgccccccgcggg                                                                                       >  1:1240389/1‑92 (MQ=255)             gCCAGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGGCGCGCGGATGGATTT.CGCCGCGGCTTTGCCTGCGGCGACCCGCC.cggggggggggc                                                                                   >  2:3092251/1‑94 (MQ=255)                aGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGGCGCGCGGATGGATTT.CGCCGCGGCTTTGCCTGCGGCGACCCGCC.cgggggggggggggg                                                                                >  1:3409251/1‑91 (MQ=255)                aGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGGCGCGCGGATGGATTT.CGCCGCGGCTTTGCCTGCGGCGACCCGCC.cgggggggggggcgc                                                                                >  2:235425/1‑91 (MQ=255)                      tgaggccccccAAGCTTTCCGGCGGGGGGGGGGGGGCGCGGGGTTGGTTTC.CGCCGCGGTTTTGCCGGCGGCGCCCCGCCGGGGGGCGCTGCCCCCCGGAc                                                                           <  1:1498140/92‑1 (MQ=255)                       gacgcccGCCAAGCTTTCCGGCGGGGGGGGGGGGGCGCGCGGATGGTTTT.CGCCGCGGTTTTCCCTGGGGCCCCCCGCCGGGGGGCGCTGCCCCCCGGAcc                                                                          <  2:1333442/95‑1 (MQ=255)                       gaaccccGCCAAGCTGTGCCGCGGGGGGGGGGGGGGGCGCGGTTGGTTTC.CGCCGGGGCTTTCCCGGGGGCCCCCCGCCGGGGGGCGCTGCCCCCCGGAcc                                                                          <  2:596199/95‑1 (MQ=255)                       aacgtccGCAAAGCTTTCCCGGGGGGGGAGGGGGGCGCGGGGATGTATTC.CCCCTCGCTTTTCCCGGGGGCGCCCCGCCGGGGGGCGCTGCCCCCCGGAcc                                                                          <  2:391066/95‑1 (MQ=255)                          ccccGCCAAGCTTTCCGGGGGGGGGGGGGGGGCGCGGGGGTGGTTTC.CGCCGCGTTTTTGCCGGCGGCCCCCCGCCGGGGGGCGCTGCCCCCCGGAccccc                                                                       <  1:2058975/98‑1 (MQ=12)                          ccTCCCCAACCTTTCCCGCGGGGGGTGGGGGGCCCGC.GTTGGTTTTCCGCCCCGGTTTTCCCGGCGCCGCCCCGCCGGGGGGCGCTGCCCCCCGGAccccc                                                                       <  1:1334402/100‑1 (MQ=9)                                  aGCTTTCCCGGGGGGGGGGGGGGGCGCGCGGATGGTTTC.CGCCGCGGTTTTCCCGGCGGCCCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGatat                                                               <  2:120073/101‑1 (MQ=16)                                          gggggggggGGGGGGGCGGGGGGTTGGTTTC.CCCCGCGGTTTGCCCTGGGGCCCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCa                                                       <  1:1339921/96‑1 (MQ=9)                                          cgggtggggGGGGGGGCGCGGGGTTGGATTC.CCCGTCTGCTCTCCCTGCGCCCCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCa                                                       <  2:4472893/96‑1 (MQ=9)                                            gggggggGGGGGGGCGCGGGGGGGGTTTC.CCCCGCGCCTTGGCCTGCGCCCCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAg                                                     <  2:306237/98‑1 (MQ=9)                                                 gggggggggCGCGGGGATGGTTTC.CGCCGCGGCTTCCCCCGCGGCCCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGaa                                                <  1:2412220/96‑1 (MQ=18)                                                                ggtggttTC.CGCCGCGGCTGCGCCGGCGCCCCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGttt                                 <  2:4642334/99‑1 (MQ=18)                                                                        tccgccgcGGTTTGGCGGGCGGCCCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTa                          <  2:2510614/99‑1 (MQ=18)                                                                        tcccccgcGGCTTTCCCGGCGCCCCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTa                          <  1:479369/97‑1 (MQ=21)                                                                              gcgGCTTTGCCCGCGGCCCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGa                    <  1:1084083/101‑1 (MQ=25)                                                                                ggCTTTGCCTGCGGCGACCCGCCGGGGGGGGCGGCCCCCCGGCCCCCCCCGCCCATTTTCGGCAAGAGGAAGGGGAAGCGGGGTTGGGGGGTACGCAGaaa                  >  1:611057/1‑101 (MQ=255)                                                                                       cgcgcgcccccccGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAGGGcc           <  1:2879283/91‑1 (MQ=25)                                                                                        ccgcgccgCCCCGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAGGGCCg          <  2:3409045/99‑1 (MQ=25)                                                                                         cgcgcccccccGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAGGGCCgg         <  2:4138687/93‑1 (MQ=25)                                                                                           cgcccccccGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAGGGCCgggg       <  1:27549/95‑1 (MQ=25)                                                                                             gccccccGCCGGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAGGGCCGGGGtt     <  2:1485700/97‑1 (MQ=255)                                                                                                cgccgccgGGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAGGGCCGGGGTTTcc  <  2:4674385/99‑1 (MQ=255)                                                                                                       |                                                                                              CGAACTATCAGGGCCAGGGGCTGAAGCTCGCCAAGCATTTCCGCGCGGGGTGAGGGGCGCGCGGATGGATTT.CGCCGCGGCTTTGCCTGCGGCGACCCGCC.GGGGGCGCTGCCCCCCGGACCCCCCGGGATATTTTCGGCAAGATGAAGGGGAAGCGGTGTTTGCCGGTACACAGAAAAAGGGCCGGGGTTTCC  >  NC_014034/2070533‑2070726 Base quality scores:  ATCG  < 3 ≤  ATCG  < 23 ≤  ATCG  < 31 ≤  ATCG  < 35 ≤  ATCG  < 41 ≤  ATCG