Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,130,656 C→G 11.1% L83L (CTG→CTC) ‡ mtr ← tryptophan transporter of high affinity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,130,6560CG11.1% 287.2 / 22.3 108L83L (CTG→CTC) ‡mtrtryptophan transporter of high affinity
Reads supporting (aligned to +/- strand):  ref base C (44/51);  new base G (5/7);  total (50/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

Only 100 of 115 total aligned reads displayed.

TATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCAGGTTAGCTTCCAGAATCATC  >  minE/2130586‑2130725
                                                                      |                                                                     
tataGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGccc                                                                       <  1:2750719/71‑1 (MQ=255)
 ataGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGccc                                                                       <  1:3340983/70‑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)
           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)
                gATATAGAGAACGAAGGCAACGGTAATACCGCTGATAATGTTCCAGGTGTTACCGATTAAATCTTTGGTGa                                                       >  1:3970433/1‑71 (MQ=255)
                 atatAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGAt                                                      <  1:2089740/71‑1 (MQ=255)
                  tataGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATg                                                     >  1:1344534/1‑71 (MQ=255)
                    taGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTg                                                        <  1:841671/66‑1 (MQ=255)
                    taGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTg                                                        >  1:1008746/1‑66 (MQ=255)
                    taGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGa                                                       <  1:3855934/67‑1 (MQ=255)
                    taGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGa                                                       <  1:5029024/67‑1 (MQ=255)
                    taGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGa                                                       <  1:5022463/67‑1 (MQ=255)
                      gagCACATAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGa                                                       <  1:2824126/65‑1 (MQ=255)
                      gagCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCATATCttt                                                            <  1:11791/60‑1 (MQ=255)
                      gagCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCttt                                                            <  1:4053417/60‑1 (MQ=255)
                      gagCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCttt                                                            <  1:1146724/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:4613688/71‑1 (MQ=255)
                         cacaAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCaa                                              <  1:3111314/71‑1 (MQ=255)
                          acgaaGGCAACGGTAATACCGCTGATAATGTTCCAGGTGTTACCGATTAAATCTTTGGTGATGGTGTTaa                                              <  1:4131665/67‑1 (MQ=255)
                          acgaaGGCAACGGTAATACCGCTGATAATGTTCCAGGTGTTACCGATTAAATCTTTGGTGATGGTGTTaa                                              <  1:5535236/67‑1 (MQ=255)
                          acaAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCaa                                              <  1:3226932/70‑1 (MQ=255)
                           cgaaGGCAACGGTAATACCGCTGATAATGTTCCAGGTGTTACCGATTAAATCTTTGGTGATGGTGTTaaaa                                            >  1:190574/3‑71 (MQ=255)
                           cgaaGGCAACGGTAATACCGCTGATAATGTTCCAGGTGTTACCGATTAAATCTTTGGTGATGGTGTTaaaa                                            >  1:1372428/3‑71 (MQ=255)
                           cgaaGGCAACGGTAATACCGCTGATAATGTTCCAGGTGTTACCGATTAAATCTTTGGTGATGGTGTTaaaa                                            <  1:223424/69‑1 (MQ=255)
                           caAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATg                                                     <  1:3569765/62‑1 (MQ=255)
                           caAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATg                                                     <  1:307056/62‑1 (MQ=255)
                           caAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCaaaa                                            <  1:4942605/71‑1 (MQ=255)
                             aaGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCaaaa                                            >  1:2756128/1‑69 (MQ=255)
                              aGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAAt                                                                >  1:5260916/1‑48 (MQ=255)
                              aGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAAt                                                                >  1:5668755/1‑48 (MQ=255)
                              aGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAAt                                                                >  1:3630603/1‑48 (MQ=255)
                              aGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATc                                                               <  1:4633698/49‑1 (MQ=255)
                              aGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATc                                                               <  1:1220998/49‑1 (MQ=255)
                              aGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATc                                                               <  1:721850/49‑1 (MQ=255)
                              aGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTc                                         >  1:4489128/1‑71 (MQ=255)
                              aGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTc                                         <  1:2801628/71‑1 (MQ=255)
                              aGGCAATGGAAATGCCGTTAACCACGTTCCAGCCTTTGCCCAGCAAAt                                                                >  1:2848161/1‑48 (MQ=255)
                              aGGCAACGGTAATACCGCTGATAATGTTCCAGGTGTTACCGATTAAATCTTTGGTGATGgtgtta                                               >  1:4137638/1‑63 (MQ=255)
                              aGGCAACGGTAATACCGCTGATAATGTTCCAGGTGTTACCGATTAAATCTTTGGTGATGGTGTTAAAACtg                                         <  1:3774747/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:5610506/1‑51 (MQ=255)
                                  aaTGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGt                                                         >  1:5439009/1‑51 (MQ=255)
                                  aaTGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGt                                                         >  1:4753208/1‑51 (MQ=255)
                                  aaTGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGt                                                         >  1:940966/1‑51 (MQ=255)
                                   aTGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCaaa                                                                 >  1:1236361/1‑42 (MQ=255)
                                   aTGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCaaa                                                                 >  1:5349269/1‑42 (MQ=255)
                                   aTGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCt                                                              <  1:4840014/45‑1 (MQ=255)
                                   aTGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCt                                                              <  1:2854481/45‑1 (MQ=255)
                                   aTGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTgg                                                          <  1:1796894/49‑1 (MQ=255)
                                   aTGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTgg                                                          <  1:1938430/49‑1 (MQ=255)
                                   aTGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTgg                                                          <  1:1941406/49‑1 (MQ=255)
                                    tGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCtt                                                             <  1:1311217/45‑1 (MQ=255)
                                    tGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTg                                                           >  1:3599880/1‑47 (MQ=255)
                                    tGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGt                                                         <  1:1357031/49‑1 (MQ=255)
                                     ggAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGc                                                                    >  1:3105797/1‑37 (MQ=255)
                                     ggAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGc                                                                    >  1:883130/1‑37 (MQ=255)
                                     ggAAATGCCGTTGACCACGTTCCAGCCTTTGCCAAGCAAATCttt                                                            >  1:1344116/1‑45 (MQ=255)
                                       aaaTGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGAtt                                >  1:2521635/1‑71 (MQ=255)
                                       aaaTGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGAtt                                <  1:4854646/71‑1 (MQ=255)
                                        aaTGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCg                                        <  1:2759106/62‑1 (MQ=255)
                                          tGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTc                               <  1:1520045/69‑1 (MQ=255)
                                          tGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTc                               <  1:2190586/69‑1 (MQ=255)
                                           gCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTg                             >  1:4191754/1‑70 (MQ=255)
                                            ccGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGt                            <  1:3230345/70‑1 (MQ=255)
                                            ccGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGt                            >  1:5087710/1‑70 (MQ=255)
                                             cGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCaa                                              <  1:268737/51‑1 (MQ=255)
                                              gTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCaaaa                                            >  1:3957565/1‑52 (MQ=255)
                                              gTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAAcc                                    <  1:2726522/60‑1 (MQ=255)
                                                 gACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTc                       <  1:126705/70‑1 (MQ=255)
                                                    cACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCAGGt                   >  1:869825/1‑71 (MQ=255)
                                                      cGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCAGGTTa                 <  1:3566885/71‑1 (MQ=255)
                                                      cGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCAGGTTa                 <  1:1545702/71‑1 (MQ=255)
                                                             gCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCAGGTTAGCTTCCa          <  1:481581/71‑1 (MQ=255)
                                                              ccTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTc                       >  1:3424217/1‑57 (MQ=255)
                                                              ccTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCa                      >  1:2062936/1‑58 (MQ=255)
                                                              ccTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCa                      >  1:4726282/1‑58 (MQ=255)
                                                                tttGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTg                             <  1:2077361/49‑1 (MQ=255)
                                                                  tGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCAGGTTAGCTTCCAGAat      <  1:5644314/70‑1 (MQ=255)
                                                                  tGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCAGGTTAGCTTCCAGAat      <  1:2934140/70‑1 (MQ=255)
                                                                    cccAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGAtt                                <  1:3203028/42‑1 (MQ=255)
                                                                    cccAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCAGGTTAGCTTCCAGAatcat   <  1:560616/71‑1 (MQ=255)
                                                                    cccAGCAAATCCTTGGTGATGGTGTCAAAACTCGAACCGAtt                                <  1:5203672/42‑1 (MQ=255)
                                                                     ccAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCAGGTTAGCTTCCAGAatcatc  <  1:3729583/71‑1 (MQ=255)
                                                                     ccAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCAGGTTAGCTTCCAGAatcat   >  1:3842728/1‑70 (MQ=255)
                                                                      |                                                                     
TATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCAGGTTAGCTTCCAGAATCATC  >  minE/2130586‑2130725

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: