Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,388,186	A→G	41.2%	N76N (AAT→AAC)	yigL ←	predicted hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,388,186	0	A	G	41.2%	19.3 / 44.8	63	N76N (AAT→AAC)	yigL	predicted hydrolase
Reads supporting (aligned to +/- strand):  ref base A (11/26);  new base G (22/4);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 2.62e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

TTGTCGTTGACTACGCCAAACAGATCGCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGACTTAATCTCCAGATT  >  minE/2388121‑2388251                                                                  |                                                                  ttGTCGTTGACTACGCCAAACAGATCGCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTa                                                                <  1:666444/69‑1 (MQ=255) ttGTCGTTGACTACGCCAAACAGATCGCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTa                                                                <  1:3412243/69‑1 (MQ=255) ttGTCGTTGACTACGCCAAACAGATCGCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTa                                                                <  1:3066533/69‑1 (MQ=255) ttGTCGTTGACTACGCCAAACAGATCGCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTa                                                                <  1:2332429/69‑1 (MQ=255) ttGTCGTTGACTACGCCAAACAGATCGCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTa                                                                <  1:2247117/69‑1 (MQ=255) ttGTCGTTGACTACGCCAAACAGATCGCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTa                                                                <  1:1713551/69‑1 (MQ=255) ttGTCGTTGACTACGCCAAACAGATCGCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTa                                                                <  1:14611/69‑1 (MQ=255)                          cGCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCAc                                       <  1:971682/69‑1 (MQ=255)                          cGCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCAc                                       <  1:2372789/69‑1 (MQ=255)                          cGCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCAc                                       <  1:1735913/69‑1 (MQ=255)                          cGCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCAc                                       <  1:1571577/69‑1 (MQ=255)                          cGCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCAc                                       <  1:1455187/69‑1 (MQ=255)                           gCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCAcc                                      <  1:1436329/69‑1 (MQ=255)                           gCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCAcc                                      <  1:3008474/69‑1 (MQ=255)                           gCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCAc                                       <  1:1896685/68‑1 (MQ=255)                                      cgcgATCCAGGTTATGATCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTa                            >  1:2578235/1‑68 (MQ=255)                                      cgcgATCCAGGTTATGAGCAAATATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTa                            >  1:1190348/1‑68 (MQ=255)                                      cgcgATCCAGGTTATGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTa                            >  1:1134455/1‑68 (MQ=255)                                      cgcgATCCAGGTTATGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTa                            >  1:1880586/1‑68 (MQ=255)                                      cgcgATCAAGGTTATGAGCAAAAATCAGGTTGCCATCAAGATCGTGAACGCGCGCGCCATTGGAGGTa                            >  1:276227/1‑68 (MQ=255)                                       gcgATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAAt                          >  1:1017582/1‑69 (MQ=255)                                       gcgATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAAt                          >  1:41172/1‑69 (MQ=255)                                       gcgATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAAt                          >  1:2995231/1‑69 (MQ=255)                                       gcgATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAAt                          <  1:1945759/69‑1 (MQ=255)                                       gcgATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAAt                          <  1:1874949/69‑1 (MQ=255)                                            ccAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCAt                       >  1:3142476/1‑67 (MQ=255)                                               ggTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTaa                   <  1:865891/68‑1 (MQ=255)                                               ggTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTaa                   <  1:3064321/68‑1 (MQ=255)                                               ggTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTaa                   <  1:1722473/68‑1 (MQ=255)                                               ggTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTaa                   <  1:2985614/68‑1 (MQ=255)                                               ggTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTaa                   <  1:2906638/68‑1 (MQ=255)                                               ggTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTaa                   <  1:1355763/68‑1 (MQ=255)                                               ggTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTaa                   <  1:2640553/68‑1 (MQ=255)                                               ggTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTaa                   <  1:2484541/68‑1 (MQ=255)                                                 ttATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGa                 >  1:1484006/1‑68 (MQ=255)                                                 ttATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGa                 >  1:1732665/1‑68 (MQ=255)                                                 ttATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGa                 >  1:2604691/1‑68 (MQ=255)                                                  tATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCa                        >  1:234902/1‑60 (MQ=255)                                                  tATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGAc                >  1:1853521/1‑68 (MQ=255)                                                  tATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGAc                >  1:1152082/1‑68 (MQ=255)                                                   aTGAGCAAGAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                        >  1:1341656/1‑59 (MQ=255)                                                   aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAAt                          >  1:1419630/1‑57 (MQ=255)                                                   aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATc                         >  1:1795305/1‑58 (MQ=255)                                                   aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                        >  1:1948842/1‑59 (MQ=255)                                                   aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                        >  1:1901317/1‑59 (MQ=255)                                                   aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                        >  1:1115152/1‑59 (MQ=255)                                                   aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                        >  1:521874/1‑59 (MQ=255)                                                   aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                        >  1:429727/1‑59 (MQ=255)                                                   aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                        >  1:3143561/1‑59 (MQ=255)                                                   aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                        >  1:137106/1‑59 (MQ=255)                                                   aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                        >  1:2418287/1‑59 (MQ=255)                                                   aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                        >  1:2968155/1‑59 (MQ=255)                                                   aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                        >  1:2253128/1‑59 (MQ=255)                                                   aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                        >  1:1761502/1‑59 (MQ=255)                                                   aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                        >  1:2750752/1‑59 (MQ=255)                                                   aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                        >  1:1738687/1‑59 (MQ=255)                                                   aTGAGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCa                        >  1:255217/1‑59 (MQ=255)                                                   aTGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGACtt              <  1:1985291/69‑1 (MQ=255)                                                      aGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCATGTAAGACTTAAt           <  1:3063921/69‑1 (MQ=255)                                                      aGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCATGTAAGACTTAAt           <  1:3297449/69‑1 (MQ=255)                                                      aGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCATGTAAGACTTAAt           <  1:3309369/69‑1 (MQ=255)                                                      aGCAAAAATCAGGTTGCCATCCAGATCGTGAACGCGCGCGCCATTGGAGGTAATCATGTAAGACTTAAt           <  1:1119772/69‑1 (MQ=255)                                                               cAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGACTTAATCTCCAGAtt  >  1:2107349/1‑69 (MQ=255)                                                                  |                                                                  TTGTCGTTGACTACGCCAAACAGATCGCTGGCAATGTCGCGATCCAGGTTATGAGCAAAAATCAGATTACCATCCAGATCGTGAACGCGCGCACCATTGGAGGTAATCATGTAAGACTTAATCTCCAGATT  >  minE/2388121‑2388251

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