Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,904,690	C→T	29.3%	Q468Q (CAG→CAA)	treC ←	trehalose‑6‑P hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,904,690	0	C	T	29.3%	67.8 / 28.0	68	Q468Q (CAG→CAA)	treC	trehalose‑6‑P hydrolase
Reads supporting (aligned to +/- strand):  ref base C (28/20);  new base T (13/7);  total (41/27)
Fisher's exact test for biased strand distribution p-value = 7.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.36e-01

AGCACCACAATACAGGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGAATCATCGGCCAGC  >  minE/2904626‑2904750                                                                 |                                                             aGCACCACAATACAGGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCt                                                          >  1:2051733/1‑69 (MQ=255) aGCACCACAATACAGGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCt                                                          >  1:905975/1‑69 (MQ=255)     ccacAATACAGGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACg                                                      <  1:232992/69‑1 (MQ=255)      cacaATACAGGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGc                                                     <  1:2296277/69‑1 (MQ=255)      cacaATACAGGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGc                                                     <  1:2857933/69‑1 (MQ=255)      cacaATACAGGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGc                                                     <  1:1252453/69‑1 (MQ=255)       acaATACAGGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCa                                                    >  1:35044/1‑69 (MQ=255)       acaATACAGGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCa                                                    >  1:1867018/1‑69 (MQ=255)       acaATACAGGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCa                                                    >  1:1011222/1‑69 (MQ=255)       acaATACAGGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCa                                                    >  1:1060936/1‑69 (MQ=255)          aTACAGGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGt                                                  >  1:1778498/1‑68 (MQ=255)          aTACAGGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGt                                                  >  1:743705/1‑68 (MQ=255)             cAGGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTgaa                                               >  1:2115705/1‑66 (MQ=255)                 gCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTaa                                          >  1:740331/1‑69 (MQ=255)                       tGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACtttt                                     >  1:3360435/1‑68 (MQ=255)                       tGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACtttt                                     >  1:3065785/1‑68 (MQ=255)                            gcagATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAgg                               >  1:131365/1‑69 (MQ=255)                            gcagATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAgg                               >  1:424273/1‑69 (MQ=255)                            gcagATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCT‑CTTACGCAGTGCGATTAACTTTTGGTAgg                               >  1:1296103/1‑68 (MQ=255)                            gcagATCCTGGTAATTGCCCCATGTCAGGATGGCTTACTGCTTACGCAGTGCGATTAACTTTTGGTAgg                               >  1:564875/1‑69 (MQ=255)                             cagATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTgg                                   <  1:986331/64‑1 (MQ=255)                             cagATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTgg                                   <  1:2478085/64‑1 (MQ=255)                             cagATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTgg                                   <  1:2353486/64‑1 (MQ=255)                             cagATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTgg                                   <  1:1615826/64‑1 (MQ=255)                             cagATCCTGGTAATTGCCCCATGTCAGGACGGCTTCCTGCTTACGCAGTGCGTTTAACTTTTgg                                   <  1:2948535/64‑1 (MQ=255)                                  ccTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTaa                          >  1:1983799/1‑68 (MQ=255)                                  ccTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTaa                          >  1:2653615/1‑68 (MQ=255)                                   cTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTaaaa                        >  1:2868556/1‑69 (MQ=255)                                        aaTTGCCCCATGTCAGGATGGCTTCCTGCTTACGc                                                     <  1:1978627/35‑1 (MQ=255)                                           cgCCCCATGTCAGGACGGCTTCTTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACAGAgg                <  1:917058/68‑1 (MQ=255)                                           cgCCCCATGTCAGGACGGCTTCTTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACAGAgg                <  1:1825451/68‑1 (MQ=255)                                           cgCCCCATGTCAGGACGGCTTCTTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACAGAgg                <  1:2668238/68‑1 (MQ=255)                                           cgCCCCATGTCAGGACGGCTTCTTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACAGAgg                <  1:901110/68‑1 (MQ=255)                                           cgCCCCATGTCAGGACGGCTTCTTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACAGAgg                <  1:2540014/68‑1 (MQ=255)                                           cgCCCCATGTCAGGACGGCTTCTTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACAGAgg                <  1:1184488/68‑1 (MQ=255)                                           cgCCCCATGTCAGGACGGCTTCTTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACAGAgg                <  1:1580513/68‑1 (MQ=255)                                               ccATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGAatc            <  1:87016/69‑1 (MQ=255)                                               ccATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGAatc            <  1:1863581/69‑1 (MQ=255)                                               ccATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGAatc            <  1:2441048/69‑1 (MQ=255)                                               ccATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGAatc            <  1:1229636/69‑1 (MQ=255)                                               ccATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGAatc            <  1:211636/69‑1 (MQ=255)                                               ccATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGAatc            <  1:426217/69‑1 (MQ=255)                                                    tCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGAATCATCg        >  1:3416036/1‑68 (MQ=255)                                                    tCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGAATCATCg        >  1:1370549/1‑68 (MQ=255)                                                    tCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGAATCATCg        >  1:1194302/1‑68 (MQ=255)                                                    tCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGAATCATCg        >  1:2817619/1‑68 (MQ=255)                                                    tCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGAATCATCg        >  1:2813752/1‑68 (MQ=255)                                                    tCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGAATCATCg        >  1:54364/1‑68 (MQ=255)                                                    tCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGAATCATCg        >  1:2000140/1‑68 (MQ=255)                                                    tCAGGATGGCTTCCTGCTTACGCAGTGCGACTAACTTTTGGTAGGTGTAAAACACCGAGGAATCATCg        >  1:3391686/1‑68 (MQ=255)                                                     cAGGACGGCTTCTTGCTTACGCAGTGCGATTAACTTTTGGTAGgt                              >  1:1146487/1‑45 (MQ=255)                                                     cAGGACGGCTTCTTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGACTCATCg        >  1:683474/1‑67 (MQ=255)                                                     cAGGACGGCTTCTTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGACTCATCg        >  1:195109/1‑67 (MQ=255)                                                     cAGGACGGCTTCTTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGACTCATCg        >  1:1873509/1‑67 (MQ=255)                                                     cAGGACGGCTTCTTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGACTCATCg        >  1:2368626/1‑67 (MQ=255)                                                     cAGGACGGCTTCTTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGACTCATCg        >  1:2636178/1‑67 (MQ=255)                                                     cAGGACGGCTTCTTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGACTCATCg        >  1:2668579/1‑67 (MQ=255)                                                     cAGGACGGCTTCTTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGACTCATCg        >  1:737279/1‑67 (MQ=255)                                                     cAGGACGGCTTCTTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGACTCATCg        >  1:1249579/1‑67 (MQ=255)                                                     cAGGACGGCTTCTTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGACTCATCg        >  1:451272/1‑67 (MQ=255)                                                     cAGGACGGCTTCTTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGACTCATCg        >  1:3055097/1‑67 (MQ=255)                                                     cAGGACGGCTTCTTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGACTCATCg        >  1:3280289/1‑67 (MQ=255)                                                     cAGGACGGCTTCTTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGACTCATCg        >  1:1503914/1‑67 (MQ=255)                                                         aTGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGAATCATCGGCCAgc  >  1:3296692/1‑69 (MQ=255)                                                            gCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGATGAatc            <  1:2699540/56‑1 (MQ=255)                                                            gCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGAatc            <  1:2787574/56‑1 (MQ=255)                                                            gCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGAatc            <  1:2586995/56‑1 (MQ=255)                                                            gCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGAatc            <  1:3426282/56‑1 (MQ=255)                                                                 |                                                             AGCACCACAATACAGGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGAATCATCGGCCAGC  >  minE/2904626‑2904750

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