Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	1,210,013	T→C	L102P (CTG→CCG)	ymfS →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,210,013	0	T	C	100.0%	160.8 / NA	47	L102P (CTG→CCG)	ymfS	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (27/20);  total (27/20)

GATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGC  >  W3110S.gb/1209949‑1210078                                                                 |                                                                  gATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCCGAAAc                                                              >  1:1577002/1‑70 (MQ=255)     ttATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCCGAAACAata                                                          >  1:4485260/1‑70 (MQ=255)              gCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTgg                                                >  1:4225885/1‑71 (MQ=255)                 gACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTg                                             <  1:3640493/71‑1 (MQ=255)                     gTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGAtt                                         <  1:3793665/71‑1 (MQ=255)                       ggCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTa                                        >  1:2874022/1‑70 (MQ=255)                             gTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTgg                                   >  1:4088673/1‑69 (MQ=255)                             gTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTgg                                   >  1:3127358/1‑69 (MQ=255)                              tAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAg                                >  1:1643773/1‑71 (MQ=255)                              tAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAg                                >  1:2628502/1‑71 (MQ=255)                              tAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAg                                >  1:4040110/1‑71 (MQ=255)                              tAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAg                                >  1:1266635/1‑71 (MQ=255)                               aaGGCTGCATTGGGAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAg                                >  1:1137816/1‑70 (MQ=255)                                aGGCTGCATTGGGAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTg                                             >  1:2293133/1‑56 (MQ=255)                                aGGCTGCATTGGGAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGc                               <  1:2736042/70‑1 (MQ=255)                                aGGCTGCATTGGGAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGc                               <  1:2437774/70‑1 (MQ=255)                                 ggCTGCATTGGGAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGa                                  <  1:1378594/66‑1 (MQ=255)                                  gCTGCATTGGGAAGACTT‑AAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACt                            <  1:1436843/70‑1 (MQ=25)                                   cTGCATTGGGAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCa                              >  1:2532634/1‑68 (MQ=255)                                   cTGCATTGGGAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCa                              >  1:2492587/1‑68 (MQ=255)                                      cATTGGGAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTgg                          <  1:2533996/69‑1 (MQ=255)                                      cATTGGGAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTgg                          <  1:3696711/69‑1 (MQ=255)                                      cATTGGGAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTgg                          <  1:4592924/69‑1 (MQ=255)                                      cATTGGGAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTgg                          <  1:3663730/69‑1 (MQ=255)                                        ttGGGAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCa                              <  1:2312720/63‑1 (MQ=255)                                          gggAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTgg                    >  1:4344313/1‑71 (MQ=255)                                          gggAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTgg                    >  1:1299676/1‑71 (MQ=255)                                          gggAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTgg                    >  1:1372694/1‑71 (MQ=255)                                          gggAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTgg                    >  1:1811054/1‑71 (MQ=255)                                          gggAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTgg                    >  1:2197127/1‑71 (MQ=255)                                          gggAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTgg                    >  1:2202432/1‑71 (MQ=255)                                          gggAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTg                     >  1:2569325/1‑70 (MQ=255)                                            gAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGtt                  <  1:534280/71‑1 (MQ=255)                                             aaGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGtt                  <  1:812299/70‑1 (MQ=255)                                               gACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGc                               <  1:4131062/55‑1 (MQ=255)                                                aCTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGa                                  <  1:3540652/51‑1 (MQ=255)                                                 cTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATa              >  1:270397/1‑70 (MQ=255)                                                 cTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATa              >  1:492882/1‑70 (MQ=255)                                                 cTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATa              >  1:1958381/1‑70 (MQ=255)                                                  ttAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACa            >  1:167038/1‑71 (MQ=255)                                                     aagaagATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATcc         <  1:3035860/71‑1 (MQ=255)                                                     aagaagATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATcc         <  1:861719/71‑1 (MQ=255)                                                      agaagaTGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATcc         <  1:1350686/70‑1 (MQ=255)                                                          gaTGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGc                               <  1:544317/44‑1 (MQ=255)                                                          gaTGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTgc    >  1:3640373/1‑71 (MQ=255)                                                          gaTGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTgc    >  1:1431297/1‑71 (MQ=255)                                                            tGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTgcgc  <  1:3462538/71‑1 (MQ=255)                                                                 |                                                                  GATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGC  >  W3110S.gb/1209949‑1210078

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