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%	207.0 / NA	48	coding (3047/3363 nt)	kefA	fused mechanosensitive channel proteins
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (39/9);  total (39/9)

CTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGA  >  W3110S.gb/488749‑488863                                                          |                                                          ctgacACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGaaaagtg                                                     <  1:2176232/64‑4 (MQ=255)             cgcgTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTa                                        >  1:1101875/1‑64 (MQ=255)             cgcgTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTa                                        >  1:2805255/1‑64 (MQ=255)             cgcgTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTa                                        >  1:1811799/1‑64 (MQ=255)             cgcgTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTa                                        >  1:2910077/1‑64 (MQ=255)             cgcgTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTa                                        >  1:2925115/1‑64 (MQ=255)             cgcgTCTGGTGATCCGTATCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTa                                        >  1:1422632/1‑64 (MQ=255)                 tCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGaa                                   >  1:2085638/1‑65 (MQ=255)                 tCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGaa                                   >  1:2995853/1‑65 (MQ=255)                 tCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGaa                                   >  1:477236/1‑65 (MQ=255)                         tCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGa                           >  1:1663064/1‑65 (MQ=255)                         tCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGa                           >  1:600724/1‑65 (MQ=255)                         tCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGa                           >  1:846996/1‑65 (MQ=255)                         tCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGa                           >  1:2016666/1‑65 (MQ=255)                                     tgcccTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCaa               <  1:1007591/62‑1 (MQ=255)                                     tGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCaa               <  1:3103158/65‑1 (MQ=255)                                     tGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCaa               <  1:787590/65‑1 (MQ=255)                                     tGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCaa               <  1:2776014/65‑1 (MQ=255)                                     tGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCaa               <  1:233455/65‑1 (MQ=255)                                     tGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCaa               <  1:937019/65‑1 (MQ=255)                                     tGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCaa               <  1:1974950/65‑1 (MQ=255)                                     tGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCaa               <  1:1215348/65‑1 (MQ=255)                                     tGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCaa               <  1:1338817/65‑1 (MQ=255)                                                tacGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa     >  1:2257823/3‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGc      >  1:1853486/1‑63 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGc      >  1:318145/1‑63 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa     >  1:1623610/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa     >  1:1126228/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa     >  1:1152926/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa     >  1:523121/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa     >  1:47833/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa     >  1:1338015/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa     >  1:405034/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa     >  1:3105838/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa     >  1:211066/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa     >  1:2936396/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa     >  1:2799956/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa     >  1:2682854/1‑64 (MQ=255)                                                tCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCa     >  1:2177868/1‑64 (MQ=255)                                                  cGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGa  >  1:2939642/1‑65 (MQ=255)                                                  cGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGa  >  1:1665316/1‑65 (MQ=255)                                                  cGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGa  >  1:1709591/1‑65 (MQ=255)                                                  cGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGa  >  1:446637/1‑65 (MQ=255)                                                  cGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGa  >  1:2809418/1‑65 (MQ=255)                                                  cGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGa  >  1:185161/1‑65 (MQ=255)                                                  cGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGa  >  1:519701/1‑65 (MQ=255)                                                  cGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGa  >  1:1314563/1‑65 (MQ=255)                                                  cGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGa  >  1:831473/1‑65 (MQ=255)                                                  cGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGa  >  1:2026888/1‑65 (MQ=255)                                                          |                                                          CTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGA  >  W3110S.gb/488749‑488863

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