Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,524,816	A→G	50.7%	V186A (GTA→GCA)	yfcJ ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,524,816	0	A	G	50.7%	‑1.5 / 62.8	71	V186A (GTA→GCA)	yfcJ	predicted transporter
Reads supporting (aligned to +/- strand):  ref base A (9/26);  new base G (13/23);  total (22/49)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.04e-01

GCCAGGGCCGGTACTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCA  >  minE/1524759‑1524872                                                          |                                                         gCCAGGGCCGGTACTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTgtggtg                                               <  1:2327818/69‑1 (MQ=255) gCCAGGGCCGGTACTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTgtggtg                                               <  1:456956/69‑1 (MQ=255) gCCAGGGCCGGTACTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTAAGGGTAATGCCATTgtggtg                                               >  1:2740386/1‑69 (MQ=255) gCCAGGGCCGGTACTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTAAGGGTAATGCCATTgtggtg                                               >  1:523422/1‑69 (MQ=255) gCCAGGGCCGGTACTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTAAGGGTAATGCCATTgtggtg                                               >  1:2225111/1‑69 (MQ=255) gCCAGGGCCGGTACTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTAAGGGTAATGCCATTgtggtg                                               >  1:2294517/1‑69 (MQ=255) gCCAGGGCCGGTACTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTAAGGGTAATGCCATTgtggtg                                               >  1:643693/1‑69 (MQ=255) gCCAGGGCCGGTACTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTAAGGGTAATGCCATTgtggtg                                               >  1:1955249/1‑69 (MQ=255) gCCAGGGCCGGTACTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTAAGGGTAATGCCATTgtggtg                                               >  1:74986/1‑69 (MQ=255) gCCAGGGCCGGTACTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTAAGGGTAATGCCATTgtggtg                                               >  1:1467661/1‑69 (MQ=255) gCCAGGGCCGGTACTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTAAGGGTAATGCCATTgtggtg                                               >  1:818854/1‑69 (MQ=255) gCCAGGGCCGGTACTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTAAGGGTAATGCCATTgtggtg                                               >  1:1376443/1‑69 (MQ=255) gCCAGGGCCGGTACTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTAAGGGTAATGCCATTgtggtg                                               >  1:329928/1‑69 (MQ=255) gCCAGGGCCGGTACTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTAAGGGTAATGCCATTgtggtg                                               >  1:3407431/1‑69 (MQ=255) gCCAGGGCCGGTACTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTAAGGGTAATGCCATTgtggtg                                               >  1:876530/1‑69 (MQ=255)             aCTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGc                                    >  1:543424/1‑68 (MQ=255)             aCTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGc                                    >  1:573995/1‑68 (MQ=255)             aCTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGc                                    >  1:1719329/1‑68 (MQ=255)             aCTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGc                                    >  1:2736595/1‑68 (MQ=255)                 tGCGCACTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTgtggt                                                >  1:868366/1‑52 (MQ=255)                 tGCGCACTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCaa                               >  1:315549/1‑69 (MQ=255)                 tGCGCACTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCaa                               >  1:886128/1‑69 (MQ=255)                 tGCGCACTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCaa                               >  1:1159716/1‑69 (MQ=255)                 tGCGCACTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCaa                               >  1:1217630/1‑69 (MQ=255)                       cTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGt                         <  1:464076/69‑1 (MQ=255)                       cTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGt                         <  1:690444/69‑1 (MQ=255)                       cTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGt                         <  1:2404191/69‑1 (MQ=255)                       cTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGt                         <  1:2313088/69‑1 (MQ=255)                       cTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGt                         <  1:2054488/69‑1 (MQ=255)                       cTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGt                         <  1:301375/69‑1 (MQ=255)                       cTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGt                         <  1:2162214/69‑1 (MQ=255)                       cTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGt                         <  1:1115368/69‑1 (MQ=255)                       cTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGt                         <  1:3072062/69‑1 (MQ=255)                       cTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGt                         <  1:2271000/69‑1 (MQ=255)                       cTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGt                         <  1:3237713/69‑1 (MQ=255)                       cTGTGCCGTAACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGt                         <  1:2038910/69‑1 (MQ=255)                         gtgCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCaaa                              <  1:193489/62‑1 (MQ=255)                         gtgCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCaaa                              <  1:122541/62‑1 (MQ=255)                              gTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGTAATGGc                   <  1:2931143/68‑1 (MQ=255)                              gTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGTAATGGc                   <  1:36994/68‑1 (MQ=255)                              gTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGTAATGGc                   <  1:1168355/68‑1 (MQ=255)                              gTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGTAATGGc                   <  1:858895/68‑1 (MQ=255)                              gTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGTAATGGc                   <  1:1268736/68‑1 (MQ=255)                              gTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGTAATGGc                   <  1:1496220/68‑1 (MQ=255)                              gTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGTAATGGc                   <  1:1673372/68‑1 (MQ=255)                              gTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGTAATGGc                   <  1:1703169/68‑1 (MQ=255)                              gTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGTAATGGc                   <  1:1783030/68‑1 (MQ=255)                               ttACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGTAATGGc                   <  1:2539593/67‑1 (MQ=255)                                              aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:745721/69‑1 (MQ=255)                                              aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:557377/69‑1 (MQ=255)                                              aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:706657/69‑1 (MQ=255)                                              aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:658395/69‑1 (MQ=255)                                              aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:2274082/69‑1 (MQ=255)                                              aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:2964725/69‑1 (MQ=255)                                              aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:1175901/69‑1 (MQ=255)                                              aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:1210673/69‑1 (MQ=255)                                              aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:1224998/69‑1 (MQ=255)                                              aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:145531/69‑1 (MQ=255)                                              aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:1552615/69‑1 (MQ=255)                                              aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:1694601/69‑1 (MQ=255)                                              aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:1755998/69‑1 (MQ=255)                                              aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:185659/69‑1 (MQ=255)                                              aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:1977183/69‑1 (MQ=255)                                              aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:2062857/69‑1 (MQ=255)                                              aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:2091464/69‑1 (MQ=255)                                              aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:2105658/69‑1 (MQ=255)                                              aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:226710/69‑1 (MQ=255)                                              aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:2379799/69‑1 (MQ=255)                                              aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:45183/69‑1 (MQ=255)                                               gTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:1012181/68‑1 (MQ=255)                                               gTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:1104699/68‑1 (MQ=255)                                                          |                                                         GCCAGGGCCGGTACTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCA  >  minE/1524759‑1524872

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