Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	388,228	C→T	35.9%	A49A (GCG→GCA)	ybdH ←	predicted oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	388,228	0	C	T	35.9%	58.1 / 49.8	86	A49A (GCG→GCA)	ybdH	predicted oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (28/27);  new base T (11/20);  total (39/47)
Fisher's exact test for biased strand distribution p-value = 1.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGAAA  >  minE/388167‑388291                                                              |                                                                gACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCgca                                                           >  1:1928372/1‑68 (MQ=255) gACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCgca                                                           >  1:3185159/1‑68 (MQ=255)  aCCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCgcagc                                                         >  1:2741851/1‑69 (MQ=255)  aCCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCgcagc                                                         <  1:163722/69‑1 (MQ=255)  aCCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCgcagc                                                         >  1:521918/1‑69 (MQ=255)  aCCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCgcagc                                                         >  1:525870/1‑69 (MQ=255)  aCCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCgcagc                                                         >  1:2136198/1‑69 (MQ=255)  aCCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCgcagc                                                         <  1:3038406/69‑1 (MQ=255)  aCCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCgcagc                                                         >  1:2486192/1‑69 (MQ=255)  aCCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCgcagc                                                         >  1:2520576/1‑69 (MQ=255)  aCCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCgcagc                                                         >  1:2640010/1‑69 (MQ=255)  aCCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCgcagc                                                         <  1:1226774/69‑1 (MQ=255)  aCCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCgcagc                                                         >  1:85434/1‑69 (MQ=255)  aCCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCgcagc                                                         >  1:897042/1‑69 (MQ=255)  aCCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCgcagc                                                         >  1:270759/1‑69 (MQ=255)  aCCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCgcagc                                                         >  1:1091459/1‑69 (MQ=255)  aCCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCgcagc                                                         >  1:10509/1‑69 (MQ=255)                  tGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGgcgc                                                 >  1:2472442/1‑61 (MQ=255)                  tGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGcgcg                                                <  1:81345/62‑1 (MQ=255)                  tGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGcgcg                                                <  1:2139727/62‑1 (MQ=255)                  tGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGcgcg                                                >  1:156362/1‑62 (MQ=255)                  tGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGcgcg                                                <  1:491254/62‑1 (MQ=255)                  tGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGcgcg                                                <  1:484020/62‑1 (MQ=255)                  tGCTTTGCCCCTGGAAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGcgcg                                                >  1:486383/1‑62 (MQ=255)                                   tCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGcgcg                                                <  1:778118/45‑1 (MQ=255)                                    tcAAACGCAGGCGGCAGTTTGGTTTGTGCCGCAGCAATGGCGCTTTCGCCGTAGATCCACACCGCGCga                       <  1:1812832/68‑1 (MQ=255)                                    tcAAACGCAGGCGGCAGTTTGGTTTGTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCga                       <  1:766766/68‑1 (MQ=255)                                    tcAAACGCAGGCGGCAGTTTGGTTTGTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCga                       <  1:1099262/68‑1 (MQ=255)                                    tcAAACGCAGGCGGCAGTTTGGTTTGTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCga                       <  1:125156/68‑1 (MQ=255)                                    tcAAACGCAGGCGGCAGTTTGGTTTGTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCga                       <  1:2565690/68‑1 (MQ=255)                                    tcAAACGCAGGCGGCAGTTTGGTTTGTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCga                       <  1:828505/68‑1 (MQ=255)                                    tcAAACGCAGGCGGCAGTTTGGTTTGTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCga                       <  1:1280377/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:2152685/68‑1 (MQ=255)                                    tcAAACGCAGGCGGCAGTTTGGTTTGTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCga                       <  1:691657/68‑1 (MQ=255)                                    tcAAACGCAGGCGGCAGTTTGGTTTGTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCga                       <  1:2114552/68‑1 (MQ=255)                                    tcAAACGCAGGCGGCAGTTTGGTTTGTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCga                       <  1:1786186/68‑1 (MQ=255)                                    tcAAACGCAGGCGGCAGTTTGGTTTGTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCga                       <  1:344117/68‑1 (MQ=255)                                    tcAAACGCAGGCGGCAGTTTGGTTTGTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCga                       <  1:369670/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:1927049/1‑68 (MQ=255)                                    ccAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACcgcgcg                        >  1:1691867/1‑68 (MQ=255)                                     cAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCga                       >  1:2422431/1‑68 (MQ=255)                                     cAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCga                       >  1:295100/1‑68 (MQ=255)                                     cAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCga                       >  1:1105902/1‑68 (MQ=255)                                     cAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCga                       >  1:894977/1‑68 (MQ=255)                                     cAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCga                       >  1:3250395/1‑68 (MQ=255)                                     cAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCga                       >  1:1447303/1‑68 (MQ=255)                                     cAAACGCAGGCGGCAGTTTGGTTTGTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCga                       <  1:2215046/68‑1 (MQ=255)                                         cGCAGGCGGCAGTTTGGTTTGTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCga                       <  1:139316/64‑1 (MQ=255)                                           cTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGcgcg                                                <  1:3435776/37‑1 (MQ=255)                                            tGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAgttg               <  1:1955220/69‑1 (MQ=255)                                            tGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAgttg               <  1:284573/69‑1 (MQ=255)                                            tGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAgttg               <  1:53470/69‑1 (MQ=255)                                            tGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAgttg               <  1:70221/69‑1 (MQ=255)                                            tGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAgttg               <  1:1292113/69‑1 (MQ=255)                                            tGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAgttg               <  1:859324/69‑1 (MQ=255)                                            tGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAgttg               <  1:1039591/69‑1 (MQ=255)                                            tGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAgttg               <  1:2035608/69‑1 (MQ=255)                                            tGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAgttg               <  1:2235212/69‑1 (MQ=255)                                            tGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAgttg               <  1:433340/69‑1 (MQ=255)                                            tGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAgttg               <  1:3248426/69‑1 (MQ=255)                                            tGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAgttg               <  1:3350934/69‑1 (MQ=255)                                                  aaGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACAC‑‑CGCGAGAAAGTTGTTCATCAg       <  1:2320482/69‑1 (MQ=255)                                                    gtttggtttgCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAg       <  1:2086838/69‑1 (MQ=255)                                                    gtttggtttgCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAg       <  1:1133288/69‑1 (MQ=255)                                                    gtttggtttgCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAg       <  1:1484959/69‑1 (MQ=255)                                                    gtttggtttgCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAg       <  1:1623670/69‑1 (MQ=255)                                                    gtttggtttgCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAg       <  1:485786/69‑1 (MQ=255)                                                           ttgTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGaaa  >  1:3287685/1‑67 (MQ=255)                                                           ttgTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGaaa  >  1:1436084/1‑67 (MQ=255)                                                           ttgTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGaaa  >  1:2054610/1‑67 (MQ=255)                                                           ttgTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGaaa  >  1:1270570/1‑67 (MQ=255)                                                           ttgTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGaaa  >  1:367740/1‑67 (MQ=255)                                                           ttgTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGaaa  >  1:2253766/1‑67 (MQ=255)                                                           ttgTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGaaa  >  1:872369/1‑67 (MQ=255)                                                           ttgTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGaaa  >  1:2336932/1‑67 (MQ=255)                                                           ttgTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGaaa  >  1:909928/1‑67 (MQ=255)                                                           ttgTGCCGCAGCAATGGCGCGTTCGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGaaa  >  1:2709382/1‑67 (MQ=255)                                                           ttgTGCCGCAGCAATGGCGCGTTCGCCGTAGATCAACACCGCGCGAGAAAGTTGTTCATCAGTGaaa  >  1:130927/1‑67 (MQ=255)                                                              |                                                                GACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGAAA  >  minE/388167‑388291

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