Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	4,588,402	C→T	L1032L (CTG→CTA)	hsdR ←	endonuclease R

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,588,402	0	C	T	100.0%	176.9 / NA	53	L1032L (CTG→CTA)	hsdR	endonuclease R
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (34/19);  total (34/19)

GCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGC  >  W3110S.gb/4588347‑4588456                                                        |                                                       gcgGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGa                                                   >  1:1302791/1‑61 (MQ=255) gcgGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGa                                                   >  1:620914/1‑61 (MQ=255)     gCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCaa                                              >  1:343186/1‑62 (MQ=255)     gCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCaa                                              >  1:2016311/1‑62 (MQ=255)     gCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCaa                                              >  1:3034663/1‑62 (MQ=255)     gCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCaa                                              >  1:1724847/1‑62 (MQ=255)     gCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCaa                                              >  1:2064039/1‑62 (MQ=255)     gCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCaa                                              >  1:2620882/1‑62 (MQ=255)     gCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCaa                                              >  1:3315556/1‑62 (MQ=255)     gCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCaa                                              >  1:161794/1‑62 (MQ=255)     gCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCaa                                              >  1:1554690/1‑62 (MQ=255)     gCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCaa                                              >  1:543701/1‑62 (MQ=255)     gCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCaa                                              >  1:1000075/1‑62 (MQ=255)     gCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCaa                                              >  1:3764333/1‑62 (MQ=255)     gCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCaa                                              >  1:4072921/1‑62 (MQ=255)     gCGATTAATAACTGCCTGAAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCaa                                              >  1:2016974/1‑62 (MQ=255)        aTTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCAAAgg                                           >  1:1706262/1‑62 (MQ=255)        aTTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCAAAgg                                           >  1:2115700/1‑62 (MQ=255)        aTTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCAAAgg                                           >  1:3042326/1‑62 (MQ=255)        aTTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCAAAgg                                           >  1:1701067/1‑62 (MQ=255)        aTTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCAAAgg                                           >  1:1357840/1‑62 (MQ=255)             taaCTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCAAAGGCTTc                                       <  1:2112380/61‑1 (MQ=255)             taaCTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCAAAGGCTTc                                       <  1:511592/61‑1 (MQ=255)             taaCTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCAAAGGCTTc                                       <  1:805050/61‑1 (MQ=255)             taaCTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCAAAGGCTTc                                       <  1:2353969/61‑1 (MQ=255)             taaCTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCAAAGGCTTc                                       <  1:4215492/61‑1 (MQ=255)             taaCTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCAAAGGCTTc                                       <  1:1732233/61‑1 (MQ=255)              aaCTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCAAAGGCTTc                                       <  1:3765602/60‑1 (MQ=255)               aCTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCAAAGGCTTc                                       <  1:2716124/59‑1 (MQ=255)                 tgcctgcAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCAAAGGCTTCGAGGa                                  >  1:2889842/1‑62 (MQ=255)                 tgcctgcAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCAAAGGCTTCGAGGa                                  >  1:1731238/1‑62 (MQ=255)                       cAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCAATGGCTTCGAGGAAATCCt                            >  1:424956/1‑62 (MQ=255)                       cAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATCCt                            >  1:2183537/1‑62 (MQ=255)                       cAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATCCt                            >  1:4186446/1‑62 (MQ=255)                       cAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATCCt                            >  1:1235144/1‑62 (MQ=255)                       cAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATCCt                            >  1:1548336/1‑62 (MQ=255)                            ccGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATCCTgcggc                       <  1:4618418/62‑1 (MQ=255)                            ccGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATCCTgcggc                       <  1:3782654/62‑1 (MQ=255)                            ccGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATCCTgcggc                       <  1:1285352/62‑1 (MQ=255)                            ccGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATCCTgcggc                       <  1:1638577/62‑1 (MQ=255)                            ccGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATCCTgcggc                       <  1:1718015/62‑1 (MQ=255)                                  gcgcGTTCGGGGAACGTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTc                  >  1:4709085/1‑61 (MQ=255)                                  gcgcGTTCGGGGAACGTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTc                  >  1:4049935/1‑61 (MQ=255)                                  gcgcGTTCGGGGAACGTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTc                  >  1:873276/1‑61 (MQ=255)                                  gcgcGTTCGGGGAACGTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTc                  >  1:1625071/1‑61 (MQ=255)                                  gcgcGTTCGGGGAACGTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTc                  >  1:1082372/1‑61 (MQ=255)                                             gAACGTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATCCTgcggc                       <  1:2614523/45‑1 (MQ=255)                                                cGTTGCACTAGCCAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACa    >  1:348817/1‑61 (MQ=255)                                                 gTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACa    <  1:3186668/60‑1 (MQ=255)                                                 gTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGc  <  1:3819296/62‑1 (MQ=255)                                                 gTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGc  <  1:2527573/62‑1 (MQ=255)                                                 gTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGc  <  1:873410/62‑1 (MQ=255)                                                 gTTGCACTAGCGAGTCAAAGACTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGc  <  1:2093568/62‑1 (MQ=255)                                                        |                                                       GCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGC  >  W3110S.gb/4588347‑4588456

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