Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,628,496 A→G 21.6% E322E (GAA→GAG livJ → leucine/isoleucine/valine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,628,4960AG21.6% 92.4 / 18.6 69E322E (GAA→GAGlivJleucine/isoleucine/valine transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (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

CCTACGCCGCGCTGCAATCTTTGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGTAATGGGACCGCTGACCTGGGATGAGAAAG  >  minE/2628438‑2628563
                                                          |                                                                   
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:2910943/61‑1 (MQ=255)
                   tttGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCg                                                <  1:269386/61‑1 (MQ=255)
                   tttGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCg                                                >  1:1399961/1‑61 (MQ=255)
                   tttGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCg                                                <  1:2754822/61‑1 (MQ=255)
                   tttGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCg                                                >  1:3072970/1‑61 (MQ=255)
                   tttGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCg                                                <  1:2358321/61‑1 (MQ=255)
                   tttGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCg                                                >  1:625852/1‑61 (MQ=255)
                   tttGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCg                                                >  1:217487/1‑61 (MQ=255)
                   tttGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCg                                                >  1:810366/1‑61 (MQ=255)
                   tttGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCg                                                <  1:3012695/61‑1 (MQ=255)
                   tttGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCg                                                >  1:943045/1‑61 (MQ=255)
                       cAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTg                                       <  1:2800593/66‑1 (MQ=255)
                        aGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGAt                                    <  1:1469901/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:177519/68‑1 (MQ=255)
                         ggCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATAc                                  <  1:2158145/69‑1 (MQ=255)
                         ggCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATAc                                  <  1:683476/69‑1 (MQ=255)
                         ggCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATAc                                  <  1:2026585/69‑1 (MQ=255)
                            gggCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGt                               <  1:3356279/69‑1 (MQ=255)
                            gggCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGt                               <  1:2695261/69‑1 (MQ=255)
                            gggCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGt                               <  1:2187145/69‑1 (MQ=255)
                            gggCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGt                               <  1:91666/69‑1 (MQ=255)
                             ggCCTGAATCAGTCTGACGATCCGGCGGAGATCGCTAAATACCTGAAAGCGAACTCGGTTGATACGGtg                              <  1:3115631/69‑2 (MQ=255)
                             ggCCTGAATCAGTCTGACGATCCGGCGGAGATCGCTAAATACCTGAAAGCGAACTCGGTTGATACGGtg                              <  1:3250649/69‑2 (MQ=255)
                             ggCCTGAATCAGTCTGACGATCCGGCGGAGATCGCTAAATACCTGAAAGCGAACTCGGTTGATACGGtg                              <  1:3286455/69‑2 (MQ=255)
                             ggCCTGAATCAGTCTGACGATCCGGCGGAGATCGCTAAATACCTGAAAGCGAACTCGGTTGATACGGtg                              <  1:2936418/69‑2 (MQ=255)
                             ggCCTGAATCAGTCTGACGATCCGGCGGAGATCGCTAAATACCTGAAAGCGAACTCGGTTGATACGGtg                              <  1:3033380/69‑2 (MQ=255)
                             ggCCTGAATCAGTCTGACGATCCGGCGGAGATCGCTAAATACCTGAAAGCGAACTCGGTTGATACGGtg                              <  1:2991134/69‑2 (MQ=255)
                             ggCCTGAATCAGTCTGACGATCCGGCGGAGATCGCTAAATACCTGAAAGCGAACTCGGTTGATACGGtg                              <  1:1702533/69‑2 (MQ=255)
                             ggCCTGAATCAGTCTGACGATCCGGCGGAGATCGCTAAATACCTGAAAGCGAACTCGGTTGATACGGtg                              <  1:1945972/69‑2 (MQ=255)
                             ggCCTGAATCAGTCTGACGATCCGGCGGAGATCGCTAAATACCTGAAAGCGAACTCGGTTGATACGGtg                              <  1:266942/69‑2 (MQ=255)
                                   aaTCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGTAATgg                          <  1:2110076/67‑1 (MQ=255)
                                       aGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCg                                                <  1:2953149/41‑1 (MQ=255)
                                       aGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCg                                                <  1:891105/41‑1 (MQ=255)
                                       aGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCg                                                <  1:2203202/41‑1 (MQ=255)
                                       aGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCg                                                <  1:2474467/41‑1 (MQ=255)
                                       aGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCg                                                <  1:2773431/41‑1 (MQ=255)
                                       aGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCg                                                <  1:2959647/41‑1 (MQ=255)
                                          cTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGTAATGGGACCGCTGa                 <  1:2955784/69‑1 (MQ=255)
                                          cTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGTAATGGGACCGCTGa                 <  1:718167/69‑1 (MQ=255)
                                          cTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGTAATGGGACCGCTGa                 <  1:1583574/69‑1 (MQ=255)
                                          cTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGTAATGGGACCGCTGa                 <  1:3397302/69‑1 (MQ=255)
                                          cTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGTAATGGGACCGCTGa                 <  1:1440106/69‑1 (MQ=255)
                                           tGACGATCCGGCGGAGATCGCTAAATACCTGAAAGCGAACTCGGTTGATACGGTGATGGGGCCGCTGAc                >  1:1188045/1‑69 (MQ=255)
                                           tGACGATCCGGCGGAGATCGCTAAATACCTGAAAGCGAACTCGGTTGATACGGTGATGGGGCCGCTGAc                >  1:868661/1‑69 (MQ=255)
                                           tGACGATCCGGCGGAGATCGCTAAATACCTGAAAGCGAACTCGGTTGATACGGTGATGGGGCCGCTGAc                >  1:345662/1‑69 (MQ=255)
                                           tGACGATCCGGCGGAGATCGCTAAATACCTGAAAGCGAACTCGGTTGATACGGTGATGGGGCCGCTGAc                >  1:1332954/1‑69 (MQ=255)
                                           tGACGATCCGGCGGAGATCGCTAAATACCTGAAAGCGAACTCGGTTGATACGGTGATGGGGCCGCTGAc                >  1:1400855/1‑69 (MQ=255)
                                           tGACGATCCGGCGGAGATCGCTAAATACCTGAAAGCGAACTCGGTTGATACGGTGATGGGGCCGCTGAc                >  1:1344134/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:597803/61‑1 (MQ=255)
                                                  ccGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGTAATGGGACCGCTGa                 <  1:3012547/61‑1 (MQ=255)
                                                         aaaTCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGt                               <  1:1407718/40‑1 (MQ=255)
                                                         aaaTCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGTAATGGGACCGCTGACCTGGGATGAGAAAg  >  1:3408827/1‑69 (MQ=255)
                                                          |                                                                   
CCTACGCCGCGCTGCAATCTTTGCAGGCGGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGTAATGGGACCGCTGACCTGGGATGAGAAAG  >  minE/2628438‑2628563

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: