Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	454,239	G→T	54.7%	A282A (GCC→GCA)	ybgH ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	454,239	0	G	T	54.7%	‑1.2 / 45.3	42	A282A (GCC→GCA)	ybgH	predicted transporter
Reads supporting (aligned to +/- strand):  ref base G (10/9);  new base T (20/3);  total (30/12)
Fisher's exact test for biased strand distribution p-value = 1.95e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

GCGGTTAACGAAGCGGTCGATATAAAGGCTAATCGAGCTACCGCCCTGTTGTGCGAAGGCCCAGAACAACATACTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTG  >  minE/454182‑454298                                                          |                                                            gcgGTTAACGAAGCGGTCGATATAAAGGCTAATCGAGCTACCGCCCTGTTGTGCGAAGGCCCAGaacaa                                                  <  1:1440498/69‑1 (MQ=255) gcgGTTAACGAAGCGGTCGATATAAAGGCTAATCGAGCTACCGCCCTGTTGTGCGAAGGCCCAGaacaa                                                  <  1:3268128/69‑1 (MQ=255)                   gATATAAAGGCTAATCGAGCTACCCTCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                >  1:3045005/1‑69 (MQ=255)                   gATATAAAGGCTAATCGAGCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                >  1:281568/1‑69 (MQ=255)                   gATATAAAGGCTAATCGAGCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                >  1:913099/1‑69 (MQ=255)                   gATATAAAGGCTAATCGAGCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                >  1:747532/1‑69 (MQ=255)                   gATATAAAGGCTAATCGAGCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                >  1:707075/1‑69 (MQ=255)                   gATATAAAGGCTAATCGAGCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                >  1:670355/1‑69 (MQ=255)                   gATATAAAGGCTAATCGAGCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                >  1:645792/1‑69 (MQ=255)                   gATATAAAGGCTAATCGAGCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                >  1:522320/1‑69 (MQ=255)                   gATATAAAGGCTAATCGAGCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                >  1:488157/1‑69 (MQ=255)                   gATATAAAGGCTAATCGAGCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                >  1:3027596/1‑69 (MQ=255)                   gATATAAAGGCTAATCGAGCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                >  1:1135104/1‑69 (MQ=255)                   gATATAAAGGCTAATCGAGCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                >  1:2783203/1‑69 (MQ=255)                   gATATAAAGGCTAATCGAGCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                >  1:2656722/1‑69 (MQ=255)                   gATATAAAGGCTAATCGAGCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                >  1:2622523/1‑69 (MQ=255)                   gATATAAAGGCTAATCGAGCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                >  1:1648013/1‑69 (MQ=255)                   gATATAAAGGCTAATCGAGCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                >  1:11688/1‑69 (MQ=255)                   gATATAAAGGCTAATCGAGCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                >  1:129256/1‑69 (MQ=255)                   gATATAAAGGCTAATCGAGCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                >  1:1534273/1‑69 (MQ=255)                   gATATAAAGGCTAATCGAGCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                >  1:1692809/1‑69 (MQ=255)                   gATATAAAGGCTAATCGAGATACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAg                                >  1:3081202/1‑69 (MQ=255)                         aaGGCTAATCGAGCTACCGCCCTGTTGTGCGAAGGCCCAGAACAACATACTGAAGAAGGTGAGCGTcac                          >  1:2105492/1‑69 (MQ=255)                         aaGGCTAATCGAGCTACCGCCCTGTTGTGCGAAGGCCCAGAACAACATACTGAAGAAGGTGAGCGTcac                          >  1:720240/1‑69 (MQ=255)                         aaGGCTAATCGAGCTACCGCCCTGTTGTGCGAAGGCCCAGAACAACATACTGAAGAAGGTGAGCGTcac                          >  1:265833/1‑69 (MQ=255)                                  cGAGCTACCGCCCTGTTGTGCGAAGGCCCAGAACAACATACTGAAGAAGGTg                                  >  1:2933455/1‑52 (MQ=255)                                  cGAGCTACCGCCCTGTTGTGCGAAGGCCCAGAACAACATACTGAAGAAGGTg                                  >  1:1965174/1‑52 (MQ=255)                                  cGAGCTACCGCCCTGTTGTGCGAAGGCCCAGAACAACATACTGAAGAAGGTg                                  >  1:825206/1‑52 (MQ=255)                                  cGAGCTACCGCCCTGTTGTGCGAAGGCCCAGAACAACATACTGAAGAAGGTg                                  >  1:32057/1‑52 (MQ=255)                                  cGAGCTACCGCCCTGTTGTGCGAAGGCCCAGAACAACATACTGAAGAAGGTg                                  >  1:2335684/1‑52 (MQ=255)                                  cGAGCTACCGCCCTGTTGTGCGAAGGCCCAGAACAACATACTGAAGAAGGTGAGCGTCACAATCAGc                   <  1:182124/67‑1 (MQ=255)                                  cGAGCTACCGCCCTGTTGTGCGAAGGCCCAGAACAACATACTGAAGAAGGTGAGCGTCACAATCAGc                   <  1:1222627/67‑1 (MQ=255)                                  cGAGCTACCGCCCTGTTGTGCGAAGGCCCAGAACAACATACTGAAGAAGGTGAGCGTCACAATCAGc                   <  1:1414855/67‑1 (MQ=255)                                     gCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCa               <  1:286840/68‑1 (MQ=255)                                     gCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCa               <  1:3367426/68‑1 (MQ=255)                                     gCTACCCCCCTGTTGTGCGAATGCCCAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCa               <  1:3187811/68‑1 (MQ=255)                                       tACCGCCCTGTTGTGCGAAGGCCCAGAACAACATACTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCt            <  1:1993122/69‑1 (MQ=255)                                       tACCGCCCTGTTGTGCGAAGGCCCAGAACAACATACTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCt            <  1:1251521/69‑1 (MQ=255)                                       tACCGCCCTGTTGTGCGAAGGCCCAGAACAACATACTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCt            <  1:1581935/69‑1 (MQ=255)                                       tACCGCCCTGTTGTGCGAAGGCCCAGAACAACATACTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCt            <  1:811201/69‑1 (MQ=255)                                                 ttgtGCGAAGGCCCAGAACAACATACTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTg  >  1:58892/1‑69 (MQ=255)                                                 ttgtGCGAAGGCCCAGAACAACATACTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTg  >  1:1747834/1‑69 (MQ=255)                                                          |                                                            GCGGTTAACGAAGCGGTCGATATAAAGGCTAATCGAGCTACCGCCCTGTTGTGCGAAGGCCCAGAACAACATACTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTG  >  minE/454182‑454298

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