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%	71.3 / NA	22	pseudogene (367/393 nt)	yhiK	ECK3474:JW5674:b3489; hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (2/20);  total (2/20)

TATGGGGCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAGGGCTTAGTTGTTTTCAACTACATTAAAACTCCATGAATACG  >  W3110S.gb/4007478‑4007580                                                              |                                          tatGGGGCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAagg                                         >  1:1352631/1‑64 (MQ=255) tatGGGGCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAagg                                         >  1:497872/1‑64 (MQ=255) tatGGGGCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAagg                                         >  1:2638764/1‑64 (MQ=255)                  tGAAGGAAGCACGTGATAGAAATAATAAAATAATTTATTCCTCAAGGCTTAGTTGTTTTCAACt                        <  1:1455495/64‑1 (MQ=255)                  tGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAGTTGTTTTCAACt                        <  1:2627102/64‑1 (MQ=255)                  tGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAGTTGTTTTCAACt                        <  1:985790/64‑1 (MQ=255)                  tGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAGTTGTTTTCAACt                        <  1:962117/64‑1 (MQ=255)                  tGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAGTTGTTTTCAACt                        <  1:854432/64‑1 (MQ=255)                  tGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAGTTGTTTTCAACt                        <  1:2709501/64‑1 (MQ=255)                  tGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAGTTGTTTTCAACt                        <  1:2655322/64‑1 (MQ=255)                  tGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAGTTGTTTTCAACt                        <  1:2628149/64‑1 (MQ=255)                  tGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAGTTGTTTTCAACt                        <  1:1020379/64‑1 (MQ=255)                  tGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAGTTGTTTTCAACt                        <  1:2383067/64‑1 (MQ=255)                  tGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAGTTGTTTTCAACt                        <  1:1816166/64‑1 (MQ=255)                  tGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAGTTGTTTTCAACt                        <  1:17527/64‑1 (MQ=255)                  tGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAGTTGTTTTCAACt                        <  1:1039172/64‑1 (MQ=255)                               tGATAAAAATAATAAAATAATTTATTCCTCAATGATTATTTGTTTTCAACt                        <  1:1497575/51‑1 (MQ=255)                                       ataataaaataatTTATTCCTCAAGGCTTAGTTGTTTTCAACTACATTAAAACTCCATGAATACg  <  1:2607419/65‑1 (MQ=255)                                       ataataaaataatTTATTCCTCAAGGCTTAGTTGTTTTCAACTACATTAAAACTCCATGAATACg  >  1:2668170/1‑65 (MQ=255)                                       ataataaaataatTTATTCCTCAAGGCTTAGTTGTTTTCAACTACATTAAAACTCCATGAATACg  <  1:1549034/65‑1 (MQ=255)                                       ataataaaataatTTATTCCTCAAGGCTTAGTTGTTTTCAACTACATTAAAACTCCATGAATACg  <  1:2931113/65‑1 (MQ=255)                                       ataataaaataatTTATTCCTCAAGGCTTAGTTGTTTTCAACTACATTAAAACTCCATGAATACg  <  1:557278/65‑1 (MQ=255)                                       ataataaaataatTTATTCCTCAAGGCTTAGTTGTTTTCAACTACATTAAAACTCCATGAATACg  <  1:621569/65‑1 (MQ=255)                                       ataataaaataatTTATTCCTCAAGGCTTAGTTGTTTTCAACTACATTAAAACTCCATGAATACg  <  1:1175850/65‑1 (MQ=255)                                       ataataaaataatTTATTCCTCAAGGCTTAGTTGTTTTCAACTACATTAAAACTCCATGAATACg  >  1:1131632/1‑65 (MQ=255)                                                              |                                          TATGGGGCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAGGGCTTAGTTGTTTTCAACTACATTAAAACTCCATGAATACG  >  W3110S.gb/4007478‑4007580

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