Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,426,862	T→G	100%	G164G (GGA→GGC)	rffD ←	UDP‑N‑acetyl‑D‑mannosaminuronic acid dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,426,862	0	T	G	98.2%	181.9 / ‑4.4	54	G164G (GGA→GGC)	rffD	UDP‑N‑acetyl‑D‑mannosaminuronic acid dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (53/0);  total (53/1)
Fisher's exact test for biased strand distribution p-value = 1.85e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCGACCATTACCTGTCCTGGTAACACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCTGCTCGCCCACCTG  >  minE/2426848‑2426921               |                                                            tCGACCATTACCTGTCCTGGTAACACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCTGCTCGc         <  1:255711/67‑1 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:534273/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:2088008/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:2278354/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:2291763/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:2303546/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:23099/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:2313769/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:271444/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:293357/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:368050/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:373946/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:402163/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:476407/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:207727/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:633229/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:638751/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:702388/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:733748/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:799972/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:82244/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:822885/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:852974/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:864766/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:910921/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:964032/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:995702/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:1464358/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:1071168/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:1090643/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:1120678/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:1175369/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:1202150/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:1272724/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:1290180/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:1383577/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:1400035/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:140616/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:1428564/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:1439971/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:1062440/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:1475520/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:1481638/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:1499115/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:1514342/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:1636258/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:1837548/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:1881661/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:1884292/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:1976853/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACctg  >  1:2042286/1‑65 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACct   >  1:1952766/1‑64 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCAACTCGCCCACct   >  1:206551/1‑64 (MQ=255)          aCCTGGCCCGGTAATACGCGTTCCGGGCACTAAGCa                               >  1:1148027/1‑36 (MQ=19)               |                                                            TCGACCATTACCTGTCCTGGTAACACGCGTTCCGGGCAGTAAGCAATGTTGACGTCCGCCTGCTCGCCCACCTG  >  minE/2426848‑2426921

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