Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	3,824,280	T→C	pseudogene (413/687 nt)	rph →	ECK3633:JW3618:b3643; ribonuclease PH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,824,280	0	T	C	100.0%	123.9 / NA	35	pseudogene (413/687 nt)	rph	ECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (0/35);  total (0/35)

GCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGA  >  W3110S.gb/3824228‑3824333                                                     |                                                      tcTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACa                                           <  1:2894797/64‑1 (MQ=255) gCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTTGCGCTGGCAGATGCGCTACa                                           <  1:2416159/65‑1 (MQ=255) gCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACa                                           <  1:1044980/65‑1 (MQ=255) gCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACa                                           <  1:877657/65‑1 (MQ=255) gCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACa                                           <  1:1159334/65‑1 (MQ=255) gCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACa                                           <  1:339084/65‑1 (MQ=255) gCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACa                                           <  1:313714/65‑1 (MQ=255) gCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACa                                           <  1:1341532/65‑1 (MQ=255) gCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACa                                           <  1:2722861/65‑1 (MQ=255) gCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACa                                           <  1:1844884/65‑1 (MQ=255) gCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACa                                           <  1:1982073/65‑1 (MQ=255) gCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACa                                           <  1:2592042/65‑1 (MQ=255) gCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACa                                           <  1:200715/65‑1 (MQ=255) gCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACa                                           <  1:1029787/65‑1 (MQ=255)  cTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACa                                           <  1:1672458/64‑1 (MQ=255)                  tACCGCGTCGATTACGGGGGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACg                          <  1:2313966/65‑1 (MQ=255)                                            tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa  <  1:325287/63‑1 (MQ=255)                                            tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa  <  1:952684/63‑1 (MQ=255)                                            tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa  <  1:1053904/63‑1 (MQ=255)                                            tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa  <  1:728067/63‑1 (MQ=255)                                            tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa  <  1:711572/63‑1 (MQ=255)                                            tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa  <  1:698524/63‑1 (MQ=255)                                            tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa  <  1:603452/63‑1 (MQ=255)                                            tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa  <  1:415021/63‑1 (MQ=255)                                            tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa  <  1:394476/63‑1 (MQ=255)                                            tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa  <  1:1164281/63‑1 (MQ=255)                                            tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa  <  1:2258561/63‑1 (MQ=255)                                            tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa  <  1:1306874/63‑1 (MQ=255)                                            tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa  <  1:1537984/63‑1 (MQ=255)                                            tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa  <  1:2879001/63‑1 (MQ=255)                                            tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa  <  1:1715456/63‑1 (MQ=255)                                            tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa  <  1:2006018/63‑1 (MQ=255)                                            tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa  <  1:2532200/63‑1 (MQ=255)                                            tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa  <  1:2172633/63‑1 (MQ=255)                                            tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa  <  1:2255533/63‑1 (MQ=255)                                                     |                                                      GCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGA  >  W3110S.gb/3824228‑3824333

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