Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	3,960,751	C→T	51.1%	R83H (CGT→CAT)	kdgK ←	ketodeoxygluconokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,960,751	0	C	T	51.1%	‑5.5 / 66.5	45	R83H (CGT→CAT)	kdgK	ketodeoxygluconokinase
Reads supporting (aligned to +/- strand):  ref base C (17/5);  new base T (17/6);  total (34/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCCAGTAGTAGAACGTACGCTCGCCGGTGCTGTCGGTTTCAATGTAGTAAAGGCCCGGCAGACGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACGTTCTCGCCGTGCCAGGCGTCCAGCAT  >  W3110S.gb/3960687‑3960812                                                                 |                                                              ccGCCAGTAGTAGAACGTACGCTCGCCGGTGCTGTCGGTTTCAATGTAGTAAAGGCCCGGCAGACGGtttt                                                         >  1:1468436/1‑71 (MQ=255)            agAACGTACGCTCGCCGGTGCTGTCGGTTTCAATGTAGTAAAGGCCCGGCAGACGGTTTTCCATCCGTTgg                                              >  1:3619842/1‑71 (MQ=255)             gAACGTACGCTCGCCGGTGCTGTCGGTTTCAATGTAGTAAAGGCCCGGCAGACGGTTTTCCATCCGTTggg                                             >  1:2375730/1‑71 (MQ=255)             gAACGTACGCTCGCCGGTGCTGTCGGTTTCAATGTAGTAAAGGCCCGGCAGACGGTTTTCCATCCGTTgg                                              >  1:1842361/1‑70 (MQ=255)               aCGTACGCTCGCCGGTGCTGTCGGTTTCAATGTAGTAAAGGCCCGGCAGATGGTTTTCCATCCGTTGGGTc                                           <  1:2354520/71‑1 (MQ=255)                cGTACGCTCGCCGGTGCTGTCGGTTTCAATGTAGTAAAGGCCCGGCAGATGGTTTTCAATCCGTTGGGTCa                                          <  1:2973574/71‑1 (MQ=255)                  tACGCTCGCCGGTGCTGTCGGTTTCAATGTAGTAAAGGCCCGGCAGATGGTTTTCCATCCGTTggg                                             >  1:3992360/1‑66 (MQ=255)                  tACGCTCGCCGGTGCTGTCGGTTTCAATGTAGTAAAGGCCCGGCAGATGGTTTTCCATCCGTTGGGTCAgg                                        >  1:112272/1‑71 (MQ=255)                  tACGCTCGCCGGTGCTGTCGGTTTCAATGTAGTAAAGGCCCGGCAGATGGTTTTCCATCCGTTGGGTCAgg                                        >  1:2520209/1‑71 (MQ=255)                  tACGCTCGCCGGTGCTGTCGGTTTCAATGTAGTAAAGGCCCGGCAGATGGTTTTCCATCCGTTGGGTCAgg                                        >  1:1942236/1‑71 (MQ=255)                    cGCTCGCCTGTGCTGTCGGTTTCAATGTAGTAAAGGCCCGGCAGATGGTTTTCCATCCGTTGGGTCAGGGa                                      <  1:1278158/71‑1 (MQ=255)                    cGCTCGCCGGTGCTGTCGGTTTCAATGTAGTAAAGGCCCGGCAGATGGTTTTCCATCCGTTGGGTCAgg                                        >  1:67973/1‑69 (MQ=255)                    cGCTCGCCGGTGCTGTCGGTTTCAATGTAGTAAAGGCCCGGCAGATGGTTTTCCATCCGTTGGGTCAgg                                        >  1:1932058/1‑69 (MQ=255)                    cGCTCGCCGGTGCTGTCGGTTTCAATGTAGTAAAGGCCCGGCAGACGGTTTTCCATCCGTTGGGTCAGGGa                                      <  1:319106/71‑1 (MQ=255)                            ggTGCTGTCGGTTTCAATGTAGTAAAGGCCCGGCAGATGGTTTTCCATCCGTTGGGTCAGGGAAGTATc                                <  1:425251/69‑1 (MQ=255)                                cTGTCGGTTTCAATGTAGTAAAGGCCCGGCAGATGGtt                                                           >  1:95890/1‑38 (MQ=255)                                cTGTCGGTTTCAATGTAGTAAAGGCCCGGCAGATGGtt                                                           >  1:1209707/1‑38 (MQ=255)                                   tCGGTTTCAATGCAGTAAAGGCCCGGCAGATGGTTCTCCATCCGTTGGGTCAGGGAAGTATc                                <  1:3866325/62‑1 (MQ=255)                                    cGGTTTCAATGTAGTAAAGGCCCGGCAGATGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACGtt                          >  1:1604570/1‑67 (MQ=255)                                    cGGTTTCAATGTAGTAAAGGCCCGGCAGATGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACGTTCTCg                      >  1:1912017/1‑71 (MQ=255)                                     ggTTTCAATGTAGTAAAGGCCCGGCAGACGGTTTTCCATCCGTTGGGTCAGGGaa                                     >  1:422477/1‑55 (MQ=255)                                     ggTTTCAATGTAGTAAAGGCCCGGCAGACGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACGTTCTCGc                     >  1:3812284/1‑71 (MQ=255)                                       tttCAATGTAGTAAAGGCCCGGCAGACGGTTTTCCATCCGTTGGGTCAGGGa                                      <  1:3666085/52‑1 (MQ=255)                                        ttCAATGTAGTAAAGGCCCGGCAGACGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACGTTCTCGCCGt                  >  1:3006253/1‑71 (MQ=255)                                          cAATGTAGTAAAGGCCCGGCAGACGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACGTtct                        >  1:1494774/1‑63 (MQ=255)                                              gtagtaAAGGCCCGGCAGATGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACGTTCTCGCCGTGCCAgg            >  1:2464846/1‑71 (MQ=255)                                              gtagtaAAGGCCCGGCAGATGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACGTTCTCGCCGTGCCAgg            >  1:2220030/1‑71 (MQ=255)                                              gtagtaAAGGCCCGGCAGATGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACGTTCTCGCCGTGCCAgg            >  1:1939551/1‑71 (MQ=255)                                              gtagtaAAGGCCCGGCAGATGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACGTTCTCGCCGTGCCAg             >  1:3659675/1‑70 (MQ=255)                                               tagtaAAGGCCCGGCAGACGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACGTTCTCGCCGTGCCAGGc           >  1:3692145/1‑71 (MQ=255)                                                    aaGGCCCGGCAGATGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACGTTCTCGCCGTGCCAGGCGTc        >  1:1384819/1‑69 (MQ=255)                                                    aaGGCCCGGCAGATGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACGTTCTCGCCGTGCCAGGCGTc        >  1:3709738/1‑69 (MQ=255)                                                    aaGGCCCGGCAGATGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACCTtctc                       >  1:1695777/1‑54 (MQ=255)                                                    aaGGCCCGGCAGACGGTTTTCCATCCGTTGGGTCa                                          >  1:2552520/1‑35 (MQ=255)                                                    aaGGCCCGGCAGACGGTTTTCCATCCGTTGGGTCAGGGAAg                                    >  1:1006840/1‑41 (MQ=255)                                                    aaGGCCCGGCAGACGGTTTTCCATCCGTTGGGTCAGGGAAGTATCaa                              >  1:1592138/1‑47 (MQ=255)                                                     aGGCCCGGCAGATGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACGt                           <  1:2716258/49‑1 (MQ=255)                                                      ggCCCGGCAGACGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACGTTCTCGCCGTGCCAGGCGTCCAGc    <  1:3552639/71‑1 (MQ=255)                                                       gCCCGGCAGACGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACGTTCTCGCCGTGCCAGGCGTCCAGCa   >  1:3569791/1‑71 (MQ=255)                                                       gCCCGGCAGACGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACGTTCTCGCCGTGCCAGGCGTCCAGCa   >  1:182703/1‑71 (MQ=255)                                                             cAGACGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACGTTCTCGCCGTGCCAGGCGTCCAGc    <  1:3706377/64‑1 (MQ=255)                                                             cAGACGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACGTTCTCGCCGTGCCAGGCGTCCAGc    <  1:2906163/64‑1 (MQ=255)                                                               gACGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACGTTCTCGCCGTGCCAGGCGTCCAGCAt  >  1:1667692/1‑64 (MQ=255)                                                               gACGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACGTTCTCGCCGTGCCAGGCGTCCAGCAt  >  1:1352530/1‑64 (MQ=255)                                                               gACGGTTTTACATCCGTTGGGTCAGGGAAGTATCAACGTTCTCGCCGTGCCAGGCGTCCAGCAt  >  1:2007870/1‑64 (MQ=255)                                                                 |                                                              CCGCCAGTAGTAGAACGTACGCTCGCCGGTGCTGTCGGTTTCAATGTAGTAAAGGCCCGGCAGACGGTTTTCCATCCGTTGGGTCAGGGAAGTATCAACGTTCTCGCCGTGCCAGGCGTCCAGCAT  >  W3110S.gb/3960687‑3960812

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