Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	2,355,666	(T)5→4	coding (16/1077 nt)	glpQ ←	periplasmic glycerophosphodiester phosphodiesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,355,662	0	T	.	100.0%	132.3 / NA	31	coding (20/1077 nt)	glpQ	periplasmic glycerophosphodiester phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (24/7);  total (24/7)

GCCGCCATTGCACTGCTTCCCATGACTATAGTGCTCATCATGATCGCCATGCTAAGGTTTTTCAGCGTCAATTTCATGCCATTAGCCTCCGTTGCGTTCTTGCAGT  >  W3110S.gb/2355605‑2355710                                                          |                                                 gccgccATTGCACTGCTTCCCATGACTATAGTGCTCATCATGATCGCCATGCTAAGGttttc                                              >  1:2897317/1‑61 (MQ=255) gccgccATTGCACTGCTTCCCATGACTATAGTGCTCATCATGATCGCCATGCTAAGGtttt                                               >  1:1174161/1‑61 (MQ=255)              tGCTTCCCATGACTATAGTGCTCATCATGATCGCCATGCTAAGG‑TTTTCAGCGTCAATTTCa                                >  1:1299696/1‑62 (MQ=255)              tGCTTCCCATGACTATAGTGCTCATCATGATCGCCATGCTAAGG‑TTTTCAGCGTCAATTTCa                                >  1:924498/1‑62 (MQ=255)                                tGCTCATCATGATCGCCATGCTATGG‑TTTTCAGCGTCAATTTCATGCCATTAGCCTCCGTTg              >  1:924002/1‑62 (MQ=255)                                tGCTCATCATGATCGCCATGCTAAGG‑TTTTCAGCGTCAATTTCATGCCATTAGCCTCCGTTg              >  1:2499190/1‑62 (MQ=255)                                tGCTCATCATGATCGCCATGCTAAGG‑TTTTCAGCGTCAATTTCATGCCATTAGCCTCCGTTg              >  1:3105026/1‑62 (MQ=255)                                tGCTCATCATGATCGCCATGCTAAGG‑TTTTCAGCGTCAATTTCATGCCATTAGCCTCCGTTg              >  1:58504/1‑62 (MQ=255)                                tGCTCATCATGATCGCCATGCTAAGG‑TTTTCAGCGTCAATTTCATGCCATTAGCCTCCGTTg              >  1:2595881/1‑62 (MQ=255)                                tGCTCATCATGATCGCCATGCTAAGG‑TTTTCAGCGTCAATTTCATGCCATTAGCCTCCGTTg              >  1:2535635/1‑62 (MQ=255)                                tGCTCATCATGATCGCCATGCTAAGG‑TTTTCAGCGTCAATTTCATGCCATTAGCCTCCGTTg              >  1:2200766/1‑62 (MQ=255)                                tGCTCATCATGATCGCCATGCTAAGG‑TTTTCAGCGTCAATTTCATGCCATTAGCCTCCGTTg              >  1:2135143/1‑62 (MQ=255)                                tGCTCATCATGATCGCCATGCTAAGG‑TTTTCAGCGTCAATTTCATGCCATTAGCCTCCGTTg              >  1:1673208/1‑62 (MQ=255)                                        atGATCGCCATGCTAAGG‑TTTTCAGCGTCAATTTCATGCCATTAGCCTCCGtt               <  1:655071/53‑1 (MQ=255)                                        atGATCGCCATGCTAAGG‑TTTTCAGCGTCAATTTCATGCCATTAGCCTCCGtt               <  1:2741970/53‑1 (MQ=255)                                        atGATCGCCATGCTAAGG‑TTTTCAGCGTCAATTTCATGCCATTAGCCTCCGtt               <  1:793272/53‑1 (MQ=255)                                        atGATCGCCATGCTAAGG‑TTTTCAGCGTCAATTTCATGCCATTAGCCTCCGTTGCGTTc         <  1:2916630/59‑1 (MQ=255)                                        atGATCGCCATGCTAAGG‑TTTTCAGCGTCAATTTCATGCCATTAGCCTCCGTTGCGTTc         <  1:2971216/59‑1 (MQ=255)                                        atGATCGCCATGCTAAGG‑TTTTCAGCGTCAATTTCATGCCATTAGCCTCCGTTGCGTTc         <  1:481451/59‑1 (MQ=255)                                        atGATCGCCATGCTAAGG‑TTTTCAGCGTCAATTTCATGCCATAAGCCTCCGtt               <  1:2523685/53‑1 (MQ=38)                                               ccATGCTAAGG‑TTTTCAGCGTCAATTTCATGCCATTAGCCTCCGTTGCGTTCTTGCAGt  >  1:589016/1‑59 (MQ=255)                                               ccATGCTAAGG‑TTTTCAGCGTCAATTTCATGCCATTAGCCTCCGTTGCGTTCTTGCAGt  >  1:1256716/1‑59 (MQ=255)                                               ccATGCTAAGG‑TTTTCAGCGTCAATTTCATGCCATTAGCCTCCGTTGCGTTCTTGCAGt  >  1:654357/1‑59 (MQ=255)                                               ccATGCTAAGG‑TTTTCAGCGTCAATTTCATGCCATTAGCCTCCGTTGCGTTCTTGCAGt  >  1:579467/1‑59 (MQ=255)                                               ccATGCTAAGG‑TTTTCAGCGTCAATTTCATGCCATTAGCCTCCGTTGCGTTCTTGCAGt  >  1:3020541/1‑59 (MQ=255)                                               ccATGCTAAGG‑TTTTCAGCGTCAATTTCATGCCATTAGCCTCCGTTGCGTTCTTGCAGt  >  1:2645349/1‑59 (MQ=255)                                               ccATGCTAAGG‑TTTTCAGCGTCAATTTCATGCCATTAGCCTCCGTTGCGTTCTTGCAGt  >  1:2273538/1‑59 (MQ=255)                                               ccATGCTAAGG‑TTTTCAGCGTCAATTTCATGCCATTAGCCTCCGTTGCGTTCTTGCAGt  >  1:2190761/1‑59 (MQ=255)                                               ccATGCTAAGG‑TTTTCAGCGTCAATTTCATGCCATTAGCCTCCGTTGCGTTCTTGCAGt  >  1:1976807/1‑59 (MQ=255)                                               ccATGCTAAGG‑TTTTCAGCGTCAATTTCATGCCATTAGCCTCCGTTGCGTTCTTGCAGt  >  1:1942064/1‑59 (MQ=255)                                               ccATGCTAAGG‑TTTTCAGCGTCAATTTCATGCCATTAGCCTCCGTTGCGTTCTTGCAGt  >  1:1932006/1‑59 (MQ=255)                                               ccATGCTAAGG‑TTTTCAGCGTCAATTTCATGCCATTAGCCTCCGTTGCGTTCTTGCAGt  >  1:1728405/1‑59 (MQ=255)                                               ccATGCTAAGG‑TTTTCAGCGTCAATTTCAAGCCATTAGACTCCGt                >  1:931868/1‑45 (MQ=255)                                                          |                                                 GCCGCCATTGCACTGCTTCCCATGACTATAGTGCTCATCATGATCGCCATGCTAAGGTTTTTCAGCGTCAATTTCATGCCATTAGCCTCCGTTGCGTTCTTGCAGT  >  W3110S.gb/2355605‑2355710

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