Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,610,182	A→C	100%	V58V (GTT→GTG)	eutC ←	ethanolamine ammonia‑lyase, small subunit (light chain)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,610,182	0	A	C	75.0%	19.8 / 3.2	16	V58V (GTT→GTG)	eutC	ethanolamine ammonia‑lyase, small subunit (light chain)
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base C (12/0);  total (14/2)
Fisher's exact test for biased strand distribution p-value = 5.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGCGCGATGCGGATTTTCAACACCAATCCACGCTTTTGCTTCAGCGGAACCTAAATCCAGCGCACAGCTTTCCGAGGTCACCGGTG  >  minE/1610164‑1610249                   |                                                                    tGCGCGATGCGGATTTTCAACACCAATCCACGCTTTTGCTTCAGCGGAACCTAAATCCAGCGCACAGCt                   <  1:3036200/69‑1 (MQ=255)    gcgATGTGGGTTTTCCACACCAATCCACGCTTTTGCTTCCGCAGAACCTAAATCCAGCGCGCAGCttt                 >  1:1370827/1‑68 (MQ=255)    gcgATGTGGGTTTTCCACACCAATCCACGCTTTTGCTTCCGCAGAACCTAAATCCAGCGCGCAGCttt                 >  1:1389853/1‑68 (MQ=255)    gcgATGTGGGTTTTCCACACCAATCCACGCTTTTGCTTCCGCAGAACCTAAATCCAGCGCGCAGCttt                 >  1:1676638/1‑68 (MQ=255)    gcgATGTGGGTTTTCCACACCAATCCACGCTTTTGCTTCCGCAGAACCTAAATCCAGCGCGCAGCttt                 >  1:1927811/1‑68 (MQ=255)    gcgATGTGGGTTTTCCACACCAATCCACGCTTTTGCTTCCGCAGAACCTAAATCCAGCGCGCAGCttt                 >  1:239512/1‑68 (MQ=255)    gcgATGTGGGTTTTCCACACCAATCCACGCTTTTGCTTCCGCAGAACCTAAATCCAGCGCGCAGCttt                 >  1:2402759/1‑68 (MQ=255)    gcgATGTGGGTTTTCCACACCAATCCACGCTTTTGCTTCCGCAGAACCTAAATCCAGCGCGCAGCttt                 >  1:2558734/1‑68 (MQ=255)    gcgATGTGGGTTTTCCACACCAATCCACGCTTTTGCTTCCGCAGAACCTAAATCCAGCGCGCAGCttt                 >  1:2984407/1‑68 (MQ=255)    gcgATGTGGGTTTTCCACACCAATCCACGCTTTTGCTTCCGCAGAACCTAAATCCAGCGCGCAGCttt                 >  1:3123474/1‑68 (MQ=255)    gcgATGTGGGTTTTCCACACCAATCCACGCTTTTGCTTCCGCAGAACCTAAATCCAGCGCGCAGCttt                 >  1:484170/1‑68 (MQ=255)    gcgATGTGGGTTTTCCACACCAATCCACGCTTTTGCTTCCGCAGAACCTAAATCCAGCGCGCAGCttt                 >  1:804166/1‑68 (MQ=255)    gcgATGTGGGTTTTCCACACCAATCCACGCTTTTGCTTCCGCAGAACCTAAATCCAGCGCGCAGCttt                 >  1:987017/1‑68 (MQ=255)        tGCGGATTTTCAACACCAATCCACGCTTTTGCTTCAGCGGAACCTAAATCCAGCGCACAGCTTTCCGAg            <  1:2471890/69‑1 (MQ=255)                  cAACACCAATCCACGCTTTTGCTTCAGCGGAACCTAAATCCAGCGCACAGCTTTCCGAGGTCACCGGTg  >  1:1857865/1‑69 (MQ=255)                  cAACACCAATCCACGCTTTTGCTTCAGCGGAACCTAAATCCAGCGCACAGCTTTCCGAGGTCACCGGTg  >  1:961173/1‑69 (MQ=255)                   |                                                                    TGCGCGATGCGGATTTTCAACACCAATCCACGCTTTTGCTTCAGCGGAACCTAAATCCAGCGCACAGCTTTCCGAGGTCACCGGTG  >  minE/1610164‑1610249

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