Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	4,007,539	G→A	pseudogene (367/393 nt)	yhiK →	ECK3474:JW5674:b3489; hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,007,539	0	G	A	100.0%	191.1 / NA	57	pseudogene (367/393 nt)	yhiK	ECK3474:JW5674:b3489; hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (50/7);  total (50/7)

GGCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAGGGCTTAGTTGTTTTCAACTACATTAAAACTCCATGAATACGGTGGTATGTCTTT  >  W3110S.gb/4007483‑4007593                                                         |                                                       ggCTAAATTCCCTGAAGGAAGCGCGTGATAAAAATAGTAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:1899049/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATATAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:364390/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAGAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:1085104/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:255111/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:1913358/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:1030489/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:2075689/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:2090441/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:2124349/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:2133707/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:214322/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:2173591/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:2330512/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:2334798/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:2441261/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:2506864/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:1915978/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:275177/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:40400/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:434874/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:677492/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:681156/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:692927/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:699086/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:788594/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:793386/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:894633/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:1467031/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:1031478/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:1121173/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:1153891/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:1167791/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:1205936/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:1228084/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:1326718/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:1346784/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:1363077/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:1399893/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:1405829/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:183610/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:1475088/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:1478257/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:1483079/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:149635/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:1539332/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:159179/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:1637161/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAgt                                                >  1:1824584/1‑65 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAg                                                 >  1:516595/1‑64 (MQ=255) ggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAg                                                 >  1:168919/1‑64 (MQ=255)                                               ttATTCCTCAAGGCTTAGTTGTTTTCAACTACATTAAAACTCCATGAATACGGTGGTATGTCttt  <  1:2398770/65‑1 (MQ=255)                                               ttATTCCTCAAGGCTTAGTTGTTTTCAACTACATTAAAACTCCATGAATACGGTGGTATGTCttt  <  1:2324724/65‑1 (MQ=255)                                               ttATTCCTCAAGGCTTAGTTGTTTTCAACTACATTAAAACTCCATGAATACGGTGGTATGTCttt  <  1:2297196/65‑1 (MQ=255)                                               ttATTCCTCAAGGCTTAGTTGTTTTCAACTACATTAAAACTCCATGAATACGGTGGTATGTCttt  <  1:668967/65‑1 (MQ=255)                                               ttATTCCTCAAGGCTTAGTTGTTTTCAACTACATTAAAACTCCATGAATACGGTGGTATGTCttt  <  1:1199946/65‑1 (MQ=255)                                               ttATTCCTCAAGGCTTAGTTGTTTTCAACTACATTAAAACTCCATGAATACGGT‑GTATGTCttt  <  1:2278419/64‑1 (MQ=255)                                                tATTCCTCAAGGCTTAGTTGTTTTCAACTACATTAAAACTCCATGAATACGGTGGTATGTCttt  <  1:1743071/64‑1 (MQ=255)                                                         |                                                       GGCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAGGGCTTAGTTGTTTTCAACTACATTAAAACTCCATGAATACGGTGGTATGTCTTT  >  W3110S.gb/4007483‑4007593

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