Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	412,090	A→G	22.9%	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	22.9%	60.7 / 17.5	35	V89A (GTG→GCG)	leuS	leucyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (16/11);  new base G (6/2);  total (22/13)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

GGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCT  >  minE/412025‑412159                                                                  |                                                                      ggTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCt                                                                  >  1:3759542/1‑71 (MQ=255)    ttttCATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTcc                                                               >  1:2638367/1‑71 (MQ=255)     tttCATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCg                                                              >  1:3820476/1‑71 (MQ=255)        cATATACGCGTTGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCgca                                                            <  1:2791936/70‑1 (MQ=255)        cATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCgca                                                            <  1:901488/70‑1 (MQ=255)            taCGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCgcagg                                                          >  1:377990/1‑68 (MQ=255)              cgcgATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCgcaggcagcc                                                     >  1:986798/1‑69 (MQ=255)                 gATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCGCAGGCAGACCa                                                   <  1:3694827/70‑1 (MQ=255)                  aTGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGcc                                                                   <  1:544225/53‑1 (MQ=255)                  aTGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCaa                                                  >  1:364127/1‑70 (MQ=255)                   tgttgtCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCGCAGGCAGACCaaa                                                 >  1:2622561/1‑70 (MQ=255)                   tgttgtCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAAc                                                <  1:226503/71‑1 (MQ=255)                    gttgtCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCGCAGGCAGACCaa                                                  >  1:1585193/1‑68 (MQ=255)                     ttgtCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCa                                                                       >  1:960508/1‑46 (MQ=255)                       gtCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGt                                            >  1:4079000/1‑71 (MQ=255)                          gTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTcc                                          >  1:1307477/1‑70 (MQ=255)                           tACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCa                                        <  1:3754659/71‑1 (MQ=255)                             cGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAg                                       <  1:485444/70‑1 (MQ=255)                                 cACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGAt                                  >  1:3753283/1‑71 (MQ=255)                                      tGCCGGAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCt                             >  1:2765148/1‑71 (MQ=255)                                               ggTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGGTCGGCTGCAGGACGt                    <  1:3129167/71‑1 (MQ=255)                                                gtgtTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGtt                   <  1:2549838/71‑1 (MQ=255)                                                  gttgttTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGtttt                 >  1:1386500/1‑71 (MQ=255)                                                   ttgttTTTCACCGCCACGCCTTCCGCAGGCAGACCAAAc                                                <  1:889248/39‑1 (MQ=255)                                                      tttttCACCGCCGCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGcc             <  1:2568458/71‑1 (MQ=255)                                                       ttttCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGcc             >  1:2185482/1‑70 (MQ=255)                                                       ttttCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGcc             >  1:2789775/1‑70 (MQ=255)                                                       ttttCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGcc             >  1:2388905/1‑70 (MQ=255)                                                        tttCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGccc            <  1:1919531/70‑1 (MQ=255)                                                        tttCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGcc             >  1:315636/1‑69 (MQ=255)                                                        tttCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGcc             >  1:2117940/1‑69 (MQ=255)                                                        tttCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGcc             >  1:1076850/1‑69 (MQ=255)                                                         ttCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGAt                                  >  1:1762420/1‑47 (MQ=255)                                                         ttCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCa           >  1:2758207/1‑70 (MQ=255)                                                          tCACCGCCGCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGCCGTTTTTGCCCAg          >  1:196089/1‑70 (MQ=255)                                                                 cACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCt  <  1:392691/71‑1 (MQ=255)                                                                  |                                                                      GGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCT  >  minE/412025‑412159

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