Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,157,473	G→T	64.3%	A328A (GCC→GCA)	astA ←	arginine succinyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,157,473	0	G	T	64.3%	33.2 / 45.4	56	A328A (GCC→GCA)	astA	arginine succinyltransferase
Reads supporting (aligned to +/- strand):  ref base G (16/4);  new base T (10/26);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 2.40e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGG  >  minE/1157417‑1157537                                                         |                                                                 tAAAGTGATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCTGCGTGGCAttt                                                       <  1:356580/68‑1 (MQ=255) tAAAGTCATGTTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCTGCGTGGCAttt                                                       <  1:2045234/68‑1 (MQ=255) tAAAGTCATGCTGTTTTCTCCTCTGCGCTCAGGCGCACCAGACGAACGCGATCCCCTGCGTGGCAttt                                                       <  1:2624593/68‑1 (MQ=255) tAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAGCGCGATCCCCGGCGTGGCATTTg                                                      >  1:478597/1‑69 (MQ=255) tAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCTGCGTGGTAttt                                                       <  1:3265016/68‑1 (MQ=255) tAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCTGCGTGGCAttt                                                       <  1:3239964/68‑1 (MQ=255) tAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCTGCGTGGCAttt                                                       <  1:2191596/68‑1 (MQ=255) tAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCTGCGTGGCAttt                                                       <  1:237683/68‑1 (MQ=255) tAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCTGCGTGGCAttt                                                       <  1:2879339/68‑1 (MQ=255) tAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCTGCGTGGCAttt                                                       <  1:2908834/68‑1 (MQ=255) tAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCTGCGTGGCAttt                                                       <  1:3032131/68‑1 (MQ=255) tAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCTGCGTGGCAttt                                                       <  1:3196947/68‑1 (MQ=255) tAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCTGCGTGGCAttt                                                       <  1:2614296/68‑1 (MQ=255) tAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCTGCGTGGCAttt                                                       <  1:3290100/68‑1 (MQ=255) tAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCTGCGTGGCAttt                                                       <  1:724531/68‑1 (MQ=255) tAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCTGCGTGGCAttt                                                       <  1:930614/68‑1 (MQ=255) tAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCTGCGTGGCAttt                                                       <  1:954928/68‑1 (MQ=255) tAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCTGCGTGGCAttt                                                       <  1:2036183/68‑1 (MQ=255) tAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCTGCGTGGCAttt                                                       <  1:1091132/68‑1 (MQ=255) tAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCTGCGTGGCAttt                                                       <  1:1479423/68‑1 (MQ=255) tAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCTGCGTGGCAttt                                                       <  1:1969525/68‑1 (MQ=255) tAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCTGCGTGGCAttt                                                       <  1:1982042/68‑1 (MQ=255) tAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCTGCGTGGCAttt                                                       <  1:2174175/68‑1 (MQ=255) tAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTg                                                      >  1:178293/1‑69 (MQ=255) tAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTg                                                      >  1:3312945/1‑69 (MQ=255) tAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTg                                                      >  1:1227510/1‑69 (MQ=255) tAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTg                                                      >  1:2240970/1‑69 (MQ=255) tAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTg                                                      >  1:3123438/1‑69 (MQ=255) tAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTg                                                      >  1:835688/1‑69 (MQ=255) tAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTg                                                      >  1:2884089/1‑69 (MQ=255) tAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTg                                                      >  1:2442069/1‑69 (MQ=255)  aaaGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCTGCGTGGCAttt                                                       <  1:2274314/67‑1 (MQ=255)  aaaGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCTGCGTGGCAttt                                                       <  1:540475/67‑1 (MQ=255)   aaGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTTAg                                                    <  1:2536449/69‑1 (MQ=255)   aaGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAg                                                    <  1:567567/69‑1 (MQ=255)   aaGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAg                                                    >  1:541908/1‑69 (MQ=255)   aaGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAg                                                    >  1:457720/1‑69 (MQ=255)   aaGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAg                                                    >  1:3311987/1‑69 (MQ=255)   aaGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAg                                                    <  1:2294984/69‑1 (MQ=255)   aaGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACAAGACGAACGCGATCCCCGGCGTGGCATTTGAg                                                    >  1:2536439/1‑69 (MQ=255)      tCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCTGCGTGGCAttt                                                       <  1:1068686/63‑1 (MQ=255)        aTGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCTGCGTGGCATTTGAGGGCa                                                >  1:1124495/1‑68 (MQ=255)        aTGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCTGCGTGGCATTTGAGGGCa                                                >  1:3278509/1‑68 (MQ=255)        aTGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCTGCGTGGCATTTGAGGGCa                                                >  1:41991/1‑68 (MQ=255)        aTGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCTGCGTGGCATTTGAGGGCa                                                >  1:1404043/1‑68 (MQ=255)        aTGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCTGCGTGGCATTTGAGGGCa                                                >  1:2340843/1‑68 (MQ=255)        aTGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCTGCGTGGCATTTGAGGGCa                                                >  1:1577187/1‑68 (MQ=255)        aTGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCTGCGTGGCATTTGAGGGCa                                                >  1:1611620/1‑68 (MQ=255)        aTGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCTGCGTGGCATTTGAGGGCa                                                >  1:1909068/1‑68 (MQ=255)        aTGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCTGCGTGGCATTTGAGGGCa                                                >  1:2371967/1‑68 (MQ=255)        aTGCTGTTTTCTCCTCTGCGCACAGGCGCAACAGACGAACGCGATCCCCTGCGTGGCATTTGAGGGc                                                 >  1:3340388/1‑67 (MQ=255)           cTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCa                                            >  1:44908/1‑69 (MQ=255)           cTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCa                                            >  1:3053823/1‑69 (MQ=255)                ttCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTtgt                                       >  1:1364296/1‑69 (MQ=255)                        tGCGCACAGGCGCACCAGACGAACGCGATCCCCTGCGTGGCAttt                                                       <  1:1835690/45‑1 (MQ=255)                                                     ccccGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCgg  <  1:1218562/69‑1 (MQ=255)                                                         |                                                                 TAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGG  >  minE/1157417‑1157537

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