Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	410,370	A→G	60.9%	G662G (GGT→GGC)	leuS ←	leucyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	410,370	0	A	G	60.9%	21.1 / 27.4	46	G662G (GGT→GGC)	leuS	leucyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (10/8);  new base G (20/8);  total (30/16)
Fisher's exact test for biased strand distribution p-value = 3.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.37e-01

GTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCAAACATCATAAACA  >  minE/410314‑410425                                                         |                                                        gtgtgCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTCCTGc                                             <  1:322762/69‑1 (MQ=255)    tgCTCGTAAACCAGTTTCCAGACACGTTTCAGGAGGCGGTTAGCCCCTTCCACACCGGATTCCTGCCAt                                          >  1:2182725/1‑69 (MQ=255)    tgCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTCCTGCCAt                                          >  1:1386713/1‑69 (MQ=255)    tgCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTCCTGCCAt                                          >  1:2942226/1‑69 (MQ=255)    tgCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTCCTGCCAt                                          >  1:909814/1‑69 (MQ=255)                    tCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGAt                                                   >  1:3251271/1‑44 (MQ=255)                    tCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGAt                                                   >  1:805770/1‑44 (MQ=255)                    tCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGAt                                                   <  1:800046/44‑1 (MQ=255)                    tCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGAt                                                   <  1:1787496/44‑1 (MQ=255)                    tCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGAt                                                   <  1:1882417/44‑1 (MQ=255)                    tCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGAt                                                   >  1:500270/1‑44 (MQ=255)                        gACACGTTTCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCgg                      <  1:2273883/69‑1 (MQ=255)                        gACACGTTTCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCgg                      <  1:1001407/69‑1 (MQ=255)                        gACACGTTTCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCgg                      <  1:630420/69‑1 (MQ=255)                        gACACGTTTCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCgg                      <  1:776119/69‑1 (MQ=255)                        gACACGTTTCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCgg                      <  1:2367564/69‑1 (MQ=255)                        gACACGTTTCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCgg                      <  1:1068114/69‑1 (MQ=255)                        gACACGTTTCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCgg                      <  1:1886955/69‑1 (MQ=255)                        gACACGTTTCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCgg                      <  1:1281444/69‑1 (MQ=255)                               ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGc                >  1:2814417/1‑68 (MQ=255)                               ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGc                >  1:989278/1‑68 (MQ=255)                               ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCa               >  1:1035246/1‑69 (MQ=255)                               ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCa               >  1:605115/1‑69 (MQ=255)                               ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCa               >  1:499764/1‑69 (MQ=255)                               ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCa               >  1:949381/1‑69 (MQ=255)                               ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCa               >  1:3381766/1‑69 (MQ=255)                               ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCa               >  1:332800/1‑69 (MQ=255)                               ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCa               >  1:2776140/1‑69 (MQ=255)                               ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCa               >  1:1273758/1‑69 (MQ=255)                               ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCa               >  1:1401039/1‑69 (MQ=255)                               ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCa               >  1:1633610/1‑69 (MQ=255)                               ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCa               >  1:1734053/1‑69 (MQ=255)                               ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCa               >  1:247641/1‑69 (MQ=255)                               ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCa               >  1:2696075/1‑69 (MQ=255)                               ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCa               >  1:3160128/1‑69 (MQ=255)                               ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCa               >  1:279146/1‑69 (MQ=255)                               ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCa               >  1:3065416/1‑69 (MQ=255)                               ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATAAGCCGGAGAAGCa               >  1:485440/1‑69 (MQ=255)                               atCAGGAAGCGGTTAGCCCCTTCAACGCCGGATTCCTGCCATTAGAGAGTCATATCAGCCGgaga                   >  1:1393398/2‑65 (MQ=255)                                      aGCGGTTAGCCCCTTCCACACCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCAAAcatc         >  1:3051133/1‑68 (MQ=255)                                            tAGCCCCTTCCACACCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCAAACATCATAAACa  <  1:2814055/69‑1 (MQ=255)                                            tAGCCCCTTCCACACCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCAAACATCATAAACa  <  1:3049389/69‑1 (MQ=255)                                            tAGCCCCTTCCACACCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCAAACATCATAAACa  <  1:809209/69‑1 (MQ=255)                                             aGCCCCTTCCACACCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCAAACATCATAAACa  <  1:3137525/68‑1 (MQ=255)                                                 ccTTCCACACCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCAAACATCATa      >  1:716907/1‑60 (MQ=255)                                                 ccTTCCACACCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCAAACATCATa      >  1:3353661/1‑60 (MQ=255)                                                         |                                                        GTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCAAACATCATAAACA  >  minE/410314‑410425

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