Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	1,538,270	A→G	F21S (TTT→TCT)	narV ←	nitrate reductase 2 (NRZ), gamma subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,538,270	0	A	G	100.0%	147.8 / NA	40	F21S (TTT→TCT)	narV	nitrate reductase 2 (NRZ), gamma subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (40/0);  total (40/0)

TAGCGCAGCCAGCTGCCGAGGAAAAACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGAC  >  W3110S.gb/1538246‑1538313                         |                                            tAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTaaaaaaa     >  1:1028260/1‑65 (MQ=255) tAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTaaaaaaa     >  1:1015290/1‑65 (MQ=255) tAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTaaaaaaa     >  1:1287544/1‑65 (MQ=255) tAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTaaaaaaa     >  1:2443530/1‑65 (MQ=255) tAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTaaaaaaa     >  1:137518/1‑65 (MQ=255) tAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTaaaaaaa     >  1:745415/1‑65 (MQ=255) tAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTaaaaaaa     >  1:409656/1‑65 (MQ=255) tAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTaaaaaaa     >  1:537159/1‑65 (MQ=255)    cgcAGCCAGCTGCCGAGTAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGAc  >  1:793941/1‑65 (MQ=255)    cgcAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGAc  >  1:655046/1‑65 (MQ=255)    cgcAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGAc  >  1:2412985/1‑65 (MQ=255)    cgcAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGAc  >  1:2439537/1‑65 (MQ=255)    cgcAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGAc  >  1:431060/1‑65 (MQ=255)    cgcAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGAc  >  1:432015/1‑65 (MQ=255)    cgcAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGAc  >  1:480980/1‑65 (MQ=255)    cgcAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGAc  >  1:639470/1‑65 (MQ=255)    cgcAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGAc  >  1:2026501/1‑65 (MQ=255)    cgcAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGAc  >  1:711031/1‑65 (MQ=255)    cgcAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGAc  >  1:911096/1‑65 (MQ=255)    cgcAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGAc  >  1:946772/1‑65 (MQ=255)    cgcAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGAc  >  1:969288/1‑65 (MQ=255)    cgcAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGAc  >  1:979065/1‑65 (MQ=255)    cgcAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGAc  >  1:2190447/1‑65 (MQ=255)    cgcAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGAc  >  1:1263540/1‑65 (MQ=255)    cgcAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGAc  >  1:1267596/1‑65 (MQ=255)    cgcAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGAc  >  1:1328687/1‑65 (MQ=255)    cgcAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGAc  >  1:1419252/1‑65 (MQ=255)    cgcAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGAc  >  1:1433946/1‑65 (MQ=255)    cgcAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGAc  >  1:1451839/1‑65 (MQ=255)    cgcAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGAc  >  1:1540528/1‑65 (MQ=255)    cgcAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGAc  >  1:1637523/1‑65 (MQ=255)    cgcAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGAc  >  1:1650372/1‑65 (MQ=255)    cgcAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGAc  >  1:1678221/1‑65 (MQ=255)    cgcAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGAc  >  1:1749225/1‑65 (MQ=255)    cgcAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGAc  >  1:1831237/1‑65 (MQ=255)    cgcAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGAc  >  1:1838375/1‑65 (MQ=255)    cgcAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGAc  >  1:1840207/1‑65 (MQ=255)    cgcAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGAc  >  1:1974071/1‑65 (MQ=255)    cgcAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGAc  >  1:2361986/1‑65 (MQ=255)    cgcAGCCAGCTGACGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGAc  >  1:1018521/1‑65 (MQ=255)                         |                                            TAGCGCAGCCAGCTGCCGAGGAAAAACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGAC  >  W3110S.gb/1538246‑1538313

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