Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	874,871	T→C	37.3%	A29A (GCT→GCC)	acnA →	aconitate hydratase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	874,871	0	T	C	37.3%	39.9 / 47.5	75	A29A (GCT→GCC)	acnA	aconitate hydratase 1
Reads supporting (aligned to +/- strand):  ref base T (28/19);  new base C (14/14);  total (42/33)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGCTCGAAAACC  >  minE/874808‑874929                                                                |                                                           cAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATc                                                       <  1:1111319/69‑1 (MQ=255) cAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATc                                                       <  1:2153597/69‑1 (MQ=255) cAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATc                                                       <  1:1356457/69‑1 (MQ=255) cAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATc                                                       <  1:2228367/69‑1 (MQ=255) cAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATc                                                       <  1:1266785/69‑1 (MQ=255) cAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATc                                                       <  1:3433129/69‑1 (MQ=255) cAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATc                                                       <  1:3285907/69‑1 (MQ=255) cAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATc                                                       <  1:267164/69‑1 (MQ=255) cAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATc                                                       <  1:268574/69‑1 (MQ=255) cAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATc                                                       <  1:851306/69‑1 (MQ=255) cAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATc                                                       <  1:2808874/69‑1 (MQ=255) cAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATc                                                       <  1:1090310/69‑1 (MQ=255) cAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATc                                                       <  1:2884213/69‑1 (MQ=255) cAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATc                                                       <  1:3054151/69‑1 (MQ=255) cAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATc                                                       <  1:1607271/69‑1 (MQ=255)    tAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACt                                                    >  1:1580006/1‑69 (MQ=255)      aagACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCCAAATCACTgg                                                  >  1:1648174/3‑69 (MQ=255)      aagACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCCAAATCACTgg                                                  >  1:1016194/3‑69 (MQ=255)      aagACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCCAAATCACTgg                                                  >  1:2128274/3‑69 (MQ=255)      aagACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCCAAATCACTgg                                                  >  1:2371878/3‑69 (MQ=255)      aagACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCCAAATCACTgg                                                  >  1:2722034/3‑69 (MQ=255)      aagACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCCAAATCACTgg                                                  >  1:3192616/3‑69 (MQ=255)      aagACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCCAAATCACTgg                                                  >  1:3147830/3‑69 (MQ=255)      aagACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCCAAATCACTgg                                                  >  1:1270289/3‑69 (MQ=255)      aagACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCCAAATCACTgg                                                  >  1:1045856/3‑69 (MQ=255)      aagACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCCAAATCACTgg                                                  >  1:1266219/3‑69 (MQ=255)      aagACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCCAAATCACTgg                                                  >  1:164131/3‑69 (MQ=255)      aagACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCCAAATCACTgg                                                  >  1:1095203/3‑69 (MQ=255)      aagACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCCAAATCACTg                                                   >  1:3108414/3‑68 (MQ=255)      aagACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCCAAATCACTg                                                   >  1:3254693/3‑68 (MQ=255)        gACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGc                                                <  1:523012/69‑1 (MQ=255)         acacGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGc                                                <  1:2290602/68‑1 (MQ=255)           acGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATc                                                       <  1:2013301/59‑1 (MQ=255)               tGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCa                                         <  1:746859/69‑1 (MQ=255)                gCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCAc                                        >  1:721705/1‑69 (MQ=255)                gCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCAc                                        >  1:2249877/1‑69 (MQ=255)                gCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCAc                                        >  1:1351522/1‑69 (MQ=255)                gCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCAc                                        >  1:860175/1‑69 (MQ=255)                gCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCAc                                        >  1:109102/1‑69 (MQ=255)                      cAAATATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCCAAATCACTGGGCGATATCACCCGTCt                                  <  1:2028190/69‑1 (MQ=255)                      cAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCCAAATCACTGGGCGATATCACCCGTCt                                  <  1:1240167/69‑1 (MQ=255)                      cAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCCAAATCACTGGGCGATATCACCCGTCt                                  <  1:2495516/69‑1 (MQ=255)                      cAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCCAAATCACTGGGCGATATCACCCGTCt                                  <  1:2658436/69‑1 (MQ=255)                      cAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCCAAATCACTGGGCGATATCACCCGTCt                                  <  1:10852/69‑1 (MQ=255)                      cAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCCAAATCACTGGGCGATATCACCCGTCt                                  <  1:317678/69‑1 (MQ=255)                      cAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCCAAATCACTGGGCGATATCACCCGTCt                                  <  1:2324520/69‑1 (MQ=255)                      cAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCCAAATCACTGGGCGATATCACCCGTCt                                  <  1:2540437/69‑1 (MQ=255)                      cAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCCAAATCACTGGGCGATATCACCCGTCt                                  <  1:2194270/69‑1 (MQ=255)                      cAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCCAAATCACTGGGCGATATCACCCGTCt                                  <  1:1430511/69‑1 (MQ=255)                       aaaGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCCAACTCACTGGGCGATATCACCCGTCt                                  <  1:106515/68‑1 (MQ=255)                       aaaGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCCAAATCACTGGGCGATATCACCCGTCt                                  <  1:2936129/68‑1 (MQ=255)                                         tactacAGCCTGCCTCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGtt               >  1:1956159/1‑69 (MQ=255)                                         tactacAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGtt               >  1:1645841/1‑69 (MQ=255)                                         tactacAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGtt               >  1:1623748/1‑69 (MQ=255)                                         tactacAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGtt               >  1:308143/1‑69 (MQ=255)                                         tactacAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGtt               >  1:856597/1‑69 (MQ=255)                                         tactacAGCCTGCCGCTTGCTGCTAAATAACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGtt               >  1:1893572/1‑69 (MQ=255)                                          actacAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGttt              >  1:2696473/1‑69 (MQ=255)                                          actacAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGttt              >  1:2692687/1‑69 (MQ=255)                                          actacAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGttt              >  1:224492/1‑69 (MQ=255)                                          actacAGCCTGCCGCTTGCTGCCAAATCACTGGGCGATATCACCCGTCt                                  <  1:1158615/49‑1 (MQ=255)                                           ctacAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGtttt             >  1:1039927/1‑69 (MQ=255)                                           ctacAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGtttt             >  1:3071539/1‑69 (MQ=255)                                           ctacAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGtttt             >  1:2528694/1‑69 (MQ=255)                                              cAGCCTGCCGCTTGCTGCTAAATCACTGGGCGatat                                           >  1:1452638/1‑36 (MQ=255)                                              cAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTg            >  1:1096253/1‑67 (MQ=255)                                              cAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGc           >  1:322847/1‑68 (MQ=255)                                              cAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGc           >  1:2673246/1‑68 (MQ=255)                                              cAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGc           >  1:2283284/1‑68 (MQ=255)                                              cAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGc           >  1:2141598/1‑68 (MQ=255)                                              cAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGc           >  1:1954624/1‑68 (MQ=255)                                              cAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGc           >  1:1598379/1‑68 (MQ=255)                                                  cTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGCTCGaa      >  1:2556311/1‑69 (MQ=255)                                                   tGCCGCTTGCTGCCAAATCACTGGGCGATATCACCCGTCt                                  <  1:448467/40‑1 (MQ=255)                                                      cGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGCTCGAAAAcc  >  1:895838/1‑69 (MQ=255)                                                                |                                                           CAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGCTCGAAAACC  >  minE/874808‑874929

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