Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,814,577	G→A	60.8%	S200S (AGC→AGT)	alaS ←	alanyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,814,577	0	G	A	60.8%	7.8 / 15.5	28	S200S (AGC→AGT)	alaS	alanyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (4/7);  new base A (10/7);  total (14/14)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

GGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGGCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTG  >  minE/1814516‑1814631                                                              |                                                       ggTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCg                                                 <  1:1304516/69‑1 (MQ=255)    tGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGACTTCCCGGAg                                              <  1:1486315/69‑1 (MQ=255)    tGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGACTTCCCGGAg                                              <  1:830567/69‑1 (MQ=255)    tGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGACTTCCCGGAg                                              <  1:667629/69‑1 (MQ=255)    tGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGACTTCCCGGAg                                              <  1:2155763/69‑1 (MQ=255)    tGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGACTTCCCGGAg                                              <  1:2304017/69‑1 (MQ=255)    tGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGACTTCCCGGAg                                              <  1:2455300/69‑1 (MQ=255)    tGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGACTTCCCGGAg                                              <  1:2587281/69‑1 (MQ=255)              gAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGACTTCCCGGAGGGCCGCCCCa                                    >  1:119489/1‑69 (MQ=255)              gAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGACTTCCCGGAGGGCCGCCCCa                                    >  1:842899/1‑69 (MQ=255)              gAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGACTTCCCGGAGGGCCGCCCCa                                    >  1:507903/1‑69 (MQ=255)              gAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGACTTCCCGGAGGGCCGCCCCa                                    >  1:369809/1‑69 (MQ=255)              gAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGACTTCCCGGAGGGCCGCCCCa                                    >  1:2926632/1‑69 (MQ=255)              gAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGACTTCCCGGAGGGCCGCCCCa                                    >  1:2815415/1‑69 (MQ=255)              gAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGACTTCCCGGAGGGCCGCCCCa                                    >  1:1842692/1‑69 (MQ=255)              gAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGACTTCCCGGAGGGCCGCCCCa                                    >  1:1829359/1‑69 (MQ=255)              gAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGACTTCCCGGAGGGCCGCCCCa                                    >  1:1680388/1‑69 (MQ=255)              gAAGACGATGTTCCAGATCTCAATGTAGAGGTCGCCGTCTTCTTCCGGACTTCCCGGAGGGCCGCCCCa                                    >  1:2846734/1‑69 (MQ=255)                    gATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGGCCCCCCCAAAtgtg                              <  1:1152706/69‑1 (MQ=255)                    gATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGGCCCCCCCAAAtgtg                              <  1:960707/69‑1 (MQ=255)                               tctcAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGGCCCCCCCAAATGTGGTCGCCGTGAt                   >  1:168606/1‑69 (MQ=255)                               tctcAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGGCCCCCCCAAATGTGGTCGCCGTGAt                   >  1:1667361/1‑69 (MQ=255)                                ctcAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGGCCCCCCCAAATGTGGTCGCCGt                      >  1:764268/1‑65 (MQ=255)                                ctcAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGGCCCCCCCAAATGTGGTCGCCGTGAt                   >  1:2865041/1‑68 (MQ=255)                                       tAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGGCCCCCCCAAATGTGGTCGGCGTGATCGTagaag           <  1:2817651/69‑1 (MQ=255)                                          cGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGGCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGAt         <  1:2234810/68‑1 (MQ=255)                                                 cGTCTTCTTCCGGGCTTCCCGGAGGGCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTg  <  1:3065896/68‑1 (MQ=255)                                                 cGTCTTCTTCCGGGCTTCCCGGAGGGCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTg  <  1:62496/68‑1 (MQ=255)                                                              |                                                       GGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGGCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTG  >  minE/1814516‑1814631

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