Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,367,123	C→T	31.8%	P110P (CCG→CCA)	pepQ ←	proline dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,367,123	0	C	T	31.8%	84.5 / 46.0	94	P110P (CCG→CCA)	pepQ	proline dipeptidase
Reads supporting (aligned to +/- strand):  ref base C (50/14);  new base T (16/14);  total (66/28)
Fisher's exact test for biased strand distribution p-value = 2.80e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

ACCGGACCGATATAACCGATATTGCCGCGCGCAGCAGGCAGCAGGCTACCAATGCCATCGGCTTTCGGCAGCGCGATCACTTCCACATCTTCAGTCCAGAAAGAGGTCGGCAGCGGTTCGACGTTGT  >  minE/2367058‑2367184                                                                  |                                                              aCCGGACCGATATAACCGATATTGCCGCGCGCAGCAGGCAGCAGGCTACCAATGCCATCGGCTTTCGGc                                                            >  1:2880708/1‑69 (MQ=255) aCCGGACCGATATAACCGATATTGCCGCGCGCAGCAGGCAGCAGGCTACCAATGCCATCGGCTTTCGGc                                                            >  1:2302900/1‑69 (MQ=255) aCCGGACCGATATAACCGATATTGCCGCGCGCAGCAGGCACCAGGCTACCAATGCCATCGGCTTTCGGc                                                            >  1:3409155/1‑69 (MQ=255)         gATATAACCGATATTGCCGCGCGCAGCAGGCAGCAGGCTACCAATGCCATCGGCTTTCGGCAGCGCGAt                                                    <  1:2663082/69‑1 (MQ=255)         gATATAACCGATATTGCCGCGCGCAGCAGGCAGCAGGCTACCAATGCCATCGGCTTTCGGCAGCGCGAt                                                    <  1:1955001/69‑1 (MQ=255)         gATATAACCGATATTGCCGCGCGCAGCAGGCAGCAGGCTACCAATGCCATCGGCTTTCGGCAGCGCGAt                                                    <  1:3291510/69‑1 (MQ=255)         gATATAACCGATATTGCCGCGCGCAGCAGGCAGCAGGCTACCAATGCCATCGGCTTTCGGCAGCGCGAt                                                    <  1:3296449/69‑1 (MQ=255)         gATATAACCGATATTGCCGCGCGCAGCAGGCAGCAGGCTACCAATGCCATCGGCTTTCGGCAGCGCGAt                                                    <  1:2107497/69‑1 (MQ=255)         gATATAACCGATATTGCCGCGCGCAGCAGGCAGCAGGCTACCAATGCCATCGGCTTTCGGCAGCGCGAt                                                    <  1:532880/69‑1 (MQ=255)          atatAACCGATATTGCCGCGCGCAGCAGGCAGCAGGCTACCAATGCCATCGGCTTTCGGCAGCGCGATc                                                   <  1:861201/69‑1 (MQ=255)          atatAACCGATATTGCCGCGCGCAGCAGGCAGCAGGCTACCAATGCCATCGGCTTTCGGCAGCGCGATc                                                   <  1:1123845/69‑1 (MQ=255)          atatAACCGATATTGCCGCGCGCAGCAGGCAGCAGGCTACCAATGCCATCGGCTTTCGGCAGCGCGATc                                                   <  1:2802706/69‑1 (MQ=255)          atatAACCGATATTGCCGCGCGCAGCAGGCAGCAGGCTACCAATGCCATCGGCTTTCGGCAGCGCGATc                                                   <  1:1400332/69‑1 (MQ=255)          atatAACCGATATTGCCGCGCGCAGCAGGCAGCAGGCTACCAATGCCATCGGCTTTCGGCAGCGCGATc                                                   <  1:748409/69‑1 (MQ=255)          atatAACCGATATTGCCGCGCGCAGCAGGCAGCAGGCTACCAATGCCATCGGCTTTCGGCAGCGCGATc                                                   <  1:1602306/69‑1 (MQ=255)             taACCGATATTGCCGCGCGCAGCAGGCAGCAGGCTACCAATGCCATCGGCTTTCGGCAGCGCGATCACt                                                >  1:2428201/1‑69 (MQ=255)             taACCGATATTGCCGCGCGCAGCAGGCAGCAGGCTACCAATGCCATCGGCTTTCGGCAGCGCGATCACt                                                >  1:2281616/1‑69 (MQ=255)             taACCGATATTGCCGCGCGCAGCAGGCAGCAGGCTACCAATGCCATCGGCTTTCGGCAGCGCGATCACt                                                >  1:3034945/1‑69 (MQ=255)             taACCGATATTGCCGCGCGCAGCAGGCAGCAGGCTACCAATGCCATCGGCTTTCGGCAGCGCGATCACt                                                >  1:140115/1‑69 (MQ=255)             taACCGATATTGCCGCGCGCAGCAGGCAGCAGGCTACCAATGCCATCGGCTTTCGGCAGCGCGATCACt                                                >  1:2004627/1‑69 (MQ=255)             taACCGATATTGCCGCGCGCAGCAGGCAGCAGGCTACCAATGCCATCGGCTTTCGGCAGCGCGATCACt                                                >  1:1547235/1‑69 (MQ=255)              aaCCGATATTGCCTCGCGCAGCAGGCAACAGACTACCAATGCCATCGGCTTTTGGCAGCGCGATCACt                                                >  1:3236146/1‑68 (MQ=255)              aaCCGATATTGCCGCGCGCAGCAGGCAACAGACTACCAATGCCATCGGCTTTTGGCAGCGCGATCAc                                                 >  1:606781/1‑67 (MQ=255)              aaCCGATATTGCCGCGCGCAGCAGGCAACAGACTACCAATGCCATCGGCTTTTGGCAGCGCGATCACt                                                >  1:1027365/1‑68 (MQ=255)              aaCCGATATTGCCGCGCGCAGCAGGCAACAGACTACCAATGCCATCGGCTTTTGGCAGCGCGATCACt                                                >  1:2036154/1‑68 (MQ=255)              aaCCGATATTGCCGCGCGCAGCAGGCAACAGACTACCAATGCCATCGGCTTTTGGCAGCGCGATCACt                                                >  1:2309225/1‑68 (MQ=255)              aaCCGATATTGCCGCGCGCAGCAGGCAACAGACTACCAATGCCATCGGCTTTTGGCAGCGCGATCACt                                                >  1:2162283/1‑68 (MQ=255)              aaCCGATATTGCCGCGCGCAGCAGGCAACAGACTACCAATGCCATCGGCTTTTGGCAGCGCGATCACt                                                >  1:677097/1‑68 (MQ=255)              aaCCGATATTGCCGCGCGCAGCAGGCAACAGACTACCAATGCCATCGGCTTTTGGCAGCGCGATCACt                                                >  1:1813551/1‑68 (MQ=255)              aaCCGATATTGCCGCGCGCAGCAGGCAACAGACTACCAATGCCATCGGCTTTTGGCAGCGCGATCACt                                                >  1:1618798/1‑68 (MQ=255)              aaCCGATATTGCCGCGCGCAGCAGGCAACAGACTACCAATGCCATCGGCTTTTGGCAGCGCGATCACt                                                >  1:1465651/1‑68 (MQ=255)              aaCCGATATTGCCGCGCGCAGCAGGCAACAGACTACCAATGCCATCGGCTTTTGGCAGCGCGATCACt                                                >  1:136058/1‑68 (MQ=255)              aaCCGATATTGCCGCGCGCAGCAGGCAACAGACTACCAATGCCATCGGCTTTTGGCAGCGCGATCACt                                                >  1:1176506/1‑68 (MQ=255)              aaCCGATATTGCCGCGCGCAGCAGGCAACAGACTACCAATGCCATCGGCTTTTGGCAGCGCGATCACt                                                >  1:861017/1‑68 (MQ=255)              aaCCGATATTGCCGCGCGCAGCAGGCAACAGACTACCAATGCCATCGGCTTTTGGCAGCGCGATCACt                                                >  1:3406599/1‑68 (MQ=255)              aaCCGATATTGCCGCGCGCAGCAGGCAACAGACTACCAATGCCATCGGCTTTTGGCAGCGCGATCACt                                                >  1:1065063/1‑68 (MQ=255)              aaCCGATATTGCCGCGCGCAGCAGGCAACAGACTACCAATGCCATCGGCTTTTGGCAGCGCGATCACt                                                >  1:1043066/1‑68 (MQ=255)               aCCGATATTGCCGCGCGCAGCAGGCAGCAGGCTACCAATGCCATCGGCTTTCGGCAGCGCGATc                                                   >  1:2312270/1‑64 (MQ=255)               aCCGATATTGCCGCGCGCAGCAGGCAGCAGGCTACCAATGCCATCGGCTTTCGGCAGCGCGATCACTTc                                              >  1:3109239/1‑69 (MQ=255)               aCCGATATTGCCGCGCGCAGCAGGCAGCAGGCTACCAATGCCATCGGCTTTCGGCAGCGCGATCACTTc                                              >  1:520906/1‑69 (MQ=255)               aCCGATATTGCCGCGCGCAGCAGGCAGCAGGCTACCAATGCCATCGGCTTTCGGCAGCGCGATCACTTc                                              >  1:2225780/1‑69 (MQ=255)               aCCGATATTGCCGCGCGCAGCAGGCAGCAGGCTACCAATGCCATCGGCTTTCGGCAGCGCGATCACTTc                                              >  1:1882941/1‑69 (MQ=255)                ccGATATTGCCGCGCGCAGCAGGCAGCAGGCTACCAATGCCATCGGCTTTCGGCAGCGCGATCACTTc                                              >  1:923766/1‑68 (MQ=255)                ccGATATTGCCGCGCGCAGCAGGCAGCAGGCTACCAATGCCATCGGCTTTCGGCAGCGCGATCACTTc                                              >  1:212863/1‑68 (MQ=255)                ccGATATTGCCGCGCGCAGCAGGCAGCAGGCTACCAATGCCATCGGCTTTCGGCAGCGCGATCACTTc                                              >  1:1898673/1‑68 (MQ=255)                ccGATATTGCCGCGCGCAGCAGGCAGCAGGCTACCAATGCCATCGGCTTTCGGCAGCGCGATCACTTc                                              >  1:3111143/1‑68 (MQ=255)                ccGATATTGCCGCGCGCAGCAGGCAGCAGGCTACCAATGCCATCGGCTTTCGGCAGCGCGATCACTTc                                              >  1:2958927/1‑68 (MQ=255)                        gCCGCGCGCAGCAGGCAGCAGGCTACCAATGCCATCGGCTTTCGGCAGcgcg                                                      >  1:613900/1‑52 (MQ=255)                        gCCGCGCGCAGCAGGCAGCAGGCTACCAATGCCATCGGCTTTCGGCAGCGCGATCAc                                                 >  1:2051034/1‑57 (MQ=255)                        gCCGCGCGCAGCAGGCAGCAGGCTACCAATGCCATCGGCTTTCGGCAGCGCGATCAc                                                 >  1:1113603/1‑57 (MQ=255)                        gCCGCGCGCAGCAGGCAGCAGGCTACCAATGCCATCGGCTTTCGGCAGCGCGATCAc                                                 >  1:1622973/1‑57 (MQ=255)                        gCCGCGCGCAGCAGGCAGCAGGCTACCAAAGCCATCGGCTTTCGGCAGCGCGATctc                                                 >  1:157265/1‑55 (MQ=255)                                     ggcagcaggcTACCAATGCCATCGGCTTTCGGCAGCGCGATCACTTCCACATCTTCAGTCCAGAAagag                        <  1:1209995/69‑1 (MQ=255)                                       cagcaggcTACCAATGCCATCGGCTTTCGGCAGCGCGATCACTTCCACATCTTCAGTCCAGAAAGAGGt                      >  1:3433950/1‑69 (MQ=255)                                       cagcaggcTACCAATGCCATCGGCTTTCGGCAGCGCGATCACTTCCACATCTTCAGTCCAGAAAGAGGt                      >  1:1284697/1‑69 (MQ=255)                                       cagcaggcTACCAATGCCATCGGCTTTCGGCAGCGCGATCACTTCCACATCTTCAGTCCAGAAAGAGGt                      >  1:855512/1‑69 (MQ=255)                                       cagcaggcTACCAATGCCATCGGCTTTCGGCAGCGCGATCACTTCCACATCTTCAGTCCAGAAAGAGGt                      >  1:1783094/1‑69 (MQ=255)                                       cagcaggcTACCAATGCCATCGGCTTTCGGCAGCGCGATCACTTCCACATCTTCAGTCCAGAAAGAGGt                      >  1:2446151/1‑69 (MQ=255)                                       cagcaggcTACCAATGCCATCGGCTTTCGGCAGCGCGATCACTTCCACATCTTCAGTCCAGAAAGAGGt                      >  1:657872/1‑69 (MQ=255)                                        agcaggcTACCAATGCCATCGGCTTTCGGCAGCGCGATCACTTCcac                                           >  1:2737315/1‑47 (MQ=255)                                        agcaggcTACCAATGCCATCGGCTTTCGGCAGCGCGATCACTTCCACATCTTCAGTCCAGAAAGATGTc                     >  1:2908953/1‑69 (MQ=255)                                        agcaggcTACCAATGCCATCGGCTTTCGGCAGCGCGATCACTTCCACATCTTCAGTCCAGAAAGAGGTc                     >  1:763381/1‑69 (MQ=255)                                        agcaggcTACCAATGCCATCGGCTTTCGGCAGCGCGATCACTTCCACATCTTCAGTCCAGAAAGAGGTc                     >  1:2745477/1‑69 (MQ=255)                                        agcaggcTACCAATGCCATCGGCTTTCGGCAGCGCGATCACTTCCACATCTTCAGTCCAGAAAGAGGTc                     >  1:780475/1‑69 (MQ=255)                                        agcaggcTACCAATGCCATCGGCTTTCGGCAGCGCGATCACTTCCACATCTTCAGTCCAGAAAGAGGTc                     >  1:2452058/1‑69 (MQ=255)                                        agcaggcTACCAATGCCATCGGCTTTCGGCAGCGCGATCACTTCCACATCTTCAGTCCAGAAAGAGGTc                     >  1:3405218/1‑69 (MQ=255)                                        agcaggcTACCAATGCCATCGGCTTTCGGCAGCGCGATCACTTCCACATCTTCAGTCCAGAAAGAGGTc                     >  1:1952649/1‑69 (MQ=255)                                           aggcTACCAATGCCATCGGCTTTCGGCAGCGCGATCACTTCCACATCTTCAGTCCAGAAAGAGGTCGGc                  <  1:1098388/69‑1 (MQ=255)                                                    aTGCCATCGGCTTTTGGCAGCGCGATCACTTCTACATCTTCAGTCCAGAAGGAGTTCGGCa                 <  1:540176/61‑1 (MQ=255)                                                    aTGCCATCGGCTTTTGGCAGCGCGATCACTTCTACATCTTCAGTCCAGAAGGAGTTCGGCa                 <  1:990263/61‑1 (MQ=255)                                                    aTGCCATCGGCTTTTGGCAGCGCGATCACTTCTACATCTTCAGTCCAGAAGGAGTTCGGCa                 <  1:623410/61‑1 (MQ=255)                                                    aTGCCATCGGCTTTTGGCAGCGCGATCACTTCTACATCTTCAGTCCAGAAGGAGTTCGGCa                 <  1:687537/61‑1 (MQ=255)                                                    aTGCCATCGGCTTTTGGCAGCGCGATCACTTCTACATCTTCAGTCCAGAAGGAGTTCGGCa                 <  1:665288/61‑1 (MQ=255)                                                    aTGCCATCGGCTTTTGGCAGCGCGATCACTTCTACATCTTCAGTCCAGAAGGAGTTCGGCa                 <  1:2201199/61‑1 (MQ=255)                                                    aTGCCATCGGCTTTTGGCAGCGCGATCACTTCTACATCTTCAGTCCAGAAGGAGTTCGGCa                 <  1:1036097/61‑1 (MQ=255)                                                    aTGCCATCGGCTTTTGGCAGCGCGATCACTTCTACATCTTCAGTCCAGAAGGAGTTCGGCa                 <  1:1105156/61‑1 (MQ=255)                                                    aTGCCATCGGCTTTTGGCAGCGCGATCACTTCTACATCTTCAGTCCAGAAGGAGTTCGGCa                 <  1:1976741/61‑1 (MQ=255)                                                    aTGCCATCGGCTTTTGGCAGCGCGATCACTTCTACATCTTCAGTCCAGAAGGAGTTCGGCa                 <  1:509248/61‑1 (MQ=255)                                                    aTGCCATCGGCTTTTGGCAGCGCGATCACTTCTACATCTTCAGTCCAGAAGGAGTTCGGCa                 <  1:2232057/61‑1 (MQ=255)                                                    aTGCCATCGGCTTTTGGCAGCGCGATCACTTCTACATCTTCAGTCCAGAAGGAGTTCGGCa                 <  1:3071883/61‑1 (MQ=255)                                                    aTGCCATCGGCTTTTGGCAGCGCGATCACTTCTACATCTTCAGTCCAGAAGGAGTTCGGCa                 <  1:3204760/61‑1 (MQ=255)                                                    aTGCCATCGGCTTTTGGCAGCGCGATCACTTCTACATCTTCAGTCCAGAAGGAGTTCGGCa                 <  1:460689/61‑1 (MQ=255)                                                            ggCTTTCGGCAGCGCGATCACTTCCACATCTTCAGTCCAGAAAGAGGTCGGCAGCGGTTCGACGTtgt  >  1:2705022/1‑68 (MQ=255)                                                            ggCTTTCGGCAGCGCGATCACTTCCACATCTTCAGTCCAGAAAGAGGTCGGCAGCGGTTCGACGTtgt  >  1:6684/1‑68 (MQ=255)                                                            ggCTTTCGGCAGCGCGATCACTTCCACATCTTCAGTCCAGAAAGAGGTCGGCAGCGGTTCGACGTtgt  >  1:2164795/1‑68 (MQ=255)                                                            ggCTTTCGGCAGCGCGATCACTTCCACATCTTCAGTCCAGAAAGAGGTCGGCAGCGGTTCGACGTtgt  >  1:581292/1‑68 (MQ=255)                                                            ggCTTTCGGCAGCGCGATCACTTCCACATCTTCAGTCCAGAAAGAGGTCGGCAGCGGTTCGACGTtgt  >  1:1853157/1‑68 (MQ=255)                                                            ggCTTTCGGCAGCGCGATCACTTCCACATCTTCAGTCCAGAAAGAGGTCGGCAGCGGTTCGACGTtgt  >  1:1599723/1‑68 (MQ=255)                                                            ggCTTTCGGCAGCGCGATCACTTCCACATCTTCAGTCCAGAAAGAGGTCGGCAGCGGTTCGACGTtgt  >  1:83524/1‑68 (MQ=255)                                                            ggCTTTCGGCAGCGCGATCACTTCCACATCTTCAGTCCAGAAAGAGGTCGGCAGCGGTTCGACGTtgt  >  1:1366143/1‑68 (MQ=255)                                                            ggCTTTCGGCAGCGCGATCACTTCCACATCTTCAGTCCAGAAAGAGGTCGGCAGCGGTTCGACGTtgt  >  1:1114092/1‑68 (MQ=255)                                                            ggCTTTCGGCAGCGCGATCACTTCCACATCTTCAGTCCAGAAAGAGGTCGGCAGCGGTTCGACGTtgt  >  1:408110/1‑68 (MQ=255)                                                            ggCTTTCGGCAGCGCGATCACTTCCACATCTTCAGTCCAGAAAGAGGTCGGCAGCGGTTCGACGTtgt  >  1:1092055/1‑68 (MQ=255)                                                            ggCTTTCGGCAGCGCGATCACTTCCACATCTTCAGTCCAGAAAGAGGTCGGCAGCGGGTCGACGt     >  1:2348645/1‑65 (MQ=255)                                                                  |                                                              ACCGGACCGATATAACCGATATTGCCGCGCGCAGCAGGCAGCAGGCTACCAATGCCATCGGCTTTCGGCAGCGCGATCACTTCCACATCTTCAGTCCAGAAAGAGGTCGGCAGCGGTTCGACGTTGT  >  minE/2367058‑2367184

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