Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	909,573	C→T	46.5%	L584L (CTG→TTG)	ycjT →	predicted hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	909,573	0	C	T	46.5%	19.0 / 98.7	71	L584L (CTG→TTG)	ycjT	predicted hydrolase
Reads supporting (aligned to +/- strand):  ref base C (27/11);  new base T (21/12);  total (48/23)
Fisher's exact test for biased strand distribution p-value = 6.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTCACGCGCAGAAGTGAACGAGATGCAGATCCTCAAACAAGCTGATGTGGTGATGCTCAATTACATGCTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAATTTTATGAACCGCGCACTATTCACGA  >  minE/909505‑909643                                                                     |                                                                       attCACGCGCAGAAGTGAACGAGATGCAGATCCTCAAACAAGCTGATGTGGTGATGCTCAATTACAtgttg                                                                      >  1:1318099/1‑71 (MQ=255)     aCGCGCAGAAGTGAACGAGATGCAGATCCTCAAACAAGCTGATGTGGTGATGCTCAATTACATGCTGCCgg                                                                  <  1:1589800/71‑1 (MQ=255)             aaGTGAACGAGATGCAGATCCTCAAACAAGCTGATGTGGTGATGCTCAATTACATGTTGCCGGAGCAGTtc                                                          <  1:493809/71‑1 (MQ=255)               gTGAACGAGATGCAGATCCTCAAACAAGCTGATGTGGTGATGCTCAATTACATGCTGCCGGAGCAGTtct                                                         >  1:737976/1‑70 (MQ=255)               gTGAACGAGATGCAGATCCTCAAACAAGCTGATGTGGTGATGCTCAATTACATGCTGCCGGAGCAGTtct                                                         >  1:3662900/1‑70 (MQ=255)                 gAACGAGATGCAGATCCTCAAACAAGCTGATGTGGTGATGCTCAATTACAtgttg                                                                      <  1:2309305/55‑1 (MQ=255)                  aaCGAGATGCAGATCCTCAAACAAGCTGATGTGGTGATGCTCAATTACATGTTGCCGGAGCAGTtctg                                                        >  1:1074262/1‑67 (MQ=255)                  aaCGAGATGCAGATCCTCAAACAAGCTGATGTGGTGATGCTCAATTACATGTTGCCGGAGCAGTtctg                                                        >  1:1932310/1‑67 (MQ=255)                  aaCGAGATGCAGATCCTCAAACAAGCTGATGTGGTGATGCTCAATTACATGTTGCCGGAGCAGTtctg                                                        >  1:2508839/1‑67 (MQ=255)                  aaCGAGATGCAGATCCTCAAACAAGCTGATGTGGTGATGCTCAATTACATGTTGCCGGAGCAGTtctg                                                        >  1:2466786/1‑67 (MQ=255)                    cGAGATGCAGATCCTCAAACAAGCTGATGTGGTGATGCTCAATTACATGCTGCCGGAGCAg                                                             >  1:1368994/1‑61 (MQ=255)                    cGAGATGCAGATCCTCAAACAAGCTGATGTGGTGATGCTCAATTACATGCTGCCGGAGCAg                                                             >  1:3096618/1‑61 (MQ=255)                    cGAGATGCAGATCCTCAAACAAGCTGATGTGGTGATGCTCAATTACATGCTGCCGGAGCAGTTCTCAg                                                      >  1:3729438/1‑68 (MQ=255)                    cGAGATGCAGATCCTCAAACAAGCTGATGTGGTGATGCTCAATTACATGCTGCCGGAGCAGTTCTCAg                                                      >  1:1687992/1‑68 (MQ=255)                    cGAGATGCAGATCCTCAAACAAGCTGATGTGGTGATGCTCAATTACATGCTGCCGGAGCAGTTCTCAg                                                      >  1:481781/1‑68 (MQ=255)                    cGAGATGCAGATCCTCAAACAAGCTGATGTGGTGATGCTCAATTACATGCTGCCGGAGCAGTTCTCAg                                                      >  1:2834494/1‑68 (MQ=255)                      agaTGCAGATCCTCAAACAAGCTGATGTGGTGATGCTCAATTACATGTTGCCGGAGCAGTTCTCAGCGGCa                                                 >  1:2153408/1‑71 (MQ=255)                      agaTGCAGATCCTCAAACAAGCTGATGTGGTGATGCTCAATTACATGTTGCCGGAGCAGTTCTCAGCGGCa                                                 >  1:2103938/1‑71 (MQ=255)                      agaTGCAGATCCTCAAACAAGCTGATGTGGTGATGCTCAATTACATGTTGCCGGAGCAGTTCTCAGCGGCa                                                 >  1:2305910/1‑71 (MQ=255)                      agaTGCAGATCCTCAAACAAGCTGATGTGGTGATGCTCAATTACATGTTGCCGGAGCAGTTCTCAGCGGCa                                                 >  1:4196262/1‑71 (MQ=255)                        aTGCAGATCCTCAAACAAGCTGATGTGGTGATGCTCAATTACATGCTGCCGGAGCAGTTCTCAGCGGCATc                                               <  1:1592902/71‑1 (MQ=255)                          gCAGATCCTCAAACAAGCTGATGTGGTGATGCTCAATTACATGCTGCCGGAGCAGTTCTCAGCGGCATCgt                                             >  1:2442007/1‑71 (MQ=255)                              aTCCTCAAACAAGCTGATGTGGTGATGCTCAATTACAtgttgc                                                                     <  1:1832431/43‑1 (MQ=255)                              aTCCTCAAACAAGCTGATGTGGTGATGCTCAATTACAtgttgc                                                                     <  1:3024570/43‑1 (MQ=255)                              aTCCTCAAACAAGCTGATGTGGTGATGCTCAATTACAtgttgc                                                                     <  1:3529537/43‑1 (MQ=255)                              aTCCTCAAACAAGCTGATGTGGTGATGCTCAATTACAtgttgc                                                                     <  1:1506827/43‑1 (MQ=255)                              aTCCTCAAACAAGCTGATGTGGTGATGCTCAATTACAtgttgc                                                                     <  1:3963849/43‑1 (MQ=255)                              aTCCTCAAACAAGCTGATGTGGTGATGCTCAATTACAtgttgc                                                                     <  1:229073/43‑1 (MQ=255)                                ccTCAAACAAGCTGATGTGGTGATGCTCAATTACATGCTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTg                                       >  1:1824808/1‑71 (MQ=255)                                 cTCAAACAAGCTGATGTGGTGATGCTCAATTACATGCTGCCGGAGCAGTTCTCAGCGGCAt                                                >  1:812144/1‑61 (MQ=255)                                 cTCAAACAAGCTGATGTGGTGATGCTCAATTACATGCTGCCGGAGCAGTTCTCAGCGGCAt                                                >  1:3147625/1‑61 (MQ=255)                                 cTCAAACAAGCTGATGTGGTGATGCTCAATTACATGCTGCCGGAGCAGTTCTCAGCGGCAt                                                >  1:2013774/1‑61 (MQ=255)                                   cAAACAAGCTGATGTGGTGATGCTCAATTACATGTTGCCGGAGCAGTTCTCAGCGGCATCg                                              >  1:2763795/1‑61 (MQ=255)                                   cAAACAAGCTGATGTGGTGATGCTCAATTACATGTTGCCGGAGCAGTTCTCAGCGGCATCg                                              >  1:2991898/1‑61 (MQ=255)                                    aaaCAAGCTGATGTGGTGATGCTCAATTACATGTTGCCGTAGCAGTTCTCAGCGGCATCGTGTCTTGCCa                                    <  1:2241264/70‑1 (MQ=255)                                    aaaCAAGCTGATGTGGTGATGCTCAATTACATGTTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCa                                    <  1:1955172/70‑1 (MQ=255)                                       cAAGCTGATGTGGTGATGCTCAATTACATGCTGCCGGAGCAg                                                             >  1:1061278/1‑42 (MQ=255)                                         aGCTGATGTGGTGATGCTCAATTACATGTTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGc                              <  1:202770/71‑1 (MQ=255)                                          gCTGATGTGGTGATGCTCAATTACATGCTGCCGGAGCAGTTCTCAGCGGCATc                                               >  1:3128856/1‑53 (MQ=255)                                          gCTGATGTGGTGATGCTCAATTACATGCTGCCGGAGCAGTTCTCAGCGGCATc                                               >  1:593119/1‑53 (MQ=255)                                             gATGTGGTGATGCTCAATTACATGTTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAAtt                          >  1:2505532/1‑71 (MQ=255)                                             gATGTGGTGATGCTCAATTACATGTTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAAtt                          >  1:428427/1‑71 (MQ=255)                                              aTGTGGTGATGCTCAATTACATGTTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAAttt                         <  1:3543464/71‑1 (MQ=255)                                              aTGTGGTGATGCTCAATTACATGTTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAAttt                         >  1:4045593/1‑71 (MQ=255)                                              aTGTGGTGATGCTCAATTACATGTTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAAttt                         <  1:1326553/71‑1 (MQ=255)                                              aTGTGGTGATGCTCAATTACATGCTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAAttt                         <  1:304872/71‑1 (MQ=255)                                              aTGTGGTGATGCTCAATTACATGCTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAAttt                         <  1:29972/71‑1 (MQ=255)                                              aTGTGGTGATGCTCAATTACATGCTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAAttt                         <  1:1504786/71‑1 (MQ=255)                                              aTGTGGTGACGCTCAATTACATGCTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAAttt                         <  1:3285551/71‑1 (MQ=255)                                                gtggtgATGCTCAATTACATGCTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAATTTTa                       >  1:314073/1‑71 (MQ=255)                                                gtggtgATGCTCAATTACATGCTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAATTTTa                       >  1:794523/1‑71 (MQ=255)                                                gtggtgATGCTCAATTACATGCTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAATTTTa                       >  1:1897410/1‑71 (MQ=255)                                                gtggtgATGCTCAATTACATGCTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAATTTTa                       >  1:4139049/1‑71 (MQ=255)                                                gtggtgATGCTCAATTACATGCTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAATTTTa                       >  1:2783535/1‑71 (MQ=255)                                                gtggtgATGCTCAATTACATGCTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAATTTTa                       >  1:3575977/1‑71 (MQ=255)                                                 tggtgATGCTCAATTACATGTTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAATTTTAt                      >  1:535593/1‑71 (MQ=255)                                                 tggtgATGCTCAATTACATGTTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAATTTTAt                      >  1:497762/1‑71 (MQ=255)                                                 tggtgATGCTCAATTACATGTTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAATTTTAt                      >  1:4171506/1‑71 (MQ=255)                                                 tggtgATGCTCAATTACATGCTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAATTTTa                       <  1:2230734/70‑1 (MQ=255)                                                 tggtgATGCTCAATTACATGCTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAATTTTAt                      <  1:120410/71‑1 (MQ=255)                                                 tggtgATGCTCAATTACATGCTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAATTTTAt                      <  1:3409213/71‑1 (MQ=255)                                                  ggtgATGCTCAATTACATGTTGCCGGAGCAGTTCTCAGCgg                                                   <  1:3439720/41‑1 (MQ=255)                                                   gtgATGCTCAATTACATGTTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAATTTTATg                     <  1:1937630/70‑1 (MQ=255)                                                        gCTCAATTACATGTTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAAttt                         <  1:2856918/61‑1 (MQ=255)                                                          tCAATTACATGCTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAATTTTATg                     <  1:1283070/63‑1 (MQ=255)                                                              ttACATGTTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAATTTTa                       >  1:1627155/1‑57 (MQ=255)                                                              ttACATGTTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAATTTTATGAACCGCGCACTa         <  1:3285824/71‑1 (MQ=255)                                                              ttACATGTTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAATTTTATGAACCGCGCACTa         <  1:3266271/71‑1 (MQ=255)                                                              ttACATGTTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAATTTTATGAACCGCGCACTa         <  1:3995777/71‑1 (MQ=255)                                                              ttACATGCTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAAt                                  <  1:171379/46‑1 (MQ=255)                                                               tACATGCTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAATTTTATGAACCGCGCa            >  1:3732993/1‑67 (MQ=255)                                                                aCATGTTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAATTTTATGAACCGCGCACTAtt       >  1:4055085/1‑71 (MQ=255)                                                                aCATGTTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAATTTTATGAACCGCGCACTAtt       >  1:1361804/1‑71 (MQ=255)                                                                aCATGTTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAATTTTATGAACCGCGCACTAt        >  1:1356739/1‑70 (MQ=255)                                                                aCATGTTGCCGGAGCAGTTCTCAGCGGAATCGTGTCTTGCCAATCTGCAAtttt                        >  1:797866/1‑54 (MQ=255)                                                                aCATGCTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATc                                 >  1:3487357/1‑45 (MQ=255)                                                                aCATGCTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCt                                >  1:3646176/1‑46 (MQ=255)                                                                aCATGCTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAATTTTATGAACCGCGCACTAt        >  1:154170/1‑70 (MQ=255)                                                                aCATGCTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAATTTTATGAACCGCGCACTAt        >  1:312668/1‑70 (MQ=255)                                                                   tgctgcCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAATTTTATGAACCGCGCACTATTCAc    <  1:153151/71‑1 (MQ=255)                                                                     ctgcCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAATTTTATGAACCGCGCACTATTCACGa  >  1:379539/1‑71 (MQ=255)                                                                     ctgcCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAATTTTATGAACCGCGCACTATTCACGa  >  1:1412375/1‑71 (MQ=255)                                                                     ctgcCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAATTTTATGAACCGCGCACTATTCACGa  >  1:3584814/1‑71 (MQ=255)                                                                     |                                                                       ATTCACGCGCAGAAGTGAACGAGATGCAGATCCTCAAACAAGCTGATGTGGTGATGCTCAATTACATGCTGCCGGAGCAGTTCTCAGCGGCATCGTGTCTTGCCAATCTGCAATTTTATGAACCGCGCACTATTCACGA  >  minE/909505‑909643

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