Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	4,544,265	(T)7→6	intergenic (‑82/‑1374)	yjhA ← / → fimB	N‑acetylnuraminic acid outer membrane channel protein/tyrosine recombinase/inversion of on/off regulator of fimA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,544,259	0	T	.	100.0%	352.2 / NA	82	intergenic (‑76/‑1380)	yjhA/fimB	N‑acetylnuraminic acid outer membrane channel protein/tyrosine recombinase/inversion of on/off regulator of fimA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (0/82);  total (0/82)

ATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGC  >  W3110S.gb/4544236‑4544298                        |                                        aTTATCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:836796/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGTCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:833135/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:432287/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:619222/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:612744/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:605398/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:598059/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:582556/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:538768/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:482309/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:465517/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:636435/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:425326/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:421881/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:408133/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:374739/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:346379/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:334879/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:308582/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:302266/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:283740/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:83912/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:988452/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:983450/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:958525/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:95768/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:952138/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:936809/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:905753/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:886350/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:862146/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:635029/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:793116/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:782680/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:760870/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:758226/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:736346/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:653308/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:649690/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:1053483/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:1151749/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:279011/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:1373427/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:1328908/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:1292154/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:1257386/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:119588/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:1188644/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:1186958/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:1179038/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:1178306/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:1388408/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:1112944/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:1096806/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:1096745/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:1094841/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:1084919/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:1079977/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:1065949/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:1054017/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:1033527/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:1397908/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:1509060/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:1513029/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:226631/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:1578052/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:1588740/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:1597287/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:163031/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:1641703/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:1645204/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:201769/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:193513/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:19274/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGCTAGc  <  1:1439297/62‑1 (MQ=255) aTTAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGGCTTTGGCAACGGGTAGc  <  1:173243/62‑1 (MQ=255) aTTAGCATTTATGTTGAGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:260606/62‑1 (MQ=255) aTTAGCATCTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:887989/62‑1 (MQ=255)  ttAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:323315/61‑1 (MQ=255)   tAGCATTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:166079/60‑1 (MQ=255)       aTTTATGTTGTGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:287014/56‑1 (MQ=255)            tgttgtGTAATA‑TTTTTTGCCAGGCTTATAGTGTCTTTGGCAACGGGTAGc  <  1:1548664/51‑1 (MQ=255)                        |                                        ATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGC  >  W3110S.gb/4544236‑4544298

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