Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	1,006,848	A→T	K159* (AAA→TAA)	ycbW →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,006,848	0	A	T	100.0%	171.6 / NA	57	K159* (AAA→TAA)	ycbW	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (30/27);  total (30/27)

CTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTCT  >  W3110S.gb/1006782‑1006914                                                                   |                                                                   cTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAAt                                                                >  1:482585/1‑71 (MQ=255) cTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAAt                                                                >  1:3141428/1‑71 (MQ=255)   aTGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATc                                                               >  1:4454633/1‑70 (MQ=255)   aTGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCa                                                              >  1:813251/1‑71 (MQ=255)   aTGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCa                                                              >  1:4333943/1‑71 (MQ=255)   aTGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCa                                                              >  1:3801721/1‑71 (MQ=255)   aTGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCa                                                              >  1:2203918/1‑71 (MQ=255)    tGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATc                                                               <  1:698934/69‑1 (MQ=255)    tGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATc                                                               <  1:3073639/69‑1 (MQ=255)    tGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCAt                                                             >  1:4175812/1‑71 (MQ=255)    tGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCAt                                                             >  1:2818568/1‑71 (MQ=255)    tGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCAt                                                             >  1:3349742/1‑71 (MQ=255)    tGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCAt                                                             >  1:2930273/1‑71 (MQ=255)      ccTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCATg                                                            >  1:3794605/1‑70 (MQ=255)      ccTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCATg                                                            <  1:2673917/70‑1 (MQ=255)      ccTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCATg                                                            >  1:1963405/1‑70 (MQ=255)          ctgGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCATGAACGa                                                       <  1:887678/71‑1 (MQ=255)          ctgGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCATGAACGa                                                       >  1:2785010/1‑71 (MQ=255)          ctgGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCATGAACGa                                                       >  1:2760937/1‑71 (MQ=255)          ctgGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCATGAACGa                                                       <  1:148968/71‑1 (MQ=255)                caACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCATGAACGACAGGCt                                                 >  1:3717526/1‑71 (MQ=255)                caACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCATGAACGACAGGCt                                                 >  1:1361649/1‑71 (MQ=255)                  aCCCTGCGTTGTTATTGCGGTTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCATGAACGACAGGCTGa                                               <  1:812277/71‑1 (MQ=255)                  aCCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCATGAACGa                                                       >  1:3841113/1‑63 (MQ=255)                  aCCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCATGAACGACAGGCTg                                                >  1:4286149/1‑70 (MQ=255)                  aCCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCATGAACGACAGGCTg                                                >  1:1701809/1‑70 (MQ=255)                  aCCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCATGAACGACAGGCTg                                                >  1:198664/1‑70 (MQ=255)                  aCCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCATGAACGACAGGCTGa                                               <  1:2174638/71‑1 (MQ=255)                  aCCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTGGGGCGACGCCATTTAAATCATGAACGACAGGCTGa                                               <  1:2269652/71‑1 (MQ=255)                      tGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCa                                                              >  1:2639593/1‑52 (MQ=255)                       gCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCATGAACGACAGGCTGAAACCg                                          >  1:2900355/1‑71 (MQ=255)                        cGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCATGAACGACAGGCTGAAACCg                                          <  1:2446901/70‑1 (MQ=255)                            gttATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCATGAACGACAGGCTGAAACCGCAGGt                                     >  1:1677043/1‑71 (MQ=255)                             ttATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCATGAACGACAGGCTGAAACCGCAGGtt                                    <  1:2020985/71‑1 (MQ=255)                             ttATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCATGAACGACAGGCTGAAACCGCAGGtt                                    <  1:833495/71‑1 (MQ=255)                             ttATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCATGAACGACAGGCTGAAACCGCAGGtt                                    <  1:1629535/71‑1 (MQ=255)                             ttATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCATGAACGACAGGCTGAAACCGCAGGtt                                    <  1:4388782/71‑1 (MQ=255)                             ttATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCATGAACGACAGGCTGAAACCGCAGGtt                                    <  1:3469989/71‑1 (MQ=255)                             ttATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCATGAACGACAGGCTGAAACCGCAGGtt                                    <  1:2779549/71‑1 (MQ=255)                                ttGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCATGAACGACAGGCTGAAACCGCAGGTTAAt                                 <  1:644963/71‑1 (MQ=255)                                   cGGGTCGCGCGATGCAGTTGGGCGAGGCCATTTAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGtt                              <  1:3174142/71‑1 (MQ=255)                                   cGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGtt                              <  1:3314463/71‑1 (MQ=255)                                   cGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGtt                              <  1:1030322/71‑1 (MQ=255)                                    gggTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGtt                              >  1:2914783/1‑70 (MQ=255)                                                gCAGTTGGGCGACGCCATTTAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGcc                 >  1:1397208/1‑71 (MQ=255)                                                   gTTGGGCGACGCCATTTAAATCATGAACGACAGGCTGAAACCGCAGGtt                                    <  1:2977755/49‑1 (MQ=255)                                                   gTTGGGCGACGCCATTTAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCg              <  1:58151/71‑1 (MQ=255)                                                     tGGGCGACGCCATTTAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGaa            >  1:3075809/1‑71 (MQ=255)                                                      gggCGACGCCATTTAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGaa            <  1:3141146/70‑1 (MQ=255)                                                       ggCGACGCCATTTAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGtt                              <  1:3502413/51‑1 (MQ=255)                                                          gACGCCATTTAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAgg        <  1:4298114/70‑1 (MQ=255)                                                             gCCATTTAAATCATGAACGACAGGCTGAAACCG‑AGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTc   >  1:3203063/1‑71 (MQ=255)                                                               cATTTAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTCt  <  1:4253321/71‑1 (MQ=255)                                                               cATTTAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTCt  <  1:2127401/71‑1 (MQ=255)                                                               cATTTAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCAAACAGGCAGTCt  <  1:3376675/71‑1 (MQ=255)                                                                aTTTAAATCATGAACGACAGGCTGAAACCGCAGGTTaa                                  >  1:3418395/1‑38 (MQ=255)                                                                aTTTAAATCATGAACGACAGGCTGAAACCGCAGGTTaa                                  >  1:3786776/1‑38 (MQ=255)                                                                   |                                                                   CTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTCT  >  W3110S.gb/1006782‑1006914

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