Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	504,079	G→A	33.2%	H270Y (CAC→TAC)	ybhA ←	predicted hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	504,079	0	G	A	33.2%	71.2 / 44.9	87	H270Y (CAC→TAC)	ybhA	predicted hydrolase
Reads supporting (aligned to +/- strand):  ref base G (45/13);  new base A (10/19);  total (55/32)
Fisher's exact test for biased strand distribution p-value = 1.36e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGC  >  minE/504020‑504138                                                            |                                                            ggCCTGTGGCTGTATGCGCAAATTAAAAGTGTGTCGATACCCGCCTGATTAAATCAGGTAGc                                                           >  1:1832403/1‑62 (MQ=255) ggCCTGTGGCTGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATaaa                                                     >  1:3007367/1‑68 (MQ=255) ggCCTGTGGCTGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATaaa                                                     >  1:1578078/1‑68 (MQ=255) ggCCTGTGGCTGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATaaa                                                     >  1:2625911/1‑68 (MQ=255) ggCCTGTGGCTGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATaaa                                                     >  1:3407090/1‑68 (MQ=255) ggCCTGTGGCTGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATaaa                                                     >  1:769909/1‑68 (MQ=255) ggCCTGTGGCTGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATaaa                                                     >  1:2900669/1‑68 (MQ=255) ggCCTGTGGCTGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATaaa                                                     >  1:2742835/1‑68 (MQ=255) ggCCTGTGGCTGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATaaa                                                     >  1:1097661/1‑68 (MQ=255) ggCCTGTGGCTGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATaa                                                      >  1:2426430/1‑67 (MQ=255)    cTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATg                                                   <  1:1750990/67‑1 (MQ=255)    cTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATg                                                   <  1:730289/67‑1 (MQ=255)           tGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATAAATGAAGTGGGCa                                          <  1:2193938/69‑1 (MQ=255)           tGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATAAATGAACTggga                                           <  1:79775/68‑2 (MQ=255)           tGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATAAATGAACTGGGCa                                          <  1:2718810/69‑1 (MQ=255)           tGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATAAATGAACTGGGCa                                          <  1:680359/69‑1 (MQ=255)           tGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATAAATGAACTGGGCa                                          <  1:3150022/69‑1 (MQ=255)           tGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATAAATGAACTGGGCa                                          <  1:3226460/69‑1 (MQ=255)           tGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATAAATGAACTGGGCa                                          <  1:3260968/69‑1 (MQ=255)           tGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATAAATGAACTGGGCa                                          <  1:1965996/69‑1 (MQ=255)           tGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATAAATGAACTGGGCa                                          <  1:677404/69‑1 (MQ=255)           tGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATAAATGAACTGGGCa                                          <  1:1521803/69‑1 (MQ=255)           tGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATAAATGAACTGGGCa                                          <  1:783028/69‑1 (MQ=255)           tGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATAAATGAACTGGGCa                                          <  1:1425753/69‑1 (MQ=255)           tGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATAAATGAACTGGGCa                                          <  1:2512398/69‑1 (MQ=255)           tGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATAAATGAACTGGGCa                                          <  1:869912/69‑1 (MQ=255)           tGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATAAATGAACTGGGCa                                          <  1:924128/69‑1 (MQ=255)           tGTATGCGCAAATTAAAAGTGTGACGATAACCGCCTGATTAAATCAGGTAGCTATAAATGAACTGGGCa                                          <  1:2074594/69‑1 (MQ=255)           tGTATGCGCAAATTAAAAGTATGTCGATAACCGCCTGATTAAATCAGGTAGCTATAAATGAACTGGGCa                                          <  1:308231/69‑1 (MQ=255)           tGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGc                                           <  1:1656917/68‑1 (MQ=255)           tGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGc                                           <  1:1493329/68‑1 (MQ=255)           tGTACGCGCAAATTAAAAGGGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGc                                           <  1:1139105/68‑1 (MQ=255)               cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc                                      >  1:56132/1‑69 (MQ=255)               cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc                                      >  1:2245751/1‑69 (MQ=255)               cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc                                      >  1:2816304/1‑69 (MQ=255)               cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc                                      >  1:1088767/1‑69 (MQ=255)               cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc                                      >  1:777737/1‑69 (MQ=255)               cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc                                      >  1:791869/1‑69 (MQ=255)               cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc                                      >  1:845355/1‑69 (MQ=255)               cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc                                      >  1:938763/1‑69 (MQ=255)               cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc                                      >  1:980224/1‑69 (MQ=255)               cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc                                      >  1:1657840/1‑69 (MQ=255)               cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc                                      >  1:1076296/1‑69 (MQ=255)               cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc                                      >  1:1144374/1‑69 (MQ=255)               cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc                                      >  1:1365890/1‑69 (MQ=255)               cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc                                      >  1:1406161/1‑69 (MQ=255)               cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc                                      >  1:1442991/1‑69 (MQ=255)               cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc                                      >  1:1613608/1‑69 (MQ=255)               cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc                                      >  1:2239018/1‑69 (MQ=255)               cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc                                      >  1:1693942/1‑69 (MQ=255)               cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc                                      >  1:1942610/1‑69 (MQ=255)               cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc                                      >  1:2045327/1‑69 (MQ=255)                   cAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTc                                  <  1:1470469/69‑1 (MQ=255)                        tAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGc                                           >  1:1695398/1‑55 (MQ=255)                           aaGTGTGTCGATAACCGGCTGATTAAATCAGTTAGCTTTAAATGAACTGGGCa                                          <  1:3409871/53‑1 (MQ=255)                              tgtgtCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTAt                       <  1:1989824/69‑1 (MQ=255)                              tgtgtCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTAt                       <  1:1503120/69‑1 (MQ=255)                              tgtgtCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTAt                       <  1:1568845/69‑1 (MQ=255)                              tgtgtCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTAt                       <  1:2195060/69‑1 (MQ=255)                              tgtgtCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTAt                       <  1:1280579/69‑1 (MQ=255)                              tgtgtCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTAt                       <  1:2992091/69‑1 (MQ=255)                                   cGATAACCGCCTGATTAAATCAGGTAGCTATAAATGAACTGGGCa                                          <  1:1917343/45‑1 (MQ=255)                                    gATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGgtgt                          >  1:1441148/1‑60 (MQ=255)                                    gATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGtt                         >  1:768112/1‑61 (MQ=255)                                    gATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGtt                         >  1:706560/1‑61 (MQ=255)                                    gATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGtt                         >  1:678437/1‑61 (MQ=255)                                    gATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGtt                         >  1:2872723/1‑61 (MQ=255)                                     aTAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCa                  <  1:500215/67‑1 (MQ=255)                                                  tAAATCAGGTGGCTATATATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGgcg   >  1:1789061/1‑69 (MQ=255)                                                  tAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGgcg   >  1:2455191/1‑69 (MQ=255)                                                  tAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGgcg   >  1:216084/1‑69 (MQ=255)                                                  tAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGgcg   >  1:2376014/1‑69 (MQ=255)                                                  tAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGgcg   >  1:2594520/1‑69 (MQ=255)                                                  tAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGgcg   >  1:3419514/1‑69 (MQ=255)                                                  tAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGgcg   >  1:1272619/1‑69 (MQ=255)                                                  tAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGgcg   >  1:1673753/1‑69 (MQ=255)                                                  tAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGgcg   >  1:2683718/1‑69 (MQ=255)                                                  tAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGgc    >  1:1726887/1‑68 (MQ=255)                                                   aaaTCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGgcgc  >  1:837880/1‑69 (MQ=255)                                                   aaaTCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGgcgc  >  1:2803289/1‑69 (MQ=255)                                                   aaaTCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGgcgc  >  1:564390/1‑69 (MQ=255)                                                   aaaTCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGgcgc  >  1:3290877/1‑69 (MQ=255)                                                   aaaTCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGgcgc  >  1:2943816/1‑69 (MQ=255)                                                   aaaTCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGgcgc  >  1:1912777/1‑69 (MQ=255)                                                   aaaTCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGgcg   >  1:1399608/1‑68 (MQ=255)                                                   aaaTCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGgcg   >  1:1553489/1‑68 (MQ=255)                                                   aaaTCAAGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGgcgc  >  1:230/1‑69 (MQ=255)                                                            |                                                            GGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGC  >  minE/504020‑504138

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