Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,681,918	A→G	38.9%	R598R (CGT→CGC)	yhgF ←	predicted transcriptional accessory protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,681,918	0	A	G	38.9%	8.1 / 17.0	18	R598R (CGT→CGC)	yhgF	predicted transcriptional accessory protein
Reads supporting (aligned to +/- strand):  ref base A (4/7);  new base G (3/4);  total (7/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 = 9.32e-01

CACACCGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAGTGCCTGCTGTGTTGCTGCCAGAA  >  minE/2681876‑2681976                                           |                                                           cacaCCGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTGCGCAGTTCGCTGCTGTTa                                          <  1:88194/61‑1 (MQ=255) cacaCCGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTGCGCAGTTCGCTGCTGTTa                                          <  1:449960/61‑1 (MQ=255) cacaCCGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTGCGCAGTTCGCTGCTGTTa                                          <  1:38956/61‑1 (MQ=255)            tctcATCAGTAAAGTCAGACGCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCt                              <  1:544391/62‑1 (MQ=255)                     tAAAGTCAGACGCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAGTGc                      <  1:126733/61‑1 (MQ=255)                     tAAAGTCAGACGCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAGTGc                      <  1:453247/61‑1 (MQ=255)                     tAAAGTCAGACGCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAGTGc                      <  1:348733/61‑1 (MQ=255)                     tAAAGTCAGACGCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAGTGc                      <  1:133780/61‑1 (MQ=255)                       aaGTCAGACGCTTTCAGGTTGCGCAGTTCGCTGCTGTTACCCATCAGATc                               <  1:213701/50‑1 (MQ=255)                              aCGCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAGTGCCTGCTGTGtt            >  1:59017/1‑62 (MQ=255)                                gCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAGTGCCTGCTGTGTtgc          >  1:72562/1‑62 (MQ=255)                                gCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAGTGCCTGCTGTGTtgc          >  1:97232/1‑62 (MQ=255)                                gCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAGTGCCTGCTGTGTtga          >  1:131834/1‑61 (MQ=255)                                       ggTTGCGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAGTGCCTGCTGTGTTGCTGCCAGa   >  1:296381/1‑62 (MQ=255)                                       ggTTGCGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAGTGCCTGCTGTGTTGCTGCCAGa   >  1:317518/1‑62 (MQ=255)                                       ggTTGCGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAGTGCCTGCTGTGTTGCTGCCAGa   >  1:116981/1‑62 (MQ=255)                                        gTTACGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAGTGCCTGCTGTGTTGCTGCCAGaa  <  1:25823/62‑1 (MQ=255)                                        gTTACGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAGTGCCTGCTGTGTTGCTGCCAGaa  <  1:138409/62‑1 (MQ=255)                                           |                                                           CACACCGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAGTGCCTGCTGTGTTGCTGCCAGAA  >  minE/2681876‑2681976

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