Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	1,930,964	T→A	D153V (GAT→GTT)	yebE ←	conserved hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,930,964	0	T	A	100.0%	108.7 / NA	35	D153V (GAT→GTT)	yebE	conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (13/22);  total (13/22)

GAAATAGATTTCCAGCGCCTCTTCTTCATTGCGGACCCCGGTAGCCAGGCGTTGTGGATCCAGCGGTTGTTCGATTGCCTGCTCAATGAGTACACGCCCCT  >  W3110S.gb/1930906‑1931006                                                           |                                           gAAATAGATTTCCAGCGCCTCTTCTTCATTGCGGACCCCGGTAGCCAGGCGTTGTGGAACCa                                         <  1:3743784/62‑1 (MQ=255) gAAATAGATTTCCAGCGCCTCTTCTTCATTGCGGACCCCGGTAGCCAGGCGTTGTGGAACCa                                         <  1:2787014/62‑1 (MQ=255) gAAATAGATTTCCAGCGCCTCTTCTTCATTGCGGACCCCGGTAGCCAGGCGTTGTGGAACCa                                         <  1:3931172/62‑1 (MQ=255)  aaaTAGATTTCCAGCGCCTCTTCTTCATTGCGGACCCCGGTAGCCAGGCGTTGTGGAACCa                                         <  1:2088366/61‑1 (MQ=255)             cAGCGCCTCTTCTTCATTGCGGACCCCGGTAGCCAGGCGTTGTGGAACCAGCGGTTGTTc                               >  1:137554/1‑60 (MQ=255)             cAGCGCCTCTTCTTCATTGCGGACCCCGGTAGCCAGGCGTTGTGGAACCAGCGGTTGTTCg                              >  1:4509305/1‑61 (MQ=255)             cAGCGCCTCTTCTTCATTGCGGACCCCGGTAGCCAGGCGTTGTGGAACCAGCGGTTGTTCg                              >  1:4148292/1‑61 (MQ=255)             cAGCGCCTCTTCTTCATTGCGGACCCCGGTAGCCAGGCGTTGTGGAACCAGCGGTTGTTCg                              >  1:4313840/1‑61 (MQ=255)             cAGCGCCTCTTCTTCATTGCGGACCCCGGTAGCCAGGCGTTGTGGAACCAGCGGTTGTTCGa                             >  1:1604955/1‑62 (MQ=255)             cAGCGCCTCTTCTTCATTGCGGACCCCGGTAGCCAGGCGTTGTGGAACCAGCGGTTGTTCGa                             >  1:2320532/1‑62 (MQ=255)             cAGCGCCTCTTCTTCATTGCGGACCCCGGTAGCCAGGCGTTGTGGAACCAGCGGTTGTTCGa                             >  1:1056947/1‑62 (MQ=255)                cgcCTCTTCTTCATTGCGGACCCCGGTAGCCAGGCGTTGTGGAACCAGCGGTTGTTCGATTg                          <  1:156140/62‑1 (MQ=255)                cgcCTCTTCTTCATTGCGGACCCCGGTAGCCAGGCGTTGTGGAACCAGCGGTTGTTCGATTg                          <  1:4150045/62‑1 (MQ=255)                cgcCTCTTCTTCATTGCGGACCCCGGTAGCCAGGCGTTGTGGAACCAGCGGTTGTTCGATTg                          <  1:2032789/62‑1 (MQ=255)                cgcCTCTTCTTCATTGCGGACCCCGGTAGCCAGGCGTTGTGGAACCAGCGGTTGTTCGATTg                          <  1:13946/62‑1 (MQ=255)                cgcCTCTTCTTCATTGCGGACCCCGGTAGCCAGGCGTTGTGGAACCAGCGGTTGTTCGATTg                          <  1:656050/62‑1 (MQ=255)                    tcttctTCATTGCGGACCCCGGTAGCCAGGCGTTGTGGAACCAGCTGTTGTTCGATTGCCTg                      <  1:1009289/62‑1 (MQ=255)                          tcATTGCGGACCCCGGTAGCCAGGCGTTGTGGAACCAGCGgt                                    >  1:1048254/1‑42 (MQ=255)                          tcATTGCGGACCCCGGTAGCCAGGCGTTGTGGAACCAGCGgt                                    >  1:4592181/1‑42 (MQ=255)                          tcATTGCGGACCCCGGTAGCCAGGCGTTGTGGAACCAGCGgt                                    >  1:3552660/1‑42 (MQ=255)                          tcATTGCGGACCCCGGTAGCCAGGCGTTGTGGAACCAGCGgt                                    >  1:2594301/1‑42 (MQ=255)                          tcATTGCGGACCCCGGTAGCCAGGCGTTGTGGAACCAGCGgt                                    >  1:160914/1‑42 (MQ=255)                          tcATTGCGGACCCCGGTAGCCAGGCGTTGTGGAACCAGCGgt                                    >  1:1218697/1‑42 (MQ=255)                              tGCGGACCCCGGTAGCCAGGCGTTGTGGAACCAGCGGTTGTTCGATTGCCTGCTCAATGAGt            <  1:4636112/62‑1 (MQ=255)                               gCGGACCCCGGTAGCCAGGCGTTGTGGAACCAGCGGTTGTTCGATTg                          <  1:2008734/47‑1 (MQ=255)                               gCGGACCCCGGTAGCCAGGCGTTGTGGAACCAGCGGTTGTTCGATTg                          <  1:4378886/47‑1 (MQ=255)                               gCGGACCCCGGTAGCCAGGCGTTGTGGAACCAGCGGTTGTTCGATTg                          <  1:485767/47‑1 (MQ=255)                                     cccGGTAGCCAGGCGTTGTGGAACCAGCGGTTGTTCGATTg                          <  1:3182406/41‑1 (MQ=255)                                     cccGGTAGCCAGGCGTTGTGGAACCAGCGGTTGTTCGATTg                          <  1:3839382/41‑1 (MQ=255)                                     cccGGTAGCCAGGCGTTGTGGAACCAGCGGTTGTTCGATTg                          <  1:2187096/41‑1 (MQ=255)                                     cccGGTAGCCAGGCGTTGTGGAACCAGCGGTTGTTCGATTg                          <  1:1938381/41‑1 (MQ=255)                                     cccGGTAGCCAGGCGTTGTGGAACCAGCGGTTGTTCGATTg                          <  1:699419/41‑1 (MQ=255)                                     cccGGTAGCCAGGCGTTGTGGAACCAGCGGTTGTTCGATTg                          <  1:919007/41‑1 (MQ=255)                                        ggTAGCCAGGCGTTGTGGAACCAGCGGTTGTTCGATTGCCTGCTCAATGAGTACACGCCCCt  <  1:3884365/62‑1 (MQ=255)                                        ggTAGCCAGGCGTTGTGGAACCAGCGGTTGTTCGATTGCCTGCTCAATGAGTACACGCCCCt  <  1:1320371/62‑1 (MQ=255)                                                           |                                           GAAATAGATTTCCAGCGCCTCTTCTTCATTGCGGACCCCGGTAGCCAGGCGTTGTGGATCCAGCGGTTGTTCGATTGCCTGCTCAATGAGTACACGCCCCT  >  W3110S.gb/1930906‑1931006

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