Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,130,649 C→T 11.2% G86S (GGC→AGC) ‡ mtr ← tryptophan transporter of high affinity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,130,6490CT11.2% 285.4 / 19.5 107G86S (GGC→AGC) ‡mtrtryptophan transporter of high affinity
Reads supporting (aligned to +/- strand):  ref base C (41/54);  new base T (5/7);  total (46/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

GCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCAGGTTAGCTTCCA  >  minE/2130579‑2130717
                                                                      |                                                                    
gCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGcc                                                                      >  1:5184729/1‑71 (MQ=255)
      atatAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGcc                                                                >  1:1800618/1‑71 (MQ=255)
      atatAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGcc                                                                <  1:3590723/71‑1 (MQ=255)
       tataGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGccc                                                               <  1:2750719/71‑1 (MQ=255)
        ataGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGccc                                                               <  1:3340983/70‑1 (MQ=255)
        ataGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGcc                                                                <  1:5230412/69‑1 (MQ=255)
         taGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAg                                                             >  1:17965/1‑71 (MQ=255)
            gCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCa                                                           >  1:2994315/1‑70 (MQ=255)
            gCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCa                                                           >  1:3476297/1‑70 (MQ=255)
              aTAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAAt                                                        <  1:5670630/71‑1 (MQ=255)
               tAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATc                                                       >  1:389185/1‑71 (MQ=255)
                aGGTCAGGATATAGAGCACAAAGGCGATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCt                                                      >  1:2905092/1‑71 (MQ=255)
                 ggTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCttt                                                                   <  1:1562379/57‑1 (MQ=255)
                  gTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCtt                                                     <  1:4258585/70‑1 (MQ=255)
                    cAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCCGCAAATCTTTg                                                   >  1:2796383/1‑70 (MQ=255)
                     aGGATATAGAGAACGAAGGCAACGGTAATACCGCTGATAATGTTCCAGGTGTTACCGATTAAATCTTTGGt                                                 >  1:3723559/1‑71 (MQ=255)
                      ggATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTg                                                <  1:4305190/71‑1 (MQ=255)
                       gATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGa                                               >  1:5392089/1‑71 (MQ=255)
                       gATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGa                                               <  1:2813722/71‑1 (MQ=255)
                       gATATAGAGAACGAAGGCAACGGTAATACCGCTGATAATGTTCCAGGTGTTACCGATTAAATCTTTGGTGa                                               >  1:3970433/1‑71 (MQ=255)
                        atatAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGAt                                              <  1:2089740/71‑1 (MQ=255)
                         tataGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATg                                             >  1:1344534/1‑71 (MQ=255)
                           taGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTg                                                >  1:2056723/1‑66 (MQ=255)
                           taGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTg                                                <  1:841671/66‑1 (MQ=255)
                           taGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTg                                                >  1:1008746/1‑66 (MQ=255)
                           taGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGa                                               <  1:5022463/67‑1 (MQ=255)
                           taGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGa                                               <  1:5029024/67‑1 (MQ=255)
                             gagCACATAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGa                                               <  1:2824126/65‑1 (MQ=255)
                             gagCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCATATCttt                                                    <  1:11791/60‑1 (MQ=255)
                             gagCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCttt                                                    <  1:1146724/60‑1 (MQ=255)
                             gagCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCttt                                                    <  1:4053417/60‑1 (MQ=255)
                             gagCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCttt                                                    <  1:1267016/60‑1 (MQ=255)
                             gagCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGgtg                                          >  1:836607/1‑70 (MQ=255)
                               gCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTgg                                                  >  1:4798413/1‑60 (MQ=255)
                                cacaAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCaa                                      <  1:2050087/71‑1 (MQ=255)
                                cacaAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCaa                                      <  1:3111314/71‑1 (MQ=255)
                                cacaAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCaa                                      <  1:4613688/71‑1 (MQ=255)
                                 acgaaGGCAACGGTAATACCGCTGATAATGTTCCAGGTGTTACCGATTAAATCTTTGGTGATGGTGTTaa                                      <  1:4131665/67‑1 (MQ=255)
                                 acgaaGGCAACGGTAATACCGCTGATAATGTTCCAGGTGTTACCGATTAAATCTTTGGTGATGGTGTTaa                                      <  1:5010047/67‑1 (MQ=255)
                                 acgaaGGCAACGGTAATACCGCTGATAATGTTCCAGGTGTTACCGATTAAATCTTTGGTGATGGTGTTaa                                      <  1:5535236/67‑1 (MQ=255)
                                 acaAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCttt                                                                   <  1:2357421/41‑1 (MQ=255)
                                 acaAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCaa                                      <  1:3226932/70‑1 (MQ=255)
                                  cgaaGGCAACGGTAATACCGCTGATAATGTTCCAGGTGTTACCGATTAAATCTTTGGTGATGGTGTTaaaa                                    <  1:223424/69‑1 (MQ=255)
                                  cgaaGGCAACGGTAATACCGCTGATAATGTTCCAGGTGTTACCGATTAAATCTTTGGTGATGGTGTTaaaa                                    >  1:190574/3‑71 (MQ=255)
                                  cgaaGGCAACGGTAATACCGCTGATAATGTTCCAGGTGTTACCGATTAAATCTTTGGTGATGGTGTTaaaa                                    >  1:1372428/3‑71 (MQ=255)
                                  caAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATg                                             <  1:307056/62‑1 (MQ=255)
                                  caAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATg                                             <  1:3569765/62‑1 (MQ=255)
                                  caAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCaaaa                                    <  1:4942605/71‑1 (MQ=255)
                                    aaGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCaaaa                                    >  1:2756128/1‑69 (MQ=255)
                                     aGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAAt                                                        >  1:3630603/1‑48 (MQ=255)
                                     aGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAAt                                                        >  1:5260916/1‑48 (MQ=255)
                                     aGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATc                                                       <  1:1220998/49‑1 (MQ=255)
                                     aGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATc                                                       <  1:721850/49‑1 (MQ=255)
                                     aGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATc                                                       <  1:4633698/49‑1 (MQ=255)
                                     aGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATc                                                       <  1:4152227/49‑1 (MQ=255)
                                     aGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTc                                 <  1:2801628/71‑1 (MQ=255)
                                     aGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTc                                 >  1:4489128/1‑71 (MQ=255)
                                     aGGCAATGGAAATGCCGTTAACCACGTTCCAGCCTTTGCCCAGCAAAt                                                        >  1:2848161/1‑48 (MQ=255)
                                     aGGCAACGGTAATACCGCTGATAATGTTCCAGGTGTTACCGATTAAATCTTTGGTGATGGTGTTAAAACtg                                 <  1:3774747/71‑2 (MQ=255)
                                     aGGCAACGGTAATACCGCTGATAATGTTCCAGGTGTTACCGATTAAATCTTTGGTGATGGTGTTAAAACtg                                 <  1:2681888/71‑2 (MQ=255)
                                      ggCAACGGTAATACCGCTGATAATGTTCCAGGTGTTACCGATTAAATCTTTGGTGATGGTGTTAAAACtgg                                <  1:2852568/71‑3 (MQ=255)
                                       gCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTg                                                   >  1:2372887/1‑51 (MQ=255)
                                       gCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGa                                               >  1:3868890/1‑55 (MQ=255)
                                         aaTGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCtt                                                     <  1:3966389/47‑1 (MQ=255)
                                         aaTGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGt                                                 >  1:5439009/1‑51 (MQ=255)
                                         aaTGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGt                                                 >  1:4753208/1‑51 (MQ=255)
                                         aaTGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGt                                                 >  1:940966/1‑51 (MQ=255)
                                         aaTGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGt                                                 >  1:5610506/1‑51 (MQ=255)
                                          aTGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCaaa                                                         >  1:1236361/1‑42 (MQ=255)
                                          aTGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCaaa                                                         >  1:5349269/1‑42 (MQ=255)
                                          aTGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCt                                                      <  1:3890396/45‑1 (MQ=255)
                                          aTGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCt                                                      <  1:4840014/45‑1 (MQ=255)
                                          aTGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCt                                                      <  1:2854481/45‑1 (MQ=255)
                                          aTGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTgg                                                  <  1:1938430/49‑1 (MQ=255)
                                          aTGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTgg                                                  <  1:1941406/49‑1 (MQ=255)
                                           tGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCtt                                                     <  1:1311217/45‑1 (MQ=255)
                                           tGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCtt                                                     >  1:644446/1‑45 (MQ=255)
                                           tGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTg                                                   >  1:3599880/1‑47 (MQ=255)
                                           tGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGt                                                 <  1:1357031/49‑1 (MQ=255)
                                            ggAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGc                                                            >  1:883130/1‑37 (MQ=255)
                                            ggAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGc                                                            >  1:3105797/1‑37 (MQ=255)
                                              aaaTGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGAtt                        <  1:4854646/71‑1 (MQ=255)
                                              aaaTGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGAtt                        >  1:2521635/1‑71 (MQ=255)
                                                 tGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCaa                                      <  1:1296079/54‑1 (MQ=255)
                                                 tGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTc                       <  1:1520045/69‑1 (MQ=255)
                                                  gCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTg                     >  1:4191754/1‑70 (MQ=255)
                                                   ccGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGt                    <  1:3230345/70‑1 (MQ=255)
                                                    cGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCaa                                      <  1:268737/51‑1 (MQ=255)
                                                     gTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCaaaa                                    >  1:3957565/1‑52 (MQ=255)
                                                     gTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAAcc                            <  1:3655996/60‑1 (MQ=255)
                                                     gTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAAcc                            <  1:2726522/60‑1 (MQ=255)
                                                        gACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTc               <  1:126705/70‑1 (MQ=255)
                                                           cACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCAGGt           >  1:869825/1‑71 (MQ=255)
                                                             cGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCAGGt           >  1:1774955/1‑69 (MQ=255)
                                                             cGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCAGGTTa         <  1:1545702/71‑1 (MQ=255)
                                                             cGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCAGGTTa         <  1:3566885/71‑1 (MQ=255)
                                                                    gCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCAGGTTAGCTTCCa  <  1:481581/71‑1 (MQ=255)
                                                                     ccTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCa              >  1:4726282/1‑58 (MQ=255)
                                                                     ccTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCa              >  1:2062936/1‑58 (MQ=255)
                                                                     ccTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCa              >  1:3231600/1‑58 (MQ=255)
                                                                      |                                                                    
GCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCAGGTTAGCTTCCA  >  minE/2130579‑2130717

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: