Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,083,454	C→T	37.8%	S160S (AGC→AGT)	yqjA →	conserved inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,083,454	0	C	T	37.8%	13.1 / 14.5	29	S160S (AGC→AGT)	yqjA	conserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (10/8);  new base T (6/5);  total (16/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

TGCTGCCGACGATTGCCGGGTTATCAGGGCTGAATAACGCGCGCTTTCAGTTTTTCAACTGGATGAGCGGTCTGCTGTGGGTATTGATCCTGACAACTCTGGGTTACATGCTCGGCAAAACGCCGGT  >  minE/2083387‑2083513                                                                    |                                                            tgctgcCGACGATTGCCGGGTTATCAGGGCTGAATAACGCGCGCTTTCAGTTTTTCAACTGGATGAGCg                                                            <  1:2135760/69‑1 (MQ=255)  gctgaCGACGATTGCCGGGGTATCAGGGCTGAATAACGCGCGCTTTCAGTTTTTCAACTGGATGAGCg                                                            <  1:58734/68‑1 (MQ=255)        gacgaTTGCCGGGTTATCAGGGCTGAATAACGCGCGCTTTCAGTTCTTCAACTGGATGAGTGGTctg                                                       <  1:592795/67‑1 (MQ=255)        gacgaTTGCCGGGTTATCAGGGCTGAATAACGCGCGCTTTCAGTTCTTCAACTGGATGAGTGGTctg                                                       <  1:2948837/67‑1 (MQ=255)        gacgaTTGCCGGGTTATCAGGGCTGAATAACGCGCGCTTTCAGTTCTTCAACTGGATGAGTGGTctg                                                       <  1:2379453/67‑1 (MQ=255)        gacgaTTGCCGGGTTATCAGGGCTGAATAACGCGCGCTTTCAGTTCTTCAACTGGATGAGTGGTctg                                                       <  1:1322676/67‑1 (MQ=255)        gacgaTTGCCGGGTTATCAGGGCTGAATAACGCGCGCTTTCAGTTCTTCAACTGGATGAGTGGTctg                                                       <  1:1477871/67‑1 (MQ=255)           gaTTGCCGGGTTATCAGGGCTGAATAACGCGCGCTTTCAGTTTTTCAACTGGATGAGCGGTCTGCtgtg                                                  >  1:1138842/1‑69 (MQ=255)                   ggTTATCAGGGCTGAATAACGCGCGCTTTCAGTTTTTCAACTGGATGAGCGGTCTGCTGTGGGTATTGa                                          <  1:1048410/69‑1 (MQ=255)                   ggTTATCAGGGCTGAATAACGCGCGCTTTCAGTTTTTCAACTGGATGAGCGGTCTGCTGTGGGTATTGa                                          <  1:1199535/69‑1 (MQ=255)                          aGGGCTGAATAACGCGCGCTTTCAGTTTTTCAACTGGATGAGCGGTCTGCTGTGGGTATTGATCCTGAc                                   >  1:2993335/1‑69 (MQ=255)                          aGGGCTGAATAACGCGCGCTTTCAGTTTTTCAACTGGATGAGCGGTCTGCTGTGGGTATTGATCCTGAc                                   <  1:1486725/69‑1 (MQ=255)                          aGGGCTGAATAACGCGCGCTTTCAGTTTTTCAACTGGATGAGCGGTCTGCTGTGGGTATTGATCCTGAc                                   >  1:1781330/1‑69 (MQ=255)                          aGGGCTGAATAACGCGCGCTTTCAGTTTTTCAACTGGATGAGCGGTCTGCTGTGGGTATTGATCCTGAc                                   >  1:2343031/1‑69 (MQ=255)                          aGGGCTGAATAACGCGCGCTTTCAGTTTTTCAACTGGATGAGCGGTCTGCTGTGGGTATTGATCCTGAc                                   >  1:138414/1‑69 (MQ=255)                          aGGGCTGAATAACGCGCGCTTTCAGTTTTTCAACTGGATGAGCGGTCTGCTGTGGGTATTGATCCTGAc                                   >  1:2164644/1‑69 (MQ=255)                          aGGGCTGAATAACGCGCGCTTTCAGTTTTTCAACTGGATGAGCGGTCTGCTGTGGGTATTGATCCTGAc                                   >  1:2168527/1‑69 (MQ=255)                                      cgcgcgCTTTCAGTTTTTCAACTGGATGAGCGGGCTGCTGTGGGTATTGATCCTGAc                                   <  1:179384/57‑1 (MQ=255)                                       gcgcgcTTTCAGTTTTTCAACTGGATGAGCGGTCTGCTGTGGGTATTGATCCTGACAACTCTGGGTTa                       <  1:1034105/68‑1 (MQ=255)                                       gcgcgcTTTCAGTTTTTCAACTGGATGAGCGGTCTGCAGTGGGTATTGATCCTGACAACTCTGGGTTa                       <  1:340818/68‑1 (MQ=255)                                                     cttCAACTGGATGAGTGGTCTGCTGTGGGTGTTGATCCTGACAACTCTGGGTTACATGCTCGGCAAAAc        >  1:2119447/2‑69 (MQ=255)                                                     cttCAACTGGATGAGTGGTCTGCTGTGGGTGTTGATCCTGACAACTCTGGGTTACATGCTCGGCAAAAc        >  1:1211061/2‑69 (MQ=255)                                                     cttCAACTGGATGAGTGGTCTGCTGTGGGTGTTGATCCTGACAACTCTGGGTTACATGCTCGGCAAAAc        >  1:2697407/2‑69 (MQ=255)                                                     cttCAACTGGATGAGTGGTCTGCTGTGGGTGTTGATCCTGACAACTCTGGGTTACATGCTCGGCAAAAc        >  1:192206/2‑69 (MQ=255)                                                     cttCAACTGGATGAGTGGTCTGCTGTGGGTGTTGATCCTGACAACTCTGGGTTACATGCTCGGCAAAAc        >  1:1391008/2‑69 (MQ=255)                                                     cttCAACTGGATGAGTGGTCTGCTGTGGGTGTTGATCCTGACAACTCTGGGTTACATGCTCGGCAAAAc        >  1:5449/2‑69 (MQ=255)                                                           cTGGATGAGCGGTCTGCTGTGGGTATTGATCCTGACAACTCTGGGTTACATGCTCGGCAAAACGCCGGt  >  1:274335/1‑69 (MQ=255)                                                           cTGGATGAGCGGTCTGCTGTGGGTATTGATCCTGACAACTCTGGGTTACATGCTCGGCAAAACGCCGGt  >  1:1114075/1‑69 (MQ=255)                                                           cTGGATGAGCGGTCTGCTGTGGGTATTGATCCTGACAACTCTGGGTTACATGCTCGGCAAAACGCCGGt  >  1:1407992/1‑69 (MQ=255)                                                                    |                                                            TGCTGCCGACGATTGCCGGGTTATCAGGGCTGAATAACGCGCGCTTTCAGTTTTTCAACTGGATGAGCGGTCTGCTGTGGGTATTGATCCTGACAACTCTGGGTTACATGCTCGGCAAAACGCCGGT  >  minE/2083387‑2083513

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