Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,654,915	T→C	44.4%	L134L (CTT→CTC)	glgP →	glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,654,915	0	T	C	44.4%	2.7 / 25.9	36	L134L (CTT→CTC)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (15/5);  new base C (6/10);  total (21/15)
Fisher's exact test for biased strand distribution p-value = 4.09e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGATTGATGAAGAAAATGACCCAGGCCTCGGTAACGGTGGCCTGGGACGTCTGGCGGCTTGCTTCCTTGATTCTCTGGCGACGTTAGGGTTGCCGGGGCGCGGTTACGGCATCCGCTATGACTACGG  >  minE/2654847‑2654974                                                                     |                                                            cTGATTGATGAAGAAAATGACCCAGGCCTCGGTAACGGTGGCCTGGGACGTCTGGCGGCTTGCTTCCtt                                                             >  1:2222656/1‑69 (MQ=255)   gattgatGAAGAAAATGACCCAGGCCTCGGTAACGGTGGCCTGGGACGTCTGGCGGCTTGCTTCCTTGa                                                           >  1:3098344/1‑69 (MQ=255)           aagaaAATGACCCAGGCCTCGGTAACGGTGGCCTGGGACGTCTGGCGGCTTGCTTCCTTGGTtc                                                        >  1:1883332/1‑64 (MQ=255)           aagaaAATGACCCAGGCCTCGGTAACGGTGGCCTGGGACGTCTGGCGGCTTGCTTCCTTGATtc                                                        >  1:1087493/1‑64 (MQ=255)           aagaaAATGACCCAGGCCTCGGTAACGGTGGCCTGGGACGTCTGGCGGCTTGCTTCCTTGATtc                                                        >  1:259334/1‑64 (MQ=255)           aagaaAATGACCCAGGCCTCGGTAACGGTGGCCTGGGACGTCTGGCGGCTTGCTTCCTTGATtc                                                        >  1:731000/1‑64 (MQ=255)           aagaaAATGACCCAGGCCTCGGTAACGGTGGCCTGGGACGTCTGGCGGCTTGCTTCCTTGATtc                                                        >  1:658187/1‑64 (MQ=255)           aagaaAATGACCCAGGCCTCGGTAACGGTGGCCTGGGACGTCTGGCGGCTTGCTTCCTTGATtc                                                        >  1:2864750/1‑64 (MQ=255)           aagaaAATGACCCAGGCCTCGGTAACGGTGGCCTGGGACGTATGGCGGCTTGCTTCCTTGATtc                                                        >  1:3033408/1‑64 (MQ=255)                                       tGGCCTGGGCCGTGTGGCGGCTTGCTTCCTCGATTCTCTGGCGACGTTAGGGTTGCCGGGGCGCGGTTa                       <  1:1646214/69‑1 (MQ=38)                                       tGGCCTGGGCCGTCTGGCGGCTTGCTTCCTCGATTCTCTGGCGACGTTAGGGTTGCCGGGGCGCGGTTa                       <  1:1398148/69‑1 (MQ=255)                                       tGGCCTGGGCCGTCTGGCGGCTTGCTTCCTCGATTCTCTGGCGACGTTAGGGTTGCCGGGGCGCGGTTa                       <  1:2120900/69‑1 (MQ=255)                                       tGGCCTGGGCCGTCTGGCGGCTTGCTTCCTCGATTCTCTGGCGACGTTAGGGTTGCCGGGGCGCGGTTa                       <  1:2192441/69‑1 (MQ=255)                                       tGGCCTGGGCCGTCTGGCGGCTTGCTTCCTCGATTCTCTGGCGACGTTAGGGTTGCCGGGGCGCGGTTa                       <  1:411679/69‑1 (MQ=255)                                       tGGCCTGGGCCGTCTGGCGGCTTGCTTCCTCGATTCTCTGGCGACGTTAGGGTTGCCGGGGCGCGGTTa                       <  1:2552258/69‑1 (MQ=255)                                       tGGCCTGGGCCGTCTGGCGGCTTGCTTCCTCGATTCTCTGGCGACGTTAGGGTTGCCGGGGCGCGGTTa                       <  1:2525236/69‑1 (MQ=255)                                       tGGCCTGGGCCGTCTGGCGGCTTGCTTCATCGATTCTCTGGCGACGTTAGGGTTGCCGGGGCGCGGTTa                       <  1:848624/69‑1 (MQ=38)                                        ggCCTGGGCCGTCTGGCGGCTTGCTTCCTCGATTCTCTGGCGACGTTAggg                                        >  1:1163988/1‑51 (MQ=255)                                        ggCCTGGGCCGTCTGGCGGCTTGCTTCCTCGATTCTCTGGCGACGTTAggg                                        >  1:3359581/1‑51 (MQ=255)                                        ggCCTGGGCCGTCTGGCGGCTTGCTTCCTCGATTCTCTGGCGACGTTAggg                                        >  1:2432437/1‑51 (MQ=255)                                        ggCCTGGGCCGTCTGGCGGCTTGCTTCCTCGATTCTCTGGCGACGTTAggg                                        >  1:1992017/1‑51 (MQ=255)                                        ggCCTGGGCCGTCTGGCGGCTTGCTTCCTCGATTCTCTGGCGACGTTAggg                                        >  1:1500034/1‑51 (MQ=255)                                        ggCCTGGGCCGTCTGGCGGCTTGCTTCCTCGATTCTCTGGCGACGTTAggg                                        >  1:1369959/1‑51 (MQ=255)                                        ggCCTGGGCCGTCTGGCGGCTTGCTTCCTCGATTCTCTGGCGACGTTAGGGTTGCCGGGGCGCGGTTa                       <  1:1207158/68‑1 (MQ=255)                                        ggCCTGGGCCGTCTGGCGGCTTGCTTCCTCGATTCTCTGGCGACGTTAGGGTTGCCGGGGCGCGGTTa                       <  1:1863303/68‑1 (MQ=255)                                                aCGTCTGGCGGCTTGCTTCCTTGATTCTCTGGCGACGTTAGGGTTGCCGGGGCGCGGTTACGGCATcc               >  1:2713950/1‑68 (MQ=255)                                                     tGGCGGCTTGCTTCCTTGATTCTCTGGCGACGTTAGGGTTGCCGGGGCGCGGTTACGGCATCCGc             <  1:943367/65‑1 (MQ=255)                                                     tGGCGGCTTGCTTCCTTGATTCTCTGGCGACGTTAGGGTTGCCGGGGCGCGGTTACGGCATCCGc             <  1:2510321/65‑1 (MQ=255)                                                     tGGCGGCTTGCTTCCTTGATTCTCTGGCGACGTTAGGGTTGCCGGGGCGCGGTTACGGCATCCGc             <  1:1229887/65‑1 (MQ=255)                                                     tGGCGGCTTGCTTCCTTGATTCTCTGGCGACGTTAGGGTTGCCGGGGCGCGGTTACGGCATCCGc             <  1:1338458/65‑1 (MQ=255)                                                     tGGCGGCTTGCTTCCTTGATTCTCTGGCGACGTTAGGGTTGCCGGGGCGCGGTTACGGCATCCGc             <  1:2972292/65‑1 (MQ=255)                                                             tgcttCCTTGATTCTCTGGCGACGTTAGGGTTGCCGGGGCGCGGTTACGGCATCCGCTATGACTACgg  >  1:442558/1‑68 (MQ=255)                                                             tgcttCCTTGATTCTCTGGCGACGTTAGGGTTGCCGGGGCGCGGTTACGGCATCCGCTATGACTACgg  >  1:446767/1‑68 (MQ=255)                                                             tgcttCCTTGATTCTCTGGCGACGTTAGGGTTGCCGGGGCGCGGTTACGGCATCCGCTATGACTACgg  >  1:596334/1‑68 (MQ=255)                                                             tgcttCCTTGATTCTCTGGCGACGTTAGGGTTGCCGGGGCGCGGTTACGGCATCCGCTATGACTACgg  >  1:339273/1‑68 (MQ=255)                                                             tgcttCCTTGATTCTCTGGCGACGTTAGGGTTGCCGGGGCGCGGTTACGGCATCCGCTATGACTACgg  >  1:1601169/1‑68 (MQ=255)                                                             tgcttCCTTGATTCTCTGGCGACGTTAGGGTTGCCGGGGCGCGGTTACGGCATCCGCTATGACTACgg  >  1:2935308/1‑68 (MQ=255)                                                                     |                                                            CTGATTGATGAAGAAAATGACCCAGGCCTCGGTAACGGTGGCCTGGGACGTCTGGCGGCTTGCTTCCTTGATTCTCTGGCGACGTTAGGGTTGCCGGGGCGCGGTTACGGCATCCGCTATGACTACGG  >  minE/2654847‑2654974

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