#include <bits/stdc++.h>

using namespace std;

#define rep(i, a, n) for (int i=(a); i < (n); i++)
#define per(i, a, n) for (int i=(n)-1; i >= (a); i--)
#define ll long long
ll N, K;
ll res;
ll MOD;
vector<vector<ll> > pod;
vector<vector<ll> > v;
vector<vector<ll> > solve(ll a){
	vector<vector<ll> >puvodni;
	vector<vector<ll> > novy;
	puvodni.resize(N+1); novy.resize(N+1);
	ll j,k,l;
	for(j=0;j<=N;j++){
		puvodni[j]={0,0,0};
		novy[j]={0,0,0};
	}
	puvodni[0]={1,0,0};
	puvodni[1]={0,0,1};
	if(pod[a].size()==0) return puvodni;
	for(l=0;l<pod[a].size();l++){
		v=solve(pod[a][l]);
		for(j=0;j<=K;j++){
			for(k=0;k<3;k++){
				novy[j][k]=0;
			}
		}
		for(j=0;j<=K;j++){
			for(k=0;k<=K-j;k++){
				novy[j+k][0]+=((v[j][0]+v[j][1])*puvodni[k][0]);
				novy[j+k][0]%=MOD;
				novy[j+k][2]+=(v[j][0]*puvodni[k][2]);
				novy[j+k][2]%=MOD;
				novy[j+k][1]+=((v[j][0]+v[j][1]+v[j][2])*(puvodni[k][1]+puvodni[k][2]));
				novy[j+k][1]-=(v[j][0]*puvodni[k][2]);
				novy[j+k][1]%=MOD;
			}
		}
		/*for(j=0;j<=K;j++){
			for(k=0;k<3;k++){
				printf("%lld ", novy[j][k]);
			}
			printf("\n");
		}
		printf("\n");*/
		puvodni=novy;
	}
	/*for(j=0;j<=K;j++){
		puvodni[j][1]-=puvodni[j][2];
		puvodni[j][1]+=MOD;
		puvodni[j][1]%=MOD;
	}*/
	return puvodni;
}
int main(void){
	MOD=1000000007;
	while(scanf("%lld %lld", &N, &K)==2){
		vector<vector<ll> > sousede;
		ll x, y;
		sousede.resize(N);
		ll flag[N];
		pod.clear(); pod.resize(N);
		int i;
		for(i=0;i<N-1;i++){
			scanf("%lld %lld", &x, &y);
			x--; y--;
			sousede[x].push_back(y);
			sousede[y].push_back(x);
			flag[i]=0;
		}	
		flag[N-1]=0;
		queue<ll> fronta;
		fronta.push(0);
		while(!fronta.empty()){
			x=fronta.front();
			fronta.pop();
			flag[x]=1;
			for(i=0;i<sousede[x].size();i++){
				if(flag[sousede[x][i]]==0){
					pod[x].push_back(sousede[x][i]);
					flag[sousede[x][i]]=1;
					fronta.push(sousede[x][i]);
				}
			}
		}
		vector<vector<ll> > w;
		w=solve(0);
		printf("%lld\n", (w[K][0]+w[K][1])%MOD);
	}
	
	return 0;
}