Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	3,824,537	+G	pseudogene (670/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,533	1	.	G	100.0%	124.3 / NA	34	pseudogene (666/687 nt)	rph	ECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (9/25);  total (9/25)

GAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGA‑GGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAA  >  W3110S.gb/3824486‑3824579                                                 |                                               gAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                   <  1:3018423/62‑1 (MQ=255) gAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                   <  1:2983140/62‑1 (MQ=255) gAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                   <  1:2724983/62‑1 (MQ=255) gAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                   <  1:2141778/62‑1 (MQ=255) gAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                   <  1:594926/62‑1 (MQ=255)    ccGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAt                                >  1:3609359/1‑62 (MQ=255)    ccGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAt                                >  1:2032571/1‑62 (MQ=255)        tCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTa                            >  1:840357/1‑62 (MQ=255)        tCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTa                            >  1:3504438/1‑62 (MQ=255)        tCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTa                            >  1:1246880/1‑62 (MQ=255)        tCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTa                            >  1:1529904/1‑62 (MQ=255)        tCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTa                            >  1:1467798/1‑62 (MQ=255)                        aCTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcgg            >  1:1582201/1‑62 (MQ=255)                           caccTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcggc           <  1:1921350/57‑1 (MQ=255)                           caccTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcggc           <  1:668311/57‑1 (MQ=255)                           caccTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcggc           <  1:2411822/57‑1 (MQ=255)                           caccTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcggc           <  1:2431006/57‑1 (MQ=255)                           caccTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcggc           <  1:2454561/57‑1 (MQ=255)                           caccTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcggc           <  1:3575317/57‑1 (MQ=255)                           caccTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcggc           <  1:3364885/57‑1 (MQ=255)                            accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGa                        <  1:3744210/44‑1 (MQ=255)                            accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGa                        <  1:695376/44‑1 (MQ=255)                            accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGa                        <  1:789792/44‑1 (MQ=255)                            accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGa                        <  1:2751228/44‑1 (MQ=255)                            accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGa                        <  1:3521266/44‑1 (MQ=255)                            accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGa                        <  1:3050785/44‑1 (MQ=255)                            accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGa                        <  1:2229431/44‑1 (MQ=255)                            accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGa                        <  1:1757263/44‑1 (MQ=255)                            accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGa                        <  1:1638313/44‑1 (MQ=255)                            accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGa                        <  1:1516484/44‑1 (MQ=255)                            accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGa                        <  1:137915/44‑1 (MQ=255)                            accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGa                        <  1:1315158/44‑1 (MQ=255)                              cTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGAGCTg       <  1:764308/61‑1 (MQ=255)                                  ttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCaaa  >  1:3249756/1‑62 (MQ=255)                                                 |                                               GAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGA‑GGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAA  >  W3110S.gb/3824486‑3824579

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