Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	388,248	T→C	47.1%	K43E (AAA→GAA)	ybdH ←	predicted oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	388,248	0	T	C	47.1%	1.7 / 58.0	68	K43E (AAA→GAA)	ybdH	predicted oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (11/25);  new base C (11/21);  total (22/46)
Fisher's exact test for biased strand distribution p-value = 7.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGAAAAAATCGTGCAGGTG  >  minE/388202‑388305                                               |                                                          tcAAACGCAGGCGGCAGTTTGGTTTGTGCCGCAGCAATGGCGCTTTCGCCGTAGATCCACACCGCGCga                                     <  1:1812832/68‑1 (MQ=255) tcAAACGCAGGCGGCAGTTTGGTTTGTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCga                                     <  1:344117/68‑1 (MQ=255) tcAAACGCAGGCGGCAGTTTGGTTTGTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCga                                     <  1:2114552/68‑1 (MQ=255) tcAAACGCAGGCGGCAGTTTGGTTTGTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCga                                     <  1:2152685/68‑1 (MQ=255) tcAAACGCAGGCGGCAGTTTGGTTTGTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCga                                     <  1:2565690/68‑1 (MQ=255) tcAAACGCAGGCGGCAGTTTGGTTTGTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCga                                     <  1:1786186/68‑1 (MQ=255) tcAAACGCAGGCGGCAGTTTGGTTTGTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCga                                     <  1:2924843/68‑1 (MQ=255) tcAAACGCAGGCGGCAGTTTGGTTTGTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCga                                     <  1:2975295/68‑1 (MQ=255) tcAAACGCAGGCGGCAGTTTGGTTTGTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCga                                     <  1:3015583/68‑1 (MQ=255) tcAAACGCAGGCGGCAGTTTGGTTTGTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCga                                     <  1:304726/68‑1 (MQ=255) tcAAACGCAGGCGGCAGTTTGGTTTGTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCga                                     <  1:316721/68‑1 (MQ=255) tcAAACGCAGGCGGCAGTTTGGTTTGTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCga                                     <  1:369670/68‑1 (MQ=255) tcAAACGCAGGCGGCAGTTTGGTTTGTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCga                                     <  1:1280377/68‑1 (MQ=255) tcAAACGCAGGCGGCAGTTTGGTTTGTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCga                                     <  1:691657/68‑1 (MQ=255) tcAAACGCAGGCGGCAGTTTGGTTTGTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCga                                     <  1:125156/68‑1 (MQ=255) tcAAACGCAGGCGGCAGTTTGGTTTGTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCga                                     <  1:766766/68‑1 (MQ=255) tcAAACGCAGGCGGCAGTTTGGTTTGTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCga                                     <  1:828505/68‑1 (MQ=255) tcAAACGCAGGCGGCAGTTTGGTTTGTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCga                                     <  1:1099262/68‑1 (MQ=255) ccAACCGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACcgcgcg                                      >  1:1020956/1‑68 (MQ=255) ccAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACcgcgcg                                      >  1:985242/1‑68 (MQ=255) ccAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACcgcgcg                                      >  1:617181/1‑68 (MQ=255) ccAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACcgcgcg                                      >  1:1691867/1‑68 (MQ=255) ccAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACcgcgcg                                      >  1:1927049/1‑68 (MQ=255)  cAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCga                                     >  1:3250395/1‑68 (MQ=255)  cAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCga                                     >  1:1447303/1‑68 (MQ=255)  cAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCga                                     >  1:295100/1‑68 (MQ=255)  cAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCga                                     >  1:2422431/1‑68 (MQ=255)  cAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCga                                     >  1:1105902/1‑68 (MQ=255)  cAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCga                                     >  1:894977/1‑68 (MQ=255)  cAAACGCAGGCGGCAGTTTGGTTTGTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCga                                     <  1:2215046/68‑1 (MQ=255)      cGCAGGCGGCAGTTTGGTTTGTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCga                                     <  1:139316/64‑1 (MQ=255)         tGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAgttg                             <  1:1955220/69‑1 (MQ=255)         tGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAgttg                             <  1:3248426/69‑1 (MQ=255)         tGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAgttg                             <  1:3350934/69‑1 (MQ=255)         tGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAgttg                             <  1:1039591/69‑1 (MQ=255)         tGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAgttg                             <  1:433340/69‑1 (MQ=255)         tGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAgttg                             <  1:53470/69‑1 (MQ=255)         tGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAgttg                             <  1:2035608/69‑1 (MQ=255)         tGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAgttg                             <  1:2235212/69‑1 (MQ=255)         tGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAgttg                             <  1:284573/69‑1 (MQ=255)         tGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAgttg                             <  1:1292113/69‑1 (MQ=255)         tGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAgttg                             <  1:70221/69‑1 (MQ=255)         tGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAgttg                             <  1:859324/69‑1 (MQ=255)               aaGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACAC‑‑CGCGAGAAAGTTGTTCATCAg                     <  1:2320482/69‑1 (MQ=255)                 gtttggtttgCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAg                     <  1:485786/69‑1 (MQ=255)                 gtttggtttgCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAg                     <  1:1133288/69‑1 (MQ=255)                 gtttggtttgCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAg                     <  1:2086838/69‑1 (MQ=255)                 gtttggtttgCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAg                     <  1:1484959/69‑1 (MQ=255)                 gtttggtttgCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAg                     <  1:1623670/69‑1 (MQ=255)                        ttgTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGaaa                >  1:3287685/1‑67 (MQ=255)                        ttgTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGaaa                >  1:2054610/1‑67 (MQ=255)                        ttgTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGaaa                >  1:909928/1‑67 (MQ=255)                        ttgTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGaaa                >  1:872369/1‑67 (MQ=255)                        ttgTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGaaa                >  1:2253766/1‑67 (MQ=255)                        ttgTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGaaa                >  1:2336932/1‑67 (MQ=255)                        ttgTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGaaa                >  1:1270570/1‑67 (MQ=255)                        ttgTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGaaa                >  1:2709382/1‑67 (MQ=255)                        ttgTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGaaa                >  1:367740/1‑67 (MQ=255)                        ttgTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGaaa                >  1:1436084/1‑67 (MQ=255)                        ttgTGCCGCAGCAATGGCGCGTTCGCCGTAGATCAACACCGCGCGAGAAAGTTGTTCATCAGTGaaa                >  1:130927/1‑67 (MQ=255)                               gcagcaATGGCGCGTTCGCCGTAGATCCACACCGCGCga                                     <  1:3254748/39‑1 (MQ=255)                                   caATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGAAAAAATCGTGCAGGt   <  1:3038155/69‑1 (MQ=255)                                   caATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGAAAAAATCGTGCAGGt   <  1:2212735/69‑1 (MQ=255)                                   caATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGAAAAAATCGTGCAGGt   <  1:2191291/69‑1 (MQ=255)                                   caATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGAAAAAATCGTGCAGGt   <  1:3217356/69‑1 (MQ=255)                                        gcgcGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGAAAAAATCGTGCAGGTg  <  1:658180/65‑1 (MQ=255)                                        gcgcGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGAAAAAATCGTGCAGGTg  <  1:2436802/65‑1 (MQ=255)                                        gcgcGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGAAAAAATCGTGCAGGTg  <  1:954082/65‑1 (MQ=255)                                               |                                                          CCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGAAAAAATCGTGCAGGTG  >  minE/388202‑388305

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