Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	2,654,752	(T)5→6	coding (239/2448 nt)	glgP →	glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,654,747	1	.	T	100.0%	67.0 / NA	22	E78D (GAG→GAT)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (0/22);  total (0/22)

TCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑TTTTTGATTGGCCGTACGCTCTCCAA  >  minE/2654691‑2654773                                                          |                           tCAAACCGTGTCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:1118589/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:357277/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:2407065/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:2448951/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:2498931/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:2815755/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:2873845/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:3172878/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:3207114/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:3230875/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:3234198/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:3282594/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:3349315/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:3411126/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:3477297/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:3526559/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:3528410/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:2537043/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:3611837/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:3637693/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:402451/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:428225/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:563305/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:657061/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:685113/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:69052/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:746622/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:879921/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:970620/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:978077/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:996071/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:2399534/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:1082967/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:1095563/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:1171956/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:1181016/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:1333233/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:1356758/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:1437952/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:1672236/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:1757659/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:1951109/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:1987406/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:2024262/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:2404129/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:2030057/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:2034056/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:2102313/62‑1 (MQ=255) tCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCAGTCGATGGAG‑ttttt                       <  1:2470979/62‑1 (MQ=255) tCAAACCGTGCCCAGGTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:1563288/62‑1 (MQ=255)        gTGCGCAGTTGTCGCAAGAAACTCGTCAGGTTTACTTCCTGTCGATGGAG‑ttttt                       <  1:1443524/55‑1 (MQ=255)                ttGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:2654625/47‑1 (MQ=255)                 tGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑ttttt                       <  1:250167/46‑1 (MQ=255)                       aaGAAACTCGTCAGGTTTACTACCTGTCGATGGAGTTTTTTGATTGGCCGTACGCTCTCaaa  <  1:1603/62‑1 (MQ=255)                       aaGAAACTCGTCAGGTTTACTACCTGTCGATGGAGTTTTTTGATTGGCCGTACGCTCTCCaa  <  1:535871/62‑1 (MQ=255)                       aaGAAACTCGTCAGGTTTACTACCTGTCGATGGAGTTTTTTGATTGGCCGTACGCTCTCCaa  <  1:2644003/62‑1 (MQ=255)                       aaGAAACTCGTCAGGTTTACTACCTGTCGATGGAGTTTTTTGATTGGCCGTACGCTCTCCaa  <  1:1002943/62‑1 (MQ=255)                       aaGAAACTCGTCAGGTTTACTACCTGTCGATGGAGTTTTTTGATTGGCCGTACGCTCTCCaa  <  1:1288275/62‑1 (MQ=255)                       aaGAAACTCGTCAGGTTTACTACCTGTCGATGGAGTTTTTTGATTGGCCGTACGCTCTCCaa  <  1:1248427/62‑1 (MQ=255)                       aaGAAACTCGTCAGGTTTACTACCTGTCGATGGAGTTTTTTGATTGGCCGTACGCTCTCCaa  <  1:1196935/62‑1 (MQ=255)                       aaGAAACTCGTCAGGTTTACTACCTGTCGATGGAGTTTTTTGATTGGCCGTACGCTCTCCaa  <  1:794415/62‑1 (MQ=255)                       aaGAAACTCGTCAGGTTTACTACCTGTCGATGGAGTTTTTTGATTGGCCGTACGCTCTCCaa  <  1:892566/62‑1 (MQ=255)                       aaGAAACTCGTCAGGTTTACTACCTGTCGATGGAGTTTTTTGATTGGCCGTACGCTCTCCaa  <  1:2454199/62‑1 (MQ=255)                       aaGAAACTCGTCAGGTTTACTACCTGTCGATGGAGTTTTTTGATTGGCCGTACGCTCTCCaa  <  1:466462/62‑1 (MQ=255)                       aaGAAACTCGTCAGGTTTACTACCTGTCGATGGAGTTTTTTGATTGGCCGTACGCTCTCCaa  <  1:456430/62‑1 (MQ=255)                       aaGAAACTCGTCAGGTTTACTACCTGTCGATGGAGTTTTTTGATTGGCCGTACGCTCTCCaa  <  1:3101313/62‑1 (MQ=255)                       aaGAAACTCGTCAGGTTTACTACCTGTCGATGGAGTTTTTTGATTGGCCGTACGCTCTCCaa  <  1:1475673/62‑1 (MQ=255)                       aaGAAACTCGTCAGGTTTACTACCTGTCGATGGAGTTTTTTGATTGGCCGTACGCTCTCCaa  <  1:1536362/62‑1 (MQ=255)                       aaGAAACTCGTCAGGTTTACTACCTGTCGATGGAGTTTTTTGATTGGCCGTACGCTCTCCaa  <  1:1603999/62‑1 (MQ=255)                       aaGAAACTCGTCAGGTTTACTACCTGTCGATGGAGTTTTTTGATTGGCCGTACGCTCTCCaa  <  1:3475306/62‑1 (MQ=255)                       aaGAAACTCGTCAGGTTTACTACCTGTCGATGGAGTTTTTTGATTGGCCGTACGCTCTCCaa  <  1:1741226/62‑1 (MQ=255)                       aaGAAACTCGTCAGGTTTACTACCTGTCGATGGAGTTTTTTGATTGGCCGTACGCTCTCCaa  <  1:2956048/62‑1 (MQ=255)                       aaGAAACTCGTCAGGTTTACTACCTGTCGATGGAGTTTTTTGATTGGCCGTACGCTCTCCaa  <  1:1987853/62‑1 (MQ=255)                       aaGAAACTCGTCAGGGTTACTACCTGTCGATGGAGTTTTTTGATTGGCCGTACGCTCTCCaa  <  1:1698256/62‑1 (MQ=255)                        aGAAACTCGTCAGGTTTACTACCTGTCGATGGAGTTTTTTGATTGGCCGTACGCTCTCCaa  <  1:4469/61‑1 (MQ=255)                                                          |                           TCAAACCGTGCCCAGTTGTCGCAAGAAACTCGTCAGGTTTACTACCTGTCGATGGAG‑TTTTTGATTGGCCGTACGCTCTCCAA  >  minE/2654691‑2654773

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