Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	412,090	A→G	59.5%	V89A (GTG→GCG)	leuS ←	leucyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	412,090	0	A	G	59.5%	22.6 / 40.9	42	V89A (GTG→GCG)	leuS	leucyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (9/8);  new base G (14/11);  total (23/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

TTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTG  >  minE/412028‑412160                                                               |                                                                       ttttCAGATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTc                                                                 >  1:737738/1‑70 (MQ=255)   ttCATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCg                                                               <  1:498557/70‑1 (MQ=255)    tCATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCg                                                               >  1:2713103/1‑69 (MQ=255)      atatACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCgcagg                                                           <  1:4687293/71‑1 (MQ=255)         taCGCGATGTTGTCTTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTc                                                                 <  1:4714640/62‑1 (MQ=255)           cgcgATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCGCAGGCAGAc                                                      <  1:4475258/71‑1 (MQ=255)              gATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCGCAGGCAGACCaa                                                   <  1:187431/71‑1 (MQ=255)              gATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCGCAGGCAGACCa                                                    >  1:1670534/1‑70 (MQ=255)               aTGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCg                                                               <  1:1384771/58‑1 (MQ=255)                 gttgtCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCGCAGGCAGACCAAAcg                                                >  1:3148891/1‑71 (MQ=255)                 gttgtCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCGCAGGCAGACCAAAcg                                                >  1:340576/1‑71 (MQ=255)                  ttgtCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAAcgc                                               >  1:1194897/1‑71 (MQ=255)                  ttgtCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAAcg                                                <  1:3981078/70‑1 (MQ=255)                  ttgtCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAAcg                                                <  1:31688/70‑1 (MQ=255)                   tgtCGTACGTCCACGGTGCCGGAGCGGTTTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAAcg                                                <  1:3103002/69‑1 (MQ=255)                    gtCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCgcaggca                                                         >  1:1546089/1‑59 (MQ=255)                    gtCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCgcaggca                                                         >  1:4899165/1‑59 (MQ=255)                    gtCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCgcaggca                                                         >  1:3249420/1‑59 (MQ=255)                    gtCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCgcaggca                                                         >  1:3033253/1‑59 (MQ=255)                    gtCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCgcaggca                                                         >  1:3193693/1‑59 (MQ=255)                    gtCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCgcaggca                                                         >  1:1734484/1‑59 (MQ=255)                       gTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCgca                                                             <  1:4082168/52‑1 (MQ=255)                       gTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCgca                                                             <  1:4080263/52‑1 (MQ=255)                          cGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGc                                       >  1:3179455/1‑71 (MQ=255)                              cACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGAt                                   >  1:99393/1‑71 (MQ=255)                              cACGGTGCCGGAG‑GGTGTTGTTTTTCACCGCCGCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATc                                  >  1:926991/1‑71 (MQ=255)                               aCGGTGCCGGAGCGGTGTTGTTTTTCACCGCGGCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATc                                  >  1:505171/1‑71 (MQ=255)                                cGGTGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCg                                 <  1:83354/71‑1 (MQ=255)                                  gTGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGc                               >  1:4694087/1‑71 (MQ=255)                                  gTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGc                               >  1:5167527/1‑71 (MQ=255)                                   tGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTcc                                           <  1:1935059/58‑1 (MQ=255)                                     ccGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTg                             <  1:2639676/70‑1 (MQ=255)                                     ccGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTg                             <  1:4601966/70‑1 (MQ=255)                                                  gttTTTCACCGCCGCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGtt                    <  1:5463912/66‑1 (MQ=255)                                                  gttTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGc               >  1:1351614/1‑71 (MQ=255)                                                  gttTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGc               >  1:555186/1‑71 (MQ=255)                                                  gttTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGc               >  1:2437607/1‑71 (MQ=255)                                                    ttttCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCg                                     <  1:1213575/47‑1 (MQ=255)                                                       tCACCGCCGCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAg           <  1:144348/70‑1 (MQ=255)                                                       tCACCGCCGCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAg           >  1:1456907/1‑70 (MQ=255)                                                         aCCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCa         >  1:10218/1‑70 (MQ=255)                                                           cgccgcGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATa       <  1:5640122/70‑1 (MQ=255)                                                             ccACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGc    <  1:212344/71‑1 (MQ=255)                                                               aCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTg  >  1:240601/1‑71 (MQ=255)                                                               |                                                                       TTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTG  >  minE/412028‑412160

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