Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,330,311	A→G	61.9%	N425N (AAT→AAC)	cpsG ←	phosphomannomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,330,311	0	A	G	61.9%	55.7 / 62.5	97	N425N (AAT→AAC)	cpsG	phosphomannomutase
Reads supporting (aligned to +/- strand):  ref base A (10/27);  new base G (25/35);  total (35/62)
Fisher's exact test for biased strand distribution p-value = 1.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTCGGCAAAGGTCATGCTGATGC  >  minE/1330248‑1330360                                                                |                                                  tcgcgCTTCCATCAGCGGCAGATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                                              <  1:3253008/69‑1 (MQ=255) tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTGGGTGTTGGa                                              <  1:2711020/69‑1 (MQ=255) tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                                              <  1:1890067/69‑1 (MQ=255) tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                                              <  1:2205092/69‑1 (MQ=255) tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                                              <  1:787765/69‑1 (MQ=255) tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                                              <  1:2073922/69‑1 (MQ=255) tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                                              <  1:2095422/69‑1 (MQ=255) tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                                              <  1:2173541/69‑1 (MQ=255) tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                                              <  1:2182274/69‑1 (MQ=255) tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                                              <  1:2186149/69‑1 (MQ=255) tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                                              <  1:227347/69‑1 (MQ=255) tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                                              <  1:2280865/69‑1 (MQ=255) tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                                              <  1:757745/69‑1 (MQ=255) tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                                              <  1:2682390/69‑1 (MQ=255) tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                                              <  1:2697082/69‑1 (MQ=255) tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                                              <  1:2728727/69‑1 (MQ=255) tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                                              <  1:693399/69‑1 (MQ=255) tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                                              <  1:2798258/69‑1 (MQ=255) tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                                              <  1:2820640/69‑1 (MQ=255) tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                                              <  1:2953898/69‑1 (MQ=255) tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                                              <  1:2967059/69‑1 (MQ=255) tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                                              <  1:3206035/69‑1 (MQ=255) tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                                              <  1:663409/69‑1 (MQ=255) tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                                              <  1:3295237/69‑1 (MQ=255) tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                                              <  1:583703/69‑1 (MQ=255) tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                                              <  1:1267697/69‑1 (MQ=255) tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                                              <  1:1785917/69‑1 (MQ=255) tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                                              <  1:816503/69‑1 (MQ=255) tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                                              <  1:1393830/69‑1 (MQ=255) tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                                              <  1:1604972/69‑1 (MQ=255) tcgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACAGGTTCGGTGTTGGa                                              <  1:3090356/69‑1 (MQ=255)  cgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                                              <  1:3364509/68‑1 (MQ=255)  cgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                                              <  1:2184501/68‑1 (MQ=255)  cgcgCTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                                              <  1:2436783/68‑1 (MQ=255)   gcgcTTCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGa                                              <  1:1734240/67‑1 (MQ=255)        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTTGAggag                                          >  1:2504417/1‑64 (MQ=255)        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                                          >  1:2558524/1‑64 (MQ=255)        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                                          >  1:2634682/1‑64 (MQ=255)        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                                          >  1:1587933/1‑64 (MQ=255)        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                                          >  1:790005/1‑64 (MQ=255)        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                                          >  1:1562179/1‑64 (MQ=255)        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                                          >  1:2707665/1‑64 (MQ=255)        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                                          >  1:1547909/1‑64 (MQ=255)        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                                          >  1:2803474/1‑64 (MQ=255)        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                                          >  1:1377872/1‑64 (MQ=255)        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                                          >  1:2979991/1‑64 (MQ=255)        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                                          >  1:1289451/1‑64 (MQ=255)        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                                          >  1:596729/1‑64 (MQ=255)        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                                          >  1:979566/1‑64 (MQ=255)        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                                          >  1:1078622/1‑64 (MQ=255)        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                                          >  1:1834874/1‑64 (MQ=255)        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                                          >  1:1942829/1‑64 (MQ=255)        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                                          >  1:1991307/1‑64 (MQ=255)        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                                          >  1:2010533/1‑64 (MQ=255)        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                                          >  1:803040/1‑64 (MQ=255)        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                                          >  1:2276657/1‑64 (MQ=255)        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                                          >  1:760706/1‑64 (MQ=255)        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                                          >  1:2465076/1‑64 (MQ=255)        tCCATCAGCGGCACATCACCGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAGGAgcg                                        >  1:2480653/1‑68 (MQ=255)        tCCATCAGCGGCACATCACAGCGTGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAggag                                          >  1:3002623/1‑64 (MQ=255)                 ggCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAg                              >  1:2710097/1‑69 (MQ=255)                    acaTCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGa                                              <  1:2835833/50‑1 (MQ=255)                    acaTCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGa                                              <  1:1614442/50‑1 (MQ=255)                    acaTCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAgcgc                           <  1:2551808/69‑1 (MQ=255)                    acaTCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAgcgc                           <  1:824081/69‑1 (MQ=255)                    acaTCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAgcgc                           <  1:3010089/69‑1 (MQ=255)                    acaTCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAgcgc                           <  1:1483328/69‑1 (MQ=255)                     caTCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAgcgc                           <  1:2003240/68‑1 (MQ=255)                         aCCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTc                      >  1:665537/1‑69 (MQ=255)                         aCCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTc                      >  1:391454/1‑69 (MQ=255)                         aCCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTc                      >  1:1156901/1‑69 (MQ=255)                         aCCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTc                      >  1:950078/1‑69 (MQ=255)                         aCCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTc                      >  1:1764510/1‑69 (MQ=255)                         aCCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTc                      >  1:1560126/1‑69 (MQ=255)                         aCCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTc                      >  1:831511/1‑69 (MQ=255)                                  ttCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTCGGCAAAGGt             <  1:1443364/69‑1 (MQ=255)                                   tCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTCGGCAAAGGTc            <  1:1536158/69‑1 (MQ=255)                                    ccACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTCGGCAAAGGt             <  1:2517321/67‑1 (MQ=255)                                    ccACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTCGGCAAAGGt             <  1:769934/67‑1 (MQ=255)                                    ccACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTCGGCAAAGGt             <  1:2045885/67‑1 (MQ=255)                                    ccACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTCGGCAAAGGt             <  1:1096843/67‑1 (MQ=255)                                    ccACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTCGGCAAAGGt             <  1:1858743/67‑1 (MQ=255)                                    ccACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTCGGCAAAGGt             <  1:2175386/67‑1 (MQ=255)                                    ccACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTCGGCAAAGGt             <  1:2511927/67‑1 (MQ=255)                                    ccACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTCGGCAAAGGt             <  1:2562586/67‑1 (MQ=255)                                    ccACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTCGGCAAAGGt             <  1:378901/67‑1 (MQ=255)                                           cAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTCGGCAAAGGTCATGc        >  1:1696470/1‑65 (MQ=255)                                           cAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTCGGCAAAGGTCATGCTGAt    >  1:616051/1‑69 (MQ=255)                                                     accGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTCGGCAAAGGTCATGCTGATGc  <  1:2789054/61‑1 (MQ=255)                                                     accGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTCGGCAAAGGTCATGCTGATGc  <  1:3178721/61‑1 (MQ=255)                                                     accGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTCGGCAAAGGTCATGCTGATGc  <  1:3216337/61‑1 (MQ=255)                                                     accGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTCGGCAAAGGTCATGCTGATGc  <  1:550711/61‑1 (MQ=255)                                                     accGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTCGGCAAAGGTCATGCTGATGc  <  1:1755059/61‑1 (MQ=255)                                                     accGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTCGGCAAAGGTCATGCTGATGc  <  1:1612603/61‑1 (MQ=255)                                                     accGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTCGGCAAAGGTCATGCTGATGc  <  1:1404538/61‑1 (MQ=255)                                                     accGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTCGGCAAAGGTCATGCTGATGc  <  1:133450/61‑1 (MQ=255)                                                     accGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTCGGCAAAGGTCATGCTGATGc  <  1:965100/61‑1 (MQ=255)                                                                |                                                  TCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTCGGCAAAGGTCATGCTGATGC  >  minE/1330248‑1330360

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