Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	327,825	A→G	53.9%	R206R (CGT→CGC)	aes ←	acetyl esterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	327,825	0	A	G	53.9%	26.3 / 112.9	89	R206R (CGT→CGC)	aes	acetyl esterase
Reads supporting (aligned to +/- strand):  ref base A (27/14);  new base G (27/21);  total (54/35)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATATGCCTCTTCGTACATCTGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATC  >  minE/327757‑327893                                                                     |                                                                     aaaTATGCCTCTTCGTACATCTGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAgacga                                                                      <  1:777352/69‑1 (MQ=255)  aaTATGCCTCTTCGTACATCTGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAgacgac                                                                     <  1:1918086/69‑1 (MQ=255)    tatGCCTCTTCGTACATCTGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGGCGAGt                                                                 <  1:303352/71‑1 (MQ=255)      tGCCTCTTCGTACATCTGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGGCGAGTCa                                                               >  1:2701914/1‑71 (MQ=255)      tGCCTCTTCGTACATCTGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGGCGAGTCa                                                               >  1:2841328/1‑71 (MQ=255)       gCCTCTTCGTACATCTGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCAc                                                              >  1:581398/1‑71 (MQ=255)         ctctTCGTACATCTGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGGCGAGTCACgg                                                            <  1:1718305/71‑1 (MQ=255)          tctTCGTACATCTGCAAATCCTGTTGCTTTAAGCCATCCCAGACACCTCCCAACAGACGGCGAGTc                                                                <  1:2548814/66‑1 (MQ=255)            ttCGTACATCTGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGGCGAGTCACGGa                                                           <  1:1322894/69‑1 (MQ=255)            ttCGTACATCTGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGGCGAGTCACGGa                                                           <  1:651716/69‑1 (MQ=255)            ttCGTACATCTGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGGCGAGTCACGGa                                                           <  1:730575/69‑1 (MQ=255)            ttCGTACATCTGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGGCGAGTCACGGa                                                           <  1:2765418/69‑1 (MQ=255)            ttCGTACATCTGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGGCGAGTCACGGAAt                                                         >  1:1081107/1‑71 (MQ=255)            ttCGTACATCTGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGGCGAGTCACGGAAt                                                         >  1:3247054/1‑71 (MQ=255)             tCGTACATCTGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATc                                                        >  1:2990994/1‑71 (MQ=255)             tCGTACATCTGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATc                                                        >  1:2163022/1‑71 (MQ=255)              cGTACATCTGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGGCGAGTCACGGaa                                                          <  1:794939/68‑1 (MQ=255)              cGTACATCTGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGGCGAGTCACGGaa                                                          <  1:3023743/68‑1 (MQ=255)              cGTACATCTGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATcc                                                       >  1:4472878/1‑71 (MQ=255)               gTACATCTGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGGCGAGTCACGGAATccc                                                      <  1:4295497/71‑1 (MQ=255)                tACATCTGCAAATCCTGTTGCGTTAAGCCATCCCAGACACGGCCCAACAGACGGCGAGTCACGGAATCCCg                                                     <  1:4068120/71‑1 (MQ=255)                tACATCTGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCg                                                     >  1:2612751/1‑71 (MQ=255)                tACATCTGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCg                                                     >  1:1427794/1‑71 (MQ=255)                    tCTGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGGCGAGTCAc                                                              >  1:809330/1‑58 (MQ=255)                      tGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTAATcc                                               >  1:313667/1‑71 (MQ=255)                      tGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTAATcc                                               >  1:1843638/1‑71 (MQ=255)                        cAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTAAt                                                 <  1:831204/67‑1 (MQ=255)                        cAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGt                                                    >  1:243204/1‑64 (MQ=255)                        cAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAAt                                                 <  1:480374/67‑1 (MQ=255)                        cAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATc                                                >  1:4642461/1‑68 (MQ=255)                        cAAATCCTGTGGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATc                                                <  1:1173840/68‑1 (MQ=255)                            tCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTAATCCGTAAAg                                         >  1:63335/1‑71 (MQ=255)                            tCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTc                                                                <  1:3770173/48‑1 (MQ=255)                              cTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGcc                                       <  1:731974/71‑1 (MQ=255)                              cTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGcc                                       <  1:134770/71‑1 (MQ=255)                                 ttGCGTTAAGCCATCCCAGACACCGCCCAACAGACGGCGAGTCACGGaa                                                          >  1:3882625/1‑49 (MQ=255)                                 ttGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACgg                                                            <  1:518428/47‑1 (MQ=255)                                  tGCGTTAAGCCATCCCAGACACCGCCCAACAGacgacg                                                                    >  1:3709442/1‑38 (MQ=255)                                  tGCGTTAAGCCATCCCAGACACCGCCCAACAGacgacg                                                                    >  1:3356002/1‑38 (MQ=255)                                  tGCGTTAAGCCATCCCAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTAATCCGTaaa                                          >  1:4150128/1‑64 (MQ=255)                                    cGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCg                                              <  1:1125774/58‑1 (MQ=255)                                    cGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCa                                     <  1:2794016/67‑1 (MQ=255)                                    cGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCa                                     <  1:2420567/67‑1 (MQ=255)                                    cGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATAc                                  <  1:392065/70‑1 (MQ=255)                                    cGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATAc                                  <  1:2217957/70‑1 (MQ=255)                                     gTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAAt                                                 <  1:2930757/54‑1 (MQ=255)                                       tAAGCCATCCCAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTAAt                                                 <  1:4184019/52‑1 (MQ=255)                                        aaGCCATCCCAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACcaca                              <  1:3782275/70‑1 (MQ=255)                                        aaGCCATCCCAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAg                             >  1:2464369/1‑71 (MQ=255)                                        aaGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAg                             >  1:2775952/1‑71 (MQ=255)                                            cATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATc                                                        >  1:133668/1‑40 (MQ=255)                                             aTCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCaaaa                        >  1:4545508/1‑71 (MQ=255)                                             aTCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCaaaa                        >  1:3348286/1‑71 (MQ=255)                                                ccAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACca                                >  1:481241/1‑60 (MQ=255)                                                ccAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCAGAGCAAAACGc                     >  1:1473581/1‑71 (MQ=255)                                                ccAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCaaaa                        >  1:395145/1‑68 (MQ=255)                                                ccAGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGc                     >  1:512509/1‑71 (MQ=255)                                                ccAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAAccc                     >  1:2631205/1‑69 (MQ=255)                                                ccAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGc                     >  1:314254/1‑71 (MQ=255)                                                ccAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGc                     >  1:4552692/1‑71 (MQ=255)                                                ccAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGc                     >  1:2348479/1‑71 (MQ=255)                                                ccAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGc                     >  1:593969/1‑71 (MQ=255)                                                  aGACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACg                      <  1:3472782/68‑1 (MQ=255)                                                  aGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGcc                    >  1:1090385/1‑70 (MQ=255)                                                   gACACCGCCCAACAGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCg                  >  1:1861307/1‑71 (MQ=255)                                                   gACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATAc                                  <  1:3995201/55‑1 (MQ=255)                                                   gACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCg                  >  1:964640/1‑71 (MQ=255)                                                    acacCGCCCAACAGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGCCCa                                     >  1:2939774/1‑51 (MQ=255)                                                    acacCGCCCAACAGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCaa                          <  1:4005083/62‑1 (MQ=255)                                                    acacCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCg                  >  1:4670865/1‑70 (MQ=255)                                                     cacCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGc                 >  1:1643289/1‑70 (MQ=255)                                                      acCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGc                     >  1:2558084/1‑65 (MQ=255)                                                      acCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCaa               <  1:173741/71‑1 (MQ=255)                                                        cGCCCAACAGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAg                             >  1:2898672/1‑55 (MQ=255)                                                        cGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACt             >  1:966936/1‑71 (MQ=255)                                                         gCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACtt            >  1:2088433/1‑71 (MQ=255)                                                          cccAACAGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACtt            >  1:1723231/1‑70 (MQ=255)                                                          cccAACAGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACtt            >  1:3770079/1‑70 (MQ=255)                                                           ccAACAGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCaaaa                        >  1:1718931/1‑57 (MQ=255)                                                           ccAACAGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCaaaa                        >  1:361670/1‑57 (MQ=255)                                                            cAACAGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGCCCAt                                    <  1:2647472/44‑1 (MQ=255)                                                             aaCAGACGACGAGTCACGGAATCCCGTAATCCGTaaa                                          >  1:407760/1‑37 (MQ=255)                                                             aaCAGACGACGAGTCACGGAATACCGTAATCCGTaaa                                          >  1:2838334/1‑37 (MQ=255)                                                              aCAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTAc         >  1:2750943/1‑69 (MQ=255)                                                               cAGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGc                     <  1:3423612/56‑1 (MQ=255)                                                               cAGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGc      <  1:3489461/71‑1 (MQ=255)                                                               cAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACGACAGCAAAACGCCCGCAACTTTACCGc      <  1:2247858/71‑1 (MQ=255)                                                               cAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGc      <  1:3940423/71‑1 (MQ=255)                                                                aGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGCCCa                                     >  1:3491368/1‑39 (MQ=255)                                                                aGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTa          >  1:1669519/1‑66 (MQ=255)                                                                aGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGc      >  1:2480285/1‑70 (MQ=255)                                                                aGACGGCGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGc      >  1:4177176/1‑70 (MQ=255)                                                                aGACGGCGAGTCACGGAATCCCGTAATACTTAAAGCCCATACAACAGCAAAACGCCCGCAACtt            >  1:1842327/1‑64 (MQ=255)                                                                 gacggcGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCaa    >  1:3801728/1‑71 (MQ=255)                                                                   cgacgaGTCACGGAATCCCGTAATCCGTAAAGCCAATACCACAGCAAAACGCCCGCAACTTTACCGCAATc  >  1:191014/1‑71 (MQ=255)                                                                    gacgaGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGccc                   >  1:233459/1‑53 (MQ=255)                                                                     |                                                                     AAATATGCCTCTTCGTACATCTGCAAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATC  >  minE/327757‑327893

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑