Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	189,850:1	+T	100%	coding (1347/2142 nt)	ldcC →	lysine decarboxylase 2, constitutive

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	189,850	1	.	T	78.6%	117.5 / 30.2	42	coding (1347/2142 nt)	ldcC	lysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base . (0/9);  new base T (33/0);  total (33/9)
Fisher's exact test for biased strand distribution p-value = 2.24e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

CGCAAAGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGG‑TTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGC  >  minE/189809‑189888                                           |                                       cGCAAAGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTTCGATATCTGGCa                       >  1:1475109/1‑60 (MQ=255) cGCAAAGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTTCGATATCTGGCAAccgc                  >  1:319714/1‑65 (MQ=255) cGCAAAGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTTCGATATCTGGCAAccgc                  >  1:1711195/1‑65 (MQ=255) cGCAAAGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTTCGATATCTGGCAAccgc                  >  1:1811161/1‑65 (MQ=255) cGCAAAGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTTCGATATCTGGCAAccgc                  >  1:1849243/1‑65 (MQ=255) cGCAAAGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTTCGATATCTGGCAAccgc                  >  1:1009968/1‑65 (MQ=255) cGCAAAGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTTCGATATCTGGCAAccgc                  >  1:1998458/1‑65 (MQ=255) cGCAAAGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTTCGATATCTGGCAAccgc                  >  1:2027110/1‑65 (MQ=255) cGCAAAGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTTCGATATCTGGCAAccgc                  >  1:2068771/1‑65 (MQ=255) cGCAAAGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTTCGATATCTGGCAAccgc                  >  1:1886425/1‑65 (MQ=255) cGCAAAGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTTCGATATCTGGCAAccgc                  >  1:432509/1‑65 (MQ=255) cGCAAAGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTTCGATATCTGGCAAccgc                  >  1:551894/1‑65 (MQ=255) cGCAAAGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTTCGATATCTGGCAAccgc                  >  1:718765/1‑65 (MQ=255) cGCAAAGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTTCGATATCTGGCAAccgc                  >  1:744032/1‑65 (MQ=255) cGCAAAGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTTCGATATCTGGCAAccgc                  >  1:765931/1‑65 (MQ=255) cGCAAAGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTTCGATATCTGGCAAccgc                  >  1:818328/1‑65 (MQ=255) cGCAAAGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTTCGATATCTGGCAAccgc                  >  1:824308/1‑65 (MQ=255) cGCAAAGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTTCGATATCTGGCAAccgc                  >  1:925713/1‑65 (MQ=255) cGCAAAGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTTCGATATCTGGCAAccgc                  >  1:1637337/1‑65 (MQ=255) cGCAAAGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTTCGATATCTGGCAAccgc                  >  1:1017162/1‑65 (MQ=255) cGCAAAGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTTCGATATCTGGCAAccgc                  >  1:1048533/1‑65 (MQ=255) cGCAAAGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTTCGATATCTGGCAAccgc                  >  1:1096166/1‑65 (MQ=255) cGCAAAGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTTCGATATCTGGCAAccgc                  >  1:1144222/1‑65 (MQ=255) cGCAAAGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTTCGATATCTGGCAAccgc                  >  1:1169580/1‑65 (MQ=255) cGCAAAGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTTCGATATCTGGCAAccgc                  >  1:1182451/1‑65 (MQ=255) cGCAAAGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTTCGATATCTGGCAAccgc                  >  1:1312757/1‑65 (MQ=255) cGCAAAGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTTCGATATCTGGCAAccgc                  >  1:1461054/1‑65 (MQ=255) cGCAAAGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTTCGATATCTGGCAAccgc                  >  1:1539148/1‑65 (MQ=255) cGCAAAGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTTCGATATCTGGCAAccgc                  >  1:1584121/1‑65 (MQ=255) cGCAAAGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTTCGATATCTGGCAAccgc                  >  1:1658123/1‑65 (MQ=255) cGCAAAGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTTCGATATCTGGCAAccg                   >  1:648096/1‑64 (MQ=255) cGCAAAGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGGTTTTTTCGATATCTGGCAAcagc                  >  1:29143/1‑65 (MQ=255) cGCAAAGAGGTCCAGCGGCTCCGGGAAGAGTCTGACGGTTGGTTTTTTCGATATCTGGCAAccgc                  >  1:1369910/1‑65 (MQ=255)                                 tGACGGTTGG‑TTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGc  <  1:2131963/48‑1 (MQ=255)                                 tGACGGTTGG‑TTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGc  <  1:2096659/48‑1 (MQ=255)                                 tGACGGTTGG‑TTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGc  <  1:461883/48‑1 (MQ=255)                                 tGACGGTTGG‑TTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGc  <  1:482195/48‑1 (MQ=255)                                 tGACGGTTGG‑TTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGc  <  1:538528/48‑1 (MQ=255)                                 tGACGGTTGG‑TTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGc  <  1:671732/48‑1 (MQ=255)                                 tGACGGTTGG‑TTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGc  <  1:1571235/48‑1 (MQ=255)                                 tGACGGTTGG‑TTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGc  <  1:1582548/48‑1 (MQ=255)                                 tGACGGTTGG‑TTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGc  <  1:1032919/48‑1 (MQ=255)                                           |                                       CGCAAAGAGGTCCAGCGGCTGCGGGAAGAGTCTGACGGTTGG‑TTTTTCGATATCTGGCAACCGCCGCAGGTGGATGAAGC  >  minE/189809‑189888

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