Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,390,482	A→G	28.8%	N125N (AAT→AAC)	rhtC ←	threonine efflux system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,390,482	0	A	G	28.8%	88.7 / 39.7	90	N125N (AAT→AAC)	rhtC	threonine efflux system
Reads supporting (aligned to +/- strand):  ref base A (29/35);  new base G (11/15);  total (40/50)
Fisher's exact test for biased strand distribution p-value = 8.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAACGTTATCACCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTAAACCTTTCAGGAA  >  minE/2390422‑2390514                                                             |                                 gCCAACGTTATCACCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTa                               >  1:290266/1‑64 (MQ=255) gCCAACGTTATCACCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCgag                          >  1:2178145/1‑69 (MQ=255) gCCAACGTTATCACCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCgag                          >  1:3107103/1‑69 (MQ=255) gCCAACGTTATCACCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCgag                          >  1:2818823/1‑69 (MQ=255) gCCAACGTTATCACCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCgag                          >  1:2789396/1‑69 (MQ=255) gCCAACGTTATCACCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCgag                          >  1:842163/1‑69 (MQ=255)  ccAACGTTATCACCGACAAACAATGAGAAGACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCgaga                         >  1:1989167/1‑69 (MQ=255)  ccAACGTTATCACCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCgaga                         >  1:246006/1‑69 (MQ=255)  ccAACGTTATCACCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCgaga                         >  1:2155563/1‑69 (MQ=255)  ccAACGTTATCACCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCgaga                         >  1:748179/1‑69 (MQ=255)  ccAACGTTATCACCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCgaga                         >  1:3433287/1‑69 (MQ=255)  ccAACGTTATCACCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCgaga                         >  1:936741/1‑69 (MQ=255)  ccAACGTTATCACCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCgaga                         >  1:1686605/1‑69 (MQ=255)  ccAACGTTATCACCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCgaga                         >  1:1868759/1‑69 (MQ=255)  ccAACGTTATCACCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCgaga                         >  1:2328758/1‑69 (MQ=255)     aCGTTATCACCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTg                      <  1:2163289/69‑1 (MQ=255)     aCGTTATCACCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTg                      <  1:2483659/69‑1 (MQ=255)     aCGTTATCACCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTg                      <  1:2007423/69‑1 (MQ=255)      cattATCACCGACAAACAGCGAGAACACCGAACCAAAGTAGATAATCGCTTTAGGGTTCGCGAGATTg                      >  1:2167363/3‑68 (MQ=255)      cattATCACCGACAAACAGCGAGAACACCGAACCAAAGTAGATAATCGCTTTAGGGTTCGCGAGATTg                      >  1:2157348/3‑68 (MQ=255)      cattATCACCGACAAACAGCGAGAACACCGAACCAAAGTAGATAATCGCTTTAGGGTTCGCGAGATTg                      >  1:222187/3‑68 (MQ=255)      cattATCACCGACAAACAGCGAGAACACCGAACCAAAGTAGATAATCGCTTTAGGGTTCGCGAGATTg                      >  1:299280/3‑68 (MQ=255)      cattATCACCGACAAACAGCGAGAACACCGAACCAAAGTAGATAATCGCTTTAGGGTTCGCGAGATTg                      >  1:331323/3‑68 (MQ=255)      cattATCACCGACAAACAGCGAGAACACCGAACCAAAGTAGATAATCGCTTTAGGGTTCGCGAGATTg                      >  1:396130/3‑68 (MQ=255)      cattATCACCGACAAACAGCGAGAACACCGAACCAAAGTAGATAATCGCTTTAGGGTTCGCGAGATTg                      >  1:568352/3‑68 (MQ=255)      cattATCACCGACAAACAGCGAGAACACCGAACCAAAGTAGATAATCGCTTTAGGGTTCGCGAGATTg                      >  1:1404893/3‑68 (MQ=255)      cattATCACCGACAAACAGCGAGAACACCGAACCAAAGTAGATAATCGCTTTAGGGTTCGCGAGATTg                      >  1:14005/3‑68 (MQ=255)      cattATCACCGACAAACAGCGAGAACACCGAACCAAAGTAGATAATCGCTTTAGGGTTCGCGAGATTg                      >  1:745768/3‑68 (MQ=255)      cattATCACCGACAAACAGCGAGAACACCGAACCAAAGTAGATAATCGCTTTAGGGTTCGCGAGATTg                      >  1:2431544/3‑68 (MQ=255)          aTCACCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCa                  >  1:1495916/1‑68 (MQ=255)          aTCACCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCa                  >  1:2693287/1‑68 (MQ=255)          aTCACCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCa                  >  1:2209203/1‑68 (MQ=255)             aCCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTaa              >  1:2520620/1‑69 (MQ=255)             aCCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTaa              >  1:2691327/1‑69 (MQ=255)             aCCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTaa              <  1:2934904/69‑1 (MQ=255)             aCCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTaa              <  1:2865885/69‑1 (MQ=255)             aCCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTaa              <  1:2580004/69‑1 (MQ=255)             aCCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTaa              >  1:376356/1‑69 (MQ=255)             aCCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTaa              >  1:55819/1‑69 (MQ=255)             aCCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTaa              >  1:648480/1‑69 (MQ=255)             aCCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTaa              >  1:948179/1‑69 (MQ=255)             aCCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTaa              >  1:2426322/1‑69 (MQ=255)             aCCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTaa              <  1:1157855/69‑1 (MQ=255)             aCCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTaa              >  1:163077/1‑69 (MQ=255)             aCCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTaa              >  1:2278620/1‑69 (MQ=255)             aCCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTaa              >  1:139555/1‑69 (MQ=255)             aCCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTaa              <  1:1621384/69‑1 (MQ=255)             aCCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTaa              <  1:1526318/69‑1 (MQ=255)             aCCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTaa              <  1:2077599/69‑1 (MQ=255)               cGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTaaa             <  1:614660/68‑1 (MQ=255)               cGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTaaa             <  1:412373/68‑1 (MQ=255)                      caaTGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTaaa             <  1:3039774/61‑1 (MQ=255)                      caaTGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTaaa             <  1:1111910/61‑1 (MQ=255)                      caaTGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTaaa             <  1:3063400/61‑1 (MQ=255)                      caaTGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTaaa             <  1:1435267/61‑1 (MQ=255)                      caaTGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTaaa             <  1:1333372/61‑1 (MQ=255)                      caaTGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTaaa             <  1:848602/61‑1 (MQ=255)                      caaTGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTaaa             <  1:888807/61‑1 (MQ=255)                      caaTGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTaaa             <  1:1893924/61‑1 (MQ=255)                      caaTGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTaaa             <  1:2125723/61‑1 (MQ=255)                      caaTGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTaaa             <  1:2865905/61‑1 (MQ=255)                      caaTGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTaaa             <  1:2866166/61‑1 (MQ=255)                        aTGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTAAACCTTTCAGGa   <  1:2864446/69‑1 (MQ=255)                        aTGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTAAACCTTTCAGGa   <  1:361549/69‑1 (MQ=255)                        aTGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTAAACCTTTCAGGa   <  1:1762954/69‑1 (MQ=255)                         tGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTAAACCTTTCAGGaa  <  1:2937472/69‑1 (MQ=255)                         tGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTAAACCTTTCAGGaa  <  1:2108587/69‑1 (MQ=255)                         tGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTAAACCTTTCAGGaa  <  1:271654/69‑1 (MQ=255)                         tGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTAAACCTTTCAGGaa  <  1:2339075/69‑1 (MQ=255)                         tGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTAAACCTTTCAGGa   <  1:2423823/68‑1 (MQ=255)                         cgagaACACCGAACCAAAGTAGATAATCGCTTTAGGGTTCGCGAGATTGGTCAGTAAGccc          <  1:2144985/60‑2 (MQ=255)                         cgagaACACCGAACCAAAGTAGATAATCGCTTTAGGGTTCGCGAGATTGGTCAGTAAGccc          <  1:1090812/60‑2 (MQ=255)                         cgagaACACCGAACCAAAGTAGATAATCGCTTTAGGGTTCGCGAGATTGGTCAGTAAGccc          <  1:933080/60‑2 (MQ=255)                         cgagaACACCGAACCAAAGTAGATAATCGCTTTAGGGTTCGCGAGATTGGTCAGTAAGccc          <  1:1444366/60‑2 (MQ=255)                         cgagaACACCGAACCAAAGTAGATAATCGCTTTAGGGTTCGCGAGATTGGTCAGTAAGccc          <  1:2428116/60‑2 (MQ=255)                         cgagaACACCGAACCAAAGTAGATAATCGCTTTAGGGTTCGCGAGATTGGTCAGTAAGccc          <  1:524260/60‑2 (MQ=255)                         cgagaACACCGAACCAAAGTAGATAATCGCTTTAGGGTTCGCGAGATTGGTCAGTAAGccc          <  1:155134/60‑2 (MQ=255)                         cgagaACACCGAACCAAAGTAGATAATCGCTTTAGGGTTCGCGAGATTGGTCAGTAAGccc          <  1:389243/60‑2 (MQ=255)                         cgagaACACCGAACCAAAGTAGATAATCGCTTTAGGGTTCGCGAGATTGGTCAGTAAGccc          <  1:3861/60‑2 (MQ=255)                         cgagaACACCGAACCAAAGTAGATAATCGCTTTAGGGTTCGCGAGATTGGTCAGTAAGccc          <  1:1552070/60‑2 (MQ=255)                         cgagaACACCGAACCAAAGTAGATAATCGCTTTAGGGTTCGCGAGATTGGTCAGTAAGccc          <  1:3099059/60‑2 (MQ=255)                         cgagaACACCGAACCAAAGTAGATAATCGCTTTAGGGTTCGCGAGATTGGTCAGTAAGccc          <  1:3086301/60‑2 (MQ=255)                         cgagaACACCGAACCAAAGTAGATAATCGCTTTAGGGTTCGCGAGATTGGTCAGTAAGccc          <  1:2005650/60‑2 (MQ=255)                         cgagaACACCGAACCAAAGTAGATAATCGCTTTAGGGTTCGCGAGATTGGTCAGTAAGccc          <  1:305094/60‑2 (MQ=255)                          gagaACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTAAACCTTTCAGGaa  <  1:694269/68‑1 (MQ=255)                          gagaACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTAAACCTTTCAGGaa  <  1:82638/68‑1 (MQ=255)                          gagaACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTAAACCTTTCAGGaa  <  1:1853780/68‑1 (MQ=255)                          gagaACACCGAACCAAAGTAGATAATCGCTTTAGGGTTCGCGAGATTGGTCAGTAAGccc          <  1:2658729/60‑2 (MQ=255)                                      ccAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTAAACCTTTCAGCaa  <  1:132935/56‑1 (MQ=255)                                                  aTCGCTTTCGGATTAGCGAGATTGGTCAGTAAACCTTTCAGGa   >  1:3287679/1‑43 (MQ=255)                                                             |                                 GCCAACGTTATCACCGACAAACAATGAGAACACCGAGCCAAAGTAGATAATCGCTTTCGGATTAGCGAGATTGGTCAGTAAACCTTTCAGGAA  >  minE/2390422‑2390514

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