#include<bits/stdc++.h>
#define int long long
#define X first
#define Y second
#define lson (root<<1)
#define rson (root<<1|1)
#define ALL(x) x.begin(),x.end()
#define CLR(x,y) memset(x,y,sizeof(x))
//#define endl '\n'
using namespace std;
typedef pair<int,int> pii;
int n;
#define MAX_LOG 10
int pa[500][MAX_LOG+1][500],dep[500][500];
vector<int>edge[500];
void dfs(int root,int x,int p)
{
	pa[root][0][x]=p;
	for(int i=0;i+1<11;++i){
		pa[root][i+1][x]=pa[root][i][pa[root][i][x]];
	}
	for(int v:edge[x])
	{
		if(v==p)continue;
		dep[root][v]=dep[root][x]+1;
		dfs(root,v,x);
	}
}
inline int find_lca(int root,int a,int b){
	if(dep[root][a]>dep[root][b])swap(a,b);
	for(int i=dep[root][b]-dep[root][a],k=0;i;i/=2){
		if(i&1)b=pa[root][k][b];
		++k;
	}
	if(a==b)return a;
	for(int i=MAX_LOG;i>=0;--i){
		if(pa[root][i][a]!=pa[root][i][b]){
			a=pa[root][i][a];
			b=pa[root][i][b];
		}
	}
	return pa[root][0][a];
}
int MOD=1000000007;
int qp(int a,int b)
{
	int ans=1;
	while(b)
	{
		if(b&1)ans=ans*a%MOD;
		a=a*a%MOD;
		b>>=1;
	}
	return ans%MOD;
}
int inv(int x)
{
	return qp(x,MOD-2);
}
int dp[500][500];
int dfs2(int x,int y)
{
	if(~dp[x][y])return dp[x][y];
	if(x==0)return dp[x][y]=1;
	else if(y==0)return dp[x][y]=0;
	else return dp[x][y]=inv(2)*(dfs2(x-1,y)+dfs2(x,y-1))%MOD;
}
void solve()
{
	CLR(dp,-1);
	cin>>n;
	for(int i=0;i<n-1;i++)
	{
		int u,v;
		cin>>u>>v;
		edge[u].emplace_back(v);
		edge[v].emplace_back(u);
	}
	for(int i=1;i<=n;i++)dfs(i,i,-1);
	int E=0;
	int invn=inv(n);
	for(int i=1;i<=n;i++)
	{
		for(int j=i+1;j<=n;j++)
		{
			//j in front of i
			for(int root=1;root<=n;root++)
			{
				int lca=find_lca(root,i,j);
				if(lca==j)E+=invn;
				else if(lca==i)E+=0;
				else
				{
					int a=dep[root][j]-dep[root][lca],b=dep[root][i]-dep[root][lca];
					E+=(dfs2(a,b))*invn%MOD;
				}
				E%=MOD;
			}
		}
	}
	cout<<E<<endl;
}
signed main()
{
	ios::sync_with_stdio(0);
	cin.tie(0);
	solve();
	return 0;
}
