Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	209,475	2 bp→TC	coding (855‑856/960 nt)	accA →	acetylCoA carboxylase, carboxytransferase, alpha subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	209,475	0	G	T	76.3%	54.4 / 19.0	38	S285S (TCG→TCT)	accA	acetylCoA carboxylase, carboxytransferase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (0/9);  new base T (29/0);  total (29/9)
Fisher's exact test for biased strand distribution p-value = 6.13e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01
*	W3110S.gb	209,476	0	T	C	76.3%	52.2 / 10.9	38	L286L (TTG→CTG)	accA	acetylCoA carboxylase, carboxytransferase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (0/9);  new base C (29/0);  total (29/9)
Fisher's exact test for biased strand distribution p-value = 6.13e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCACCGTAACCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACG  >  W3110S.gb/209443‑209515                                 ||                                        gCTCACCGTAACCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCTGGc             <  1:1038672/62‑1 (MQ=255) gCTCACCGTAACCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCTGGc             <  1:98441/62‑1 (MQ=255) gCTCACCGTAACCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCTGGc             <  1:1060150/62‑1 (MQ=255) gCTCACCGTAACCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCTGGc             <  1:1061804/62‑1 (MQ=255) gCTCACCGTAACCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCTGGc             <  1:950256/62‑1 (MQ=255) gCTCACCGTAACCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCTGGc             <  1:2528/62‑1 (MQ=255) gCTCACCGTAACCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCTGGc             <  1:880158/62‑1 (MQ=255) gCTCACCGTAACCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCTGGc             <  1:75951/62‑1 (MQ=255) gCTCACCGTAACCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCTGGc             <  1:790528/62‑1 (MQ=255)            cccGGAAGCGATGGCGGCGTCTCTGAAAGCGCAActg                           >  1:1212361/1‑37 (MQ=25)            cccGGAAGCGATGGCGGCGTCTCTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGAc   >  1:256024/1‑61 (MQ=255)            cccGGAAGCGATGGCGGCGTCTCTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGAc   >  1:959022/1‑61 (MQ=255)            cccGGAAGCGATGGCGGCGTCTCTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGAc   >  1:248031/1‑61 (MQ=255)            cccGGAAGCGATGGCGGCGTCTCTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGAc   >  1:807646/1‑61 (MQ=255)            cccGGAAGCGATGGCGGCGTCTCTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGAc   >  1:216529/1‑61 (MQ=255)            cccGGAAGCGATGGCGGCGTCTCTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACg  >  1:753221/1‑62 (MQ=255)            cccGGAAGCGATGGCGGCGTCTCTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACg  >  1:760605/1‑62 (MQ=255)            cccGGAAGCGATGGCGGCGTCTCTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACg  >  1:831357/1‑62 (MQ=255)            cccGGAAGCGATGGCGGCGTCTCTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACg  >  1:736806/1‑62 (MQ=255)            cccGGAAGCGATGGCGGCGTCTCTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACg  >  1:829624/1‑62 (MQ=255)            cccGGAAGCGATGGCGGCGTCTCTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACg  >  1:302331/1‑62 (MQ=255)            cccGGAAGCGATGGCGGCGTCTCTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACg  >  1:896613/1‑62 (MQ=255)            cccGGAAGCGATGGCGGCGTCTCTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACg  >  1:993526/1‑62 (MQ=255)            cccGGAAGCGATGGCGGCGTCTCTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACg  >  1:609129/1‑62 (MQ=255)            cccGGAAGCGATGGCGGCGTCTCTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACg  >  1:533099/1‑62 (MQ=255)            cccGGAAGCGATGGCGGCGTCTCTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACg  >  1:433064/1‑62 (MQ=255)            cccGGAAGCGATGGCGGCGTCTCTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACg  >  1:365055/1‑62 (MQ=255)            cccGGAAGCGATGGCGGCGTCTCTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACg  >  1:1036049/1‑62 (MQ=255)            cccGGAAGCGATGGCGGCGTCTCTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACg  >  1:208099/1‑62 (MQ=255)            cccGGAAGCGATGGCGGCGTCTCTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACg  >  1:195502/1‑62 (MQ=255)            cccGGAAGCGATGGCGGCGTCTCTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACg  >  1:175008/1‑62 (MQ=255)            cccGGAAGCGATGGCGGCGTCTCTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACg  >  1:1472481/1‑62 (MQ=255)            cccGGAAGCGATGGCGGCGTCTCTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACg  >  1:1402542/1‑62 (MQ=255)            cccGGAAGCGATGGCGGCGTCTCTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACg  >  1:1284412/1‑62 (MQ=255)            cccGGAAGCGATGGCGGCGTCTCTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACg  >  1:1188419/1‑62 (MQ=255)            cccGGAAGCGATGGCGGCGTCTCTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACg  >  1:1108229/1‑62 (MQ=255)            cccGGAAGCGATGGCGGCGTCTCTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACg  >  1:1054361/1‑62 (MQ=255)            cccGGAAGCGATGGCGGCGTCTCTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACg  >  1:1036567/1‑62 (MQ=255)                                 ||                                        GCTCACCGTAACCCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACG  >  W3110S.gb/209443‑209515

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