Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	3,224,816	C→T	H147Y (CAT→TAT)	ebgC →	cryptic beta‑D‑galactosidase, beta subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,224,816	0	C	T	100.0%	104.9 / NA	33	H147Y (CAT→TAT)	ebgC	cryptic beta‑D‑galactosidase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/17);  total (16/17)

ATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCAT  >  W3110S.gb/3224763‑3224874                                                      |                                                           aTAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCTataac                                                       <  1:3871554/59‑1 (MQ=255) aTAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCTataac                                                       <  1:3078160/59‑1 (MQ=255) aTAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCTataac                                                       <  1:698070/59‑1 (MQ=255) aTAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCTataac                                                       <  1:3072568/59‑1 (MQ=255) aTAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCTataac                                                       <  1:2712270/59‑1 (MQ=255) aTAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCTataac                                                       <  1:2175661/59‑1 (MQ=255) aTAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCTataac                                                       <  1:2581607/59‑1 (MQ=255)     cgcgGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCTataacaaataac                                                >  1:4444969/1‑62 (MQ=255)           cAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCTataacaaataac                                                >  1:1671710/1‑56 (MQ=255)           cAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCTataacaaataac                                                >  1:2833111/1‑56 (MQ=255)           cAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCTataacaaataac                                                >  1:4644282/1‑56 (MQ=255)           cAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCTataacaaataac                                                >  1:569369/1‑56 (MQ=255)            aaaaaaGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCTataacaaa                                                    <  1:584475/51‑1 (MQ=255)            aaaaaaGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCTataacaaa                                                    <  1:315611/51‑1 (MQ=255)            aaaaaaGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCTataacaaa                                                    <  1:4311277/51‑1 (MQ=255)             aaaaaGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCTATAACAAATAACAACTAcgg                                        >  1:609243/1‑62 (MQ=255)             aaaaaGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCTATAACAAATAACAACTAcgg                                        >  1:1376048/1‑62 (MQ=255)             aaaaaGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCTATAACAAATAACAACTAcgg                                        >  1:977765/1‑62 (MQ=255)             aaaaaGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCTATAACAAATAACAACTAcgg                                        >  1:3832054/1‑62 (MQ=255)                     gTTCTCAAAGTCACCATCGAAGATGGTTATTTCTATAACAAATAACAACTACGGCGGCaaaa                                <  1:3720793/62‑1 (MQ=255)                        ctcAAAGTCACCATCGAAGATGGTTATTTCTATAACAAATAACAACTACGGCGGCaaaa                                <  1:2413857/59‑1 (MQ=255)                                 aCCATCGAAGATGGTTATTTCTATAACAAATAACAACTACGGCGGCAAAAGGAGTTTgc                       <  1:3673062/59‑1 (MQ=255)                                 aCCATCGAAGATGGTTATTTCTATAACAAATAACAACTACGGCGGCAAAAGGAGTTTgc                       <  1:4603506/59‑1 (MQ=255)                                 aCCATCGAAGATGGTTATTTCTATAACAAATAACAACTACGGCGGCAAAAGGAGTTTgc                       <  1:629114/59‑1 (MQ=255)                                 aCCATCGAAGATGGTTATTTCTATAACAAATAACAACTACGGCGGCAAAAGGAGTTTgc                       <  1:790026/59‑1 (MQ=255)                                       gAAGATGGTTATTTCTATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGc              >  1:2165656/1‑62 (MQ=255)                                           aTGGTTATTTCTATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTAc           >  1:1882611/1‑61 (MQ=255)                                              gTTATTTCTATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTa       >  1:3431491/1‑62 (MQ=255)                                              gTTATTTCTATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTa       >  1:2594914/1‑62 (MQ=255)                                              gTTATTTCTATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTa       >  1:1786935/1‑62 (MQ=255)                                              gTTATTTCTATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTa       >  1:1784862/1‑62 (MQ=255)                                               ttATTTCTATAACAAATAACAACTACGGCGGCAAAAGGAGTTTgc                       <  1:3306108/45‑1 (MQ=255)                                                   ttCTATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCAt  >  1:1948770/1‑62 (MQ=255)                                                      |                                                           ATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCAT  >  W3110S.gb/3224763‑3224874

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