Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	570,909	A→G	43.7%	S216S (TCT→TCC)	ybjE ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	570,909	0	A	G	43.7%	3.1 / 22.0	32	S216S (TCT→TCC)	ybjE	predicted transporter
Reads supporting (aligned to +/- strand):  ref base A (3/15);  new base G (9/5);  total (12/20)
Fisher's exact test for biased strand distribution p-value = 1.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

CAGTTCACGGGCCAGATCATTAAAAAACGCCGCGCTCCCGATTACCGGACCAAAAGATTCGGTCAATAAAATACCGGAAAGAGAATACCAGCCGAAACCGGAGGCC  >  minE/570855‑570960                                                       |                                                    cAGTTCACGGGCCAGATCATTAAAAAACGCCGCGCTCCCGATTACCGGACCAAAAGATTCGGTCAATaa                                       <  1:2752786/69‑1 (MQ=255) cAGTTCACGGGCCAGATCATTAAAAAACGCCGCGCTCCCGATTACCGGACCAAAAGATTCGGTCAATaa                                       <  1:881620/69‑1 (MQ=255) cAGTTCACGGGCCAGATCATTAAAAAACGCCGCGCTCCCGATTACCGGACCAAAAGATTCGGTCAATaa                                       <  1:2226662/69‑1 (MQ=255)      cACGGGCCAGATCATTAAAAAACGCCGCGCTCCCGATTACCGGACCAAAAGATTCGGTCAATAAAATAc                                  <  1:3328294/69‑1 (MQ=255)      cACGGGCCAGATCATTAAAAAACGCCGCGCTCCCGATTACCGGACCAAAAGATTCGGTCAATAAAATAc                                  <  1:36965/69‑1 (MQ=255)      cACGGGCCAGATCATTAAAAAACGCCGCGCTCCCGATTACCGGACCAAAAGATTCGGTCAATAAAATAc                                  <  1:439063/69‑1 (MQ=255)       aCGGGCCAGATGATTAAAAAACGCCGCGCTCCCGATTACCGGACCAAAAGATTCGGTCAATAAAATAc                                  <  1:877428/68‑1 (MQ=255)        cGGGCCAGATCATTAAAAAACGCCGCGCTCCCGATTACCGGACCAAAGGATTCGGTCAATAAAATAcc                                 >  1:2218813/1‑68 (MQ=255)        cGGGCCAGATCATTAAAAAACGCCGCGCTCCCGATTACCGGACCAAAGGATTCGGTCAATAAAATAcc                                 >  1:3125028/1‑68 (MQ=255)        cGGGCCAGATCATTAAAAAACGCCGCGCTCCCGATTACCGGACCAAAGGATTCGGTCAATAAAATAcc                                 >  1:2365821/1‑68 (MQ=255)        cGGGCCAGATCATTAAAAAACGCCGCGCTCCCGATTACCGGACCAAAGGATTCGGTCAATAAAATAcc                                 >  1:2251699/1‑68 (MQ=255)        cGGGCCAGATCATTAAAAAACGCCGCGCTCCCGATTACCGGACCAAAGGATTCGGTCAATAAAATAcc                                 >  1:2149078/1‑68 (MQ=255)        cGGGCCAGATCATTAAAAAACGCCGCGCTCCCGATTACCGGACCAAAGGATTCGGTCAATAAAATAcc                                 >  1:1688233/1‑68 (MQ=255)        cGGGCCAGATCATTAAAAAACGCCGCGCTCCCGATTACCGGACCAAAGGATTCGGTCAATAAAATAcc                                 >  1:1644784/1‑68 (MQ=255)        cGGGCCAGATCATTAAAAAACGCCGCGCTCCCGATTACCGGACCAAAGGATTCGGTCAATAAAATAcc                                 >  1:96754/1‑68 (MQ=255)        cGGGCCAGATCATTAAAAAACGCCGCGCTCCCGATTACCGGACCAAAGGATTCGGTCAATAAAATAcc                                 >  1:1422197/1‑68 (MQ=255)          ggCCAGATCATTAAAAAACGCCGCGCTCCCGATTACCGGACCAAAAGATTCGGTCAATAAAATACCGGa                              <  1:1550803/69‑1 (MQ=255)               gagcATTAAAAAACGCCGCGCTCCCGATTACCGGACCAAAAGATTCGGTCAATAAAATACCGGAAcgag                         <  1:1491388/66‑1 (MQ=255)               gATCATTAAAAAACGCCGCGCTCCCGATTACCGGACCAAAAGATTCGGTCAATAAAATACCGGAAagag                         <  1:2219210/69‑1 (MQ=255)                 tCATTAAAAAACGCCGCGCACCCGATTACCGGACCAAAAGATTCGGTCAATAAAATACCGGAAAGAGaa                       >  1:2018493/1‑69 (MQ=255)                  cATTAAAAAACGCCGCGCTCCCGATTACCGGACCAAAAGATTCGGTCAATAAAATACCGGAAAGAGAAt                      >  1:2663393/1‑69 (MQ=255)                   aTTAAAAAACGCCGCGCTCCCGATTACCGGACCAAAAGATTCGGTCAATAAAATACCGGAAagag                         <  1:2658851/65‑1 (MQ=255)                   aTTAAAAAACGCCGCGCTCCCGATTACCGGACCAAAAGATTCGGTCAATAAAATACCGGAAagag                         <  1:2634662/65‑1 (MQ=255)                              ccgcGCTCCCGATTACCGGACCAAAAGATTCGGTCAATAAAATACCGGAAAGAGAATACCAGCCGAAAc          <  1:298577/69‑1 (MQ=255)                              ccgcGCTCCCGATTACCGGACCAAAAGATTCGGTCAATAAAATACCGGAAAGAGAATACCAGCCGAAAc          <  1:754445/69‑1 (MQ=255)                                gcgcTCCCGATTACCGGACCAAAGGATTCGGTTAATAAAATACCGGa                              <  1:2828515/47‑1 (MQ=255)                                gcgcTCCCGATTACCGGACCAAAGGATTCGGTCAATAAAATACCGGa                              <  1:2747693/47‑1 (MQ=255)                                gcgcTCCCGATTACCGGACCAAAGGATTCGGTCAATAAAATACCGGa                              <  1:3064758/47‑1 (MQ=255)                                gcgcTCCCGATTACCGGACCAAAGGATTCGGTCAATAAAATACCGGa                              <  1:2387333/47‑1 (MQ=255)                                gcgcTCCCGATTACCGGACCAAAGGATTCGGTCAATAAAATACCGGa                              <  1:2034823/47‑1 (MQ=255)                                gcgcTCCCGATTACCGGACCAAAAGATTCGGTCAATaa                                       >  1:2400857/1‑38 (MQ=255)                                        gATTACCGGACCAAAAGATTCGGTCAATAAAATACCGGAAAGAGAATACCAGCCGAAACCGGAGGcc  <  1:1569831/67‑1 (MQ=255)                                                       |                                                    CAGTTCACGGGCCAGATCATTAAAAAACGCCGCGCTCCCGATTACCGGACCAAAAGATTCGGTCAATAAAATACCGGAAAGAGAATACCAGCCGAAACCGGAGGCC  >  minE/570855‑570960

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