Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	3,824,362	A→T	pseudogene (495/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,362	0	A	T	100.0%	87.7 / NA	28	pseudogene (495/687 nt)	rph	ECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (18/10);  total (18/10)

AGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGA  >  W3110S.gb/3824313‑3824403                                                  |                                          aGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACgg                                >  1:2180372/1‑61 (MQ=255) aGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACgg                                >  1:318934/1‑61 (MQ=255) aGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACgg                                >  1:2396765/1‑61 (MQ=255) aGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACgg                                >  1:2965032/1‑61 (MQ=255) aGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACgg                                >  1:2451466/1‑61 (MQ=255)                        gggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATa        <  1:702893/62‑1 (MQ=255)                        gggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATa        <  1:2028880/62‑1 (MQ=255)                        gggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATa        <  1:2180705/62‑1 (MQ=255)                        gggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATa        <  1:443546/62‑1 (MQ=255)                        gggATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATa        <  1:2837426/62‑1 (MQ=255)                                     cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg   >  1:2557722/1‑54 (MQ=255)                                     cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg   >  1:940742/1‑54 (MQ=255)                                     cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg   >  1:908756/1‑54 (MQ=255)                                     cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg   >  1:82213/1‑54 (MQ=255)                                     cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg   >  1:55672/1‑54 (MQ=255)                                     cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg   >  1:2788960/1‑54 (MQ=255)                                     cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg   >  1:2447823/1‑54 (MQ=255)                                     cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg   >  1:2192026/1‑54 (MQ=255)                                     cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg   >  1:161145/1‑54 (MQ=255)                                     cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg   >  1:1452667/1‑54 (MQ=255)                                     cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg   >  1:1360138/1‑54 (MQ=255)                                     cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg   >  1:1269128/1‑54 (MQ=255)                                     cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATAAGTTg   >  1:462157/1‑54 (MQ=255)                                           cTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGa  <  1:1144801/49‑1 (MQ=25)                                           cTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGa  <  1:2842257/49‑1 (MQ=25)                                           cTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGa  <  1:2252996/49‑1 (MQ=25)                                           cTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGa  <  1:799110/49‑1 (MQ=25)                                           cTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGa  <  1:1150037/49‑1 (MQ=25)                                                  |                                          AGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGA  >  W3110S.gb/3824313‑3824403

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