Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	4,001,567	G→A	R403R (CGC→CGT)	pitA ←	phosphate transporter, low‑affinity

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,001,567	0	G	A	100.0%	99.6 / NA	58	R403R (CGC→CGT)	pitA	phosphate transporter, low‑affinity
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (58/0);  total (58/0)

GATTTTCTCACCGATAGTCGTTGCCACACGGCGCCAGCCAATCATCGTACC  >  W3110S.gb/4001537‑4001587                               |                     gATTTTCTCACCGATAGTCGTTGCCAGACGACGCCAGCCAATCATCGTAcc  >  1:79478/1‑51 (MQ=25) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAatca         >  1:40267/1‑44 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:71696/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:11150/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:50511/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:5134/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:52954/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:62296/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:62971/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:65961/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:66192/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:67025/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:67372/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:67918/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:68224/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:6957/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:50307/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:73090/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:75433/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:77373/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:79291/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:80101/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:81725/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:85296/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:85692/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:87453/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:9105/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:944/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:9844/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:33967/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:11191/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:12809/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:13633/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:1375/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:16244/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:17267/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:18562/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:20661/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:21008/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:24611/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:27145/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:29298/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:31166/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:3372/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:50254/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:37255/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:38319/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:38362/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:38543/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:3915/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:39788/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:43294/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:44645/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:45562/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:45976/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:46699/1‑51 (MQ=255) gATTTTCTCACCGATAGTCGTTGCCACACGACGCCAGCCAATCATCGTAc   >  1:66288/1‑50 (MQ=255) gATTTTCTCACCGATAGTAGTTGCCACACGACGCCAGCCAATCATCGTAcc  >  1:68565/1‑51 (MQ=25)                               |                     GATTTTCTCACCGATAGTCGTTGCCACACGGCGCCAGCCAATCATCGTACC  >  W3110S.gb/4001537‑4001587

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