Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,285,088	T→C	32.9%	V150V (GTA→GTG)	argH ←	argininosuccinate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,285,088	0	T	C	32.9%	74.5 / 49.2	94	V150V (GTA→GTG)	argH	argininosuccinate lyase
Reads supporting (aligned to +/- strand):  ref base T (24/39);  new base C (16/15);  total (40/54)
Fisher's exact test for biased strand distribution p-value = 2.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

CAGGCACCAGTGCGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGCCGATTGCAGCTGCCGGTTAGCCGTCAGT  >  minE/2285028‑2285154                                                             |                                                                   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:1535560/1‑69 (MQ=255)             cgcgAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTtg                                                >  1:1140033/1‑69 (MQ=255)             cgcgAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTtg                                                >  1:923778/1‑69 (MQ=255)             cgcgAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTtg                                                >  1:3321159/1‑69 (MQ=255)             cgcgAACGTCACAGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTtg                                                >  1:2991061/1‑69 (MQ=255)                 aaCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTTTTTTGTGCg                                            <  1:2711000/69‑1 (MQ=255)                 aaCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCg                                            <  1:2208342/69‑1 (MQ=255)                 aaCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCg                                            <  1:1368591/69‑1 (MQ=255)                 aaCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCg                                            <  1:2740772/69‑1 (MQ=255)                 aaCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCg                                            <  1:1533822/69‑1 (MQ=255)                 aaCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCg                                            <  1:1061800/69‑1 (MQ=255)                 aaCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCg                                            <  1:3383346/69‑1 (MQ=255)                 aaCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCg                                            <  1:2381699/69‑1 (MQ=255)                 aaCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCg                                            <  1:2270039/69‑1 (MQ=255)                 aaCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCg                                            <  1:50722/69‑1 (MQ=255)                 aaCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCg                                            <  1:1605184/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:1615582/69‑2 (MQ=255)                                      tGCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                       <  1:2073026/69‑2 (MQ=255)                                      tGCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                       <  1:2285083/69‑2 (MQ=255)                                      tGCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                       <  1:330654/69‑2 (MQ=255)                                      tGCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                       <  1:2044930/69‑2 (MQ=255)                                      tGCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                       <  1:3278141/69‑2 (MQ=255)                                      tGCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                       <  1:3230560/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:578272/68‑2 (MQ=255)                                       gCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                       <  1:711793/68‑2 (MQ=255)                                       gCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                       <  1:1413272/68‑2 (MQ=255)                                       gCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                       <  1:1469214/68‑2 (MQ=255)                                       gCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                       <  1:2581416/68‑2 (MQ=255)                                       gCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                       <  1:1355052/68‑2 (MQ=255)                                             gAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGCCGATTGCAGCTGcc                >  1:1221619/1‑69 (MQ=255)                                             gAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGCCGATTGCAGCTGcc                >  1:3129575/1‑69 (MQ=255)                                             gAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGCCGATTGCAGCTGcc                >  1:1104946/1‑69 (MQ=255)                                                  tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                             <  1:1318674/51‑1 (MQ=255)                                                  tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                             <  1:1902326/51‑1 (MQ=255)                                                  tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                             <  1:1401683/51‑1 (MQ=255)                                                  tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                             <  1:455884/51‑1 (MQ=255)                                                  tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                             <  1:42006/51‑1 (MQ=255)                                                  tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                             <  1:33434/51‑1 (MQ=255)                                                  tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                             <  1:1150550/51‑1 (MQ=255)                                                  tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                             <  1:2981337/51‑1 (MQ=255)                                                  tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                             <  1:2449793/51‑1 (MQ=255)                                                  tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                             <  1:3299290/51‑1 (MQ=255)                                                  tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                             <  1:1622076/51‑1 (MQ=255)                                                  tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                             <  1:3273357/51‑1 (MQ=255)                                                  tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                             <  1:2959996/51‑1 (MQ=255)                                                    aCCTGTCATTACAGCGTCCTGATTGTTTTGTTCGGTTTCCACCAGCGCCGATTGCAGCTGCCGGTTAGc         >  1:3431672/1‑69 (MQ=255)                                                           aTTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGCCGATTGCAGCTGCCGGTTAGCCGTCAGt  <  1:182526/69‑1 (MQ=255)                                                           aTTACCGAGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGCCGATTGCAGCTGCCGGTTAGCCGTCAGt  <  1:2563264/69‑1 (MQ=255)                                                             |                                                                   CAGGCACCAGTGCGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGCCGATTGCAGCTGCCGGTTAGCCGTCAGT  >  minE/2285028‑2285154

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