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%	77.2 / NA	24	pseudogene (657/687 nt)	rph	ECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (12/12);  total (12/12)

GGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGA‑GGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATT  >  W3110S.gb/3824470‑3824585                                                       |                                                               ggggACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑gggggaa                                               >  1:76762/1‑68 (MQ=255)      cGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑ggggg                                                 <  1:1648494/64‑2 (MQ=255)      cGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑ggggg                                                 <  1:1191242/64‑2 (MQ=255)      cGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑ggggg                                                 <  1:1803480/64‑2 (MQ=255)      cGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑ggggg                                                 <  1:3851030/64‑2 (MQ=255)      cGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑ggggg                                                 <  1:3007664/64‑2 (MQ=255)      cGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑ggggg                                                 <  1:3185453/64‑2 (MQ=255)        gCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑gggggaa                                               >  1:1792237/1‑61 (MQ=255)        gCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑gggggaa                                               >  1:3768949/1‑61 (MQ=255)        gCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑gggggaa                                               >  1:2424052/1‑61 (MQ=255)        gCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑gggggaa                                               >  1:4122760/1‑61 (MQ=255)           gAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                         >  1:4070215/1‑68 (MQ=255)           gAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                         >  1:3419799/1‑68 (MQ=255)           gAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                         >  1:3297078/1‑68 (MQ=255)           gAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                         >  1:2572014/1‑68 (MQ=255)           gAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                         >  1:969814/1‑68 (MQ=255)           gAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                         >  1:1884782/1‑68 (MQ=255)                                      cTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcg                                             <  1:714754/37‑1 (MQ=21)                                      cTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcg                                             <  1:997401/37‑1 (MQ=21)                                      cTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcg                                             <  1:2986256/37‑1 (MQ=21)                                      cTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcg                                             <  1:2034645/37‑1 (MQ=21)                                      cTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcg                                             <  1:1686837/37‑1 (MQ=21)                                            accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACt      >  1:3095601/3‑70 (MQ=255)                                                 gttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAtt  <  1:2252404/69‑1 (MQ=255)                                                       |                                                               GGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGA‑GGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATT  >  W3110S.gb/3824470‑3824585

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