#include <bits/stdc++.h>
using namespace std;

typedef long long ll;
typedef pair<int, int> pii;
typedef pair<ll, ll> pll;

const int MAXN = 3e5;

const int MOD = 1e9+7;

struct mint
{
    int x;
    mint() : x(0) {}
    mint(int x) : x(x) {}
    mint operator + (int ot) const { return x+ot>=MOD ? x+ot-MOD : x+ot; }
    mint operator - (int ot) const { return x<ot ? x+MOD-ot : x-ot; }
    mint operator - () const { return x ? MOD-x : 0; }
    mint operator * (int ot) const { return 1ll*x*ot%MOD; }
    mint operator += (int ot) { return *this = *this + ot; }
    mint operator -= (int ot) { return *this = *this - ot; }
    mint operator *= (int ot) { return *this = *this * ot; }
    operator int() const { return x; }
};

int N;
vector<int> adj[MAXN+10];
int sz[MAXN+10];

mint dp1[MAXN+10];
mint dp2[MAXN+10][21];
mint dp3[MAXN+10][21];

int clz(int x)
{
    if(!x) return 32;
    return __builtin_clz(x);
}

void dfs(int now, int bef)
{
    sz[now]=1;
    vector<int> chd;
    for(int nxt : adj[now])
    {
        if(nxt==bef) continue;
        dfs(nxt, now);
        sz[now]+=sz[nxt];
        chd.push_back(nxt);
    }

    sort(chd.begin(), chd.end(), [&](const int &p, const int &q) { return sz[p]<sz[q]; });

    vector<mint> V[2];
    V[0].push_back(1);
    V[1].push_back(0);
    for(int i=0; i<chd.size(); i++)
    {
        int nxt=chd[i];
        int k=__lg(sz[nxt])+2;
        if(i+1==chd.size()) k=min(k, __lg((int)V[0].size())+2);
        k=min(k, 20);
        vector<mint> V2[2];
        V2[0]=V2[1]=vector<mint>(1<<k);
        for(int j=0; j<(1<<k); j++)
        {
            if(j<V[0].size()) V2[0][j]+=V[0][j]*dp1[nxt];
            for(int p=0; p<k; p++) if((j&(1<<p)) && (j^(1<<p))<V[0].size()) V2[0][j]+=V[0][j^(1<<p)]*dp2[nxt][p];
        }
        for(int j=0; j<(1<<k); j++)
        {
            if(j<V[0].size()) V2[1][j]+=V[1][j]*dp1[nxt];
            for(int p=0; p<k; p++) if((j&(1<<p)) && (j^(1<<p))<V[1].size() && clz(j)==clz(j^(1<<p))) V2[1][j]+=V[1][j^(1<<p)]*dp2[nxt][p];
            for(int p=0; p<k; p++)
            {
                if(j&(1<<p) && (j^(1<<p))<V[0].size() && clz(j)!=clz(j^(1<<p))) V2[1][j]+=V[0][j^(1<<p)]*dp3[nxt][p];
            }
        }

        V[0]=V2[0];
        V[1]=V2[1];
    }

    int k=min(__lg((int)sz[now])+1, 20);
    for(int p=0; p<=k+1 && (1<<p)-1<V[0].size(); p++)
    {
        dp1[now]+=V[0][(1<<p)-1];
        dp1[now]+=V[1][(1<<p)-1];
        
        dp2[now][p]+=V[0][(1<<p)-1];
        dp3[now][p]+=V[0][(1<<p)-1];
    }
    for(int p=0; p<=k; p++)
    {
        for(int q=p+2; q<=k+1; q++) if((((1<<q)-1)^(1<<p))<V[1].size()) dp3[now][p]+=V[1][((1<<q)-1)^(1<<p)]+V[0][((1<<q)-1)^(1<<p)];
    }

    // printf("NOW %d\n", now);
    // for(int i=0; i<V[0].size(); i++) printf("%d ", (int)V[0][i]);
    // printf(" : V0\n");
    // for(int i=0; i<V[1].size(); i++) printf("%d ", (int)V[1][i]);
    // printf(" : V1\n");
    // printf("%d : DP1\n", (int)dp1[now]);
    // for(int i=0; i<=k; i++) printf("%d ", (int)dp2[now][i]);
    // printf(" : DP2\n");
    // for(int i=0; i<=k; i++) printf("%d ", (int)dp3[now][i]);
    // printf(" : DP3\n\n");
}

int main()
{
    scanf("%d", &N);
    for(int i=1; i<N; i++)
    {
        int u, v;
        scanf("%d%d", &u, &v);
        adj[u].push_back(v);
        adj[v].push_back(u);
    }
    dfs(1, 1);
    printf("%d\n", (int)dp1[1]);
}