Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	435,104	T→C	59.2%	E32E (GAA→GAG)	ybfE ←	LexA regulated protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	435,104	0	T	C	59.2%	17.8 / 32.1	49	E32E (GAA→GAG)	ybfE	LexA regulated protein
Reads supporting (aligned to +/- strand):  ref base T (18/2);  new base C (7/22);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 6.31e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCGGATTAGTTTTCGGTCGTCCCGGTCGACGCTCGTGCGCGAACAGATCT  >  minE/435045‑435167                                                            |                                                                cGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATcgcg                                                         >  1:2719175/1‑68 (MQ=255) cGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATcgcg                                                         >  1:3286782/1‑68 (MQ=255) cGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATcgcg                                                         >  1:389273/1‑68 (MQ=255)         gCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATcgc                                                          >  1:1856115/1‑59 (MQ=255)         gCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCgg                                                >  1:513178/1‑69 (MQ=255)         gCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCgg                                                >  1:397077/1‑69 (MQ=255)         gCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCgg                                                >  1:631374/1‑69 (MQ=255)         gCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCgg                                                >  1:641471/1‑69 (MQ=255)         gCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCgg                                                >  1:1009621/1‑69 (MQ=255)         gCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCgg                                                >  1:1532467/1‑69 (MQ=255)         gCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCgg                                                >  1:2576798/1‑69 (MQ=255)         gCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCgg                                                >  1:2521940/1‑69 (MQ=255)         gCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCgg                                                >  1:2145172/1‑69 (MQ=255)         gCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCgg                                                >  1:174182/1‑69 (MQ=255)         gCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCgg                                                >  1:1755769/1‑69 (MQ=255)         gCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCgc                                                >  1:2143018/1‑68 (MQ=255)                   ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAg                                                   >  1:447420/1‑56 (MQ=255)                   ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAg                                                   >  1:3009849/1‑56 (MQ=255)                   ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAg                                                   >  1:1934968/1‑56 (MQ=255)                   ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAg                                                   >  1:1010168/1‑56 (MQ=255)                   ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAg                                                   >  1:172908/1‑56 (MQ=255)                   ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAg                                                   >  1:12876/1‑56 (MQ=255)                   ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAg                                                   >  1:1398586/1‑56 (MQ=255)                   ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGtttt                                       <  1:2034053/68‑1 (MQ=255)                   ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGtttt                                       <  1:790954/68‑1 (MQ=255)                   ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGtttt                                       <  1:1140768/68‑1 (MQ=255)                   ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGtttt                                       <  1:115207/68‑1 (MQ=255)                   ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGtttt                                       <  1:1219648/68‑1 (MQ=255)                   ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGtttt                                       <  1:443191/68‑1 (MQ=255)                   ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGtttt                                       <  1:1272485/68‑1 (MQ=255)                   ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGtttt                                       <  1:3363229/68‑1 (MQ=255)                   ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGtttt                                       <  1:3174842/68‑1 (MQ=255)                   ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGtttt                                       <  1:310764/68‑1 (MQ=255)                   ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGtttt                                       <  1:1480214/68‑1 (MQ=255)                   ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGtttt                                       <  1:2701606/68‑1 (MQ=255)                   ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGtttt                                       <  1:1537600/68‑1 (MQ=255)                   ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGtttt                                       <  1:158296/68‑1 (MQ=255)                   ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGtttt                                       <  1:1725334/68‑1 (MQ=255)                        gcgTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGTTTTCGgtcg                                 <  1:2989943/69‑1 (MQ=255)                        gcgTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGTTTTCGgtcg                                 <  1:2783943/69‑1 (MQ=255)                        gcgTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGTTTTCGgtcg                                 <  1:349054/69‑1 (MQ=255)                        gcgTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGTTTTCGgtcg                                 <  1:1713886/69‑1 (MQ=255)                        gcgTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGTTTTCGgtcg                                 <  1:1060152/69‑1 (MQ=255)                           ttttAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGtttt                                       <  1:452646/60‑1 (MQ=255)                                      tGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCGGATTAGTTTTCGgtc                                  >  1:932462/1‑54 (MQ=255)                                      tGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCGGATTAGTTTTCGGTCGTCCCGGTCGACGCt                   >  1:3071678/1‑69 (MQ=255)                                          tttATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGtttt                                       <  1:1248155/45‑1 (MQ=255)                                                       aGCTGTTCATCGCGCGAAAGCGGATTAGTTTTCGGTCGTCCCGGTCGACGCTCGTGCGCGAACAGATCt  <  1:282117/69‑1 (MQ=255)                                                       aGCTGTTCATCGCGCGAAAGCGGATTAGTTTTCGGTCGTCCCGGTCGACGCTCGTGCGCGAACAGATCt  <  1:529066/69‑1 (MQ=255)                                                            |                                                                CGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCGGATTAGTTTTCGGTCGTCCCGGTCGACGCTCGTGCGCGAACAGATCT  >  minE/435045‑435167

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