Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	923,259	Δ1 bp	100%	coding (93/1614 nt)	mppA →	murein tripeptide (L‑ala‑gamma‑D‑glutamyl‑meso‑DAP) transporter subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	923,259	0	G	.	86.7%	41.7 / 3.5	15	coding (93/1614 nt)	mppA	murein tripeptide (L‑ala‑gamma‑D‑glutamyl‑meso‑DAP) transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base . (0/13);  total (2/13)
Fisher's exact test for biased strand distribution p-value = 9.52e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.99e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ATTTCTCTTTCGTATGCTGCAGAAGTTCCGAGCGGCACAGTACTGGCAGAGAAGCAGGAGCTGGTGCGCCACATTAAA  >  minE/923215‑923292                                             |                                  aTTTCTCTTTCGTATGCTGCAGAAGTTCCGAGCGGCACAGTACTGGCAGAGAAGCAGGAGCTGGTGCGCca         >  1:2249197/1‑71 (MQ=255) aTTTCTCTTTCGTATGCTGCAGAAGTTCCGAGCGGCACAGTACTGGCAGAGAAGCAGGAGCTGGTGCGCca         >  1:874906/1‑71 (MQ=255)          tCGTATGCTGCAGAAGTTCCGAGCGGCACAGTACT‑GCAGAGAAGCAGGAGCTGGTGCGCCACATTaaa  <  1:1232194/68‑1 (MQ=255)          tCGTATGCTGCAGAAGTTCCGAGCGGCACAGTACT‑GCAGAGAAGCAGGAGCTGGTGCGCCACATTaaa  <  1:1336094/68‑1 (MQ=255)          tCGTATGCTGCAGAAGTTCCGAGCGGCACAGTACT‑GCAGAGAAGCAGGAGCTGGTGCGCCACATTaaa  <  1:1510948/68‑1 (MQ=255)          tCGTATGCTGCAGAAGTTCCGAGCGGCACAGTACT‑GCAGAGAAGCAGGAGCTGGTGCGCCACATTaaa  <  1:1529607/68‑1 (MQ=255)          tCGTATGCTGCAGAAGTTCCGAGCGGCACAGTACT‑GCAGAGAAGCAGGAGCTGGTGCGCCACATTaaa  <  1:1635698/68‑1 (MQ=255)          tCGTATGCTGCAGAAGTTCCGAGCGGCACAGTACT‑GCAGAGAAGCAGGAGCTGGTGCGCCACATTaaa  <  1:1969600/68‑1 (MQ=255)          tCGTATGCTGCAGAAGTTCCGAGCGGCACAGTACT‑GCAGAGAAGCAGGAGCTGGTGCGCCACATTaaa  <  1:2312235/68‑1 (MQ=255)          tCGTATGCTGCAGAAGTTCCGAGCGGCACAGTACT‑GCAGAGAAGCAGGAGCTGGTGCGCCACATTaaa  <  1:2324688/68‑1 (MQ=255)          tCGTATGCTGCAGAAGTTCCGAGCGGCACAGTACT‑GCAGAGAAGCAGGAGCTGGTGCGCCACATTaaa  <  1:24673/68‑1 (MQ=255)          tCGTATGCTGCAGAAGTTCCGAGCGGCACAGTACT‑GCAGAGAAGCAGGAGCTGGTGCGCCACATTaaa  <  1:344449/68‑1 (MQ=255)          tCGTATGCTGCAGAAGTTCCGAGCGGCACAGTACT‑GCAGAGAAGCAGGAGCTGGTGCGCCACATTaaa  <  1:749418/68‑1 (MQ=255)           cGTATGCTGCAGAAGTTCCGAGCGGCACAGTACT‑GCAGAGAAGCAGTAGCTGGTGCGCCACATTaaa  <  1:2337624/67‑1 (MQ=255)           cGTATGCTGCAGAAGTTCCGAGCGGCACAGTACT‑GCAGAGAAGCAGGAGCTGGTGCGCCACATTaaa  <  1:56106/67‑1 (MQ=255)                                             |                                  ATTTCTCTTTCGTATGCTGCAGAAGTTCCGAGCGGCACAGTACTGGCAGAGAAGCAGGAGCTGGTGCGCCACATTAAA  >  minE/923215‑923292

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