Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	196,892	A→G	22.2%	Y116C (TAC→TGC)	yaeL →	zinc metallopeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	196,892	0	A	G	22.2%	84.8 / 22.2	45	Y116C (TAC→TGC)	yaeL	zinc metallopeptidase
Reads supporting (aligned to +/- strand):  ref base A (19/16);  new base G (6/4);  total (25/20)
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.87e-01

CGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGAAATAGCAGCCAA  >  W3110S.gb/196834‑196955                                                           |                                                                cGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGtttt                                                     <  1:3138311/71‑1 (MQ=255)       gcgATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTtattat                                               <  1:2458541/71‑1 (MQ=255)         gATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTttatta                                                >  1:2820004/1‑68 (MQ=255)          attattGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTgg                                            >  1:3823091/1‑71 (MQ=255)           ttattGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGtt                                                       <  1:2361330/59‑1 (MQ=255)            tattGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTGCTGGCTGGTTTTTATTATTg                                             >  1:4646501/1‑68 (MQ=255)               tGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTGCTGGCTGGTTTtta                                                   >  1:1163342/1‑59 (MQ=255)                gCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTGCTGGCTGGTTTTTATTATTGGTGTccc                                      >  1:2898585/1‑71 (MQ=255)                gCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGcc                                      >  1:216811/1‑71 (MQ=255)                   gCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTgg                                   <  1:2924022/71‑1 (MQ=255)                     aGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCg                                 <  1:4124391/71‑1 (MQ=255)                     aGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCg                                 <  1:2009044/71‑1 (MQ=255)                     aGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCg                                 <  1:4307346/71‑1 (MQ=255)                       ggccGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTttat                                                  <  1:3339144/50‑1 (MQ=255)                       gTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCg                                 <  1:3350672/69‑1 (MQ=255)                        tCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTAc                              <  1:2192080/71‑1 (MQ=255)                         ccGGTTGCAAACTTCATTTTTGCTATCTTTGCCTGCTGGCTGGTTTTTATTATTgg                                            >  1:1740837/1‑56 (MQ=255)                         ccGGTTGCAAACTTCATTTTTGCTATCTTTGCCTActg                                                              >  1:4622879/1‑38 (MQ=255)                         ccGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTtta                                                   >  1:1930167/1‑49 (MQ=255)                          cGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTttatta                                                >  1:3885341/1‑51 (MQ=255)                              tGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTttatta                                                >  1:1282268/1‑47 (MQ=255)                              tGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTttat                                                  <  1:3194864/45‑1 (MQ=255)                              tGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTttat                                                  >  1:78781/1‑45 (MQ=255)                               gCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTttat                                                  >  1:2214997/1‑44 (MQ=255)                               gCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTttat                                                  >  1:676226/1‑44 (MQ=255)                               gCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTttat                                                  >  1:1705704/1‑44 (MQ=255)                               gCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCggt                       >  1:3869316/1‑71 (MQ=255)                                cAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCggt                       <  1:4595209/70‑1 (MQ=255)                                  aaCTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTattatt                                              >  1:4597618/1‑45 (MQ=255)                                   aCTTCATTTTTGCTATCTTTTCCTACTGGCTGGTTTttat                                                  <  1:2909784/40‑1 (MQ=255)                                    cTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTg                  >  1:1745830/1‑71 (MQ=255)                                     ttCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTg                  >  1:4413793/1‑70 (MQ=255)                                       cATTTTTGCTATCTTTGCCTGCTGGCTGGTTTTTATTATTg                                             >  1:684549/1‑41 (MQ=255)                                       cATTTTTGCTATCTTTGCCTGCTGGCTGGTTTTTATTATTg                                             >  1:128997/1‑41 (MQ=255)                                       cATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTAc                              >  1:2425459/1‑56 (MQ=255)                                        aTTTTTGCTATCTTTGCCTGCTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGa              <  1:1385435/71‑1 (MQ=255)                                         tttttGCTATCTTTGCCTGCTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGa              <  1:1854941/70‑1 (MQ=255)                                           tttGCTATCTTTGCCTGCTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGaaa            <  1:1029856/70‑1 (MQ=255)                                           tttGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCg                                 >  1:2140735/1‑49 (MQ=255)                                           tttGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGaaa            <  1:41538/70‑1 (MQ=255)                                            ttGGTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCggt                       >  1:2536402/1‑58 (MQ=255)                                               cTATCTTTGCCTGCTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGcta              <  1:4646095/64‑3 (MQ=255)                                                 aTCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCgg                        <  1:3690137/52‑1 (MQ=255)                                                    tttGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGAAATAGCAGCCaa  >  1:3841193/1‑71 (MQ=255)                                                          tACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGa              <  1:2177471/53‑1 (MQ=255)                                                           |                                                                CGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGTTGGCGAAATAGCAGCCAA  >  W3110S.gb/196834‑196955

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑