Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	547,694	A→G	pseudogene (114/261 nt) pseudogene (114/1259 nt)	ylbE → ylbE →	ECK0512:JW0507:b4507; hypothetical protein, N‑ter fragment ECK0512:JW0508+JW0507:b4572; hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	547,694	0	A	G	100.0%	134.8 / NA	37	pseudogene (114/261 nt) pseudogene (114/1259 nt)	ylbE ylbE	ECK0512:JW0507:b4507; hypothetical protein, N‑ter fragment ECK0512:JW0508+JW0507:b4572; hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/18);  total (19/18)

CACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAAGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAG  >  W3110S.gb/547642‑547739                                                     |                                              ctcacTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAac                                       <  1:1522110/59‑1 (MQ=255) cacaCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAac                                       <  1:45156/61‑1 (MQ=255) cacaCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAac                                       <  1:3034295/61‑1 (MQ=255) cacaCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAac                                       <  1:2842028/61‑1 (MQ=255)        gctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCtt                                >  1:993061/1‑61 (MQ=255)        gctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCtt                                >  1:1545386/1‑61 (MQ=255)        gctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCtt                                >  1:3575043/1‑61 (MQ=255)        gctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCtt                                >  1:3537143/1‑61 (MQ=255)        gctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCtt                                >  1:3465067/1‑61 (MQ=255)        gctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCtt                                >  1:3336656/1‑61 (MQ=255)        gctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACGGCtt                                >  1:4589011/1‑61 (MQ=255)            gATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACg                            >  1:3073831/1‑61 (MQ=255)            gATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGc                           >  1:2691505/1‑62 (MQ=255)            gATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGc                           >  1:249307/1‑62 (MQ=255)            gATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGc                           >  1:96931/1‑62 (MQ=255)            gATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGc                           >  1:3421544/1‑62 (MQ=255)              tgtgCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGc                           <  1:4313366/60‑1 (MQ=255)              tgtgCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCg                         <  1:909167/62‑1 (MQ=255)              tgtgCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCg                         <  1:1517873/62‑1 (MQ=255)              tgtgCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCg                         <  1:4259375/62‑1 (MQ=255)              tgtgCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCg                         <  1:4195731/62‑1 (MQ=255)               gtgCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCg                         <  1:2892867/61‑1 (MQ=255)                    aCCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGgccg                    <  1:1367241/61‑1 (MQ=255)                    aCCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGgccg                    <  1:2576909/61‑1 (MQ=255)                    aCCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGgccg                    <  1:3225641/61‑1 (MQ=255)                     ccGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaaaa                                        >  1:2395573/1‑40 (MQ=255)                     ccGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaaaa                                        >  1:4501760/1‑40 (MQ=255)                     ccGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaaaa                                        >  1:475051/1‑40 (MQ=255)                     ccGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaaaa                                        >  1:613607/1‑40 (MQ=255)                     ccGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaaaa                                        >  1:620177/1‑40 (MQ=255)                              tcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCt          <  1:1821633/61‑1 (MQ=255)                              tcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCt          <  1:1728767/61‑1 (MQ=255)                              tcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCt          <  1:1254740/61‑1 (MQ=255)                                     tCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGgag  >  1:2902193/1‑62 (MQ=255)                                      cAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAgg    <  1:1775848/59‑1 (MQ=255)                                      cAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAgg    <  1:2007328/59‑1 (MQ=255)                                      cAGCGAACTAAACAAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGgag  >  1:2458213/1‑61 (MQ=255)                                                     |                                              CACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAAGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAG  >  W3110S.gb/547642‑547739

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