Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,628,493 T→G 21.7% A321A (GCT→GCG livJ → leucine/isoleucine/valine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,628,4930TG21.7% 100.5 / 29.4 69A321A (GCT→GCGlivJleucine/isoleucine/valine transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (7/47);  new base G (6/9);  total (13/56)
Fisher's exact test for biased strand distribution p-value = 2.78e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

CCTACGCCGCGCTGCAATCTTTGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGTAATGGGACCGCTGACCT  >  minE/2628438‑2628551
                                                       |                                                          
ccTACGCCGCGCTGCAATCTTTGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAAt                                               <  1:2583798/69‑1 (MQ=255)
        gcgcTGCAATCTTTGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGaaa                                       <  1:1150801/69‑1 (MQ=255)
        gcgcTGCAATCTTTGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGaaa                                       <  1:2951619/69‑1 (MQ=255)
        gcgcTGCAATCTTTGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGaaa                                       <  1:2817201/69‑1 (MQ=255)
        gcgcTGCAATCTTTGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGaaa                                       <  1:2461315/69‑1 (MQ=255)
                aTCTTTGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTc                               <  1:2925151/69‑1 (MQ=255)
                aTCTTTGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTc                               <  1:1104446/69‑1 (MQ=255)
                aTCTTTGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTc                               <  1:2351430/69‑1 (MQ=255)
                aTCTTTGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTc                               <  1:1598505/69‑1 (MQ=255)
                aTCTTTGCACGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTc                               <  1:1508132/69‑1 (MQ=255)
                  cTTTGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCg                             <  1:2112303/69‑1 (MQ=255)
                   tttGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAgat                                    >  1:1538122/1‑59 (MQ=255)
                   tttGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGc                                     <  1:2750631/60‑1 (MQ=255)
                   tttGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCg                                    >  1:943045/1‑61 (MQ=255)
                   tttGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCg                                    >  1:810366/1‑61 (MQ=255)
                   tttGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCg                                    >  1:625852/1‑61 (MQ=255)
                   tttGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCg                                    >  1:3072970/1‑61 (MQ=255)
                   tttGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCg                                    <  1:3012695/61‑1 (MQ=255)
                   tttGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCg                                    >  1:217487/1‑61 (MQ=255)
                   tttGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCg                                    <  1:2910943/61‑1 (MQ=255)
                   tttGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCg                                    <  1:2754822/61‑1 (MQ=255)
                   tttGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCg                                    <  1:2358321/61‑1 (MQ=255)
                   tttGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCg                                    >  1:1399961/1‑61 (MQ=255)
                   tttGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCg                                    <  1:269386/61‑1 (MQ=255)
                       cAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTg                           <  1:2800593/66‑1 (MQ=255)
                        aGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGAt                        <  1:177519/68‑1 (MQ=255)
                        aGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGAt                        <  1:593907/68‑1 (MQ=255)
                        aGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGAt                        <  1:3392532/68‑1 (MQ=255)
                        aGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGAt                        <  1:2491802/68‑1 (MQ=255)
                        aGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGAt                        <  1:1469901/68‑1 (MQ=255)
                         ggCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATAc                      <  1:683476/69‑1 (MQ=255)
                         ggCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATAc                      <  1:2026585/69‑1 (MQ=255)
                         ggCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATAc                      <  1:2158145/69‑1 (MQ=255)
                            gggCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGt                   <  1:91666/69‑1 (MQ=255)
                            gggCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGt                   <  1:2695261/69‑1 (MQ=255)
                            gggCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGt                   <  1:3356279/69‑1 (MQ=255)
                            gggCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGt                   <  1:2187145/69‑1 (MQ=255)
                             ggCCTGAATCAGTCTGACGATCCGGCGGAGATCGCTAAATACCTGAAAGCGAACTCGGTTGATACGGtg                  <  1:1702533/69‑2 (MQ=255)
                             ggCCTGAATCAGTCTGACGATCCGGCGGAGATCGCTAAATACCTGAAAGCGAACTCGGTTGATACGGtg                  <  1:1945972/69‑2 (MQ=255)
                             ggCCTGAATCAGTCTGACGATCCGGCGGAGATCGCTAAATACCTGAAAGCGAACTCGGTTGATACGGtg                  <  1:3286455/69‑2 (MQ=255)
                             ggCCTGAATCAGTCTGACGATCCGGCGGAGATCGCTAAATACCTGAAAGCGAACTCGGTTGATACGGtg                  <  1:3250649/69‑2 (MQ=255)
                             ggCCTGAATCAGTCTGACGATCCGGCGGAGATCGCTAAATACCTGAAAGCGAACTCGGTTGATACGGtg                  <  1:3115631/69‑2 (MQ=255)
                             ggCCTGAATCAGTCTGACGATCCGGCGGAGATCGCTAAATACCTGAAAGCGAACTCGGTTGATACGGtg                  <  1:3033380/69‑2 (MQ=255)
                             ggCCTGAATCAGTCTGACGATCCGGCGGAGATCGCTAAATACCTGAAAGCGAACTCGGTTGATACGGtg                  <  1:2991134/69‑2 (MQ=255)
                             ggCCTGAATCAGTCTGACGATCCGGCGGAGATCGCTAAATACCTGAAAGCGAACTCGGTTGATACGGtg                  <  1:266942/69‑2 (MQ=255)
                             ggCCTGAATCAGTCTGACGATCCGGCGGAGATCGCTAAATACCTGAAAGCGAACTCGGTTGATACGGtg                  <  1:2936418/69‑2 (MQ=255)
                                   aaTCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGTAATgg              <  1:2110076/67‑1 (MQ=255)
                                       aGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCg                                    <  1:2959647/41‑1 (MQ=255)
                                       aGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCg                                    <  1:891105/41‑1 (MQ=255)
                                       aGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCg                                    <  1:2474467/41‑1 (MQ=255)
                                       aGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCg                                    <  1:2773431/41‑1 (MQ=255)
                                       aGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCg                                    <  1:2203202/41‑1 (MQ=255)
                                       aGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCg                                    <  1:2953149/41‑1 (MQ=255)
                                          cTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGTAATGGGACCGCTGa     <  1:1583574/69‑1 (MQ=255)
                                          cTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGTAATGGGACCGCTGa     <  1:3397302/69‑1 (MQ=255)
                                          cTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGTAATGGGACCGCTGa     <  1:1440106/69‑1 (MQ=255)
                                          cTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGTAATGGGACCGCTGa     <  1:2955784/69‑1 (MQ=255)
                                          cTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGTAATGGGACCGCTGa     <  1:718167/69‑1 (MQ=255)
                                           tGACGATCCGGCGGAGATCGCTAAATACCTGAAAGCGAACTCGGTTGATACGGTGATGGGGCCGCTGAc    >  1:1344134/1‑69 (MQ=255)
                                           tGACGATCCGGCGGAGATCGCTAAATACCTGAAAGCGAACTCGGTTGATACGGTGATGGGGCCGCTGAc    >  1:1332954/1‑69 (MQ=255)
                                           tGACGATCCGGCGGAGATCGCTAAATACCTGAAAGCGAACTCGGTTGATACGGTGATGGGGCCGCTGAc    >  1:1188045/1‑69 (MQ=255)
                                           tGACGATCCGGCGGAGATCGCTAAATACCTGAAAGCGAACTCGGTTGATACGGTGATGGGGCCGCTGAc    >  1:345662/1‑69 (MQ=255)
                                           tGACGATCCGGCGGAGATCGCTAAATACCTGAAAGCGAACTCGGTTGATACGGTGATGGGGCCGCTGAc    >  1:1400855/1‑69 (MQ=255)
                                           tGACGATCCGGCGGAGATCGCTAAATACCTGAAAGCGAACTCGGTTGATACGGTGATGGGGCCGCTGAc    >  1:868661/1‑69 (MQ=255)
                                            gaggATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTTGATACCGTAATGGGACCGCTGa     <  1:3279910/64‑1 (MQ=255)
                                               gATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGTAATGGGACCGCTGACCt  <  1:1872837/67‑1 (MQ=255)
                                               gATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGTAATGGGACCGCTGACCt  <  1:2686121/67‑1 (MQ=255)
                                                  ccGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGTAATGGGACCGCTGa     <  1:3012547/61‑1 (MQ=255)
                                                  ccGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGTAATGGGACCGCTGa     <  1:597803/61‑1 (MQ=255)
                                                       |                                                          
CCTACGCCGCGCTGCAATCTTTGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGTAATGGGACCGCTGACCT  >  minE/2628438‑2628551

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: