Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,998,493	C→G	44.7%	A213A (GCG→GCC)	pitB ←	phosphate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,998,493	0	C	G	44.7%	11.7 / 38.9	38	A213A (GCG→GCC)	pitB	phosphate transporter
Reads supporting (aligned to +/- strand):  ref base C (7/14);  new base G (10/7);  total (17/21)
Fisher's exact test for biased strand distribution p-value = 1.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

GCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTG  >  minE/1998427‑1998561                                                                   |                                                                     gCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAgcgc                                                                    >  1:2960214/1‑69 (MQ=255) gCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAgcgc                                                                    >  1:1698060/1‑69 (MQ=255)   cGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCaa                                                                  >  1:831638/1‑69 (MQ=255)       cccTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAGACAATCAGCGCAATAcg                                                              >  1:2785876/1‑69 (MQ=255)       cccTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATAcg                                                              >  1:2179093/1‑69 (MQ=255)       cccTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATAcg                                                              >  1:244049/1‑69 (MQ=255)       cccTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATAcg                                                              >  1:2181103/1‑69 (MQ=255)       cccTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATAcg                                                              >  1:3249519/1‑69 (MQ=255)       cccTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATAcg                                                              >  1:493239/1‑69 (MQ=255)                 ccGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGGGGAAACAATGAGGGCAATACGCGTCCAGAAt                                                    >  1:1351463/1‑69 (MQ=18)                 ccGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATGAGGGCAATACGCGTCCAGAAt                                                    >  1:3394559/1‑69 (MQ=255)                 ccGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATGAGGGCAATACGCGTCCAGAAt                                                    >  1:52076/1‑69 (MQ=255)                 ccGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATGAGGGCAATACGCGTCCAGAAt                                                    >  1:704279/1‑69 (MQ=255)                 ccGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATGAGGGCAATACGCGTCCAGAAt                                                    >  1:2241868/1‑69 (MQ=255)                 ccGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATGAGGGCAATACGCGTCCAGAAt                                                    >  1:2302244/1‑69 (MQ=255)                 ccGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATGAGGGCAATACGCGTCCAGAAt                                                    >  1:2936003/1‑69 (MQ=255)                 ccGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGAAGCGGAAACAATGAGGGCAATACGCGTCCAGAAt                                                    >  1:2266275/1‑69 (MQ=18)                 ccGTCGTTAGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATGAGGGCAATACGCGTCCAGAAt                                                    >  1:237072/1‑69 (MQ=18)                 cagtcgtTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATGAGGGCAATACGCGTCCAGAAt                                                    >  1:1797109/3‑69 (MQ=18)                                               cccGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTg                      <  1:1203587/69‑1 (MQ=255)                                               cccGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTg                      <  1:615517/69‑1 (MQ=255)                                               cccGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTg                      <  1:2874076/69‑1 (MQ=255)                                               cccGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTg                      <  1:2833990/69‑1 (MQ=255)                                               cccGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTg                      <  1:701414/69‑1 (MQ=255)                                               cccGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTg                      <  1:1831949/69‑1 (MQ=255)                                               cccGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTg                      <  1:1905483/69‑1 (MQ=255)                                                ccGCAGCGGTAACAATCAGCGCAATTCGCGTCCTGAATGGCGGTTTACGTTTGGCTTTTTTCTTTTTGc                     <  1:1264966/69‑1 (MQ=255)                                                ccGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGc                     <  1:783059/69‑1 (MQ=255)                                                ccGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGc                     <  1:2601467/69‑1 (MQ=255)                                                ccGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGc                     <  1:2878336/69‑1 (MQ=255)                                                ccGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGc                     <  1:255359/69‑1 (MQ=255)                                                ccGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGc                     <  1:172770/69‑1 (MQ=255)                                                ccGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGc                     <  1:440552/69‑1 (MQ=255)                                                          aaCAATGAGGGCAATACGCGTCCAGAATGGCGGCTTACGTTTGCCTTTTTTCttttt                       <  1:2724645/57‑1 (MQ=25)                                                          aaCAATGAGGGCAATACGCGTCCAGAATGGCGGCTTACGTTTGCCTTTTTTCttttt                       <  1:460992/57‑1 (MQ=25)                                                          aaCAATGAGGGCAATACGCGTCCAGAATGGCGGCTTACGTTTGCCTTTTTTCttttt                       <  1:435874/57‑1 (MQ=25)                                                          aaCAATGAGGGCAATACGCGTCCAGAATGGCGGCTTACGTTTGCCTTTTTTCttttt                       <  1:2102705/57‑1 (MQ=25)                                                           aCAATGAGGGCAATACGCGTCCAGAATGGCGGCTTACGTTTGCCTTTTTTCttttt                       <  1:2035537/56‑1 (MQ=25)                                                           aCAATGAGGGCAATACGCGTCCAGAATGGCGGCTTACGTTTGCCTTTTTTCttttt                       <  1:1120741/56‑1 (MQ=25)                                                           aCAATGAGGGCAATACGCGTCCAGAATGGCGGCTTACGTTTGCCTTTTTTCttttt                       <  1:85764/56‑1 (MQ=25)                                                                   cgcAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTg  <  1:3265361/69‑1 (MQ=255)                                                                   cgcAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTg  <  1:3071503/69‑1 (MQ=255)                                                                   cgcAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTg  <  1:1655625/69‑1 (MQ=255)                                                                   cgcAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTg  <  1:1569634/69‑1 (MQ=255)                                                                   cgcAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTg  <  1:2811683/69‑1 (MQ=255)                                                                   cgcAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTg  <  1:872274/69‑1 (MQ=255)                                                                   cgcAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTg  <  1:883217/69‑1 (MQ=255)                                                                   |                                                                     GCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTG  >  minE/1998427‑1998561

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