Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	547,831:1	+G	100%	pseudogene (251/261 nt) pseudogene (251/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,831	1	.	G	96.3%	93.0 / ‑3.3	27	pseudogene (251/261 nt) pseudogene (251/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 . (0/1);  new base G (14/12);  total (14/13)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

TCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑GGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCC  >  W3110S.gb/547778‑547878                                                       |                                                tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTc                                  >  1:1140813/1‑70 (MQ=255) tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCa                                 >  1:600937/1‑71 (MQ=255) tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCa                                 >  1:480862/1‑71 (MQ=255) tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCa                                 >  1:418359/1‑71 (MQ=255) tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCa                                 >  1:1115275/1‑71 (MQ=255) tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCa                                 >  1:285637/1‑71 (MQ=255) tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCa                                 >  1:206899/1‑71 (MQ=255) tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑g                                                >  1:1481335/1‑55 (MQ=255)           ggTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCa                       >  1:337376/1‑71 (MQ=255)           ggTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCa                       >  1:580590/1‑71 (MQ=255)           ggTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCa                       >  1:1443131/1‑71 (MQ=255)              tGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCa                    >  1:269389/1‑71 (MQ=255)              tGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCa                    >  1:1682880/1‑71 (MQ=255)               gggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAt                   >  1:99530/1‑71 (MQ=255)               gggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAt                   >  1:822698/1‑71 (MQ=255)                  cGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                <  1:976669/71‑1 (MQ=255)                  cGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                <  1:93269/71‑1 (MQ=255)                  cGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                <  1:1228/71‑1 (MQ=255)                  cGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                <  1:904542/71‑1 (MQ=255)                  cGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                <  1:881265/71‑1 (MQ=255)                  cGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                <  1:880119/71‑1 (MQ=255)                               gcAGGCGCTGGCAATACTGGAGCA‑GGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCTGGcc  <  1:1771996/71‑1 (MQ=255)                                cAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc  <  1:713864/71‑1 (MQ=255)                                cAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc  <  1:73054/71‑1 (MQ=255)                                cAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc  <  1:190829/71‑1 (MQ=255)                                cAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc  <  1:1647799/71‑1 (MQ=255)                                cAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc  <  1:1732875/71‑1 (MQ=255)                                cAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc  <  1:255522/71‑1 (MQ=255)                                                       |                                                TCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCA‑GGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCC  >  W3110S.gb/547778‑547878

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