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

AGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGA‑GGGGAATCGAATCCATTGTAGC  >  W3110S.gb/3824469‑3824555                                                        |                                 aGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                             >  1:1255468/1‑61 (MQ=255) aGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                             >  1:883045/1‑61 (MQ=255)      aCGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑gg                      >  1:2676933/1‑62 (MQ=255)      aCGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑gg                      >  1:248491/1‑62 (MQ=255)      aCGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑gg                      >  1:579534/1‑62 (MQ=255)      aCGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑gg                      >  1:237226/1‑62 (MQ=255)      aCGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑gg                      >  1:338839/1‑62 (MQ=255)      aCGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑gg                      >  1:2306442/1‑62 (MQ=255)      aCGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑gg                      >  1:2209367/1‑62 (MQ=255)      aCGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑gg                      >  1:2191399/1‑62 (MQ=255)      aCGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑gg                      >  1:2007043/1‑62 (MQ=255)      aCGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑gg                      >  1:2829746/1‑62 (MQ=255)      aCGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑gg                      >  1:2994168/1‑62 (MQ=255)      aCGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑gg                      >  1:712348/1‑62 (MQ=255)      aCGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑gg                      >  1:769765/1‑62 (MQ=255)      aCGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑gg                      >  1:1179635/1‑62 (MQ=255)      aCGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑gg                      >  1:1398980/1‑62 (MQ=255)      aCGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑gg                      >  1:1272223/1‑62 (MQ=255)               ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                           >  1:2130065/1‑49 (MQ=255)               ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcga              >  1:3070512/1‑62 (MQ=255)               ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcga              >  1:2193010/1‑62 (MQ=255)               ggCGAGCCGTTCACCCATGAAGAGCCACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcga              >  1:3066752/1‑62 (MQ=255)                                                ttgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAga  >  1:654790/1‑40 (MQ=37)                                                ttgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc  >  1:913020/1‑41 (MQ=37)                                                ttgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc  >  1:664851/1‑41 (MQ=37)                                                ttgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc  >  1:676347/1‑41 (MQ=37)                                                ttgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc  >  1:278605/1‑41 (MQ=37)                                                ttgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc  >  1:2426702/1‑41 (MQ=37)                                                ttgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc  >  1:2365506/1‑41 (MQ=37)                                                ttgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc  >  1:2005054/1‑41 (MQ=37)                                                ttgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc  >  1:1945753/1‑41 (MQ=37)                                                ttgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc  >  1:1831913/1‑41 (MQ=37)                                                ttgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc  >  1:1599367/1‑41 (MQ=37)                                                ttgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc  >  1:1505911/1‑41 (MQ=37)                                                ttgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc  >  1:1497811/1‑41 (MQ=37)                                                   ttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc  >  1:280586/1‑38 (MQ=37)                                                        |                                 AGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGA‑GGGGAATCGAATCCATTGTAGC  >  W3110S.gb/3824469‑3824555

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