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

CATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGA‑GGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGG  >  W3110S.gb/3824458‑3824569                                                         |                                                         catTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCttgt                                                     <  1:124706/62‑1 (MQ=255) catTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCttgt                                                     <  1:4407690/62‑1 (MQ=255) catTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCttgt                                                     <  1:2820447/62‑1 (MQ=255) catTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCttgt                                                     <  1:2743519/62‑1 (MQ=255) catTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCttgt                                                     <  1:4526575/62‑1 (MQ=255) catTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCttgt                                                     <  1:4669428/62‑1 (MQ=255) catTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCttgt                                                     <  1:1266597/62‑1 (MQ=255)         tGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCtt                                                       >  1:5510/1‑52 (MQ=255)         tGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCtt                                                       >  1:4092999/1‑52 (MQ=255)         tGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCtt                                                       >  1:3948238/1‑52 (MQ=255)         tGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCtt                                                       >  1:3463073/1‑52 (MQ=255)         tGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCtt                                                       >  1:4666145/1‑52 (MQ=255)         tGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCtt                                                       >  1:2328431/1‑52 (MQ=255)         tGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCtt                                                       >  1:3345514/1‑52 (MQ=255)         tGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCtt                                                       >  1:3173795/1‑52 (MQ=255)         tGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCtt                                                       >  1:4164260/1‑52 (MQ=255)         tGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCtt                                                       >  1:607089/1‑52 (MQ=255)         tGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCtt                                                       >  1:146192/1‑52 (MQ=255)         tGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCt                                                        >  1:2737879/1‑51 (MQ=255)         tGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATCAAGAGCTACTCACCtt                                                       >  1:2704657/1‑52 (MQ=255)                     cAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCttgt                                                     <  1:2355238/42‑1 (MQ=255)                           gCGAGCCGTTCACCCATGAAGAGCTACTCACCttgt                                                     >  1:265380/1‑36 (MQ=255)                           gCGAGCCGTTCACCCATGAAGAGCTACTCACCttgt                                                     >  1:2058951/1‑36 (MQ=255)                           gCGAGCCGTTCACCCATGAAGAGCTACTCACCttgt                                                     >  1:4131208/1‑36 (MQ=255)                           gCGAGCCGTTCACCCATGAAGAGCTACTCACCttgt                                                     >  1:2026422/1‑36 (MQ=255)                           gCGAGCCGTTCACCCATGAAGAGCTACTCACCttgt                                                     >  1:3906578/1‑36 (MQ=255)                           gCGAGCCGTTCACCCATGAAGAGCTACTCACCttgt                                                     >  1:4035182/1‑36 (MQ=255)                             gAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                         <  1:3018423/62‑1 (MQ=255)                             gAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                         <  1:594926/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)                                ccGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAt                      >  1:3609359/1‑62 (MQ=255)                                ccGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAt                      >  1:2032571/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)                                    tCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTa                  >  1:840357/1‑62 (MQ=255)                                                    aCTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcgg  >  1:1582201/1‑62 (MQ=255)                                                         |                                                         CATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGA‑GGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGG  >  W3110S.gb/3824458‑3824569

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