Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	3,249,944	C→A	8.8%	P89T (CCG→ACG)	yqjG →	predicted S‑transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,249,944	0	C	A	8.8%	211.2 / 3.5	66	P89T (CCG→ACG)	yqjG	predicted S‑transferase
Reads supporting (aligned to +/- strand):  ref base C (33/27);  new base A (2/4);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTA  >  W3110S.gb/3249874‑3250011                                                                       |                                                                    ggCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAAcc                                                                     <  1:1492941/71‑1 (MQ=255)         gcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATg                                                             >  1:4429213/1‑71 (MQ=255)         gcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATg                                                             >  1:242600/1‑71 (MQ=255)         gcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATg                                                             >  1:1883038/1‑71 (MQ=255)          cacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGc                                                            <  1:291724/71‑1 (MQ=255)               tGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGa                                                        >  1:2514356/1‑70 (MQ=255)               tGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGa                                                        >  1:4312017/1‑70 (MQ=255)                 aTCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGc                                                                  >  1:1182615/1‑58 (MQ=255)                  tCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGaaaa                                                     >  1:4537664/1‑70 (MQ=255)                   cATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACg                                                   >  1:3803287/1‑71 (MQ=255)                       cgcAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTg                                               >  1:3492105/1‑71 (MQ=255)                          aaaCTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGAc                                            >  1:2660718/1‑71 (MQ=255)                              tCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGc                                                                  >  1:3942826/1‑45 (MQ=255)                               cAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACACGCTGATGCAGGAAAACGGTTg                                               >  1:4603242/1‑63 (MQ=255)                                aaaGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCttt                                        <  1:395394/69‑1 (MQ=255)                                aaaGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCttt                                        <  1:2111528/69‑1 (MQ=255)                                   ggACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCtt                                         <  1:3372697/65‑1 (MQ=255)                                     aCTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACa                                 <  1:3302563/71‑1 (MQ=255)                                       tGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCt                                                                 <  1:162070/37‑1 (MQ=255)                                       tGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACa                                 <  1:342354/69‑1 (MQ=255)                                       tGGAACCGTTTATTTCCGTTTCCGTAGTGAACACGCTGATGCTGGAAAACGGTTGGACCttt                                        <  1:502068/62‑1 (MQ=255)                                         gAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGa                                             >  1:4246781/1‑55 (MQ=255)                                         gAACCGTTTATTTCCGTTTCCGTAGTGAACACGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGttt                             <  1:912044/71‑1 (MQ=255)                                         gAACCGTTTATTTCCGTTTCCGTAGTGAACACGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGtt                              <  1:546371/70‑1 (MQ=255)                                         gAACCGTCTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGttt                             <  1:820735/71‑1 (MQ=255)                                          aaCCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGtttt                            >  1:4744153/1‑71 (MQ=255)                                            ccGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTcc                          <  1:4094344/71‑1 (MQ=255)                                            ccGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTcc                          <  1:3768825/71‑1 (MQ=255)                                            ccGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTcc                          >  1:3777806/1‑71 (MQ=255)                                              gTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCgg                        >  1:477178/1‑71 (MQ=255)                                               tttATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGaaaa                                                     <  1:3109842/41‑1 (MQ=255)                                               tttATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGaaa                                                      <  1:403251/40‑1 (MQ=255)                                               tttATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTcc                          <  1:5114037/68‑1 (MQ=255)                                               tttATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTcc                          <  1:1599671/68‑1 (MQ=255)                                                 tATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAAc                                                    >  1:728671/1‑40 (MQ=255)                                                 tATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTc                           >  1:255381/1‑65 (MQ=255)                                                 tATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTc                           >  1:1231168/1‑65 (MQ=255)                                                 tATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAg                     >  1:1309618/1‑71 (MQ=255)                                                 tATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAg                     >  1:4709467/1‑71 (MQ=255)                                                 tATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAg                     >  1:2642035/1‑71 (MQ=255)                                                 tATTTCCGTTTCCGTAGTGAACACGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAg                     >  1:2532976/1‑71 (MQ=255)                                                  aTTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACgg                                                  >  1:321463/1‑41 (MQ=255)                                                   tttccgtttccgTAGTGAACCCGCTGATGCTGGAAAACgg                                                  <  1:1617184/40‑1 (MQ=255)                                                   tttccgtttccgTAGTGAACCCGCTGATGCTGGAAAACgg                                                  <  1:1895327/40‑1 (MQ=255)                                                    ttccgtttccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCgg                        >  1:4247731/1‑65 (MQ=255)                                                    ttccgtttccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCaa                  >  1:4528357/1‑71 (MQ=255)                                                     tccgtttccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCaa                  >  1:2073506/1‑70 (MQ=255)                                                     tccgtttccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAAc                 >  1:1250337/1‑71 (MQ=255)                                                      ccgtttccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAAc                 <  1:4147254/70‑1 (MQ=255)                                                       cgtttccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTc                           >  1:3126611/1‑59 (MQ=255)                                                       cgtttccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTc                           >  1:3478669/1‑59 (MQ=255)                                                       cgtttccgtAGTGAACACGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCg               <  1:4874611/71‑1 (MQ=255)                                                         tttccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGga             <  1:1365774/71‑1 (MQ=255)                                                         tttccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGga             <  1:44427/71‑1 (MQ=255)                                                           tccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGgaga           >  1:1573755/1‑71 (MQ=255)                                                             cgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGac          <  1:2894662/70‑1 (MQ=255)                                                              gtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTtgat                                     <  1:4414532/42‑1 (MQ=255)                                                              gtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTtgat                                     <  1:838301/42‑1 (MQ=255)                                                              gtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGTTGACAGTTTTCCGGGAGCAACCgg              <  1:3439085/65‑1 (MQ=255)                                                               tAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTtg                                       <  1:1406682/39‑1 (MQ=255)                                                               tAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACg       >  1:1178518/1‑71 (MQ=255)                                                                 gTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACg       >  1:3250122/1‑69 (MQ=255)                                                                  tGAACCCGCTGATGCTGGAAAACGTTTGGACCTTTGATGACAGTTTTcc                          <  1:1917686/49‑1 (MQ=255)                                                                  tGAACCCGCTGATGCTGGAAAACGGTTGGACCttt                                        <  1:4518273/35‑1 (MQ=255)                                                                  tGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTc                           <  1:3395458/48‑1 (MQ=255)                                                                   gAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGctc    >  1:4748306/1‑70 (MQ=255)                                                                    aaCCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTa  >  1:4911798/1‑71 (MQ=255)                                                                       |                                                                    GGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTA  >  W3110S.gb/3249874‑3250011

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