Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,130,640 C→T 8.5% V89M (GTG→ATG) ‡ mtr ← tryptophan transporter of high affinity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,130,6400CT 8.5% 279.0 / 11.5 94V89M (GTG→ATG) ‡mtrtryptophan transporter of high affinity
Reads supporting (aligned to +/- strand):  ref base C (37/49);  new base T (5/3);  total (42/52)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01

GAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCAGGT  >  minE/2130570‑2130708
                                                                      |                                                                    
gAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCAc                                                                      >  1:431476/1‑71 (MQ=255)
gAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCAc                                                                      >  1:3038275/1‑71 (MQ=255)
  aCCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGt                                                                    <  1:1730270/71‑1 (MQ=255)
  aCCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGt                                                                    <  1:2963379/71‑1 (MQ=255)
       tGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAg                                                               >  1:1530889/1‑71 (MQ=255)
         gCAGAGATATAGGCATAAGTGAGGATATAGAGAACGAAGGCAACGGTAATACCGCTGATAATGTTCCAggt                                                             <  1:2258372/71‑3 (MQ=255)
          cAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCt                                                            <  1:3900250/71‑1 (MQ=255)
          cAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCt                                                            <  1:631903/71‑1 (MQ=255)
           aGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTc                                                                  <  1:1768982/64‑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:2810825/71‑1 (MQ=255)
                           gTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCttt                                           <  1:2512851/71‑1 (MQ=255)
                            tCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTg                                          <  1:2372374/71‑1 (MQ=255)
                              aGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGc                                                              <  1:596894/49‑1 (MQ=255)
                              aGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGc                                                              <  1:2374901/49‑1 (MQ=255)
                              aGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGc                                                              <  1:2156590/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:3774950/1‑71 (MQ=255)
                                 atatAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGAt                                     >  1:1591118/1‑71 (MQ=255)
                                  tataGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCaa                                                 <  1:3151871/58‑1 (MQ=255)
                                   ataGAGAACGAAGGCAACGGTAATACCGCTGATAATGTTCCAGGTGTTACCGATTAAATCTTTGGTGATgg                                   >  1:398382/1‑71 (MQ=255)
                                   ataGAGAACGAAGGCAACGGTAATACCGCTGATAATGTTCCAGGTGTTACCGATTAAATCTTTGGTGATgg                                   >  1:2353629/1‑71 (MQ=255)
                                    taGAGCACTAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGgt                                  <  1:177020/71‑1 (MQ=255)
                                    taGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATgg                                   <  1:3634079/70‑1 (MQ=255)
                                    taGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGgt                                  <  1:2526857/71‑1 (MQ=255)
                                    taGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGgt                                  >  1:1037977/1‑71 (MQ=255)
                                    taGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGgt                                  >  1:1700784/1‑71 (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:2643147/4‑70 (MQ=255)
                                          acgaaGGCAACGGTAATACCGCTGATAATGTTCCAGGTGTTACCGATTAAATCTTTGGTGATGGTGTTaa                             >  1:620774/4‑70 (MQ=255)
                                          acaAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCaa                             >  1:3648120/1‑70 (MQ=255)
                                          acaAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCaa                             <  1:2348795/70‑1 (MQ=255)
                                           caAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCttt                                           >  1:2175401/1‑55 (MQ=255)
                                           caAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCttt                                           >  1:3412681/1‑55 (MQ=255)
                                           caAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCaaaa                           <  1:1876149/71‑1 (MQ=255)
                                           caAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCaaaa                           >  1:369498/1‑71 (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:3152333/1‑71 (MQ=255)
                                              aGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTc                        <  1:2069571/71‑1 (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)
                                                   aTGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGt                                        <  1:1991982/50‑1 (MQ=255)
                                                   aTGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCaa                             <  1:467502/61‑1 (MQ=255)
                                                   aTGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAAcc                   >  1:2725876/1‑71 (MQ=255)
                                                   aTGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAAcc                   >  1:1096749/1‑71 (MQ=255)
                                                   aTGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAAcc                   <  1:2346839/71‑1 (MQ=255)
                                                   aTGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAAcc                   <  1:1703473/71‑1 (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:1531108/1‑40 (MQ=255)
                                                      gAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAAt                                               >  1:2238315/1‑40 (MQ=255)
                                                      gAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACt                         >  1:889928/1‑62 (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)
                                                               ttGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAAtt       >  1:2915736/1‑71 (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)
                                                                      |                                                                    
GAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCAGGT  >  minE/2130570‑2130708

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: