Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,285,082	T→C	33.6%	P152P (CCA→CCG)	argH ←	argininosuccinate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,285,082	0	T	C	33.6%	67.9 / 49.8	92	P152P (CCA→CCG)	argH	argininosuccinate lyase
Reads supporting (aligned to +/- strand):  ref base T (24/37);  new base C (16/15);  total (40/52)
Fisher's exact test for biased strand distribution p-value = 2.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGCACCAGTGCGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGCCGATTGCAGCTGCCGGTTAGC  >  minE/2285028‑2285147                                                       |                                                                  cAGGCACCAGTGCGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTc                                                     <  1:62728/69‑1 (MQ=255) cAGGCACCAGTGCGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTc                                                     <  1:1097154/69‑1 (MQ=255) cAGGCACCAGTGCGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTc                                                     <  1:14965/69‑1 (MQ=255) cAGGCACCAGTGCGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTc                                                     <  1:1142129/69‑1 (MQ=255) cAGGCACCAGTGCGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTc                                                     <  1:2308348/69‑1 (MQ=255) cAGGCACCAGTGCGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTc                                                     <  1:1980164/69‑1 (MQ=255) cAGGCACCAGTGCGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTc                                                     <  1:2214221/69‑1 (MQ=255) cAGGCACCAGTGCGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTc                                                     <  1:2118923/69‑1 (MQ=255)  aGGCACCAGTGCGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTc                                                     <  1:2998934/68‑1 (MQ=255)   ggCACCAGTGAGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCCGGCATCACCGCGTCCt                                                   >  1:1693620/1‑69 (MQ=255)   ggCACCAGTGAGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCCGGCATCACCGCGTCCt                                                   >  1:97007/1‑69 (MQ=255)   ggCACCAGTGAGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCCGGCATCACCGCGTCCt                                                   >  1:1727930/1‑69 (MQ=255)   ggCACCAGTGAGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCCGGCATCACCGCGTCCt                                                   >  1:1949567/1‑69 (MQ=255)   ggCACCAGTGAGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCCGGCATCACCGCGTCCt                                                   >  1:1953782/1‑69 (MQ=255)   ggCACCAGTGAGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCCGGCATCACCGCGTCCt                                                   >  1:2035367/1‑69 (MQ=255)   ggCACCAGTGAGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCCGGCATCACCGCGTCCt                                                   >  1:2060026/1‑69 (MQ=255)   ggCACCAGTGAGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCCGGCATCACCGCGTCCt                                                   >  1:2269145/1‑69 (MQ=255)   ggCACCAGTGAGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCCGGCATCACCGCGTCCt                                                   >  1:2399520/1‑69 (MQ=255)   ggCACCAGTGAGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCCGGCATCACCGCGTCCt                                                   >  1:2498868/1‑69 (MQ=255)   ggCACCAGTGAGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCCGGCATCACCGCGTCCt                                                   >  1:2569526/1‑69 (MQ=255)   ggCACCAGTGAGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCCGGCATCACCGCGTCCt                                                   >  1:2584852/1‑69 (MQ=255)   ggCACCAGTGAGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCCGGCATCACCGCGTCCt                                                   >  1:3320688/1‑69 (MQ=255)   ggCACCAGTGAGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCCGGCATCACCGCGTCCt                                                   >  1:1315801/1‑69 (MQ=255)   ggCACCAGTGAGCGAACGTCACAGGCTGGGCGCGTTGCAGGTGAGTGTAACCCGGCATCACCGCGTCCt                                                   >  1:3283648/1‑69 (MQ=255)   ggCACCAGTGAGAGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCCGGGATCACCGCGTCCt                                                   >  1:1277132/1‑69 (MQ=255)      aCCAGTGCGCGAACGTCCCCGGCTGGGCGCGTTGCTGGTGAGTGTACCCTGGCATTACCGCGTc                                                     >  1:2058512/1‑64 (MQ=255)      aCCAGTGCGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGAt                                                >  1:1581728/1‑69 (MQ=255)      aCCAGTGCGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGAt                                                >  1:1116491/1‑69 (MQ=255)      aCCAGTGCGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGAt                                                >  1:2411294/1‑69 (MQ=255)      aCCAGTGCGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGAt                                                >  1:1779627/1‑69 (MQ=255)         aGTGCGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTc                                                     <  1:198879/61‑1 (MQ=255)            gcgcGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGtttt                                          >  1:101529/1‑69 (MQ=255)            gcgcGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGtttt                                          >  1:105507/1‑69 (MQ=255)            gcgcGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGtttt                                          >  1:361989/1‑69 (MQ=255)            gcgcGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGtttt                                          >  1:2816863/1‑69 (MQ=255)            gcgcGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGtttt                                          >  1:2812569/1‑69 (MQ=255)            gcgcGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGtttt                                          >  1:1140845/1‑69 (MQ=255)            gcgcGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGtttt                                          >  1:1923777/1‑69 (MQ=255)             cgcgAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTtg                                         >  1:1570225/1‑69 (MQ=255)             cgcgAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTtg                                         >  1:3234703/1‑69 (MQ=255)             cgcgAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTtg                                         >  1:923778/1‑69 (MQ=255)             cgcgAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTtg                                         >  1:3321159/1‑69 (MQ=255)             cgcgAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTtg                                         >  1:1535560/1‑69 (MQ=255)             cgcgAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTtg                                         >  1:1140033/1‑69 (MQ=255)             cgcgAACGTCACAGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTtg                                         >  1:2991061/1‑69 (MQ=255)                 aaCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTTTTTTGTGCg                                     <  1:2711000/69‑1 (MQ=255)                 aaCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCg                                     <  1:50722/69‑1 (MQ=255)                 aaCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCg                                     <  1:1061800/69‑1 (MQ=255)                 aaCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCg                                     <  1:3383346/69‑1 (MQ=255)                 aaCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCg                                     <  1:2740772/69‑1 (MQ=255)                 aaCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCg                                     <  1:1533822/69‑1 (MQ=255)                 aaCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCg                                     <  1:2381699/69‑1 (MQ=255)                 aaCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCg                                     <  1:1368591/69‑1 (MQ=255)                 aaCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCg                                     <  1:1605184/69‑1 (MQ=255)                 aaCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCg                                     <  1:2270039/69‑1 (MQ=255)                 aaCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCg                                     <  1:2208342/69‑1 (MQ=255)                  aCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCg                                     <  1:1708219/68‑1 (MQ=255)                  aCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCg                                     <  1:691397/68‑1 (MQ=255)                  aCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCg                                     <  1:951235/68‑1 (MQ=255)                                      tGCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                <  1:3278141/69‑2 (MQ=255)                                      tGCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                <  1:3230560/69‑2 (MQ=255)                                      tGCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                <  1:330654/69‑2 (MQ=255)                                      tGCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                <  1:2044930/69‑2 (MQ=255)                                      tGCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                <  1:2073026/69‑2 (MQ=255)                                      tGCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                <  1:1615582/69‑2 (MQ=255)                                      tGCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                <  1:2285083/69‑2 (MQ=255)                                       gCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGCCGATTGCa               >  1:2875978/1‑69 (MQ=255)                                       gCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCTGTCTCGACCAGAGCGGATTgt                <  1:1679960/68‑2 (MQ=255)                                       gCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGAtttt                <  1:1737729/68‑3 (MQ=255)                                       gCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                <  1:711793/68‑2 (MQ=255)                                       gCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                <  1:1413272/68‑2 (MQ=255)                                       gCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                <  1:578272/68‑2 (MQ=255)                                       gCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                <  1:1469214/68‑2 (MQ=255)                                       gCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                <  1:1355052/68‑2 (MQ=255)                                       gCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                <  1:2581416/68‑2 (MQ=255)                                             gAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGCCGATTGCAGCTGcc         >  1:3129575/1‑69 (MQ=255)                                             gAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGCCGATTGCAGCTGcc         >  1:1104946/1‑69 (MQ=255)                                             gAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGCCGATTGCAGCTGcc         >  1:1221619/1‑69 (MQ=255)                                                  tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                      <  1:42006/51‑1 (MQ=255)                                                  tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                      <  1:1401683/51‑1 (MQ=255)                                                  tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                      <  1:1318674/51‑1 (MQ=255)                                                  tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                      <  1:1902326/51‑1 (MQ=255)                                                  tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                      <  1:1150550/51‑1 (MQ=255)                                                  tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                      <  1:2449793/51‑1 (MQ=255)                                                  tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                      <  1:455884/51‑1 (MQ=255)                                                  tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                      <  1:3273357/51‑1 (MQ=255)                                                  tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                      <  1:2959996/51‑1 (MQ=255)                                                  tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                      <  1:33434/51‑1 (MQ=255)                                                  tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                      <  1:2981337/51‑1 (MQ=255)                                                  tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                      <  1:3299290/51‑1 (MQ=255)                                                  tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                      <  1:1622076/51‑1 (MQ=255)                                                    aCCTGTCATTACAGCGTCCTGATTGTTTTGTTCGGTTTCCACCAGCGCCGATTGCAGCTGCCGGTTAGc  >  1:3431672/1‑69 (MQ=255)                                                       |                                                                  CAGGCACCAGTGCGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGCCGATTGCAGCTGCCGGTTAGC  >  minE/2285028‑2285147

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