Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	158,626	G→A	50.0%	S167S (AGC→AGT)	pcnB ←	poly(A) polymerase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	158,626	0	G	A	50.0%	‑3.5 / 44.0	32	S167S (AGC→AGT)	pcnB	poly(A) polymerase I
Reads supporting (aligned to +/- strand):  ref base G (7/9);  new base A (9/7);  total (16/16)
Fisher's exact test for biased strand distribution p-value = 7.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAA  >  W3110S.gb/158571‑158694                                                        |                                                                     ccTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGt                                                           <  1:1274666/67‑1 (MQ=255)    tCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAAt                                                      <  1:2297820/69‑1 (MQ=255)     cATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCt                                                                    <  1:2864629/54‑1 (MQ=255)           gccAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCgcggc                                                <  1:2290663/68‑1 (MQ=255)              aaCGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGTGAAATCGCGGCGCTg                                            >  1:2734631/1‑69 (MQ=255)                    aTCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATcgc                                                   >  1:2880637/1‑56 (MQ=255)                    aTCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTg                                            >  1:1206046/1‑63 (MQ=255)                    aTCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTg                                            >  1:419702/1‑63 (MQ=255)                    aTCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGTGAAATCgcgg                                                 >  1:2466653/1‑58 (MQ=255)                       acggacggTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGgcgc                                              >  1:1019848/1‑58 (MQ=255)                       acggacggTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGtcttc                                   >  1:1904963/1‑69 (MQ=255)                           acggTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCttcttc                                <  1:1858123/68‑1 (MQ=255)                            cggTAAAATCTGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGa                              <  1:1230928/69‑1 (MQ=255)                            cggTAAAATCTGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGa                              <  1:2821700/69‑1 (MQ=255)                              gTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGc                                         <  1:501410/56‑1 (MQ=255)                               tAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATgg                           <  1:524868/69‑1 (MQ=255)                                 aaaTCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAg                         >  1:1917375/1‑69 (MQ=255)                                  aaTCCGCTACTCTGTAATACAGGCTGTTGATAGTGAAATCGCGgcgc                                              <  1:3342017/47‑1 (MQ=255)                                    tCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGcc                       >  1:389071/1‑68 (MQ=255)                                              tGTAATACAGACTGTTGATAGTGAAATCGCGGCCCTGGGCGTCTTCTTCGATGGAGCCGAAAAtgttgt            >  1:1455484/1‑69 (MQ=255)                                                  aTACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCg                                        <  1:2041429/37‑1 (MQ=255)                                                  aTACAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTcg          >  1:3558087/1‑67 (MQ=255)                                                  aTACAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTcg          >  1:1678198/1‑67 (MQ=255)                                                  aTACAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTcg          >  1:1682503/1‑67 (MQ=255)                                                  aTACAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTcg          >  1:2203434/1‑67 (MQ=255)                                                  aTACAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTcg          >  1:968092/1‑67 (MQ=255)                                                  aTACAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTATTCGATGGAGCCGAAAAtg                >  1:206612/1‑61 (MQ=255)                                                     cAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAAt                 <  1:180702/57‑1 (MQ=255)                                                     cAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAAt                 <  1:3111712/57‑1 (MQ=255)                                                     cAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAAt                 <  1:1678789/57‑1 (MQ=255)                                                     cAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAAt                 <  1:185171/57‑1 (MQ=255)                                                     cAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAAt                 <  1:1696906/57‑1 (MQ=255)                                                       ggCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCACgcagca   <  1:3177561/69‑1 (MQ=255)                                                        gCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCaa  <  1:984007/69‑1 (MQ=255)                                                        |                                                                     CCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAA  >  W3110S.gb/158571‑158694

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