Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	4,431,950	A→G	V52A (GTT→GCT)	yjfZ ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,431,950	0	A	G	100.0%	325.2 / NA	86	V52A (GTT→GCT)	yjfZ	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (56/30);  total (56/30)

CTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAACAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGA  >  W3110S.gb/4431893‑4431993                                                          |                                            cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTAc                                       >  1:1462388/1‑64 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTAc                                       >  1:524817/1‑64 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:2168041/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:2435353/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:2402703/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:2322404/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:2285723/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:2238500/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:2221445/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:2209181/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:2440400/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:2165409/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:2133677/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:1962077/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:1951765/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:1908938/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:1896543/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:1874554/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:1963607/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:2440404/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:2442815/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:2446446/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:2488299/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:371699/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:372043/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:372113/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:396113/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:422265/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:516649/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:521170/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:727269/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:755351/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:828038/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:902811/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:1362064/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:1000200/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:114445/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:116165/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:1168812/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:117927/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:1201601/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:1211245/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:1310478/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:1311435/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:147572/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:149857/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:1512661/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:1532868/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:1540303/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:1794900/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:1632734/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:1697863/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:1765886/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:1783490/1‑65 (MQ=255) cTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                      >  1:1658603/1‑65 (MQ=255) cTCGTCGCT‑ACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTAc                                       >  1:1673512/1‑63 (MQ=255)                                     agaTCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGa  <  1:814239/65‑1 (MQ=255)                                     agaTCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGa  <  1:765760/65‑1 (MQ=255)                                     agaTCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGa  <  1:945574/65‑1 (MQ=255)                                     agaTCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGa  <  1:1225770/65‑1 (MQ=255)                                     agaTCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGa  <  1:1078021/65‑1 (MQ=255)                                     agaTCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGa  <  1:716216/65‑1 (MQ=255)                                     agaTCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGa  <  1:6501/65‑1 (MQ=255)                                     agaTCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGa  <  1:839821/65‑1 (MQ=255)                                     agaTCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGa  <  1:1711040/65‑1 (MQ=255)                                     agaTCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGa  <  1:1198264/65‑1 (MQ=255)                                     agaTCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGa  <  1:960006/65‑1 (MQ=255)                                     agaTCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGa  <  1:1222746/65‑1 (MQ=255)                                     agaTCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGa  <  1:2019234/65‑1 (MQ=255)                                     agaTCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGa  <  1:365412/65‑1 (MQ=255)                                     agaTCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGa  <  1:290026/65‑1 (MQ=255)                                     agaTCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGa  <  1:2092/65‑1 (MQ=255)                                     agaTCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGa  <  1:1338666/65‑1 (MQ=255)                                     agaTCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGa  <  1:1353969/65‑1 (MQ=255)                                     agaTCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGa  <  1:2441670/65‑1 (MQ=255)                                     agaTCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGa  <  1:179624/65‑1 (MQ=255)                                     agaTCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGa  <  1:1402198/65‑1 (MQ=255)                                     agaTCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGa  <  1:1688150/65‑1 (MQ=255)                                     agaTCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGa  <  1:2403327/65‑1 (MQ=255)                                     agaTCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGa  <  1:2398340/65‑1 (MQ=255)                                     agaTCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGa  <  1:2368226/65‑1 (MQ=255)                                     agaTCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGa  <  1:1649639/65‑1 (MQ=255)                                      gaTCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGa  <  1:1044295/64‑1 (MQ=255)                                      gaTCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGa  <  1:1846781/64‑1 (MQ=255)                                       aTCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGa  <  1:755835/63‑1 (MQ=255)                                                       tcAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGa  <  1:975591/47‑1 (MQ=255)                                                          |                                            CTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAACAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGA  >  W3110S.gb/4431893‑4431993

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