Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,450,298	A→G	58.8%	W78R (TGG→CGG)	yifE ←	conserved hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,450,298	0	A	G	58.8%	5.2 / 17.1	17	W78R (TGG→CGG)	yifE	conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (5/2);  new base G (4/6);  total (9/8)
Fisher's exact test for biased strand distribution p-value = 3.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

TGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCACACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGTT  >  minE/2450255‑2450357                                            |                                                            tgAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCAcgcg                                                       <  1:191757/50‑1 (MQ=255) tgAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCAcgcg                                                       <  1:49932/50‑1 (MQ=255) tgAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCAcgcg                                                       <  1:460299/50‑1 (MQ=255)        gTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCt                                                >  1:258567/1‑50 (MQ=255)        gTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCt                                                >  1:353304/1‑50 (MQ=255)                cGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTc                            <  1:354671/62‑1 (MQ=255)                  cTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGc                          <  1:339543/62‑1 (MQ=255)                  cTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGc                          <  1:375031/62‑1 (MQ=255)                  cTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGc                          >  1:39139/1‑62 (MQ=255)                  cTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGc                          >  1:246442/1‑62 (MQ=255)                                    cTTGGACCACACGCGTTCTGCTTCTTTCACTGGCTCACGTTCGCCACGGCATACTGCTACGa        <  1:552251/62‑1 (MQ=255)                                    cTTGGACCACACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGa        <  1:467077/62‑1 (MQ=255)                                         aCCACACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGt   >  1:189862/1‑62 (MQ=255)                                         aCCACACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGt   >  1:472560/1‑62 (MQ=255)                                          ccACACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGtt  >  1:64454/1‑62 (MQ=255)                                          ccACACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGACACGGCATACTGCTACGAAGAGtt  >  1:285121/1‑62 (MQ=255)                                          ccACACGCGTTCTGCTTCAGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGtt  >  1:229819/1‑62 (MQ=255)                                            |                                                            TGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCACACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGTT  >  minE/2450255‑2450357

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