Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,607,828	C→T	G161G (GGC→GGT)	tam →	trans‑aconitate 2‑methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,607,828	0	C	T	100.0%	72.5 / NA	21	G161G (GGC→GGT)	tam	trans‑aconitate 2‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (8/13);  total (8/13)

GAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCA  >  NC_000913/1607751‑1607902                                                                              |                                                                           gAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGTGtt                                                                         <  2:297631/81‑1 (MQ=255)    ccGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGTGTTCAt                                                                      <  1:297632/81‑1 (MQ=255)     cGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGTGTTCATg                                                                     <  1:297633/81‑1 (MQ=255)     cGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGTGTTCATg                                                                     >  1:297676/1‑81 (MQ=255)      gACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGTGTTCATGc                                                                    >  1:297663/1‑81 (MQ=255)            tgtgCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGTGTTCATGCTtacta                                                              >  1:297679/1‑81 (MQ=255)              tgCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGTGTTCATGCTTACTACg                                                            <  1:297634/81‑1 (MQ=255)                cTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGTGTTCATGCTTACTACGat                                                          <  2:297635/81‑1 (MQ=255)                    tGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGTGTTCATGCTTACTACGATAttt                                                      <  1:297636/81‑1 (MQ=255)                     gcgcGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGTGTTCATGCTTACTACGATAtttt                                                     >  2:297677/1‑81 (MQ=255)                      cgcgAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGTGTTCATGCTTACTACGATATTTTg                                                    <  2:297637/81‑1 (MQ=255)                           aGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGTGTTCATGCTTACTACGATATTTTGAGCGa                                               <  2:297638/81‑1 (MQ=255)                              tGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGTGTTCATGCTTACTACGATATTTTGAGCGAAGc                                            >  2:297680/1‑81 (MQ=255)                                    ggAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGTGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGAtg                                      <  2:297640/81‑1 (MQ=255)                                    ggAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGTGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGAt                                       <  2:297639/80‑1 (MQ=255)                                                   ggATCGCGGGCGTGAGCCGTTGGCTGGTGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTg                       >  2:297683/1‑81 (MQ=255)                                                     aTCGCGGGCGTGAGCCGTTGGCTGGTGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGc                     >  1:297678/1‑81 (MQ=255)                                                      tCGCGGGCGTGAGCCGTTGGCTGGTGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCg                    <  1:297641/81‑1 (MQ=255)                                                              gTGAGCCGTTGGCTGGTGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACct            <  1:297642/81‑1 (MQ=255)                                                                gAGCCGTTGGCTGGTGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACctac          >  1:297682/1‑81 (MQ=255)                                                                 aGCCGTTGGCTGGTGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACctac          <  2:297643/80‑1 (MQ=255)                                                                         gctggtGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCa  <  1:297644/80‑1 (MQ=255)                                                                          ctggtGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCAc    >  1:297695/1‑77 (MQ=255)                                                                          ctggtGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCAc    >  1:297688/1‑77 (MQ=255)                                                                              |                                                                           GAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCA  >  NC_000913/1607751‑1607902

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑