#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int, int> pii;
#define pb push_back
#define all(a) a.begin(), a.end()
#define sz(a) ((int)a.size())
#ifdef ABS
void db() {}
template <typename T>
ostream& operator << (ostream &o, vector <T> vec) {
    o << "{"; int f = 0;
    for (T i : vec) o << (f++ ? " " : "") << i;
    return o << "}";
}
template <typename T, typename ...U> void db(T i, U ...j) { cerr << i << " ", db(j...); }
#define ppr(c, x...) cerr << "\e[1;" << c << "m", db(__LINE__, "[" + string(#x) + "]", x), cerr << "\e[0m" << endl
#define bug(x...) ppr(32, x)
#define bugv(x...) ppr(36, vector(x))
#define safe ppr(33, "safe")
#else
#define bug(x...) void(0)
#define bugv(x...) void(0)
#define safe void(0)
#endif
const int mod = 998244353, N = 300007, M = 205;

int add(int x, int y){
    x+=y; if(x>=mod) x-=mod;
    return x;
}
int sub(int x, int y){
    x-=y; if(x<0) x+=mod;
    return x;
}
int mul(int x, int y){
    return ((ll)x)*y%mod;
}
int Pow(int x, ll y=mod-2){
    int res=1;
    for(; y; x=mul(x,x),y>>=1) if(y&1) res=mul(res,x);
    return res;
}

int fac[M],inv[M],ifac[M];

int C(int n, int m){
    if(m<0||m>n) return 0;
    return mul(fac[n],mul(ifac[m],ifac[n-m]));
}

int m;

vector<int> add(vector<int> a, vector<int> b){
    vector<int> res(max(sz(a),sz(b)));
    for(int i=0; i<sz(a); ++i) res[i]=a[i];
    for(int i=0; i<sz(b); ++i) res[i]=add(res[i],b[i]);
    return res;
}
vector<int> sub(vector<int> a, vector<int> b){
    vector<int> res(max(sz(a),sz(b)));
    for(int i=0; i<sz(a); ++i) res[i]=a[i];
    for(int i=0; i<sz(b); ++i) res[i]=sub(res[i],b[i]);
    return res;
}
vector<int> mul(vector<int> a, vector<int> b){
    vector<int> res(min(sz(a)+sz(b)-1,m+1));
    for(int i=0; i<sz(a); ++i) for(int j=0; j<sz(b); ++j){
        if(i+j<sz(res)){
            res[i+j]=add(res[i+j],mul(a[i],b[j]));
        }
    }
    return res;
}

int n,pw[M],val[M];
vector<int> adj[N];
vector<int> dp[N][2],res;

void dfs(int u, int fa){
    dp[u][0]={1},dp[u][1]={0};
    for(auto v: adj[u]) if(v!=fa){
        dfs(v,u);
        vector<int> nw[2];
        nw[0]={0},nw[1]={0};
        nw[1]=add(nw[1],mul(dp[u][1],add(dp[v][0],dp[v][1])));
        nw[1]=add(nw[1],mul(mul(dp[u][0],vector<int>({1,1})),dp[v][0]));
        nw[0]=add(nw[0],mul(dp[u][0],dp[v][1]));
        //nw[0]=add(nw[0],dp[v][1]);
        //nw[1]=add(nw[1],mul(vector<int>({1,1}),dp[v][0]));
        dp[u][0]=add(dp[u][0],nw[0]);
        dp[u][1]=add(dp[u][1],nw[1]);
        //bug(u,v,dp[u][0],dp[u][1]);
    }
    for(int t: {0,1}){
        for(int i=0; i<sz(dp[u][t]); ++i){
            dp[u][t][i]=mul(dp[u][t][i],2);
        }
    }
    dp[u][0][0]=sub(dp[u][0][0],1);
    res=add(res,dp[u][0]);
    res=add(res,dp[u][1]);
    for(auto v: adj[u]) if(v!=fa){
        res=sub(res,mul(vector<int>({1,1}),dp[v][0]));
        res=sub(res,dp[v][1]);
    }
    //bug(u,dp[u][0],dp[u][1]);
}

void solve() {
    cin >> n >> m;
    for(int i=0; i<n; ++i) adj[i].clear();
    for(int i=0; i<n-1; ++i){
        int u,v; cin >> u >> v; u--,v--;
        adj[u].pb(v),adj[v].pb(u);
    }
    for(int i=0; i<=m; ++i){
        pw[i]=Pow(i,m);
    }
    for(int i=0; i<=m; ++i){
        val[i]=pw[i];
        for(int j=0; j<i; ++j) val[i]=sub(val[i],mul(val[j],C(i,j)));
    }
    res={0};
    dfs(0,-1);
    int tot=0;
    for(int i=0; i<sz(res); ++i) tot=add(tot,mul(res[i],val[i]));
    cout << tot << "\n";
}

int main() {
    ios::sync_with_stdio(false), cin.tie(0);
    fac[0]=inv[1]=ifac[0]=1;
    for(int i=1; i<M; ++i) fac[i]=mul(fac[i-1],i);
    for(int i=2; i<M; ++i) inv[i]=mul(inv[mod%i],mod-mod/i);
    for(int i=1; i<M; ++i) ifac[i]=mul(ifac[i-1],inv[i]);
    int t; cin >> t;
    while(t--) solve();
}