Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,647,428	C→T	A377V (GCG→GTG)	gor →	glutathione reductase (NADPH)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,647,428	0	C	T	95.0%	47.9 / NA	20	A377V (GCG→GTG)	gor	glutathione reductase (NADPH)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (10/9);  minor base G (0/1);  total (10/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

CGATAACACGGGTGCAGTGGAGCTGACACCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCA‑GGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTA  >  NC_000913/3647204‑3647657                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                      cGATAACACGGGTGCAGTGGAGCTGACACCGGTGGCAGTTGCAGCGGGTCGACGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGCTTACAGCAACATTCCGACCGTGGTCTTCAGCCCTCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCCGTATGGCGACGATCA‑GGAGAAAGTGTATAAATCCTCTTTCACCGTGCTGTATACCGCCGTCACCACTCAcc                                                                                                                                                                                                             >  1:130627/1‑251 (MQ=255)                       cTGACACCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCCTCTGGATTACAGCAACATTCCGAACGTCGTCTTCACACATCCGCCGATTGGTACTGTTGGTTTAACGGACCCGCAGGCGCGCGAGCAGTATGGCGACGATCA‑GGTGAACGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACAGCCAGCCGTGCCGCATGCAGCt                                                                                                                                                                                       >  1:245120/1‑251 (MQ=255)                         gACACCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCA‑GGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTgg                                                                                                                                                                                     <  1:231327/251‑1 (MQ=255)                             ccGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCA‑GGTGAAAGTGTATAAAGCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGtgtg                                                                                                                                                                                 <  1:146771/251‑1 (MQ=255)                                   gCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCA‑GGTGAAAGTGTATAAAGCCTCTTTCACCGTGAGGTAGACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTgn                                                                                                                                                                           <  1:33235/251‑2 (MQ=255)                                                    ccGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCA‑GGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTg                                                                                                                                                          >  1:97883/1‑251 (MQ=255)                                                          ctctGACCGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTATTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCA‑GGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTc                                                                                                                                                        >  1:300556/1‑247 (MQ=255)                                                                 ctccTGTTTAATAACAAGCCGGATGAGCATCTGGATTAGAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGAGCCGATCA‑GGTGAAAGTGTATACATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACgg                                                                                                                                             >  1:338744/3‑251 (MQ=255)                                                                   ccTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCA‑GGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGc                                                                                                                                            >  1:169152/1‑250 (MQ=255)                                                                            aaaCAAGCCGGATGAGCATCTGGAGTACAGCAACATTCCGACCGTGGTCTTCATCCATCCGCCGATTGGTACTGTTGGTTGAACGGAACCGCAGGCGCGCGAGCAGTATGGCGATGATCA‑GGTGAAAGTGTATAAATCCGCTTTCACCGTGATGTATACCGCCGTCACCACTTACCTCCAGCCGTGCCTCATGAAGCGGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTg                                                                                                                                  <  1:255242/250‑1 (MQ=255)                                                                                                 ggATTACAGCAACATGCCGACCGTGGTCTTCTGCCATCCGCCGATTGGTGCTGTTGGTTTAACGGAACCGCAGTCGCGCGAGCAGTATGGCGACGATCA‑GGTGAAAGTGTATAAGGCCGCTGTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAATCTGGTGTGCGGTGGATCGGAAGATAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGc                                                                                                                <  1:324752/248‑1 (MQ=255)                                                                                                       cAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCA‑GGTGAAAGTATATAAATCCTCTTTCACCGTGATGTAGACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATTAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACTGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTcg                                                                                                       <  1:348310/251‑1 (MQ=255)                                                                                                       cAGCAACATCCCGACCGTGGTCGTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCA‑GGTGAAAGGATATAAAGCCTCTTTCACCGTGATGTATACCCCCGTCCCCACTCACCGCCAGCCGTGCCTCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTcg                                                                                                       <  1:359695/251‑1 (MQ=255)                                                                                                                          gTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCA‑GGTGAAAGTGTATAAATCCTCTTTCACCGGGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCgn                                                                                                                                                                      <  1:44103/169‑2 (MQ=255)                                                                                                                                      ccgccgATTGGTACTGGTGGTTGAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCA‑GGTGAAAGTGTATAAACCCTCTTTCCCCGTGATGTATACCGCCTTCACCACTCACCGCCAGCCGTGTCGCATGAAGCTGGTGTGCGTGGGATCGGAAGATAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGc                                                                                                           <  1:316707/216‑1 (MQ=255)                                                                                                                                        gccgATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGAGACCAGTATGGCGACGATCA‑GGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCaaaaa                                                                        >  1:227725/1‑249 (MQ=255)                                                                                                                                         ccgATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCA‑GGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTgg                                                                                                 >  1:225159/1‑223 (MQ=255)                                                                                                                                                        ggTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCA‑GGTGCCAGTGTATCAATCCTCTTTCACCGTGATGTATCCCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGgaaga                                                                                                                                                                  >  1:191729/1‑143 (MQ=255)                                                                                                                                                                       cAGGCGCTCGAGCAGTATGGCGACGAGCATGGT‑AAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCGTTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCTCGGCAGTGGCTTTGTTATGTACGAAATGTTGCAGGGCTTCGCGTTGGCGTTGAATATGGGGGCAACTAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAAc                                       <  1:259705/251‑1 (MQ=255)                                                                                                                                                                                                              gtATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACAACTCACCGCCAGCCGTGCCGCATGAAGATGGAGTGCGTTGGATCGGAAGAGAAGATTGTCGGAATTAAAGGCATTGGCTTTGGAATGGACGAAATGTTGAAGGGCTTCGCGGGGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGCGTTCGTGACAAAGCGTTAACTGTTa  >  1:248190/1‑250 (MQ=255)                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                      CGATAACACGGGTGCAGTGGAGCTGACACCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCA‑GGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTA  >  NC_000913/3647204‑3647657

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