Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	488,810	(A)5→4	coding (3051/3363 nt)	kefA →	fused mechanosensitive channel proteins

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	488,806	0	A	.	100.0%	193.2 / NA	45	coding (3047/3363 nt)	kefA	fused mechanosensitive channel proteins
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (45/0);  total (45/0)

CTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCA  >  W3110S.gb/488749‑488860                                                          |                                                       ctgacACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGaaaagtg                                                  <  1:1881274/64‑4 (MQ=255) ctgacACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGaaaagtg                                                  <  1:1224359/64‑4 (MQ=255) ctgacACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGaaaagtg                                                  <  1:634222/64‑4 (MQ=255) ctgacACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGaaaagtg                                                  <  1:1284038/64‑4 (MQ=255) ctgacACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGaaaagtg                                                  <  1:2088384/64‑4 (MQ=255) ctgacACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGaaaagtg                                                  <  1:2490148/64‑4 (MQ=255) ctgacACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGaaaagtg                                                  <  1:1528417/64‑4 (MQ=255) ctgacACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGaaaagtg                                                  <  1:78075/64‑4 (MQ=255) ctgacACTACTACGCGTCTGGTGATCCGTCTCGGCGGGGCCTATGGCTCCGATCTGGaaaagtg                                                  <  1:2076697/64‑4 (MQ=255)  tgacACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGaaaagtg                                                  <  1:996231/63‑4 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAg                               >  1:715228/1‑35 (MQ=38)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCAccc              >  1:1307130/1‑52 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTg       >  1:1898928/1‑59 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:2495749/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:283422/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:2483065/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:2477227/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:2130472/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:288061/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:229197/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:442191/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:506090/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:537050/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:640262/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:657773/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:724969/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:77463/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:806002/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:865070/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:966129/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:989439/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:1765471/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:1010196/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:1245239/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:1254542/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:1299930/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:1305548/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:1458239/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:1538795/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:1571401/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:1571631/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:1632120/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:1759001/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:2254154/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:1864021/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:1945092/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:196514/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:2129009/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:2129275/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:1006980/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:2171938/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:2216805/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:2245283/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAACGCGGCGACTGAGCACCCAAGGGTGATGc   >  1:2309287/1‑63 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTCAAGGCGGCGACTGAGCACCCAAGGGTGATGCa  >  1:1946135/1‑64 (MQ=255)                                                          |                                                       CTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCA  >  W3110S.gb/488749‑488860

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