Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,130,637 C→T 8.5% V90I (GTC→ATC)  mtr ← tryptophan transporter of high affinity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,130,6370CT 8.5% 278.0 / 11.8 94V90I (GTC→ATC) mtrtryptophan transporter of high affinity
Reads supporting (aligned to +/- strand):  ref base C (38/48);  new base T (5/3);  total (43/51)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

TCGAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCAGG  >  minE/2130568‑2130707
                                                                     |                                                                      
tCGAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGAc                                                                        >  1:1333389/1‑70 (MQ=255)
tCGAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGAc                                                                        >  1:2752857/1‑70 (MQ=255)
  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: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:596894/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:1591118/1‑71 (MQ=255)
                                   atatAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGAt                                    <  1:3234720/71‑1 (MQ=255)
                                   atatAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGAt                                    >  1:3774950/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: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:2346839/71‑1 (MQ=255)
                                                     aTGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAAcc                  <  1:1703473/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: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)
                                                                     |                                                                      
TCGAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTGGTGATGGTGTCAAAACTCGAACCGATTCTGTAATTCAGG  >  minE/2130568‑2130707

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: