# 初期設定
import sys
sys.setrecursionlimit(2000)
input = sys.stdin.readline
mod = 10**9+7

# 入力 + Cを圧縮
n, c = map(int, input().split())
minc = min(n * 6 - 5, c)
edge = [[] for _ in range(n)]
for _ in range(n - 1):
  a, b = map(int, input().split())
  a -= 1
  b -= 1
  edge[a].append(b)
  edge[b].append(a)

# dp初期化(初期値は1です)
dp = [minc * [1] for _ in range(n)]

# 再帰で実装します
# 返り値はdp_v
def dfs(r, v):
  # vの親がr、vの子がu
  for u in edge[v]:
    if r == u:
      continue
    # dp_uを取得
    dp_u = dfs(v, u)
    for i in range(minc):
      # dp_{u,i-3} + dp_{u,i+3}を求める
      f = 0
      if i > 2:
        f += dp_u[i - 3]
      if i < minc - 3:
        f += dp_u[i + 3]
      
      # かける
      # p_v=iにすることが不可能な場合、f=0になる子が存在するのでこのままでよいです
      dp[v][i] *= f
      dp[v][i] %= mod
  return dp[v]

# 頂点1を根としてDFS、dp_1を取得
rootdp = dfs(-1, 0)

# rootdpの総和に加え、Cを圧縮しただけansに加算
ans = sum(rootdp)
ans += rootdp[minc // 2] * (c - minc)
print(ans % mod)