Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	3,786,897	T→G	51.1%	D245A (GAC→GCC)	thlA_2 ←	Acetyl‑CoA acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	3,786,897	0	T	G	51.1%	7.6 / 15.0	15	D245A (GAC→GCC)	thlA_2	Acetyl‑CoA acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (4/3);  new base G (6/2);  total (10/5)
Fisher's exact test for biased strand distribution p-value = 6.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.15e-03

AATCGGAGTAGCGCGTGACCAAGGCTGCCGCCGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTC‑GGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCT  >  KT2440/3786781‑3787029                                                                                                                     |                                                                                                                                      aaTCGGAGTAGCGCGTGACCAAGGCTGCCGCCGCTCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTTCACGCCACCCGCCTGCACCGGGCGCACTCGGGCCAGGTCGGCCGCGGG‑GGTTGGGCGGGGGTGGc                                                                                                                    >  6:271292/1‑135 (MQ=255) aaTCGGAGTAGCGCGTGACCAAGGCTGCCGCCGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCCCCCGCGTGCCACCGGCGCAGTCGGGCCAGTGCGGCCGCCGG‑GGGAGGCCGTGGGTGGCt                                                                                                                   >  7:27528/1‑136 (MQ=255)      gAGTAGCGCGTGACCAAGGCTGCCGCCGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCCTGCCCCGGGCGCCGTCGGGGCCGTGCGGGCGCGGCGGGTGGGCCggggggggggggc                                                                                                                >  5:226361/1‑130 (MQ=255)         tAGCGCGTGACCAAGGCTGCCGCCGCGCCATCGACCACCGCGCACCTGTTGCCAGCCGTCTGCACGCCCCCGGCGTGCACCGGGCGCAGTCGGGCCCGTGCGGCCGCGGG‑CGGAGGCCGGGGGGGGGTGGCCTgg                                                                                                            >  3:129409/1‑135 (MQ=255)                        gcggccgccGCGCCATCGCCCCCCGCCCAACTGTTGCCCGCCGGCGCCCCCCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTC‑GGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTcacg                                                                                             <  3:172885/132‑1 (MQ=255)                        gCTGCCGCCGCGCCAGCGCCCCCCCCGCAAGTGGTGCCAGCGGGCTGCACGCCACGGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTC‑GGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTcacg                                                                                             <  6:266796/135‑1 (MQ=255)                                 gcgcCATCGCCCCCCGCGCAAGTGTGGCCGGCCGGCGGCCCGCCCCGGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGGC‑GGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGgca                                                                                   <  7:180005/136‑1 (MQ=255)                                                    aaCTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTC‑GGTAGGCCGTGGGTGGCTGTCCTGGCTGCCCTCACTCACGCCACCCCGCAAAGCAATGGCGCGtggctg                                                                >  7:51123/1‑136 (MQ=255)                                                    aaCTGTTGCCAGCCGTCTGCACGCCACCCGCGTGCACCGGGCGCCGTCGGGCCAGTGCGGCCGGGGC‑GGGAGGCCGGGGGGGGCGGGGCTGGGGGACCGCCCTCACCCCCCCCCGCCCAGCAAAGCCGCCtgggg                                                                 >  2:202907/1‑133 (MQ=255)                                                                                    cgcccggggCGCGGTGGGGCCCGGGGGGCGGGGTG‑GGGAGGCGGGGGGGGGCGGCCCGGGCTCCCCCCACACCCCCCAGGGGCCAATGAATTGCCGGGGGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTg                                 <  2:69355/129‑1 (MQ=255)                                                                                         cGGGCGCAGTCGGGCCAGTGCGGCCGCGTC‑GGTAGGCCGTGGGTGGCTGTCCTGGCTGCCCTCCCTCACGCCACCCCGCCAAGCAATTCCGCGGGGCTGGCCGCCGTCGACCTCCACTGCCTGCGCGGTGgcgcgg                           >  7:262377/1‑134 (MQ=255)                                                                                              gcAGTCGGGCCAGTGCGGCCGCGTC‑GGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCCGCAATGCAAAGCCGCGTGGCGGGCCGCCGTCTACCTCCAATGCCTGCGCGGTGCCGCTGGCgc                       >  1:207920/1‑134 (MQ=255)                                                                                                         aGTGCGCCCGCGGC‑GGTCGGCGGGGGGGGGCTGGCCGGGCGGACCTCACCCACGCCCCGCGGAAATGCAATGCCGGGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCtcgtcg            <  4:136393/135‑1 (MQ=255)                                                                                                              ggCCGCGGC‑GGTAGGCCGTGGGTGGCTGTCCTGGCTGCCCTCACTCACGCCCCGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGCTTTCCTAGTCGCACCa       >  2:187838/1‑135 (MQ=255)                                                                                                                   cgTC‑GGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCt  >  1:58544/1‑135 (MQ=255)                                                                                                                     |                                                                                                                                      AATCGGAGTAGCGCGTGACCAAGGCTGCCGCCGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTC‑GGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCT  >  KT2440/3786781‑3787029

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