Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	2,599,304	G→A	M57I (ATG→ATA)	bcp →	thiol peroxidase, thioredoxin‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,599,304	0	G	A	100.0%	100.1 / NA	32	M57I (ATG→ATA)	bcp	thiol peroxidase, thioredoxin‑dependent
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/20);  total (12/20)

TTATTTCTACCCGAAAGCCATGACCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACC  >  W3110S.gb/2599238‑2599349                                                                   |                                              ttATTTCTACCCGAAAGCCATGACCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATAGAt                                            >  1:1238866/1‑70 (MQ=255)      tCTACCCGAAAGCCATTACCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATAGa                                             >  1:368521/1‑64 (MQ=255)      tCTACCCGAAAGCCATGACCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATAGa                                             >  1:614889/1‑64 (MQ=255)      tCTACCCGAAAGCCATGACCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATAGa                                             >  1:3369986/1‑64 (MQ=255)      tCTACCCGAAAGCCATGACCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATAGa                                             >  1:2947873/1‑64 (MQ=255)      tCTACCCGAAAGCCATGACCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATAGa                                             >  1:2889790/1‑64 (MQ=255)                   cATGACCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATAGATGAGTTGAAAAAAGCGGGCg                         >  1:103045/1‑71 (MQ=255)                   cATGACCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATAGATGAGTTGAAAAAAGCGGGCg                         >  1:4054695/1‑71 (MQ=255)                   cATGACCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATAGATGAGTTGAAAAAAGCGGGCg                         >  1:3836292/1‑71 (MQ=255)                   cATGACCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATAGATGAGTTGAAAAAAGCGGGCg                         >  1:2032342/1‑71 (MQ=255)                        cccccGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATAGATGAGTTGAAAAAAGCGGGCGTTGa                     <  1:38709/70‑1 (MQ=255)                             ggCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATAGATGAGTTGAAAAAAGCGGGCGTTGATGTGCt               <  1:1972093/71‑1 (MQ=255)                                  tACCGTACAGGCCTGCGGCTTACGCGATAACATAGATGAGTTGaaaaa                                 <  1:787457/48‑1 (MQ=255)                                  tACCGTACAGGCCTGCGGCTTACGCGATAACATAGATGAGTTGaaaaa                                 <  1:3552513/48‑1 (MQ=255)                                  tACCGTACAGGCCTGCGGCTTACGCGATAACATAGATGAGTTGaaaaa                                 <  1:3467349/48‑1 (MQ=255)                                  tACCGTACAGGCCTGCGGCTTACGCGATAACATAGATGAGTTGaaaaa                                 <  1:1467272/48‑1 (MQ=255)                                  tACCGTACAGGCCTGCGGCTTACGCGATAACATAGATGAGTTGaaaaa                                 <  1:1575176/48‑1 (MQ=255)                                  tACCGTACAGGCCTGCGGCTTACGCGATAACATAGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTggg            <  1:120777/69‑1 (MQ=255)                                  tACCGTACAGGCCTGCGGCTTACGCGATAACATAGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTggg            <  1:2854073/69‑1 (MQ=255)                                  tACCGTACAGGCCTGCGGCTTACGCGATAACATAGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTggg            <  1:448967/69‑1 (MQ=255)                                  tACCGTACAGGCCTGCGGCTTACGCGATAACATAGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTggg            <  1:444061/69‑1 (MQ=255)                                  tACCGTACAGGCCTGCGGCTTACGCGATAACATAGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTggg            <  1:4011065/69‑1 (MQ=255)                                  tACCGTACAGGCCTGCGGCTTACGCGATAACATAGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTggg            <  1:2974526/69‑1 (MQ=255)                                  tACCGTACAGGCCTGCGGCTTACGCGATAACATAGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTggg            <  1:2836231/69‑1 (MQ=255)                                  tACCGTACAGGCCTGCGGCTTACGCGATAACATAGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGt           <  1:1828181/70‑1 (MQ=255)                                  tACCGTACAGGCCTGCGGCTTACGCGATAACATAGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGt           <  1:2777844/70‑1 (MQ=255)                                  tACCGTACAGGCCTGCGGCTTACGCGATAACATAGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGt           <  1:129611/70‑1 (MQ=255)                                    ccGTACAGGCCTGCGGCTTACGCGATAACATAGATGAGTTGAAAAAAg                               >  1:1393510/1‑48 (MQ=255)                                    ccGTACAGGCCTGCGGCTTACGCGATAACATAGATGAGTTGAAAAAAg                               >  1:1354833/1‑48 (MQ=255)                                          aGGCCTGCGGCTTACGCGATAACATAGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCAcc  <  1:277589/71‑1 (MQ=255)                                                 cGGCTTACGCGATAACATAGATTAGTTGAAAAAAGCGGGCGTTGa                     <  1:1645918/45‑1 (MQ=38)                                                 cGGCTTACGCGATAACATAGATGAGTTGAAAAAAGCGGGCGTTGa                     <  1:647685/45‑1 (MQ=255)                                                                   |                                              TTATTTCTACCCGAAAGCCATGACCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACC  >  W3110S.gb/2599238‑2599349

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