Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	3,179,239	C→T	73.7%	L55L (CTG→TTG)	ygiC →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,179,239	0	C	T	73.7%	59.2 / 24.3	38	L55L (CTG→TTG)	ygiC	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (4/6);  new base T (13/15);  total (17/21)
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.81e-01

AGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCAT  >  W3110S.gb/3179178‑3179309                                                              |                                                                       aGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGa                                                               <  1:415100/71‑1 (MQ=255)  gATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAATTGCACCAGAt                                                              >  1:1333975/1‑71 (MQ=255)   aTGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAATTGCACCAGAtg                                                             >  1:1764381/1‑71 (MQ=255)               aGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAATTGCACCAGATGTGCCTGa                                                      >  1:2708722/1‑66 (MQ=255)                 ttGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAATTGCACCAGa                                                               <  1:697360/55‑1 (MQ=255)                 ttGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAATTGCACCAGa                                                               <  1:1347863/55‑1 (MQ=255)                      ccTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCa                                                                 <  1:230516/48‑1 (MQ=255)                          gCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGt                                      <  1:397737/71‑1 (MQ=255)                           cccAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAATTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTg                                     <  1:1731778/71‑1 (MQ=255)                           cccAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAATTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTg                                     <  1:2204043/71‑1 (MQ=255)                            ccAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAATTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTg                                     <  1:1976370/70‑1 (MQ=255)                             cAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTg                                     <  1:2525293/69‑1 (MQ=255)                              aGGTTGAAAAGCTGGAAGAAGTCACCGCCGAATTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATc                                  >  1:1327701/1‑71 (MQ=255)                               ggTTGAAAAGCTGGAAGAAGTCACCGCCGAATTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCg                                 >  1:2133806/1‑71 (MQ=255)                                  tGAAAAGCTGGAAGAAGTCACCGCCGAATTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGc                                >  1:2511696/1‑69 (MQ=255)                                  tGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGc                                >  1:2647059/1‑69 (MQ=255)                                   gAAAAGCTGGAAGAAGTCACCGCCGAATTGCACCAGa                                                               >  1:703676/1‑37 (MQ=255)                                   gAAAAGCTGGAAGAAGTCACCGCCGAATTGCACCAGa                                                               >  1:706774/1‑37 (MQ=255)                                   gAAAAGCTGGAAGAAGTCACCGCCGAATTGCACCAGATGTGCCTGaaa                                                    >  1:412398/1‑48 (MQ=255)                                   gAAAAGCTGGAAGAAGTCACCGCCGAATTGCACCAGATGTGCCTGaaa                                                    >  1:2591255/1‑48 (MQ=255)                                     aaaGCTGGAAGAAGTCACCGCCGAATTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGc                            <  1:872555/70‑1 (MQ=255)                                     aaaGCTGGAAGAAGTCACCGCCGAATTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGc                            <  1:1800216/70‑1 (MQ=255)                                     aaaGCTGGAAGAAGTCACCGCCGAATTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGc                            <  1:52930/70‑1 (MQ=255)                                            gaagaaGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCt                    <  1:1038207/71‑1 (MQ=255)                                               gaagTCACCGCCGAATTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCt                    <  1:792468/68‑1 (MQ=255)                                                   tCACCGCCGAATTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtg                   <  1:1631391/65‑1 (MQ=255)                                                     aCCGCCGAATTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGa                          >  1:1154645/1‑56 (MQ=255)                                                     aCCGCCGAATTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtgat                 >  1:2016624/1‑65 (MQ=255)                                                     aCCGCCGAATTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtgat                 >  1:1592808/1‑65 (MQ=255)                                                     aCCGCCGAATTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCa            <  1:655952/70‑1 (MQ=255)                                                     aCCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCa            >  1:930985/1‑70 (MQ=255)                                                     aCCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCa            >  1:1055553/1‑70 (MQ=255)                                                     aCCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCa            >  1:2400245/1‑70 (MQ=255)                                                      ccgccgAATTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCaaa          <  1:2723313/71‑1 (MQ=255)                                                           gAATTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTTATCGCCAGCGATGAGCTGATGACCAAATTCCg     <  1:2389229/71‑1 (MQ=255)                                                           gAATTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCg     <  1:800903/71‑1 (MQ=255)                                                           gAATTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCg     <  1:90512/71‑1 (MQ=255)                                                             aCTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCa   <  1:1257907/71‑1 (MQ=255)                                                              cTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCAt  >  1:2490981/1‑71 (MQ=255)                                                              |                                                                       AGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCAT  >  W3110S.gb/3179178‑3179309

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