Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,246,341	A→T	29.4%	L53L (CTT→CTA)	znuA ←	zinc transporter subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,246,341	0	A	T	29.4%	60.7 / 30.5	51	L53L (CTT→CTA)	znuA	zinc transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (10/26);  new base T (2/13);  total (12/39)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01

GTTCTGTAAGCGTTTTACATCCGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGCAATGGCAGAAGCGATGAACCC  >  minE/1246278‑1246394                                                                |                                                      gTTCTGTAAGCGTTTTACATCCGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTaa                                                  <  1:1004677/69‑1 (MQ=255) gTTCTGTAAGCGTTTTACATCCGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTaa                                                  <  1:1666979/69‑1 (MQ=255) gTTCTGTAAGCGTTTTACATCCGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTaa                                                  <  1:2257861/69‑1 (MQ=255) gTTCTGTAAGCGTTTTACATCCGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTaa                                                  <  1:121197/69‑1 (MQ=255)             ttttACATCCGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCtctgtttc                                      <  1:1880761/69‑1 (MQ=255)                 aCATCCGATGGGCGCAGTGAATAATCATGTTTTGAAGCGCCGTCAGGAAGTAAAACCtctgtttctgt                                   <  1:2668689/68‑1 (MQ=255)                 aCATCCGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCtctgtttctgt                                   <  1:2064587/68‑1 (MQ=255)                     ccGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTAccc                              <  1:1220211/69‑1 (MQ=255)                     ccGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTAccc                              <  1:2307036/69‑1 (MQ=255)                      cGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTAcccc                             >  1:2717327/1‑69 (MQ=255)                      cGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTAcccc                             >  1:1446423/1‑69 (MQ=255)                             cgcAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCGATCAGc                       <  1:439987/68‑1 (MQ=255)                             cgcAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:2558948/68‑1 (MQ=255)                             cgcAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:2391271/68‑1 (MQ=255)                             cgcAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:2212824/68‑1 (MQ=255)                             cgcAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:328388/68‑1 (MQ=255)                             cgcAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:2194651/68‑1 (MQ=255)                             cgcAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:888729/68‑1 (MQ=255)                             cgcAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:1325376/68‑1 (MQ=255)                             cgcAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:1160504/68‑1 (MQ=255)                             cgcAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:1038728/68‑1 (MQ=255)                                   gAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:19539/62‑1 (MQ=255)                                   gAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:627760/62‑1 (MQ=255)                                    aataatCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:3394776/61‑1 (MQ=255)                                      taatCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGCa                      >  1:188292/1‑60 (MQ=255)                                      taatCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGCa                      >  1:1019831/1‑60 (MQ=255)                                      taatCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGCa                      >  1:646836/1‑60 (MQ=255)                                      taatCAAGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGCa                      >  1:1022/1‑60 (MQ=255)                                        atCATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:2045452/69‑1 (MQ=255)                                        atCATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:842386/69‑1 (MQ=255)                                        atCATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:838085/69‑1 (MQ=255)                                        atCATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:721927/69‑1 (MQ=255)                                        atCATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:691073/69‑1 (MQ=255)                                        atCATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:2875961/69‑1 (MQ=255)                                        atCATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:272318/69‑1 (MQ=255)                                        atCATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:1301558/69‑1 (MQ=255)                                        atCATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:1254785/69‑1 (MQ=255)                                        atCATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:2441270/69‑1 (MQ=255)                                        atCATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:1153876/69‑1 (MQ=255)                                        atCATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:2472727/69‑1 (MQ=255)                                          cATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:2931890/67‑1 (MQ=255)                                                cTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:1491809/49‑1 (MQ=255)                                                 tGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGCGATGAAccc  >  1:2812376/1‑69 (MQ=255)                                                 tGAAGCCCCGTCAGGTAGTAAAACCTACGTTTCTGTTACCCCATCAGCAATGGCAGAAGCGATGAAccc  >  1:1064795/1‑69 (MQ=255)                                                             aGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           >  1:2550879/1‑48 (MQ=255)                                                             aGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           >  1:451971/1‑48 (MQ=255)                                                             aGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:2282096/48‑1 (MQ=255)                                                             aGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           >  1:1815814/1‑48 (MQ=255)                                                             aGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           >  1:2141935/1‑48 (MQ=255)                                                               gAAGTAAAACCTCTGTTTCTGTTACCCCATCAGCAATGGCAGAAGCGATGAAccc  <  1:29824/55‑1 (MQ=255)                                                               gAAGTAAAACCTCTGTTTCTGTTACCCCATCAGCAATGGCAGAAGCGATGAAccc  <  1:2880252/55‑1 (MQ=255)                                                                aaGTAAAACCTCTGTTTCTGTTACCCCATCAGCAATGGCAGAAGCGATGAAccc  <  1:3100925/54‑1 (MQ=255)                                                                |                                                      GTTCTGTAAGCGTTTTACATCCGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGCAATGGCAGAAGCGATGAACCC  >  minE/1246278‑1246394

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