Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,628,532 C→G 24.2% T334T (ACC→ACG livJ → leucine/isoleucine/valine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,628,5320CG24.2% 98.3 / 29.5 62T334T (ACC→ACGlivJleucine/isoleucine/valine transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (6/41);  new base G (6/9);  total (12/50)
Fisher's exact test for biased strand distribution p-value = 5.42e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-01

GGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGTAATGGGACCGCTGACCTGGGATGAGAAAGGCGATCTGAAAGGCTTTGAGTTCGGCGTATTTGACT  >  minE/2628466‑2628599
                                                                  |                                                                   
gggCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGt                                                                   <  1:91666/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)
 ggCCTGAATCAGTCTGACGATCCGGCGGAGATCGCTAAATACCTGAAAGCGAACTCGGTTGATACGGtg                                                                  <  1:266942/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:2936418/69‑2 (MQ=255)
 ggCCTGAATCAGTCTGACGATCCGGCGGAGATCGCTAAATACCTGAAAGCGAACTCGGTTGATACGGtg                                                                  <  1:1945972/69‑2 (MQ=255)
 ggCCTGAATCAGTCTGACGATCCGGCGGAGATCGCTAAATACCTGAAAGCGAACTCGGTTGATACGGtg                                                                  <  1:1702533/69‑2 (MQ=255)
       aaTCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGTAATgg                                                              <  1:2110076/67‑1 (MQ=255)
              cTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGTAATGGGACCGCTGa                                                     <  1:718167/69‑1 (MQ=255)
              cTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGTAATGGGACCGCTGa                                                     <  1:3397302/69‑1 (MQ=255)
              cTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGTAATGGGACCGCTGa                                                     <  1:1440106/69‑1 (MQ=255)
              cTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGTAATGGGACCGCTGa                                                     <  1:1583574/69‑1 (MQ=255)
              cTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGTAATGGGACCGCTGa                                                     <  1:2955784/69‑1 (MQ=255)
               tGACGATCCGGCGGAGATCGCTAAATACCTGAAAGCGAACTCGGTTGATACGGTGATGGGGCCGCTGAc                                                    >  1:1400855/1‑69 (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: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:3012547/61‑1 (MQ=255)
                      ccGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGTAATGGGACCGCTGa                                                     <  1:597803/61‑1 (MQ=255)
                             aaaTCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGt                                                                   <  1:1407718/40‑1 (MQ=255)
                             aaaTCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGTAATGGGACCGCTGACCTGGGATGAGAAAg                                      >  1:3408827/1‑69 (MQ=255)
                                tCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGTAATGGGACCGCTGACCTGGGATGAGAAAGGCg                                   <  1:2678938/69‑1 (MQ=255)
                                tCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGTAATGGGACCGCTGACCTGGGATGAGAAAGGCg                                   <  1:2424219/69‑1 (MQ=255)
                                tCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGTAATGGGACCGCTGACCTGGGATGAGAAAGGCg                                   <  1:2295670/69‑1 (MQ=255)
                                tCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGTAATGGGACCGCTGACCTGGGATGAGAAAGGCg                                   <  1:2535934/69‑1 (MQ=255)
                                               aaGCGAACTCCGTGGATACCGTAATGGGACCGCTGa                                                     <  1:2264847/36‑1 (MQ=255)
                                                aGCGAACTCCGTGGATACCGTAATGGGACCGCTGACCTGGGATGAGAAAGGCGATCTGAAAGGCTTTGa                   >  1:669536/1‑69 (MQ=255)
                                                       tCCGTGGATACCGTAATGGGACCGCTGACCTGGGATGAGAAAGGct                                   <  1:3057489/46‑2 (MQ=255)
                                                           tGGATACCGTAATGGGACCGCTGACCTGGGATGAGAAAGGCGATCTGAAAGGCTTTGAGTTCGGCGTa         <  1:1030716/68‑1 (MQ=255)
                                                           tGGATACCGTAATGGGACCGCTGACCTGGGATGAGAAAGGCGATCTGAAAGGCTTTGAGTTCGGCGTa         <  1:126400/68‑1 (MQ=255)
                                                           tGGATACCGTAATGGGACCGCTGACCTGGGATGAGAAAGGCGATCTGAAAGGCTTTGAGTTCGGCGTa         <  1:3331807/68‑1 (MQ=255)
                                                           tGGATACCGTAATGGGACCGCTGACCTGGGATGAGAAAGGCGATCTGAAAGGCTTTGAGTTCGGCGTa         <  1:2258840/68‑1 (MQ=255)
                                                           tGGATACCGTAATGGGACCGCTGACCTGGGATGAGAAAGGCGATCTGAAAGGCTTTGAGTTCGGCGTa         <  1:498541/68‑1 (MQ=255)
                                                           tGGATACCGTAATGGGACCGCTGACCTGGGATGAGAAAGGCGATCTGAAAGGCTTTGAGTTCGGCGTa         <  1:808746/68‑1 (MQ=255)
                                                            ggATACCGTAATGGGACCGCTGACCTGGGATGAGAAAGGCGATCTGAAAGGCTTTGAGTTCGGCGTa         <  1:570866/67‑1 (MQ=255)
                                                               tACCGTAATGGGACCGCTGACCTGGGATGAGAAAGGCGATCTGAAAGGCTTTGAGTTCGGCTTATTTg     <  1:441742/68‑1 (MQ=255)
                                                               tACCGTAATGGGACCGCTGACCTGGGATGAGAAAGGCGATCTGAAAGGCTTTGAGTTCGGCGTATTTg     <  1:1503317/68‑1 (MQ=255)
                                                               tACCGTAATGGGACCGCTGACCTGGGATGAGAAAGGCGATCTGAAAGGCTTTGAGTTCGGCGTATTTg     <  1:371377/68‑1 (MQ=255)
                                                               tACCGTAATGGGACCGCTGACCTGGGATGAGAAAGGCGATCTGAAAGGCTTTGAGTTCGGCGTATTTg     <  1:2065407/68‑1 (MQ=255)
                                                               tACCGTAATGGGACCGCTGACCTGGGATGAGAAAGGCGATCTGAAAGGCTTTGAGTTCGGCGTATTTg     <  1:1307741/68‑1 (MQ=255)
                                                               tACCGTAATGGGACCGCTGACCTGGGATGAGAAAGGCGATCTGAAAGGCTTTGAGTTCGGCGTATTTg     <  1:3374488/68‑1 (MQ=255)
                                                               tACCGTAATGGGACCGCTGACCTGGGATGAGAAAGGCGATCTGAAAGGCTTTGAGTTCGGCGTATTTg     <  1:1270036/68‑1 (MQ=255)
                                                               tACCGTAATGGGACCGCTGACCTGGGATGAGAAAGGCGATCTGAAAGGCTTTGAGTTCGGCGTATTTg     <  1:2357758/68‑1 (MQ=255)
                                                               tACCGTAATGGGACCGCTGACCTGGGATGAGAAAGGCGATCTGAAAGGCTTTGAGTTCGGCGTATTTg     <  1:2921748/68‑1 (MQ=255)
                                                               tACCGTAATGGGACCGCTGACCTGGGATGAGAAAGGCGATCTGAAAGGCTTTGAGTTCGGCGTATTTg     <  1:243222/68‑1 (MQ=255)
                                                                aCCGTAATGGGACCGCTGACCTGGGATGAGAAAGGCGATCTGAAAGGCTTTGAGTTCGGCGTATTTg     <  1:2657204/67‑1 (MQ=255)
                                                                aCCGTAATGGGACCGCTGACCTGGGATGAGAAAGGCGATCTGAAAGGCTTTGAGTTCGGCGTATTTg     <  1:1969262/67‑1 (MQ=255)
                                                                aCCGTAATGGGACCGCTGACCTGGGATGAGAAAGGCGATCTGAAAGGCTTTGAGTTCGGCGTATTTGAc   >  1:2072883/1‑69 (MQ=255)
                                                                aCCGTAATGGGACCGCTGACCTGGGATGAGAAAGGCGATCTGAAAGGCTTTGAGTTCGGCGTATTTGAc   >  1:1435206/1‑69 (MQ=255)
                                                                aCCGTAATGGGACCGCTGACCTGGGATGAGAAAGGCGATCTGAAAGGCTTTGAGTTCGGCGTATTTGAc   >  1:2656872/1‑69 (MQ=255)
                                                                aCCGTAATGGGACCGCTGACCTGGGATGAGAAAGGCGATCTGAAAGGCTTTGAGTTCGGCGTATTTGAc   >  1:2510062/1‑69 (MQ=255)
                                                                 ccGTAATGGGACCGCTGACCTGGGATGAGAAAGGCGATCTGAAAGGCTTTGAGTTCGGCGTATTTGACt  >  1:3111246/1‑69 (MQ=255)
                                                                  |                                                                   
GGGCCTGAATCAGTCTGACGATCCGGCTGAAATCGCCAAATACCTGAAAGCGAACTCCGTGGATACCGTAATGGGACCGCTGACCTGGGATGAGAAAGGCGATCTGAAAGGCTTTGAGTTCGGCGTATTTGACT  >  minE/2628466‑2628599

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: