Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,911,643	A→G	57.4%	L184L (TTG→CTG)	thyA ←	thymidylate synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,911,643	0	A	G	57.4%	25.3 / 54.3	75	L184L (TTG→CTG)	thyA	thymidylate synthetase
Reads supporting (aligned to +/- strand):  ref base A (7/25);  new base G (26/17);  total (33/42)
Fisher's exact test for biased strand distribution p-value = 1.04e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.30e-01

TCGCCACCGGTCCAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTCACAGGAGCGCTG  >  minE/1911579‑1911700                                                                 |                                                          tCGCCACCGGTCCAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAAc                                                       >  1:2263060/1‑69 (MQ=255)            ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:1055930/69‑1 (MQ=255)            ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:991476/69‑1 (MQ=255)            ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:887183/69‑1 (MQ=255)            ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:855503/69‑1 (MQ=255)            ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:581154/69‑1 (MQ=255)            ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:554783/69‑1 (MQ=255)            ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:389823/69‑1 (MQ=255)            ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:3403750/69‑1 (MQ=255)            ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:329587/69‑1 (MQ=255)            ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:3186362/69‑1 (MQ=255)            ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:3059240/69‑1 (MQ=255)            ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:2996489/69‑1 (MQ=255)            ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:2439183/69‑1 (MQ=255)            ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:1933146/69‑1 (MQ=255)            ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:1186381/69‑1 (MQ=255)            ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:1856346/69‑1 (MQ=255)            ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:1089503/69‑1 (MQ=255)            ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:124583/69‑1 (MQ=255)            ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:1744828/69‑1 (MQ=255)                 cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAAcgc                                                     >  1:53653/1‑55 (MQ=255)                 cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCt                                               >  1:145855/1‑61 (MQ=255)                 cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:3038999/1‑67 (MQ=255)                 cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:2654296/1‑67 (MQ=255)                 cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:515943/1‑67 (MQ=255)                 cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:2885685/1‑67 (MQ=255)                 cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:438246/1‑67 (MQ=255)                 cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:1742349/1‑67 (MQ=255)                 cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:3003066/1‑67 (MQ=255)                 cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:3344377/1‑67 (MQ=255)                 cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:3399407/1‑67 (MQ=255)                 cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:3156390/1‑67 (MQ=255)                 cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:3170550/1‑67 (MQ=255)                 cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:3288111/1‑67 (MQ=255)                 cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:57879/1‑67 (MQ=255)                 cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:1022199/1‑67 (MQ=255)                 cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:1175127/1‑67 (MQ=255)                 cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:2490957/1‑67 (MQ=255)                 cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:7202/1‑67 (MQ=255)                 cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:2313026/1‑67 (MQ=255)                 cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:2278318/1‑67 (MQ=255)                 cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:1075516/1‑67 (MQ=255)                 cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:22576/1‑67 (MQ=255)                 cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:2123292/1‑67 (MQ=255)                 cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:2081541/1‑67 (MQ=255)                 cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:2038104/1‑67 (MQ=255)                                   gATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAgg                      <  1:260018/68‑1 (MQ=255)                                   gATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAgg                      <  1:1858705/68‑1 (MQ=255)                                   gATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAgg                      <  1:3175695/68‑1 (MQ=255)                                   gATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAgg                      <  1:1798452/68‑1 (MQ=255)                                   gATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAgg                      <  1:2748990/68‑1 (MQ=255)                                    aTCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGaa                    <  1:2328934/69‑1 (MQ=255)                                          ctgctgCGCCATCATATGCACCAGTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTc              <  1:522314/69‑1 (MQ=255)                                          ctgctgCGCCATCATATGCACCAGTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTc              <  1:3066113/69‑1 (MQ=255)                                          ctgctgCGCCATCATATGCACCAGTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTc              <  1:1426607/69‑1 (MQ=255)                                          ctgctgCGCCATCATATGCACCAGTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTc              <  1:185740/69‑1 (MQ=255)                                          ctgctgCGCCATCATATGCACCAGTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTc              <  1:181693/69‑1 (MQ=255)                                          ctgctgCGCCATCATATGCACCAGTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTc              <  1:2136197/69‑1 (MQ=255)                                          ctgctgCGCCATCATATGCACCAGTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTc              <  1:249539/69‑1 (MQ=255)                                          ctgctgCGCCATCATATGCACCAGTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTc              <  1:2610362/69‑1 (MQ=255)                                          ctgctgCGCCATCATATGCACCAGTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTc              <  1:2775021/69‑1 (MQ=255)                                          ctgctgCGCCATCATATGCACCAGTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTc              <  1:2824048/69‑1 (MQ=255)                                          ctgctgCGCCATCATATGCACCAGTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTc              <  1:2868652/69‑1 (MQ=255)                                          ctgctgCGCCATCATATGCACCAGTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTc              <  1:2936344/69‑1 (MQ=255)                                          ctgctgCGCCATCATATGCACCAGTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTc              <  1:3101642/69‑1 (MQ=255)                                          ctgctgCGCCATCATATGCACCAGTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTc              <  1:3315550/69‑1 (MQ=255)                                          ctgctgCGCCATCATATGCACCAATAACGCGTATCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGt               >  1:543863/1‑68 (MQ=255)                                           tgctgcGCCATCATATGCACCAGTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTc              <  1:526036/68‑1 (MQ=255)                                           tgctgcGCCATCATATGCACCAGTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTc              <  1:2454815/68‑1 (MQ=255)                                                       atatGCACCAGTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTc              <  1:817953/56‑1 (MQ=255)                                                       atatGCACCAATAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTCACAGGAGCGCTg  >  1:530141/1‑68 (MQ=255)                                                       atatGCACCAATAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTCACAGGAGCGCTg  >  1:496938/1‑68 (MQ=255)                                                       atatGCACCAATAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTCACAGGAGCGCTg  >  1:566015/1‑68 (MQ=255)                                                       atatGCACCAATAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTCACAGGAGCGCTg  >  1:619990/1‑68 (MQ=255)                                                       atatGCACCAATAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTCACAGGAGCGCTg  >  1:1409552/1‑68 (MQ=255)                                                                aaTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAg                   <  1:551801/42‑1 (MQ=255)                                                                aaTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAg                   <  1:2654113/42‑1 (MQ=255)                                                                aaTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAg                   <  1:1612755/42‑1 (MQ=255)                                                                 |                                                          TCGCCACCGGTCCAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTCACAGGAGCGCTG  >  minE/1911579‑1911700

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