Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 483,085 G→T 11.5% A97S (GCG→TCG)  tolA → membrane anchored protein in TolA‑TolQ‑TolR complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE483,0850GT11.5% 297.7 / 12.0 93A97S (GCG→TCG) tolAmembrane anchored protein in TolA‑TolQ‑TolR complex
Reads supporting (aligned to +/- strand):  ref base G (39/43);  new base T (4/7);  total (43/50)
Fisher's exact test for biased strand distribution p-value = 5.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGAAGCGTTCTGATGAACAGCGCAAGATGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCAACTTGAGAAAGAGCGGTTAGCGGCTCAGGAGCAGAAAAAGCAG  >  minE/483015‑483153
                                                                      |                                                                    
cgcgAAGCGTTCTGATGAACAGCGCAAGATGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAg                                                                      >  1:1497983/1‑71 (MQ=255)
    aaGCGTTCTGATGAACAGCGCAAGATGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGc                                                                  >  1:32015/1‑71 (MQ=255)
    aaGCGTTCTGATGAACAGCGCAAGATGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGc                                                                  >  1:2431597/1‑71 (MQ=255)
      gCGTTCTGATGAACAGCGCAAGATGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTg                                                                <  1:1008541/71‑1 (MQ=255)
      gCGTTCTGATGAACAGCGCAAGATGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTg                                                                <  1:2507268/71‑1 (MQ=255)
      gCGTTCTGATGAACAGCGCAAGATGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTg                                                                <  1:2217870/71‑1 (MQ=255)
        gTTCTGATGAACAGCGCAAGATGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGaa                                                              <  1:4046058/71‑1 (MQ=255)
         ttCTGATGAACAGCGCAAGATGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGc                                                                     <  1:2799106/63‑1 (MQ=255)
         ttCTGATGAACAGCGCAAGATGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGaa                                                              >  1:1867434/1‑70 (MQ=255)
           cTGATGAACAGCGCAAGATGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGaa                                                              <  1:1066599/68‑1 (MQ=255)
           cTGATGAACAGCGCAAGATGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGaa                                                              <  1:332910/68‑1 (MQ=255)
           cTGATGAACAGCGCAAGATGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGAACa                                                            <  1:406893/70‑1 (MQ=255)
           cTGATGAACAGCGCAAGATGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGAACAg                                                           >  1:2276550/1‑71 (MQ=255)
           cTGATGAACAGCGCAAGATGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGAACAg                                                           >  1:2473065/1‑71 (MQ=255)
           cTGATGAACAGCGCAAGATGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAACAAGCGGCTGAACAg                                                           >  1:120830/1‑70 (MQ=255)
             gatgaACAGCGCAAGATGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGAACAg                                                           <  1:792395/69‑1 (MQ=255)
              atgaACAGCGCAAGATGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGAACAgg                                                          <  1:3782009/69‑1 (MQ=255)
               tgaACAGCGCAAGATGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGAACAGGAAc                                                       >  1:292578/1‑71 (MQ=255)
                  aCAGCGCAAGATGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAATCGGCTGAACAGGAACGcc                                                    <  1:769777/71‑1 (MQ=255)
                   cAGCGCAAGATGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCt                                                   >  1:1907347/1‑71 (MQ=255)
                     gcgcAAGATGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGAAc                                                             <  1:2261317/59‑1 (MQ=255)
                     gcgcAAGATGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGAAc                                                             <  1:2739697/59‑1 (MQ=255)
                      cgcAAGATGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTg                                                  <  1:3086972/69‑1 (MQ=255)
                      cgcAAGATGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTg                                                  >  1:3390353/1‑69 (MQ=255)
                      cgcAAGATGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTg                                                  >  1:3801888/1‑69 (MQ=255)
                      cgcAAGATGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGaa                                                <  1:576094/71‑1 (MQ=255)
                         aaGATGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCa                                             <  1:1322597/71‑1 (MQ=255)
                           gATGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGAAc                                                             <  1:3884079/53‑1 (MQ=255)
                           gATGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGAACAGGaa                                                        >  1:2529662/1‑58 (MQ=255)
                           gATGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGAACAGGaa                                                        >  1:3632287/1‑58 (MQ=255)
                           gATGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGAACAGGAACGc                                                     >  1:968882/1‑61 (MQ=255)
                           gATGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCa                                             <  1:2305067/69‑1 (MQ=255)
                           gATGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCa                                             <  1:1772861/69‑1 (MQ=255)
                           gATGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCAAc                                           <  1:3517575/71‑1 (MQ=255)
                           gATGAAGGAACAGCAGGCTGCTCAAGAACTGAGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCaa                                            >  1:1341677/1‑70 (MQ=255)
                            aTGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGaa                                                              >  1:2096458/1‑51 (MQ=255)
                             tGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAATCGGCTGAACAGGAACGCCTGAAGCAACtt                                         <  1:4058125/71‑1 (MQ=255)
                             tGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAATCGGCTGAACAGGAACGCCTGAAGCAACtt                                         <  1:3797890/71‑1 (MQ=255)
                             tGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCAACtt                                         >  1:2323840/1‑71 (MQ=255)
                             tGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCAACtt                                         >  1:756138/1‑71 (MQ=255)
                             tGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCAACtt                                         <  1:3205433/71‑1 (MQ=255)
                             tGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCAACt                                          >  1:3644426/1‑70 (MQ=255)
                              gAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAATCGGCTGAACAGGAACGCCTGAAGCAACtt                                         >  1:468046/1‑70 (MQ=255)
                                 ggAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCAACTTgaga                                     <  1:2115760/71‑1 (MQ=255)
                                 ggAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCAACTTgaga                                     <  1:3401906/71‑1 (MQ=255)
                                   aaCAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGAACAGGAAct                                                      >  1:1855739/1‑51 (MQ=255)
                                   aaCAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCAACTTg                                        >  1:1850953/1‑66 (MQ=255)
                                   aaCAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCAACTTg                                        >  1:3186739/1‑66 (MQ=255)
                                   aaCAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCAACTTGAGaaa                                   <  1:2600902/71‑1 (MQ=255)
                                       gcagGCTGCTGAAGAACTGCGTGAGAAACAATCGGCTGAACAGGAACGCCTGAAGCAACTTGAGAAAGAGc                               <  1:2560340/71‑1 (MQ=255)
                                        cagGCTGCTGAAGAACTGCGTGAGAAACAATCGGCTGAACAGGAACGCCTGAAGCAACTTGAGAAAGAGCg                              <  1:422553/71‑1 (MQ=255)
                                        cagGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCAACTTGAGAAAGAGCg                              <  1:1728808/71‑1 (MQ=255)
                                        cagGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCAACTTGAGAAAGAGCg                              <  1:2785604/71‑1 (MQ=255)
                                           gctgctGAAGAACTGCGTGAGAAACAATCGGCTGAACAGGAACGCCTGAAGCAACTTgaga                                     <  1:294834/61‑1 (MQ=255)
                                           gctgctGAAGAACTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCAACTTGAGGAAGAGCGGtt                           <  1:1070019/71‑1 (MQ=255)
                                           gctgctGAAGAACTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCAACTTGAGAAAGAGc                               <  1:416466/67‑1 (MQ=255)
                                           gctgctGAAGAACTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCAACTTGAGAAAGAGCGGtt                           <  1:1642867/71‑1 (MQ=255)
                                           gctgctGAAGAACTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCAACTTGAGAAAGAGCGGtt                           <  1:3158267/71‑1 (MQ=255)
                                             tgctgAAGAACTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCAACTTGAGAAAGAGCgg                             >  1:3242661/1‑67 (MQ=255)
                                                tgAAGAACTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCAACTTGAGAAag                                  <  1:3409987/59‑1 (MQ=255)
                                                tgAAGAACTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCAACTTGAGAAAGAGCGGTTAGCgg                      >  1:1940276/1‑71 (MQ=255)
                                                 gaagaaCTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCAAc                                           >  1:2422596/1‑49 (MQ=255)
                                                  aagaaCTGCGTGAGAAACAATCGGCTGAACAGGAACGCCTGAAGCAACTTGAGAAAGAGCGGTTAGCGGCt                    >  1:2325543/1‑71 (MQ=255)
                                                  aagaaCTGCGTGAGAAACAAGCGGCTGAACAGGAACGc                                                     >  1:3423635/1‑38 (MQ=255)
                                                  aagaaCTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCAACTTga                                       <  1:2805216/52‑1 (MQ=255)
                                                  aagaaCTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCAACTTGAGAAAGAGCGGTTAGCGGCt                    >  1:2080774/1‑71 (MQ=255)
                                                   agaaCTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCAACTTgaga                                     <  1:3569826/53‑1 (MQ=255)
                                                    gaaCTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCAACTTGAGAAAGAGCGGTTAGCGGCt                    >  1:1405850/1‑69 (MQ=255)
                                                     aaCTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCAACTTGAGAAAGAGCGGTTAGCGGCTCAg                 <  1:353953/71‑1 (MQ=255)
                                                     aaCTGCGTGAGAAACAAGAGGCTGAACAGGAACGCCTGAAGCAACTTGAGAAAGAGCGGTTAGCGGCTCAg                 <  1:881017/71‑1 (MQ=255)
                                                      aCTGCGTGAGAAACAATCGGCTGAACAGGAACGCCTGAAGCAACTTGAGAAAGAGCGGTTAGCGGCTCAg                 >  1:957627/1‑70 (MQ=255)
                                                      aCTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCAACTTGAGAAAGAGCGGTTAGCGGCTCAg                 >  1:204216/1‑70 (MQ=255)
                                                       cTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCaa                                            <  1:49310/42‑1 (MQ=255)
                                                       cTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCaa                                            <  1:1012336/42‑1 (MQ=255)
                                                       cTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCaa                                            <  1:2822444/42‑1 (MQ=255)
                                                       cTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCAACTTGAGAAAGAGCGGTTAGCgg                      <  1:3821479/64‑1 (MQ=255)
                                                       cTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCAACTTGAGAAAGAGCGGTTAGCgg                      <  1:582112/64‑1 (MQ=255)
                                                         gCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCAACt                                          >  1:3050519/1‑42 (MQ=255)
                                                         gCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCAACt                                          >  1:2200892/1‑42 (MQ=255)
                                                         gCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCAACTTGAGAAAGAGCGGTTAg                         >  1:749969/1‑59 (MQ=255)
                                                         gCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCAACTTGAGAAAGAGCGGTTAGCGGCTCAGGAg              >  1:2200819/1‑70 (MQ=255)
                                                          cGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCAACTTGAGAAAGAGCGGTTAGCGGCTCAGGAg              >  1:1447791/1‑69 (MQ=255)
                                                              agaAACAATCGGCTGAACAGGAACGCCTGAAGCAACTTGAGAAAGAGCGGTTAGCGGCTCAGGa               >  1:4097648/1‑64 (MQ=255)
                                                              agaAACAAGCGGCTGAACAGGAACGCCTGAAGCaa                                            >  1:2441089/1‑35 (MQ=255)
                                                              agaAACAAGCGGCTGAACAGGAACGCCTGAAGCAACTTga                                       >  1:1848999/1‑40 (MQ=255)
                                                               gaAACAATCGGCTGAACAGGAACGCCTGAAGCAACTTGAGAAAGAGCGGTTAGCGGCTc                   <  1:2666777/59‑1 (MQ=255)
                                                               gaAACAAGCGGCTGAACAGGAACGCCTGAAGCAACTTGAGAAAGAGCGGTTAGCGGCTCAGGAGCAGaaa        >  1:1997579/1‑70 (MQ=255)
                                                                aaaCAAGCGGCTGAACAGGAACGCCTGAAGCAACTTGAGAAAGAGCGGTTAGCg                       >  1:2879155/1‑54 (MQ=255)
                                                                aaaCAAGCGGCTGAACAGGAACGCCTGAAGCAACTTGAGAAAGAGCGGTTAGCg                       >  1:2849317/1‑54 (MQ=255)
                                                                aaaCAAGCGGCTGAACAGGAACGCCTGAAGCAACTTGAGAAAGAGCGGTTAGCGGCTCAGGAGCAGaaaaa      <  1:3132528/71‑1 (MQ=255)
                                                                aaaCAAGCGGCTGAACAGGAACGCCTGAAGCAACTTGAGAAAGAGCGGTTAGCGGCTCAGGAGCAGaaaaa      <  1:3127725/71‑1 (MQ=255)
                                                                aaaCAAGCGGCTGAACAGGAACGCCTGAAGCAACTTGAGAAAGAGCGGTTAGCGGCTCAGGAGCAGaaaaa      <  1:1339277/71‑1 (MQ=255)
                                                                 aaCAAGCGGCTGAACAGGAACGCCTGAAGCAACTTGAGAAAGAGCGGTTAGCGGCTCAGGAGCAGaaaaa      <  1:1876215/70‑1 (MQ=255)
                                                                    aaGCGGCTGAACAGGAACGCCTGAAGCAACTTGAGAAAGAGCGGTTAGCGGCTCAGGAGCAGAAAAAGCAg  >  1:173886/1‑71 (MQ=255)
                                                                      gCGGCTGAACAGGAACGCCTGAAGCAACTTGAGAAag                                  >  1:2462753/1‑37 (MQ=255)
                                                                      |                                                                    
CGCGAAGCGTTCTGATGAACAGCGCAAGATGAAGGAACAGCAGGCTGCTGAAGAACTGCGTGAGAAACAAGCGGCTGAACAGGAACGCCTGAAGCAACTTGAGAAAGAGCGGTTAGCGGCTCAGGAGCAGAAAAAGCAG  >  minE/483015‑483153

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: