Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	1,482,053	A→G	50.0%	D401G (GAT→GGT)	ynbC →	predicted hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,482,053	0	A	G	50.0%	‑3.9 / 45.5	32	D401G (GAT→GGT)	ynbC	predicted hydrolase
Reads supporting (aligned to +/- strand):  ref base A (6/10);  new base G (11/5);  total (17/15)
Fisher's exact test for biased strand distribution p-value = 1.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGGCTGGCGCGGTATTCGCCAGCGCAAAACCCATTTACAAATACTGATTAAACAAGCCGTTGCCGATCTCCACGCCAAAGGTTTAGCCGTCCGCGTGGTTGACATTGCCGCAGGGCATGGGCGCTATG  >  W3110S.gb/1481987‑1482115                                                                   |                                                               ttGGCTGGCGCGGTATTCGCCAGCGCAAAACCCATTTACAAATACTGATTAAACAAGCCGTTGCCGAtctc                                                            >  1:1122719/1‑71 (MQ=255) ttGGCTGGCGCGGTATTCGCCAGCGCAAAACCCATTTACAAATACTGATTAAACAAGCCGTTGCCGAtctc                                                            >  1:1506836/1‑71 (MQ=255)                   gCCAGCGCAAAACCCATTTACAAATACTGATTAAACAAGCCGTTGCCGATCTCCACGCCAAAGGTTTAGcc                                          <  1:1302165/71‑1 (MQ=255)                            aaaCCCATTTACAAATACTGATTAAACAAGCCGTTGCCGGTCTCCACGCCAAAGGTTTAGCCGTCCGCGTg                                 <  1:1049588/71‑1 (MQ=255)                            aaaCCCATTTACAAATACTGATTAAACAAGCCGTTGCCGGTCTCCACGCCAAAGGTTTAGCCGTCCGCGTg                                 <  1:1065042/71‑1 (MQ=255)                            aaaCCCATTTACAAATACTGATTAAACAAGCCGTTGCCGGTCTCCACGCCAAAGGTTTAGCCGTCCGCGTg                                 <  1:322642/71‑1 (MQ=255)                            aaaCCCATTTACAAATACTGATTAAACAAGCCGTTGCCGGTCTCCACGCCAAAGGTTTAGCCGTCCGCGTg                                 <  1:27552/71‑1 (MQ=255)                            aaaCCCATTTACAAATACTGATTAAACAAGCCGTTGCCGGTCTCCACGCCAAAGGTTTAGCCGTCCGCGTg                                 <  1:1642291/71‑1 (MQ=255)                                 cATTTACAAATACTGATTAAACAAGCCGTTGCCGATCTCCACGCCAAAGGTTTAGCCGTCCGCGTGGTTGa                            <  1:897566/71‑1 (MQ=255)                                 cATTTACAAATACTGATTAAACAAGCCGTTGCCGATCTCCACGCCAAAGGTTTAGCCGTCCGCGTGGTTGa                            <  1:1018571/71‑1 (MQ=255)                                 cATTTACAAATACTGATTAAACAAGCCGTTGCCGATCTCCACGCCAAAGGTTTAGCCGTCCGCGTGGTTGa                            <  1:595724/71‑1 (MQ=255)                                 cATTTACAAATACTGATTAAACAAGCCGTTGCCGATCTCCACGCCAAAGGTTTAGCCGTCCGCGTGGTTGa                            <  1:45218/71‑1 (MQ=255)                                 cATTTACAAATACTGATTAAACAAGCCGTTGCCGATCTCCACGCCAAAGGTTTAGCCGTCCGCGTGGTTGa                            <  1:199746/71‑1 (MQ=255)                                  aTTTACAAATACTGATTAAACAAGCCGTTGCCGATCTCCACGCCAAAGGTTTAGCCGTCCGCGTGGTTGAc                           <  1:1422906/71‑1 (MQ=255)                                   tttACAAATACTGATTAAACAAGCCGTTGCCGATCTCCACGCCAAAGGTTTAGCCGTCCGCGTGGTTGAc                           <  1:873468/70‑1 (MQ=255)                                     tACAAATACTGATTAAACAAGCCGTTGCCGGTCTCCACGCCAAAGGTTTAGCCg                                         >  1:1960399/1‑54 (MQ=255)                                     tACAAATACTGATTAAACAAGCCGTTGCCGGTCTCCACGCCAAAGGTTTAGCCGTCCGCGTGGTTGACa                          >  1:1491330/1‑69 (MQ=255)                                     tACAAATACTGATTAAACAAGCCGTTGCCGGTCTCCACGCCAAAGGTTTAGCCGTCCGCGTGGTTGACa                          >  1:76692/1‑69 (MQ=255)                                     tACAAATACTGATTAAACAAGCCGTTGCCGGTCTCCACGCCAAAGGTTTAGCCGTCCGCGTGGTTGACa                          >  1:1141036/1‑69 (MQ=255)                                     tACAAATACTGATTAAACAAGCCGTTGCCGGTCTCCACGCCAAAGGTTTAGCCGTCCGCGTGGTTGACa                          >  1:580817/1‑69 (MQ=255)                                     tACAAATACTGATTAAACAAGCCGTTGCCGGTCTCCACGCCAAAGGTTTAGCCGTCCGCGTGGTTGACa                          >  1:529743/1‑69 (MQ=255)                                     tACAAATACTGATTAAACAAGCCGTTGCCGGTCTCCACGCCAAAGGTTTAGCCGTCCGCGTGGTTGACa                          >  1:1233234/1‑69 (MQ=255)                                     tACAAATACTGATTAAACAAGCCGTTGCCGGTCTCCACGCCAAAGGTTTAGCCGTCCGCGTGGTTGACa                          >  1:37786/1‑69 (MQ=255)                                     tACAAATACTGATTAAACAAGCCGTTGCCGGTCTCCACGCCAAAGGTTTAGCCGTCCGCGTGGTTGACa                          >  1:278050/1‑69 (MQ=255)                                     tACAAATACTGATTAAACAAGCCGTTGCCGGTCTCCACGCCAAAGGTTTAGCCGTCCGCGTGGTTGACa                          >  1:1796336/1‑69 (MQ=255)                                     tACAAATACTGATTAAACAAGCCGTTGCCGGTCTCCACGCCAAAGGTTTAGCCGTCCGCGTGGTTGACa                          >  1:1703347/1‑69 (MQ=255)                                           tACTGATTAAACAAGCCGTTGCCGATCTCCCCGCCAAAGGTTTAGCCGTCCGCGTGGTTGACATTGcc                     >  1:1273944/1‑68 (MQ=255)                                           tACTGATTAAACAAGCCGTTGCCGATCTCCACGCCAAAGGTTTAGCCGTCCGCGTGGTTGACATTGCCg                    >  1:1475966/1‑69 (MQ=255)                                           tACTGATTAAACAAGCCGTTGCCGATCTCCACGCCAAAGGTTTAGCCGTCCGCGTGGTTGACATTGCCg                    >  1:729071/1‑69 (MQ=255)                                           tACTGATTAAACAAGCCGTTGCCGATCTCCACGCCAAAGGTTTAGCCGTCCGCGTGGTTGACATTGCCg                    >  1:1615317/1‑69 (MQ=255)                                                           cGTTGCCGATCTCCACGCCAAAGGTTTAGCCGTCCGCGTGGTTGACATTGCCGCAGGGCATGGGCGCTATg  <  1:680929/71‑1 (MQ=255)                                                           cGTTGCCGATCTCCACGCCAAAGGTTTAGCCGTCCGCGTGGTTGACATTGCCGCAGGGCATGGGCGCTATg  <  1:1056750/71‑1 (MQ=255)                                                                   |                                                               TTGGCTGGCGCGGTATTCGCCAGCGCAAAACCCATTTACAAATACTGATTAAACAAGCCGTTGCCGATCTCCACGCCAAAGGTTTAGCCGTCCGCGTGGTTGACATTGCCGCAGGGCATGGGCGCTATG  >  W3110S.gb/1481987‑1482115

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