Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	3,249,944	C→A	21.6%	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	21.6%	100.5 / 11.9	37	P89T (CCG→ACG)	yqjG	predicted S‑transferase
Reads supporting (aligned to +/- strand):  ref base C (15/14);  new base A (3/5);  total (18/19)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

CACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGC  >  W3110S.gb/3249878‑3250007                                                                   |                                                                caccgcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCt                                                             <  1:1673136/71‑1 (MQ=255)      cacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGc                                                        >  1:544219/1‑71 (MQ=255)             aTCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGc                                                              >  1:752324/1‑58 (MQ=255)             aTCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGaaaa                                                 <  1:960718/71‑1 (MQ=255)               cATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACACGCTGATGCTGGAAAACg                                               >  1:507306/1‑71 (MQ=255)                aTGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACgg                                              <  1:284118/71‑1 (MQ=255)                aTGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACgg                                              >  1:1346957/1‑71 (MQ=255)                aTGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACg                                               <  1:1494759/70‑1 (MQ=255)                aTGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACACGCTGATGCTGGAAAACgg                                              >  1:1098645/1‑71 (MQ=255)                 tGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACACGCTGATGCTGGAAAACg                                               <  1:1711317/69‑1 (MQ=255)                    gcAAACTCAAAGGACTGGAACCGTTTATTTACGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTgg                                          <  1:1538124/71‑1 (MQ=255)                        aCTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCt                                      >  1:689749/1‑71 (MQ=255)                         cTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACACGCTGATGCTGGAAAACGGTTGGACCtt                                     <  1:1029031/71‑1 (MQ=255)                          tCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAAccc                                                                <  1:1054016/43‑1 (MQ=255)                          tCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCt                                                       >  1:260119/1‑52 (MQ=255)                            aaaGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTtg                                   <  1:124919/70‑1 (MQ=255)                               ggACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTtgatga                               <  1:95640/71‑1 (MQ=255)                               ggACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTtgatga                               <  1:12865/71‑1 (MQ=255)                                   tGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCt                                                       <  1:416113/43‑1 (MQ=255)                                   tGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGt                           >  1:770101/1‑71 (MQ=255)                                          gTTTATTTCCGTTTCCGTAGTGAACACGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCgg                    <  1:1646759/71‑1 (MQ=255)                                           tttATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCgg                    <  1:1699679/70‑1 (MQ=255)                                            ttATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGAc                              <  1:963206/59‑1 (MQ=255)                                            ttATTTCCGTTTCCGTAGTGAACACGCTGATGCTGGAAAACGGTTGGACCttt                                    >  1:204517/1‑53 (MQ=255)                                             tATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACgg                                              <  1:135290/42‑1 (MQ=255)                                             tATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAg                 >  1:397298/1‑71 (MQ=255)                                             tATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAg                 >  1:17655/1‑71 (MQ=255)                                             tATTTCCGTTTCCGTAGTGAACCCGCTGATGCCGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAg                 >  1:1703621/1‑71 (MQ=255)                                              aTTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGtt                                            >  1:260531/1‑43 (MQ=255)                                               tttccgtttccgTAGTGAACCCGCTGATGCTGGAAAACGGtt                                            >  1:1401541/1‑42 (MQ=255)                                               tttccgtttccgTAGTGAACACGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCa               <  1:556149/71‑1 (MQ=255)                                                          gtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTc                       >  1:696485/1‑52 (MQ=255)                                                          gtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGacac    >  1:1113415/1‑71 (MQ=255)                                                           tAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTc                       <  1:1727020/51‑1 (MQ=255)                                                           tAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACg   <  1:1094678/71‑1 (MQ=255)                                                            aGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAAc             >  1:1467346/1‑60 (MQ=255)                                                             gTGAACACGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGc  <  1:953855/70‑1 (MQ=255)                                                                aaCCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCg                     >  1:869299/1‑48 (MQ=255)                                                                   |                                                                CACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGC  >  W3110S.gb/3249878‑3250007

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