Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	1,959,670	C→T	S59S (TCG→TCA)	torY ←	TMAO reductase III (TorYZ), cytochrome c‑type subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,959,670	0	C	T	100.0%	166.6 / NA	50	S59S (TCG→TCA)	torY	TMAO reductase III (TorYZ), cytochrome c‑type subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (50/0);  total (50/0)

AATATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTCGAAAAGTGGACAGTTCCCT  >  W3110S.gb/1959628‑1959689                                           |                    aaTATGGCAATCGGCACATTGCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTccc   >  1:576003/1‑61 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTTTCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:1043950/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTccc   >  1:1787009/1‑61 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTccc   >  1:1899539/1‑61 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:487623/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:2521134/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:2564183/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:2572027/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:2691496/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:2760809/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:2773476/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:312777/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:384256/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:1022924/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:390181/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:471283/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:82152/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:503847/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:997193/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:537751/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:537911/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:540599/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:542023/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:954198/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:602938/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:683974/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:792498/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:1793155/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:1044247/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:1066522/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:1240951/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:1259222/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:1504318/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:1516643/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:1550852/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:1610473/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:1664117/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:1671254/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:2470532/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:1805973/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:1858517/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:1891779/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:194878/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:2073932/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:2166658/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:232462/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:2355922/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:2388804/1‑62 (MQ=255) aaTATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCt  >  1:2390517/1‑62 (MQ=255) aaTATGGCAATCGCCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAg        >  1:1784782/1‑56 (MQ=255)                                           |                    AATATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTCGAAAAGTGGACAGTTCCCT  >  W3110S.gb/1959628‑1959689

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