Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,529,257	C→T	50.0%	T55T (ACG→ACA)	yfcL ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,529,257	0	C	T	50.0%	‑5.4 / 30.3	36	T55T (ACG→ACA)	yfcL	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (8/10);  new base T (11/7);  total (19/17)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

AGCGCCTGGTCACGCGGCGACAATTCACCAGCACCAATGGCTTTTTCCAGACTCTGGCTAACGGTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTC  >  minE/1529192‑1529317                                                                  |                                                             aGCGCCTGGTCACGCGGCGACAATTCACCAGCACCAATGGCTTTTTCCAGACTCTGGCTAACGGTCGTa                                                           >  1:1717842/1‑69 (MQ=255)  gcgcCTGGTCACGCGGTGACAATTCACCCGCGCCAATGGCTTTTTCCAGACTCTGGCTAACGGTTGTa                                                           >  1:1141408/1‑68 (MQ=255)  gcgcCTGGTCACGCGGTGACAATTCACCCGCGCCAATGGCTTTTTCCAGACTCTGGCTAACGGTTGTa                                                           >  1:3301481/1‑68 (MQ=255)  gcgcCTGGTCACGCGGTGACAATTCACCCGCGCCAATGGCTTTTTCCAGACTCTGGCTAACGGTTGTa                                                           >  1:2831688/1‑68 (MQ=255)  gcgcCTGGTCACGCGGTGACAATTCACCCGCGCCAATGGCTTTTTCCAGACTCTGGCTAACGGTTGTa                                                           >  1:2197430/1‑68 (MQ=255)  gcgcCTGGTCACGCGGTGACAATTCACCCGCGCCAATGGCTTTTTCCAGACTCTGGCTAACGGTTGTa                                                           >  1:2118921/1‑68 (MQ=255)  gcgcCTGGTCACGCGGTGACAATTCACCCGCGCCAATGGCTTTTTCCAGACTCTGGCTAACGGTTGTa                                                           >  1:20817/1‑68 (MQ=255)  gcgcCTGGTCACGCGGTGACAATTCACCCGCGCCAATGGCTTTTTCCAGACTCTGGCTAACGGTTGTa                                                           >  1:2040647/1‑68 (MQ=255)  gcgcCTGGTCACGCGGTGACAATTCACCCGCGCCAATGGCTTTTTCCAGACTCTGGCTAACGGTTGTa                                                           >  1:1987156/1‑68 (MQ=255)  gcgcCTGGTCACGCGGTGACAATTCACCCGCGCCAATGGCTTTTTCCAGACTCTGGCTAACGGTTGTa                                                           >  1:1898416/1‑68 (MQ=255)  gcgcCTGGTCACGCGGTGACAATTCACCCGCGCCAATGGCTTTTTCCAGACTCTGGCTAACGGTTGTa                                                           >  1:1744965/1‑68 (MQ=255)  gcgcCTGGTCACGCGGTGACAATTCACCCGCGCCAATGGCTTTTTCCAGACTCTGGCTAACGGTTGTa                                                           >  1:1619735/1‑68 (MQ=255)   cgcCTGGTCACGCGGCGACAATTCACCAGCACCAATGGCTTTTTCCAGACTCTGGCTAACGGTCGTATg                                                         <  1:2778017/69‑1 (MQ=255)   cgcCTGGTCACGCGGCGACAATTCACCAGCACCAATGGCTTTTTCCAGACTCTGGCTAACGGTCGTATg                                                         <  1:3346363/69‑1 (MQ=255)    gcCTGGTCACGCGGCGACAATTCACCAGCACCATTGGCTTTTTCCAGAGTCTGGCTAACGGTCGTATg                                                         <  1:618649/68‑1 (MQ=255)       tGGTCACGCGGCGACAATTCACCAGCACCAATGGCTTTTTCCAGACTCTGGCTAACGGTCGTATGCAcc                                                     <  1:534212/69‑1 (MQ=255)        ggTCACGCGGCGACAATTCACCAGCACCAATGGCTTTTTCCAGACTCTGGCTAACGGTCGTATGCAcc                                                     <  1:2724270/68‑1 (MQ=255)           cACGCGGCGACAATTCACCAGCACCAATGGCTTTTTCCAGACTCTGGCTAACGGTCGTATGCACCGc                                                   >  1:2519614/1‑67 (MQ=255)                 gcgACAATTCACCAGCACCAATGGCTTTTTCCAGACTCTGGCTAACGGTCGTATGCACCGCCTGAGCgg                                           <  1:874159/69‑1 (MQ=255)                      aaTTCACCCGCGCCAATGGCTTTTTCCAGACTCTGGCTAACGGTTGTATGTACCGCCTGAGCGGAGTg                                       <  1:624053/68‑1 (MQ=255)                      aaTTCACCCGCGCCAATGGCTTTTTCCAGACTCTGGCTAACGGTTGTATGTACCGCCTGAGCGGAGTg                                       <  1:821134/68‑1 (MQ=255)                      aaTTCACCCGCGCCAATGGCTTTTTCCAGACTCTGGCTAACGGTTGTATGTACCGCCTGAGCGGAGTg                                       <  1:483136/68‑1 (MQ=255)                      aaTTCACCCGCGCCAATGGCTTTTTCCAGACTCTGGCTAACGGTTGTATGTACCGCCTGAGCGGAGTg                                       <  1:977265/68‑1 (MQ=255)                      aaTTCACCCGCGCCAATGGCTTTTTCCAGACTCTGGCTAACGGTTGTATGTACCGCCTGAGCGGAGTg                                       <  1:435577/68‑1 (MQ=255)                      aaTTCACCCGCGCCAATGGCTTTTTCCAGACTCTGGCTAACGGTTGTATGTACCGCCTGAGCGGAGTg                                       <  1:3431122/68‑1 (MQ=255)                      aaTTCACCCGCGCCAATGGCTTTTTCCAGACTCTGGCTAACGGTTGTATGTACCGCCTGAGCGGAGTg                                       <  1:1179411/68‑1 (MQ=255)                                        gCTTTTTCCAGACTCTGGCTACCGGTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGt                    >  1:2419140/1‑69 (MQ=255)                                        gCTTTTTCCAGACTCTGGCTAACGGTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAg                     <  1:1591248/68‑1 (MQ=255)                                        gCTTTTTCCAGACTCTGGCTAACGGTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGt                    >  1:3358380/1‑69 (MQ=255)                                        gCTTTTTCCAGACTCTGGCTAACGGTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGt                    >  1:2983018/1‑69 (MQ=255)                                        gCTTTTTCCAGACTCTGGCTAACGGTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGt                    >  1:2879430/1‑69 (MQ=255)                                        gCTTTTTCCAGACTCTGGCTAACGGTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGt                    >  1:548865/1‑69 (MQ=255)                                        gCTTTTTCCAGACTCTGGCTAACGGTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGt                    >  1:1763183/1‑69 (MQ=255)                                                           tAACGGTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTc  <  1:419484/68‑1 (MQ=255)                                                           tAACGGTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTc  <  1:698237/68‑1 (MQ=255)                                                           tAACGGTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTc  <  1:1672731/68‑1 (MQ=255)                                                                  |                                                             AGCGCCTGGTCACGCGGCGACAATTCACCAGCACCAATGGCTTTTTCCAGACTCTGGCTAACGGTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTC  >  minE/1529192‑1529317

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