Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,266,204	T→A	13.9%	V80V (GTA→GTT)	tufB ←	protein chain elongation factor EF‑Tu

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,266,204	0	T	A	13.9%	89.3 / 8.1	36	V80V (GTA→GTT)	tufB	protein chain elongation factor EF‑Tu
Reads supporting (aligned to +/- strand):  ref base T (9/22);  new base A (3/2);  total (12/24)
Fisher's exact test for biased strand distribution p-value = 3.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.24e-01

TCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGG  >  minE/2266137‑2266222                                                                    |                   tCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTAc                  >  1:1514593/1‑70 (MQ=35) tCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACg                 >  1:3369030/1‑71 (MQ=34) tCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACg                 >  1:4584516/1‑71 (MQ=34) tCGCGCCGTCCATCTGCGCAGCACAGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACg                 >  1:2615751/1‑71 (MQ=21)      ccGTCCATCTGCGCAGCACCGGTGATCATGTTTTCAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGc            >  1:2021762/1‑71 (MQ=21)       cggccATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCg           <  1:1234672/68‑1 (MQ=21)       cGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCg           <  1:1024165/71‑1 (MQ=34)       cGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCg           <  1:2576480/71‑1 (MQ=34)       cGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCg           <  1:3536838/71‑1 (MQ=34)       cGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCg           <  1:5083395/71‑1 (MQ=34)        gTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTTCCCCGGGCAGTCTACGTGTGCGt          <  1:116911/71‑1 (MQ=21)        gTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGt          <  1:1005779/71‑1 (MQ=34)        gTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGt          <  1:758690/71‑1 (MQ=34)         tCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTa         <  1:4593377/71‑1 (MQ=34)          ccATCTGCGCAGCACCGGTAATCATGTTCTTGATGTAGTCGGCGTGGCCTGGGCAGTCAACGTGCGCGTa         >  1:2323996/1‑70 (MQ=2)          ccATCTGCGCAGCACCGGTAATCATGTTCTTGATGTAGTCGGCGTGGCCTGGGCAGTCAACGTGCGCGTa         >  1:2374542/1‑70 (MQ=2)          ccATCTGCGCAGCACCGGTAATCATGTTCTTGATGTAGTCGGCGTGGCCTGGGCAGTCAACGTGCGCGTAg        >  1:4104371/1‑71 (MQ=2)           cATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAg        <  1:1112507/70‑1 (MQ=35)           cATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAg        <  1:2274659/70‑1 (MQ=35)           cATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGt       >  1:4327106/1‑71 (MQ=34)           cATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGt       >  1:2603834/1‑71 (MQ=34)           cATCTGCGCAGCACCGGTAATCATGTTCTTGATGTAGTCGGCGTGGCCTGGGCAGTCAACGTGAGCGTAGt       <  1:1730082/71‑1 (MQ=2)           cATCTGCGCAGCACCGGTAATCATGTTCTTGATGTAGTCGGCGTGGCCTGGGCAGTCAACGTGAGCGTAGt       <  1:4400135/71‑1 (MQ=2)             tCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGt       <  1:2374937/69‑1 (MQ=35)             tCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGt       <  1:4568808/69‑1 (MQ=35)             tCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGt       <  1:1465794/69‑1 (MQ=35)             tCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGt       <  1:1075632/69‑1 (MQ=35)             tCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGt       <  1:770768/69‑1 (MQ=35)             tCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGt       <  1:88499/69‑1 (MQ=35)             tCTGCGCAGCACCGGTGATCATGTGTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGt       <  1:738858/69‑1 (MQ=25)              cTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGt       <  1:3748934/68‑1 (MQ=35)               tGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTGTACGTGTGCGTAGTGACg   >  1:2382090/1‑71 (MQ=21)               tGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACg   >  1:2796724/1‑71 (MQ=34)                 cgcAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACgg  <  1:3849002/70‑1 (MQ=31)                 cgcAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACgg  <  1:1102974/70‑1 (MQ=31)                 cgcAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACgg  <  1:1077606/70‑1 (MQ=31)                                                                    |                   TCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGG  >  minE/2266137‑2266222

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