Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,130,648 C→G 7.1% G86A (GGC→GCC) ‡ mtr ← tryptophan transporter of high affinity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,130,6480CG 7.1% 292.5 / 8.8 98G86A (GGC→GCC) ‡mtrtryptophan transporter of high affinity
Reads supporting (aligned to +/- strand):  ref base C (44/47);  new base G (5/2);  total (49/49)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.02e-01

Only 100 of 104 total aligned reads displayed.

CAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCAGGTTAGCTTCCAG  >  minE/2130580‑2130718
                                                                    |                                                                      
cAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCt                                                                      <  1:3900250/71‑1 (MQ=255)
cAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCt                                                                      <  1:631903/71‑1 (MQ=255)
 aGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCtt                                                                     <  1:3786149/71‑1 (MQ=255)
 aGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCt                                                                      <  1:476426/70‑1 (MQ=255)
   aaaTATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTg                                                                   <  1:2108386/71‑1 (MQ=255)
      tataGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGcc                                                                 <  1:1300912/70‑1 (MQ=255)
        taGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCa                                                               <  1:205213/70‑1 (MQ=255)
           gCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCaa                                                           <  1:2968351/71‑1 (MQ=255)
             aTAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGc                                                             <  1:3927449/67‑1 (MQ=255)
             aTAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAAt                                                         <  1:572894/71‑1 (MQ=255)
              tAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATc                                                        >  1:79700/1‑71 (MQ=255)
               aGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCt                                                       >  1:342043/1‑71 (MQ=255)
               aGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCt                                                       >  1:680087/1‑71 (MQ=255)
                ggTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCtt                                                      <  1:2775292/71‑1 (MQ=255)
                ggTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCt                                                       <  1:1548219/70‑1 (MQ=255)
                 gTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCttt                                                     <  1:2512851/71‑1 (MQ=255)
                 gTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCttt                                                     <  1:2810825/71‑1 (MQ=255)
                  tCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTg                                                    <  1:2372374/71‑1 (MQ=255)
                    aGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGc                                                                        <  1:2374901/49‑1 (MQ=255)
                    aGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGc                                                                        <  1:2156590/49‑1 (MQ=255)
                    aGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGc                                                                        <  1:596894/49‑1 (MQ=255)
                     gtatatatAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCaaa                                                          <  1:769920/60‑1 (MQ=255)
                      gATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGc                                                                        <  1:1485245/47‑1 (MQ=255)
                       atatAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGAt                                               <  1:3234720/71‑1 (MQ=255)
                       atatAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGAt                                               >  1:1591118/1‑71 (MQ=255)
                        tataGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCaa                                                           <  1:3151871/58‑1 (MQ=255)
                         ataGAGAACGAAGGCAACGGTAATACCGCTGATAATGTTCCAGGTGTTACCGATTAAATCTTTGGTGATgg                                             >  1:2353629/1‑71 (MQ=255)
                         ataGAGAACGAAGGCAACGGTAATACCGCTGATAATGTTCCAGGTGTTACCGATTAAATCTTTGGTGATgg                                             >  1:398382/1‑71 (MQ=255)
                          taGAGCACTAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGgt                                            <  1:177020/71‑1 (MQ=255)
                          taGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATgg                                             <  1:3634079/70‑1 (MQ=255)
                          taGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGgt                                            >  1:1700784/1‑71 (MQ=255)
                          taGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGgt                                            >  1:1037977/1‑71 (MQ=255)
                          taGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGgt                                            <  1:2526857/71‑1 (MQ=255)
                          taGAGAACGAAGGCAACGGTAATACCGCTGATAATGTTCCAGGTGTTACCGATTAAATCTTTGGTGATGgt                                            <  1:1728274/71‑1 (MQ=255)
                           agagCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGgt                                            >  1:378976/1‑70 (MQ=255)
                            gagCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGgtg                                           <  1:2954836/70‑1 (MQ=255)
                              gCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTc                                         <  1:2107454/70‑1 (MQ=255)
                               cacaAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGAt                                               <  1:58497/63‑1 (MQ=255)
                               cacaAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCaa                                       <  1:1773326/71‑1 (MQ=255)
                               cacaAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCaa                                       <  1:507759/71‑1 (MQ=255)
                               cacaAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCaa                                       <  1:3182582/71‑1 (MQ=255)
                                acgaaGGCAACGGTAATACCGCTGATAATGTTCCAGGTGTTACCGATTAAATCTTTGGTGATGGTGTTaa                                       >  1:620774/4‑70 (MQ=255)
                                acgaaGGCAACGGTAATACCGCTGATAATGTTCCAGGTGTTACCGATTAAATCTTTGGTGATGGTGTTaa                                       >  1:2643147/4‑70 (MQ=255)
                                acaAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCaa                                       >  1:3648120/1‑70 (MQ=255)
                                acaAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCaa                                       <  1:2348795/70‑1 (MQ=255)
                                 caAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCttt                                                     >  1:3412681/1‑55 (MQ=255)
                                 caAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCttt                                                     >  1:2175401/1‑55 (MQ=255)
                                 caAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCaaaa                                     >  1:369498/1‑71 (MQ=255)
                                 caAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCaaaa                                     <  1:1876149/71‑1 (MQ=255)
                                  gaaGGCAACGGTAATACCGCTGATAATGTTCCAGGTGTTACCGATTAAATCTTTGGTGTTGGTGTTaaaa                                     <  1:3035046/69‑1 (MQ=255)
                                   aaGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTg                                                    >  1:3254630/1‑54 (MQ=255)
                                   aaGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTg                                                    >  1:1053184/1‑54 (MQ=255)
                                   aaGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACt                                   <  1:2386780/71‑1 (MQ=255)
                                    aGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTc                                  <  1:2069571/71‑1 (MQ=255)
                                    aGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTc                                  >  1:3152333/1‑71 (MQ=255)
                                    aGGCAACGGTAATACCGCTGATAATGTTCCAGGTGTTACCGATTAAATCTTTGGTGATGGTGTTAAAACtg                                  >  1:643712/1‑70 (MQ=255)
                                      gCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCtt                                                      <  1:837689/49‑1 (MQ=255)
                                      gCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGa                                <  1:1154618/71‑1 (MQ=255)
                                        aaTGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTg                                                    >  1:1236744/1‑49 (MQ=255)
                                         aTGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCaa                                       <  1:467502/61‑1 (MQ=255)
                                         aTGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAAcc                             <  1:1703473/71‑1 (MQ=255)
                                         aTGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAAcc                             <  1:2346839/71‑1 (MQ=255)
                                         aTGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAAcc                             >  1:2725876/1‑71 (MQ=255)
                                         aTGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAAcc                             >  1:1096749/1‑71 (MQ=255)
                                         aTGGAAATGACGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAAcc                             >  1:787359/1‑71 (MQ=255)
                                          tGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACt                                   >  1:696586/1‑64 (MQ=255)
                                          tGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACt                                   >  1:2888344/1‑64 (MQ=255)
                                          tGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAAcc                             >  1:2521322/1‑70 (MQ=255)
                                           ggAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGgtgt                                          >  1:3820079/1‑56 (MQ=255)
                                           ggAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGa                           >  1:3016953/1‑71 (MQ=255)
                                            gAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAAt                                                         >  1:2238315/1‑40 (MQ=255)
                                            gAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAAt                                                         >  1:1531108/1‑40 (MQ=255)
                                            gAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACt                                   >  1:3011779/1‑62 (MQ=255)
                                                 gCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGaa                               >  1:2631467/1‑61 (MQ=255)
                                                  ccGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTc                                  <  1:2119005/57‑1 (MQ=255)
                                                  ccGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGt                     >  1:3055926/1‑70 (MQ=255)
                                                   cGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTg                      <  1:3256697/68‑1 (MQ=255)
                                                   cGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTaa                   >  1:1615672/1‑71 (MQ=255)
                                                   cGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTaa                   >  1:188991/1‑71 (MQ=255)
                                                     ttGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAAtt                 <  1:1675423/71‑1 (MQ=255)
                                                      tGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAAtt                 >  1:1418809/1‑70 (MQ=255)
                                                      tGACCACGTTCCAGCCTTTGCCCAGCAAATCTTCGTGATGGTGTCAAAACTCGAACCGATTCTGTAAtt                 >  1:86946/1‑69 (MQ=39)
                                                        aCCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCAg              <  1:1114843/71‑1 (MQ=255)
                                                         ccACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCAgg             <  1:570598/71‑1 (MQ=255)
                                                           aCGTTCCAGCCTTTGCCCAGCAAATCTTAGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCACGt            >  1:3297985/1‑70 (MQ=255)
                                                            cGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGa                           <  1:3722550/54‑1 (MQ=255)
                                                            cGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTaa                   >  1:3629755/1‑62 (MQ=255)
                                                            cGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCAGGt            >  1:3764707/1‑69 (MQ=255)
                                                            cGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCAGGt            >  1:683287/1‑69 (MQ=255)
                                                              ttCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCAGGTTAGc        >  1:1519517/1‑71 (MQ=255)
                                                              ttCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCAGGTTAGc        <  1:2760230/71‑1 (MQ=255)
                                                                ccAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCAGGTTAGCtt      >  1:3097267/1‑71 (MQ=255)
                                                                ccAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCAGGTTAGCt       <  1:200495/70‑1 (MQ=255)
                                                                 cAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCa               >  1:3736095/1‑61 (MQ=255)
                                                                 cAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCAGGTTAGCTTc     <  1:3153347/71‑1 (MQ=255)
                                                                  aGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCg                                 <  1:1587678/42‑1 (MQ=255)
                                                                  aGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCAGGTTAGCTTc     >  1:369586/1‑70 (MQ=255)
                                                                   gCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGaa                               >  1:1707710/1‑43 (MQ=255)
                                                                    ccTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCAGGTTAGCTTc     >  1:713591/1‑68 (MQ=255)
                                                                    ccTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCAGGTTAGCTTCCAg  <  1:1345481/71‑1 (MQ=255)
                                                                    |                                                                      
CAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCAGGTTAGCTTCCAG  >  minE/2130580‑2130718

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: