#include<bits/stdc++.h>
#define ll long long  
#define pb push_back
using namespace std;
const int N=6e5+10;
const ll mod=1e9+7;
vector<ll>f[N],g[N];
ll s[N];
int n,h[N],l[N],t[N],ct,m;
struct node{
    ll x,y;
}bb[N];
bool cmp(node a,node b){return a.y<b.y;}
void add(int u,int v){t[++ct]=v,l[ct]=h[u],h[u]=ct;}
void dfs(int x,int fa){
    for(int i=h[x];i;i=l[i]) if(t[i]!=fa) dfs(t[i],x);
    m=0;
    for(int i=h[x];i;i=l[i]) if(t[i]!=fa) bb[++m]={t[i],(ll)f[t[i]].size()};
    sort(bb+1,bb+m+1,cmp);
    vector<ll>a,b;
    a.pb(1),b.pb(0);
    for(ll i=1,j,k;i<=m;i++){
        j=bb[i].x,k=f[j].size();
        vector<ll> c(1<<k),d(1<<k);
        for(ll o=0;o<a.size();o++){
            (c[o]+=a[o]*s[j]%mod)%=mod;
            (d[o]+=b[o]*s[j]%mod)%=mod;
            for(ll l=0;l<k;l++)if(!((o>>l)&1)){
                (c[o|(1<<l)]+=a[o]*f[j][l]%mod)%=mod;
                if(o>(1<<l)) (d[o|(1<<l)]+=b[o]*f[j][l]%mod)%=mod;
                else (d[o|(1<<l)]+=a[o]*g[j][l]%mod)%=mod;
            }
        }
        a.swap(c),b.swap(d);
    }
    ll k=f[bb[m].x].size()+1;
    f[x].resize(k),g[x].resize(k);
    for(ll i=0,j;i<k;i++){
        j=(1<<i)-1;
        (f[x][i]+=a[j])%=mod,(g[x][i]+=a[j])%=mod;
        for(ll l=0;l<=i-2;l++) (g[x][i]+=a[j^(1<<l)])%=mod;
        for(ll l=0;l<=i-2;l++) (g[x][i]+=b[j^(1<<l)])%=mod;
        (s[x]+=a[j]+b[j])%=mod;
    }
    while(f[x].back()==0&&g[x].back()==0) f[x].pop_back(),g[x].pop_back();
    
}
int main(){
    scanf("%d",&n);
    for(int i=1;i<n;i++){
        int u,v;
        scanf("%d%d",&u,&v);
        add(u,v),add(v,u);
    }
    dfs(1,0);
    printf("%lld",s[1]);
}