#include <bits/stdc++.h>
using namespace std;
#define F first
#define S second
#define pb push_back
#define ll long long
#define Mid ((L+R)/2)
#define pi pair<int,int>
#define pii pair<pi,int>
#define deb(x) cout<<#x<<"="<<x<<endl
#define go ios::sync_with_stdio(0); cin.tie(0); cout.tie(0)
const int N = 202;
const int mod=1e9+7;
int n,ans,dp[N][N],a[N][N];
int mul(int x,int y)
{
    return (ll) x*y %mod;
}
int sum(int x,int y)
{
    if((x+=y)>=mod)
        x-=mod;
    return x;
}
int po(int x,int y)
{
    if(!y)  return 1;
    if(y&1) return mul(x,po(x,y-1));
    int z=po(x,y/2);
    return mul(z,z);
}
int inv(int x)
{
    return po(x,mod-2);
}

int main()
{
    go;
    cin>>n;
    for(int i=1; i<=n; i++)
        for(int j=1; j<=n; j++)
            a[i][j]=(i==j)?0:1e9;
    for(int i=1,x,y; i<n; i++)
    {
        cin>>x>>y;
        a[x][y]=a[y][x]=1;
    }
    for(int k=1; k<=n; k++)
        for(int i=1; i<=n; i++)
            for(int j=1; j<=n; j++)
                a[i][j]=min(a[i][j],a[i][k]+a[k][j]);
    for(int i=0; i<=n; i++)
        dp[0][i]=1;
    for(int i=1; i<=n; i++)
        for(int j=1; j<=n; j++)
            dp[i][j]=mul(inv(2),sum(dp[i][j-1],dp[i-1][j]));
    for(int root=1; root<=n; root++)
        for(int i=1; i<=n; i++)
            for(int j=i+1; j<=n; j++)
            {
                int ri=a[root][i],rj=a[root][j],ij=a[i][j];
                int com=(ri+rj-ij)/2;
                ri-=com;
                rj-=com;
                ans=sum(ans,dp[rj][ri]);
            }
    cout<<mul(inv(n),ans);
    return 0;
}
