Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,062,527	A→G	61.4%	T236T (ACA→ACG)	dnaG →	DNA primase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,062,527	0	A	G	61.4%	12.4 / 16.8	31	T236T (ACA→ACG)	dnaG	DNA primase
Reads supporting (aligned to +/- strand):  ref base A (7/5);  new base G (8/11);  total (15/16)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGTTTTGGCGGGCGCGTGCTGGGCAACGATACCCCCAAATACCTGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTCTTTATGAAGCGCAGCAGGATAACGCTGAACCCA  >  minE/2062467‑2062594                                                             |                                                                    tGGTTTTGGCGGGCGGGTGCTGGGCAACGATACCCCCAAATACCTGAACTCGCCGGAAACAGACAttt                                                              <  1:766923/68‑1 (MQ=255)      ttGGCGGGCGCGTGCTGGGCAACGATACCCCCAAATACCTGAACTCGCCGGAAACAGACa                                                                 <  1:2565038/60‑1 (MQ=255)      ttGGCGGGCGCGTGCTGGGCAACGATACCCCCAAATACCTGAACTCGCCGGAAACAGACa                                                                 <  1:2276750/60‑1 (MQ=255)      ttGGCGGCCGCGTGCTGGGCAACGATACCCCCAAATACCTGAACTCGCCGGAAACAGACa                                                                 <  1:1549848/60‑1 (MQ=255)                           aagATACCCCCAAATACCTGAACTCGCCGGAAACGGACATTTTCCATAAAGGCCGCCAGCTTTa                                        <  1:368039/62‑1 (MQ=255)                           aCGATACCCCCAAATACCTGAACTCGCCGGAAACGGACATTTTCCATAAAGGCCGCCAGCTTTa                                        <  1:1164052/64‑1 (MQ=255)                           aCGATACCCCCAAATACCTGAACTCGCCGGAAACGGACATTTTCCATAAAGGCCGCCAGCTTTa                                        <  1:919470/64‑1 (MQ=255)                           aCGATACCCCCAAATACCTGAACTCGCCGGAAACGGACATTTTCCATAAAGGCCGCCAGCTTTa                                        <  1:845995/64‑1 (MQ=255)                           aCGATACCCCCAAATACCTGAACTCGCCGGAAACGGACATTTTCCATAAAGGCCGCCAGCTTTa                                        <  1:806432/64‑1 (MQ=255)                           aCGATACCCCCAAATACCTGAACTCGCCGGAAACGGACATTTTCCATAAAGGCCGCCAGCTTTa                                        <  1:504888/64‑1 (MQ=255)                           aCGATACCCCCAAATACCTGAACTCGCCGGAAACGGACATTTTCCATAAAGGCCGCCAGCTTTa                                        <  1:3079207/64‑1 (MQ=255)                           aCGATACCCCCAAATACCTGAACTCGCCGGAAACGGACATTTTCCATAAAGGCCGCCAGCTTTa                                        <  1:1752487/64‑1 (MQ=255)                           aCGATACCCCCAAATACCTGAACTCGCCGGAAACGGACATTTTCCATAAAGGCCGCCAGCTTTa                                        <  1:10581/64‑1 (MQ=255)                           aCGATACCCCCAAATACCTGAACTCGCCGGAAACGGACATTTTCCATAAAGGCCGCCAGCTTTa                                        <  1:1452537/64‑1 (MQ=255)                           aCGATACCCCCAAATACCTGAACTCGCCGGAAACGGACATTTTCCATAAAGGCCGCCAGCTTTa                                        <  1:147612/64‑1 (MQ=255)                                         tACCTGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTCTTTATGAAgcgc                       >  1:78465/1‑67 (MQ=255)                                          aCCTGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTCTTTATGAAGCgcagc                    >  1:52881/1‑69 (MQ=255)                                                          aaCGGACATTTTCCATAAAGGCCGCCAGCTTTACGGTCTTTATGAAGGGCAGCAGGATAACGCTGAAc     >  1:2537082/1‑68 (MQ=255)                                                          aaCGGACATTTTCCATAAAGGCCGCCAGCTTTACGGTCTTTATGAAGCGCAGCAGGATAACGCTGAAc     >  1:2789374/1‑68 (MQ=255)                                                          aaCGGACATTTTCCATAAAGGCCGCCAGCTTTACGGTCTTTATGAAGCGCAGCAGGATAACGCTGAAc     >  1:2811072/1‑68 (MQ=255)                                                          aaCGGACATTTTCCATAAAGGCCGCCAGCTTTACGGTCTTTATGAAGCGCAGCAGGATAACGCTGAAc     >  1:206062/1‑68 (MQ=255)                                                          aaCGGACATTTTCCATAAAGGCCGCCAGCTTTACGGTCTTTATGAAGCGCAGCAGGATAACGCTGAAc     >  1:327717/1‑68 (MQ=255)                                                          aaCGGACATTTTCCATAAAGGCCGCCAGCTTTACGGTCTTTATGAAGCGCAGCAGGATAACGCTGAAc     >  1:274099/1‑68 (MQ=255)                                                          aaCGGACATTTTCCATAAAGGCCGCCAGCTTTACGGTCTTTATGAAGCGCAGCAGGATAACGCTGAAc     >  1:902902/1‑68 (MQ=255)                                                          aaCGGACATTTTCCATAAAGGCCGCCAGCTTTACGGTCTTTATGAAGCGCAGCAGGATAACGCTGAAc     >  1:2553171/1‑68 (MQ=255)                                                            cAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTCTTTATGAAGCGCAGCAGGATAACGCTGAACCCa  >  1:2865401/1‑69 (MQ=255)                                                            cAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTCTTTATGAAGCGCAGCAGGATAACGCTGAACCCa  >  1:2994942/1‑69 (MQ=255)                                                            cAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTCTTTATGAAGCGCAGCAGGATAACGCTGAACCCa  <  1:491533/69‑1 (MQ=255)                                                            cAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTCTTTATGAAGCGCAGCAGGATAACGCTGAACCCa  >  1:1347130/1‑69 (MQ=255)                                                            cAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTCTTTATGAAGCGCAGCAGGATAACGCTGAACCCa  >  1:1258709/1‑69 (MQ=255)                                                            cAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTCTTTATGAAGCGCAGCAGCATAACGCTGAAccc   >  1:2143419/1‑68 (MQ=255)                                                             |                                                                    TGGTTTTGGCGGGCGCGTGCTGGGCAACGATACCCCCAAATACCTGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTCTTTATGAAGCGCAGCAGGATAACGCTGAACCCA  >  minE/2062467‑2062594

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