Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	1,210,013	T→C	34.2%	L102P (CTG→CCG)	ymfS →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,210,013	0	T	C	34.2%	38.7 / 34.9	38	L102P (CTG→CCG)	ymfS	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (10/15);  new base C (9/4);  total (19/19)
Fisher's exact test for biased strand distribution p-value = 1.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATC  >  W3110S.gb/1209949‑1210070                                                                 |                                                          gATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCCGaaa                                                       >  1:354501/1‑69 (MQ=255) gATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCCGAAAc                                                      >  1:292422/1‑70 (MQ=255) gATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCCGAAAc                                                      >  1:972206/1‑70 (MQ=255) gATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCCGAAAc                                                      >  1:971777/1‑70 (MQ=255) gATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCCGAAAc                                                      >  1:250353/1‑70 (MQ=255) gATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCCGAAAc                                                      >  1:1316490/1‑70 (MQ=255) gATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCCGAAAc                                                      >  1:1865741/1‑70 (MQ=255) gATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCCGAAAc                                                      >  1:1705765/1‑70 (MQ=255)    tttATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCCGAAACAata                                                  >  1:48414/1‑71 (MQ=255)              gCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTgg                                        >  1:1701328/1‑71 (MQ=255)              gCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTgg                                        >  1:1517397/1‑71 (MQ=255)              gCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTgg                                        >  1:1235032/1‑71 (MQ=255)                              tAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAg                        <  1:661269/71‑1 (MQ=255)                              tAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAg                        <  1:477678/71‑1 (MQ=255)                              tAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAg                        <  1:1453532/71‑1 (MQ=255)                              tAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAg                        <  1:203692/71‑1 (MQ=255)                                aGGCTGCATT‑GGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCAc                     <  1:422503/71‑1 (MQ=255)                                 ggCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCAc                     <  1:1123623/71‑1 (MQ=255)                                 ggCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCAc                     <  1:787574/71‑1 (MQ=255)                                 ggCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCAc                     <  1:484399/71‑1 (MQ=255)                                 ggCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCAc                     <  1:884005/71‑1 (MQ=255)                                 ggCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCAc                     <  1:985541/71‑1 (MQ=255)                                 ggCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCAc                     <  1:1718830/71‑1 (MQ=255)                                 ggCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCAc                     <  1:1529616/71‑1 (MQ=255)                                 ggCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCAc                     <  1:1527315/71‑1 (MQ=255)                                 ggCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCAc                     <  1:1202425/71‑1 (MQ=255)                                 ggCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCAc                     <  1:1049218/71‑1 (MQ=255)                                                  ttAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACa    >  1:1884859/1‑71 (MQ=255)                                                  ttAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACa    >  1:1883900/1‑71 (MQ=255)                                                  ttAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACa    >  1:577471/1‑71 (MQ=255)                                                  ttAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACa    >  1:729550/1‑71 (MQ=255)                                                  ttAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACa    >  1:960611/1‑71 (MQ=255)                                                  ttAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACa    >  1:996092/1‑71 (MQ=255)                                                     aagaagATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATc  <  1:257370/70‑1 (MQ=255)                                                     aagaagATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATc  <  1:1859440/70‑1 (MQ=255)                                                     aagaagATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATc  <  1:168426/70‑1 (MQ=255)                                                     aagaagATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATc  <  1:1204958/70‑1 (MQ=255)                                                             gAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACa    >  1:1963896/1‑60 (MQ=255)                                                                 |                                                          GATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATC  >  W3110S.gb/1209949‑1210070

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