Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	611,679	G→A	50.0%	Y551Y (TAC→TAT)	ycaO ←	conserved hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	611,679	0	G	A	50.0%	‑5.5 / 28.1	34	Y551Y (TAC→TAT)	ycaO	conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (9/8);  new base A (8/9);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01

ACGACTGATGTGCAGCAAACGCGTGCAGATCGCTATCTACCGGTTGCAGGCCGTAAAACGCCGCTTCGCCGCTCATTGCCGCACTGGCGGCTTCTACGGCATCTGCGCCGTA  >  minE/611627‑611738                                                     |                                                            acGACTGATGTGCAGCAAACGCGTGCAGATCGCTATCTACCGGTTGCAGGCCGTAAAACGCCGCTTcgc                                             >  1:1365801/1‑69 (MQ=255) acGACTGATGTGCAGCAAACGCGTGCAGATCGCTATCTACCGGTTGCAGGCCGTAAAACGCCGCTTcgc                                             >  1:2979467/1‑69 (MQ=255) acGACTGATGTGCAGCAAACGCGTGCAGATCGCTATCTACCGGTTGCAGGCCGTAAAACGCCGCTTcgc                                             >  1:287714/1‑69 (MQ=255) acGACTGATGTGCAGCAAACGCGTGCAGATCGCTATCTACCGGTTGCAGGCCGTAAAACGCCGCTTcgc                                             >  1:2544832/1‑69 (MQ=255)     cTGATGTGCAGCAAACGCGTGCAGATCGCTATCTACCGGTTGCAGGCCGTAAAACGCCGCTTCGCCGCt                                         <  1:1033378/69‑1 (MQ=255)     cTGATGTGCAGCAAACGCGTGCAGATCGCTATCTACCGGTTGCAGGCCGTAAAACGCCGCTTCGCCGCt                                         <  1:790856/69‑1 (MQ=255)     cTGATGTGCAGCAAACGCGTGCAGATCGCTATCTACCGGTTGCAGGCCGTAAAACGCCGCTTCGCCGCt                                         <  1:507658/69‑1 (MQ=255)     cTGATGTGCAGCAAACGCGTGCAGATCGCTATCTACCGGTTGCAGGCCGTAAAACGCCGCTTCGCCGCt                                         <  1:426740/69‑1 (MQ=255)     cTGATGTGCAGCAAACGCGTGCAGATCGCTATCTACCGGTTGCAGGCCGTAAAACGCCGCTTCGCCGCt                                         <  1:2340131/69‑1 (MQ=255)             cagcaAACGCGTGCAGATCGCTATCTACTGGTTGCAAGCCATAAAACGCCGCTTCGCCGCTCATTGCCg                                 <  1:1092155/69‑1 (MQ=255)             cagcaAACGCGTGCAGATCGCTATCTACTGGTTGCAAGCCATAAAACGCCGCTTCGCCGCTCATTGCCg                                 <  1:985842/69‑1 (MQ=255)             cagcaAACGCGTGCAGATCGCTATCTACTGGTTGCAAGCCATAAAACGCCGCTTCGCCGCTCATTGCCg                                 <  1:314063/69‑1 (MQ=255)             cagcaAACGCGTGCAGATCGCTATCTACTGGTTGCAAGCCATAAAACGCCGCTTCGCCGCTCATTGCCg                                 <  1:3128333/69‑1 (MQ=255)             cagcaAACGCGTGCAGATCGCTATCTACTGGTTGCAAGCCATAAAACGCCGCTTCGCCGCTCATTGCCg                                 <  1:2656571/69‑1 (MQ=255)             cagcaAACGCGTGCAGATCGCTATCTACTGGTTGCAAGCCATAAAACGCCGCTTCGCCGCTCATTGCCg                                 <  1:2256487/69‑1 (MQ=255)             cagcaAACGCGTGCAGATCGCTATCTACTGGTTGCAAGCCATAAAACGCCGCTTCGCCGCTCATTGCCg                                 <  1:1539667/69‑1 (MQ=255)             cagcaAACGCGTGCAGATCGCTATCTACTGGTTGCAAGCCATAAAACGCCGCTTCGCCGCTCATTGCCg                                 <  1:1860650/69‑1 (MQ=255)             cagcaAACGCGTGCAGATCGCTATCTACTGGTTGCAAGCCATAAAACGCCGCTTCGCCGCTCATTGCCg                                 <  1:1635585/69‑1 (MQ=255)               gcaAACGCGTGCAGATCGCTATCTACTGGTTGCAAGCCATAAAACGCCGCTTCGCCGCTCATTGCCGCa                               >  1:2524026/1‑69 (MQ=255)               gcaAACGCGTGCAGATCGCTATCTACTGGTTGCAAGCCATAAAACGCCGCTTCGCCGCTCATTGCCGCa                               >  1:1385783/1‑69 (MQ=255)               gcaAACGCGTGCAGATCGCTATCTACTGGTTGCAAGCCATAAAACGCCGCTTCGCCGCTCATTGCCGCa                               >  1:645264/1‑69 (MQ=255)               gcaAACGCGTGCAGATCGCTATCTACTGGTTGCAAGCCATAAAACGCCGCTTCGCCGCTCATTGCCGCa                               >  1:1490521/1‑69 (MQ=255)               gcaAACGCGTGCAGATCGCTATCTACTGGTTGCAAGCCATAAAACGCCGCTTCGCCGCTCATTGCCGCa                               >  1:386582/1‑69 (MQ=255)               gcaAACGCGTGCAGATCGCTATCTACTGGTTGCAAGCCATAAAACGCCGCTTCGCCGCTCATTGCCGCa                               >  1:2680260/1‑69 (MQ=255)               gcaAACGCGTGCAGATCGCTATCTACTGGTTGCAAGCCATAAAACGCCGCTTCGCCGCTCATTGCCGCa                               >  1:2399793/1‑69 (MQ=255)               gcaAACGCGTGCAGATCGCTATCTACTGGTTGCAAGCCATAAAACGCCGCTTCGCCGCTCATTGCCGCa                               >  1:2516240/1‑69 (MQ=255)                               cGCTATCTACCGGTTGCAGGCCGTAAAACGCCGCTTCGCCGCTCATTGCCGCACTGGCGGCTTCt                   >  1:886684/1‑65 (MQ=255)                               cGCTATCTACCGGTTGCAGGCCGTAAAACGCCGCTTCGCCGCTCATTGCCGCACTGGCGGCTTCTa                  >  1:3070453/1‑66 (MQ=255)                               cGCTATCTACCGGTTGCAGGCCGTAAAACGCCGCTTCGCCGCTCATTGCCGCACTGGCGGCTTCTa                  >  1:110650/1‑66 (MQ=255)                               cGCTATCTACCGGTTGCAGGCCGTAAAACGCCGCTTCGCCGCTCATTGCCGCACTGGCGGCTTCTACgg               >  1:3350109/1‑69 (MQ=255)                               cGCTATCTACCGGTTGCAGGCCGTAAAACGCCGCTTCGCCGCTCATTGCCGCACTGGCGGCTTCTACgg               >  1:2630909/1‑69 (MQ=255)                               cGCTAACTACCGGTTGCAGGCCGTAAAACGCCGCTTCGCCGCTCATTGCCGCACTGGCGGCTTCTa                  <  1:2249221/66‑1 (MQ=255)                                            ttGCAGGCCGTAAAACGCCGCTTCGCCGCTCATTGCCGCACTGGCGGCTTCTACGGCATCTGCGCCGTa  <  1:1721584/69‑1 (MQ=255)                                              gCAGGCCGTAAAACGCCGCTTCGCCGCTCATTGCCGCACTGGCGGCTTCTACGGc              <  1:3224359/55‑1 (MQ=255)                                                     |                                                            ACGACTGATGTGCAGCAAACGCGTGCAGATCGCTATCTACCGGTTGCAGGCCGTAAAACGCCGCTTCGCCGCTCATTGCCGCACTGGCGGCTTCTACGGCATCTGCGCCGTA  >  minE/611627‑611738

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