#include <bits/stdc++.h>
using namespace std;
vector<int> E[200005];
const int mod = 1e9 + 7;
    const int N = 200005;
int n , q;
struct hld {
    
        int tag[N * 4];
        int sum[N * 4];
        void build(int u,int l,int r) {
            tag[u] = 1 ; sum[u] = (r - l + 1);
            if(l == r) return ;
            build(u<<1 , l , (l + r >> 1));
            build(u<<1|1 , (l + r >> 1) + 1 , r) ;
            return ;
        }
        void apply(int u,int v) {
            tag[u] = 1LL * tag[u] * v % mod;
            sum[u] = 1LL * sum[u] * v % mod;
        }
        void pd(int u) {
            if(tag[u] != 1) {
                apply(u<<1 , tag[u]);
                apply(u<<1|1 , tag[u]);
                tag[u] = 1;
            }
        }
        void modify(int u,int l,int r,int ql,int qr,int v) {
            if(ql <= l && qr >= r) {
                apply(u , v);
                return ;
            }
            int md = (l + r>>1);
            pd(u) ;
            if(ql <= md) modify(u<<1 , l , md , ql , qr , v);
            if(qr > md) modify(u<<1|1 , md + 1 , r , ql , qr , v);
            sum[u] = (sum[u << 1] + sum[u << 1|1]) % mod;
            return ;
        }
        
    int fa[N] , sz[N] , sn[N] , tp[N] , dfn[N] , nfd[N] , dep[N] , id;
    hld() {id = 0;}
    void dfs1(int u,int f) {
        sz[u] = 1; fa[u] = f;
        for(auto v : E[u]) {
            if(v != f) {
                dep[v] = dep[u] + 1;
                dfs1(v , u) ; sz[u] += sz[v] ;
                if(sz[sn[u]] < sz[v]) sn[u] = v;
            }
        }
    }
    void dfs2(int u,int p) {
        dfn[u] = ++id ; nfd[id] = u ; tp[u] = p;
        if(sn[u]) dfs2(sn[u] , p) ;
        for(auto v : E[u]) {
            if(v != fa[u] && v != sn[u]) dfs2(v , v);
        }
    }
    int lca(int x,int y) {
        while(tp[x] ^ tp[y]) {
            if(dep[tp[x]] < dep[tp[y]]) swap(x , y);
            x = fa[tp[x]] ;
        }
        if(dep[x] > dep[y]) swap(x , y) ;
        return x;
    }
    void mdf(int x,int y,int v) {
        while(tp[x] ^ tp[y]) {
            if(dep[tp[x]] < dep[tp[y]]) swap(x , y);
            modify(1 , 1 , n , dfn[x] , dfn[y] , v);
            x = fa[tp[x]] ;
        }
        if(dep[x] > dep[y]) swap(x , y);
        modify(1 , 1 , n , dfn[x] , dfn[y] , v) ;
    }
}A , B;

int main() {
    ios::sync_with_stdio(false) ; cin.tie(0) ; cout.tie(0) ;
    cin >> n;
    for(int i = 1;i < n;i++) {
        int u , v;cin >> u >> v;
        E[u].push_back(v) ; E[v].push_back(u) ;
    }
    A.dfs1(1 , 0) ; A.dfs2(1 , 1); A.build(1 , 1 , n);
    B.dfs1(1 , 0) ; B.dfs2(1 , 1); B.build(1 , 1 ,n);
    cin >> q;
    int r2 = (mod + 1) / 2;
    while(q--) {
        char ch ; int u , v;
        cin >> ch  >> u >> v;
        if(ch == '+') {
            A.mdf(u , v , 2);
            B.mdf(u , v , 2);
            int l = B.lca(u , v);
            B.mdf(l , l , r2) ;
        }
        else {
            A.mdf(u , v , r2);
            B.mdf(u , v, r2);
            int l = B.lca(u , v);
            B.mdf(l , l , 2) ;
        }
        int ans = (A.sum[1] - B.sum[1] + mod ) % mod;
        cout << ans << '\n';
    }
    return 0;
}