Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,990,864	A→G	69.9%	Q286Q (CAA→CAG)	mutY →	adenine DNA glycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,990,864	0	A	G	69.9%	38.7 / 16.5	43	Q286Q (CAA→CAG)	mutY	adenine DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base A (5/8);  new base G (4/26);  total (9/34)
Fisher's exact test for biased strand distribution p-value = 1.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.24e-01

GTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGG  >  minE/1990812‑1990921                                                     |                                                          gTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCgg                                            <  1:1642259/68‑1 (MQ=255) gTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCgg                                            <  1:3248465/68‑1 (MQ=255)                      aaCGACAGATTGCTGCCGATAACCTGACGCAGCTGACCGCGTTTCGGCATACCTTCAGCCATTTCCa                        >  1:215474/1‑67 (MQ=255)                      aaCGACAGATTGCTGCCGATAACCTGACGCAGCTGACCGCGTTTCGGCATACCTTCAGCCATTTCCa                        >  1:2813848/1‑67 (MQ=255)                      aaCGACAGATTGCTGCCGATAACCTGACGCAGCTGACCGCGTTTCGGCATACCTTCAGCCATTTCCa                        >  1:3104759/1‑67 (MQ=255)                      aaCGACAGATTGCTGCCGATAACCTGACGCAGCTGACCGCGTTTCGGCATACCTTCAGCCATTTCCa                        >  1:1238766/1‑67 (MQ=255)                        cGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTa                    <  1:1022035/69‑1 (MQ=255)                        cGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTa                    <  1:2251640/69‑1 (MQ=255)                          gCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGa                  <  1:1350023/69‑1 (MQ=255)                          gCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGa                  <  1:2828957/69‑1 (MQ=255)                          gCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGa                  <  1:2798976/69‑1 (MQ=255)                          gCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGa                  <  1:2478291/69‑1 (MQ=255)                                    gcCGATAACCTGACGCAGCTGACAGCGTTTCGGCATACCTTCAGCCGTTTCCACTTAGATATTGTGCCt        <  1:1091677/69‑1 (MQ=255)                                    gcCGATAACCTGACGCAGCTGACAGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCt        <  1:367881/69‑1 (MQ=255)                                    gcCGATAACCTGACGCAGCTGACAGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCt        <  1:354992/69‑1 (MQ=255)                                    gcCGATAACCTGACGCAGCTGACAGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCt        <  1:3263092/69‑1 (MQ=255)                                    gcCGATAACCTGACGCAGCTGACAGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCt        <  1:3211066/69‑1 (MQ=255)                                    gcCGATAACCTGACGCAGCTGACAGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCt        <  1:728326/69‑1 (MQ=255)                                    gcCGATAACCTGACGCAGCTGACAGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCt        <  1:871167/69‑1 (MQ=255)                                    gcCGATAACCTGACGCAGCTGACAGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCt        <  1:905425/69‑1 (MQ=255)                                    gcCGATAACCTGACGCAGCTGACAGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCt        <  1:1043690/69‑1 (MQ=255)                                    gcCGATAACCTGACGCAGCTGACAGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCt        <  1:2727369/69‑1 (MQ=255)                                    gcCGATAACCTGACGCAGCTGACAGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCt        <  1:2510283/69‑1 (MQ=255)                                    gcCGATAACCTGACGCAGCTGACAGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCt        <  1:1657693/69‑1 (MQ=255)                                    gcCGATAACCTGACGCAGCTGACAGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCt        <  1:1052686/69‑1 (MQ=255)                                    gcCGATAACCTGACGCAGCTGACAGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCt        <  1:1081578/69‑1 (MQ=255)                                    gcCGATAACCTGACGCAGCTGACAGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCt        <  1:110238/69‑1 (MQ=255)                                    gcCGATAACCTGACGCAGCTGACAGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCt        <  1:1406692/69‑1 (MQ=255)                                    gcCGATAACCTGACGCAGCTGACAGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCt        <  1:1411217/69‑1 (MQ=255)                                    gcCGATAACCTGACGCAGCTGACAGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCt        <  1:1446843/69‑1 (MQ=255)                                    gcCGATAACCTGACGCAGCTGACAGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCt        <  1:1510753/69‑1 (MQ=255)                                    gcCGATAACCTGACGCAGCTGACAGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCt        <  1:2056333/69‑1 (MQ=255)                                    gcCGATAACCTGACGCAGCTGACAGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCt        <  1:213218/69‑1 (MQ=255)                                    gcCGATAACCTGACGCAGCTGACAGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCt        <  1:2299553/69‑1 (MQ=255)                                    gcCGATAACCTGACGCAGCTGACAGCGTTTCGGCATACCTTCAGCCACTTCCACTTAGATATTGTGCCt        <  1:2691967/69‑1 (MQ=255)                                     ccGATAACCTGACGCAGCTGACAGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCt        <  1:1929615/68‑1 (MQ=255)                                     ccGATAACCTGACGCAGCTGACAGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCt        <  1:741873/68‑1 (MQ=255)                                          aaCCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCAc                       >  1:1405209/1‑48 (MQ=255)                                          aaCCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTgg  >  1:2914174/1‑69 (MQ=255)                                          aaCCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTgg  >  1:2984860/1‑69 (MQ=255)                                          aaCCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTgg  >  1:2465032/1‑69 (MQ=255)                                          aaCCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTgg  >  1:2908812/1‑69 (MQ=255)                                                aCGCAGCTGACAGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCt        <  1:2415779/57‑1 (MQ=255)                                                     |                                                          GTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGG  >  minE/1990812‑1990921

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑