Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	3,249,944	C→A	11.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	11.8%	149.1 / 5.4	51	P89T (CCG→ACG)	yqjG	predicted S‑transferase
Reads supporting (aligned to +/- strand):  ref base C (19/25);  new base A (2/4);  total (22/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GGCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTA  >  W3110S.gb/3249874‑3250011                                                                       |                                                                    ggCGCACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAAcc                                                                     >  1:3354144/1‑71 (MQ=255)     caccgcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACGCGCt                                                                 >  1:1231255/1‑71 (MQ=255)     caccgcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGc                                                                  <  1:2035080/70‑1 (MQ=255)     caccgcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCt                                                                 <  1:2676425/71‑1 (MQ=255)     caccgcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCt                                                                 <  1:236287/71‑1 (MQ=255)       ccgcacGGTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACACGCTGa                                                               >  1:2265258/1‑71 (MQ=255)          cacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGc                                                            >  1:3089378/1‑71 (MQ=255)          cacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCAGTTTCCGTAGTGAACCCGCt                                                                 >  1:3719698/1‑66 (MQ=255)                 aTCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTgg                                                         <  1:529719/67‑1 (MQ=255)                   cATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGaa                                                       <  1:1203581/67‑1 (MQ=255)                   cATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACg                                                   >  1:2024866/1‑71 (MQ=255)                         cAAGCTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTgg                                              <  1:2119280/70‑1 (MQ=255)                         cAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAAccc                                                                    >  1:3065305/1‑48 (MQ=255)                          aaaCTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGAc                                            >  1:2886020/1‑71 (MQ=255)                             cTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCtt                                         >  1:1357843/1‑71 (MQ=255)                             cTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACACGCTGATGCTGGAAAACGGTTGGACCtt                                         <  1:3487917/71‑1 (MQ=255)                              tCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCg                                                                   <  1:3915435/44‑1 (MQ=255)                              tCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCttt                                        >  1:4300244/1‑71 (MQ=255)                                aaaGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGAt                                                              <  1:3797771/47‑1 (MQ=255)                                aaaGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTg                                               >  1:2469063/1‑62 (MQ=255)                                aaaGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCttt                                        <  1:3067719/69‑1 (MQ=255)                                aaaGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTtg                                       <  1:1440113/70‑1 (MQ=255)                                aaaGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACACGCTGATGCTGGAAAACGGTTGGACCttt                                        <  1:515618/69‑1 (MQ=255)                                aaaGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACACGCTGATGCTGGAAAACGGTTGGACCTTtga                                      >  1:3238400/1‑71 (MQ=255)                                  aGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTtgatg                                    <  1:1413699/71‑1 (MQ=255)                                  aGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTtgat                                     >  1:3560982/1‑70 (MQ=255)                                     aCTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACa                                 <  1:4396861/71‑1 (MQ=255)                                       tGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTtg                                       <  1:1318913/63‑1 (MQ=255)                                         gAACCGTTTATTTCCGTTTCTGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCttt                                        >  1:1933640/1‑60 (MQ=255)                                            ccGTTTTTTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTcc                          <  1:894834/71‑1 (MQ=255)                                             cGTTTTTTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTcc                          <  1:900257/70‑1 (MQ=255)                                             cGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGttt                             <  1:4550832/67‑1 (MQ=255)                                             cGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCg                         >  1:2873528/1‑71 (MQ=255)                                             cGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCg                         <  1:3670068/71‑1 (MQ=255)                                             cGTTTATTTCCGTTTCCGTAGTGAACACGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCg                         <  1:735900/71‑1 (MQ=255)                                               tttATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGt                                                 <  1:116551/45‑1 (MQ=255)                                                ttATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCgg                        >  1:412962/1‑69 (MQ=255)                                                ttATTTCCGTTTCCGTAGTGAACCAGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCgg                        >  1:3464860/1‑69 (MQ=255)                                                 tATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCtt                                         >  1:877584/1‑51 (MQ=255)                                                  aTTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAg                     <  1:3167193/70‑1 (MQ=255)                                                   tttccgtttccgTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGc                    <  1:2936511/70‑1 (MQ=255)                                                    ttccgtttccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCaa                  >  1:2700231/1‑71 (MQ=255)                                                        gtttccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCg               <  1:745485/70‑1 (MQ=255)                                                           tccgtAGTGAACCTGCTGATGCTGGAAAACGGTTGGACCt                                          >  1:4362965/1‑40 (MQ=255)                                                           tccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGgaga           >  1:177262/1‑71 (MQ=255)                                                            ccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTtg                                       <  1:160093/42‑1 (MQ=255)                                                            ccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGgaga           <  1:1241653/70‑1 (MQ=255)                                                               tAGTGAACACGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCgg                        <  1:2802282/54‑1 (MQ=255)                                                                  tGAACCCGCTGATGCTGGAAAACGGTT‑GACCTTTGATGACAGTTTTCCGGGAGCAACCGga             <  1:1863103/61‑1 (MQ=255)                                                                   gAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGt                               >  1:3194784/1‑43 (MQ=255)                                                                   gAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGga             <  1:4424615/61‑1 (MQ=255)                                                                    aaCCCGCTGATGCTGGAAAACGGTTGGACCTTtga                                      >  1:2700075/1‑35 (MQ=255)                                                                    aaCCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTa  >  1:4420163/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: ‑