Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	2,448,157	C→T	S316S (AGC→AGT)	trmC →	fused 5‑methylaminomethyl‑2‑thiouridine forming enzyme methyltransferase and FAD‑dependent demodification enzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,448,157	0	C	T	100.0%	45.9 / NA	27	S316S (AGC→AGT)	trmC	fused 5‑methylaminomethyl‑2‑thiouridine forming enzyme methyltransferase and FAD‑dependent demodification enzyme
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/0);  total (27/0)

GTGCTTCCGGCAATCGCCAGGGGGCGCTGTATCCGTTATTAAGCAAACACGATGAGGCGCTA  >  W3110S.gb/2448114‑2448175                                            |                   gTGCTTCCGGCAATCGCCAGGGGGCGCTGTATCCGCTTTTAAGTAAACACGATGaag       >  1:20268/1‑55 (MQ=255) gTGCTTCCGGCAATCGCCAGGGGGCGCTGTATCCGCTTTTAAGTAAACACGATGAAGCGCTa  >  1:41405/1‑62 (MQ=255) gTGCTTCCGGCAATCGCCAGGGGGCGCTGTATCCGCTTTTAAGTAAACACGATGAAGCGCTa  >  1:80227/1‑62 (MQ=255) gTGCTTCCGGCAATCGCCAGGGGGCGCTGTATCCGCTTTTAAGTAAACACGATGAAGCGCTa  >  1:80021/1‑62 (MQ=255) gTGCTTCCGGCAATCGCCAGGGGGCGCTGTATCCGCTTTTAAGTAAACACGATGAAGCGCTa  >  1:75940/1‑62 (MQ=255) gTGCTTCCGGCAATCGCCAGGGGGCGCTGTATCCGCTTTTAAGTAAACACGATGAAGCGCTa  >  1:70450/1‑62 (MQ=255) gTGCTTCCGGCAATCGCCAGGGGGCGCTGTATCCGCTTTTAAGTAAACACGATGAAGCGCTa  >  1:67185/1‑62 (MQ=255) gTGCTTCCGGCAATCGCCAGGGGGCGCTGTATCCGCTTTTAAGTAAACACGATGAAGCGCTa  >  1:6445/1‑62 (MQ=255) gTGCTTCCGGCAATCGCCAGGGGGCGCTGTATCCGCTTTTAAGTAAACACGATGAAGCGCTa  >  1:58588/1‑62 (MQ=255) gTGCTTCCGGCAATCGCCAGGGGGCGCTGTATCCGCTTTTAAGTAAACACGATGAAGCGCTa  >  1:53438/1‑62 (MQ=255) gTGCTTCCGGCAATCGCCAGGGGGCGCTGTATCCGCTTTTAAGTAAACACGATGAAGCGCTa  >  1:52336/1‑62 (MQ=255) gTGCTTCCGGCAATCGCCAGGGGGCGCTGTATCCGCTTTTAAGTAAACACGATGAAGCGCTa  >  1:51268/1‑62 (MQ=255) gTGCTTCCGGCAATCGCCAGGGGGCGCTGTATCCGCTTTTAAGTAAACACGATGAAGCGCTa  >  1:49460/1‑62 (MQ=255) gTGCTTCCGGCAATCGCCAGGGGGCGCTGTATCCGCTTTTAAGTAAACACGATGAAGCGCTa  >  1:43893/1‑62 (MQ=255) gTGCTTCCGGCAATCGCCAGGGGGCGCTGTATCCGCTTTTAAGTAAACACGATGAAGCGCTa  >  1:12334/1‑62 (MQ=255) gTGCTTCCGGCAATCGCCAGGGGGCGCTGTATCCGCTTTTAAGTAAACACGATGAAGCGCTa  >  1:40990/1‑62 (MQ=255) gTGCTTCCGGCAATCGCCAGGGGGCGCTGTATCCGCTTTTAAGTAAACACGATGAAGCGCTa  >  1:38069/1‑62 (MQ=255) gTGCTTCCGGCAATCGCCAGGGGGCGCTGTATCCGCTTTTAAGTAAACACGATGAAGCGCTa  >  1:34559/1‑62 (MQ=255) gTGCTTCCGGCAATCGCCAGGGGGCGCTGTATCCGCTTTTAAGTAAACACGATGAAGCGCTa  >  1:30772/1‑62 (MQ=255) gTGCTTCCGGCAATCGCCAGGGGGCGCTGTATCCGCTTTTAAGTAAACACGATGAAGCGCTa  >  1:25231/1‑62 (MQ=255) gTGCTTCCGGCAATCGCCAGGGGGCGCTGTATCCGCTTTTAAGTAAACACGATGAAGCGCTa  >  1:23488/1‑62 (MQ=255) gTGCTTCCGGCAATCGCCAGGGGGCGCTGTATCCGCTTTTAAGTAAACACGATGAAGCGCTa  >  1:23056/1‑62 (MQ=255) gTGCTTCCGGCAATCGCCAGGGGGCGCTGTATCCGCTTTTAAGTAAACACGATGAAGCGCTa  >  1:22233/1‑62 (MQ=255) gTGCTTCCGGCAATCGCCAGGGGGCGCTGTATCCGCTTTTAAGTAAACACGATGAAGCGCTa  >  1:20586/1‑62 (MQ=255) gTGCTTCCGGCAATCGCCAGGGGGCGCTGTATCCGCTTTTAAGTAAACACGATGAAGCGCTa  >  1:13746/1‑62 (MQ=255) gTGCTTCCGGCAATCGCCAGGGGGCGCTGTATCCGCTTTTAAGTAAACACGATGAAGCGCTa  >  1:12524/1‑62 (MQ=255) gTGCTTCCGGCAATCGACAGGGGGCGCTGTATCCGCTTTTAAGTAAACACGATGAAGCGCTa  >  1:68083/1‑62 (MQ=255)                                            |                   GTGCTTCCGGCAATCGCCAGGGGGCGCTGTATCCGTTATTAAGCAAACACGATGAGGCGCTA  >  W3110S.gb/2448114‑2448175

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: ‑