MOD = 10 ** 9 + 7 n = int(input()) g = [[] for _ in range(n)] deg = [0 for _ in range(n)] for _ in range(n - 1): a, b = map(int, input().split()) a -= 1 b -= 1 g[a].append(b) deg[b] += 1 src = None for v in range(n): if deg[v] == 0: src = v break dist = [None for _ in range(n)] dist[src] = 0 par = [None for _ in range(n)] dp = [0 for _ in range(n)] ans = 0 st = [~ src, src] while len(st) > 0: cur = st.pop() if cur >= 0: for nxt in g[cur]: dist[nxt] = dist[cur] + 1 par[nxt] = cur st.append(~ nxt) st.append(nxt) else: dp[~ cur] += 1 if not par[~ cur] is None: dp[par[~ cur]] += dp[~ cur] ans = (ans + dp[~ cur] * dist[~ cur]) % MOD print(ans)