Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,647,011	C→T	T238I (ACC→ATC)	gor →	glutathione reductase (NADPH)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,647,011	0	C	T	100.0%	17.9 / NA	8	T238I (ACC→ATC)	gor	glutathione reductase (NADPH)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (4/4);  total (4/4)

CGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAGCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTATTTACGCGGTGGGCGATAACACGGGTGCAGTGGAGCTGACACCGGTGGCAGTTGCAGCGGGT  >  NC_000913/3646810‑3647252                                                                                                                                                                                                          |                                                                                                                                                                                                                                                  cGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAGCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATATCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTg                                                                                                                                                                                                  <  1:2865/251‑1 (MQ=255)            gcgGGTTACATCGCCGTTGAGCTGGCGGGTGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTACACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAGCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATATCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGAt                                                                                                                                                                                       >  1:51156/1‑251 (MQ=255)               ngTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCATCGGCGCGAAAACGCATATGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGAGATGAACGCCGAAGGCCCGCAGCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATATCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTg                                                                                                                                                                                    >  1:18409/2‑251 (MQ=255)                                                                                                                                                           cccGCAGCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAANAATNTCGANGGNAGCCTGACGCTGGAGNTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGNNTNNNCNCGAGCCTg                                                                                                                                                                <  2:118264/131‑1 (MQ=255)                                                                                                                                                           cccGCAGCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATATCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGACACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTg                                                                                                                                                                >  1:118264/1‑131 (MQ=255)                                                                                                                                                                                          nnnnnnnnnnaannNNNCGANTGTAGCCGGACGCTGGAGCTGGACGATGCTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAAACTGGAAGCCGCTGGCGTTAAACCTAACGAAAAAGGCt                                                                                                          <  2:103912/144‑1 (MQ=255)                                                                                                                                                                                          ncGGTAGTGAAAAATATCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCt                                                                                                          >  1:103912/2‑154 (MQ=255)                                                                                                                                                                                             gTAGTGAAAAATATCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTATTTACGCGGTGGGCGATAACACGGGTGCAGTGGAGCTGACACCGGTGGCAGTTGCAGc      <  1:80674/251‑1 (MQ=255)                                                                                                                                                                                                  gAAAAATATCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGAGTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTATTTACGCGGTGGGCGATAACACGGGTGCAGTGGAGCTGACACCGGTGTCAGTTGCAGCgggn  <  1:48667/250‑2 (MQ=255)                                                                                                                                                                                                          |                                                                                                                                                                                                                                                  CGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAGCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTATTTACGCGGTGGGCGATAACACGGGTGCAGTGGAGCTGACACCGGTGGCAGTTGCAGCGGGT  >  NC_000913/3646810‑3647252

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑