#include <bits/stdc++.h>
#define int long long 
#define pb push_back
using namespace std;
const int Mod=1e9+7;
const int INF=1e6+5;
int n,col[INF],fav[INF],ifav[INF];string s,t;
vector <int> e[INF];
int ksm(int x,int y) {
    int ba=x%Mod,ans=1;
    while (y) {
        if (y&1) ans=(ans*ba)%Mod;
        ba=(ba*ba)%Mod;y>>=1;
    }
    return ans;
}
void init() {
    fav[0]=1;int N=1e6;
    for (int i=1;i<=N;i++) fav[i]=fav[i-1]*i%Mod;
    ifav[N]=ksm(fav[N],Mod-2);
    for (int i=N-1;~i;i--) ifav[i]=ifav[i+1]*(i+1)%Mod;
}
int C(int x,int y) {
    if (x<y) return 0;
    return fav[x]*ifav[y]%Mod*ifav[x-y]%Mod;
}
void DFS(int x,int fa) {
    for (int v:e[x]) {
        if (v==fa) continue;
        col[v]=col[x]^1;
        DFS(v,x);
    }
}
int A[5],B[5];
int sz1[INF][5],sz2[INF][5],res;
int X,Y,Z;
void DFS2(int x,int fa) {
    for (int w=0;w<3;w++) sz1[x][w]=sz2[x][w]=0;
    if (s[x]=='?') sz1[x][0]++;
    else if (s[x]=='0') sz1[x][1]++;
    else sz1[x][2]++;

    if (t[x]=='?') sz2[x][0]++;
    else if (t[x]=='0') sz2[x][1]++;
    else sz2[x][2]++;
    for (int v:e[x]) {
        if (v==fa) continue;
        DFS2(v,x);
        // W<=B[1]-A[1]+T2
        // B[1]-A[1]-T1
        if (Y>=0) {
            res-=(sz1[v][0]+sz1[v][1]-sz2[v][1])*C(Z,Y);
            res%=Mod;res+=Mod;res%=Mod;
            res+=(sz1[v][0]+sz2[v][0])*C(Z-1,Y-1)%Mod;res%=Mod;
            // for (int W=0;W<Y;W++) {
            //     res+=C(sz1[v][0]+sz2[v][0]-1,W)%Mod
            //         *C(Z-sz2[v][0]-sz1[v][0],Y-1-W)%Mod*(sz1[v][0]+sz2[v][0]);
            //     res%=Mod;
            // }
        }
//        v3[v].pb((sz1[v][0]+sz1[v][1]-sz2[v][1])*2);
//        v3[v].pb(-(sz1[v][0]+sz2[v][0])*2);
        
        for (int W=0;W<=min(Y,sz1[v][0]+sz1[v][1]-sz2[v][1]);W++) {
            res+=C(sz1[v][0]+sz2[v][0],W)%Mod
                *C(Z-sz2[v][0]-sz1[v][0],Y-W)%Mod*(sz1[v][0]+sz1[v][1]-sz2[v][1])*2;
            res%=Mod;
        }
        for (int W=0;W<=min(Y,sz1[v][0]+sz1[v][1]-sz2[v][1]);W++) {
            res+=C(sz1[v][0]+sz2[v][0]-1,W-1)%Mod
                *C(Z-sz2[v][0]-sz1[v][0],Y-W)%Mod*(sz1[v][0]+sz2[v][0])*-2;
            res%=Mod;
        }

        for (int i=0;i<3;i++)
            sz1[x][i]+=sz1[v][i],
            sz2[x][i]+=sz2[v][i];
    }
}
void solve() {
    for (int w=0;w<=n+3;w++) e[w].clear();
    cin>>n;res=0;
    for (int i=1;i<n;i++) {
        int x=0,y=0;
        cin>>x>>y;
        e[x].pb(y);
        e[y].pb(x);
    }
    for (int w=0;w<3;w++) A[w]=B[w]=0;
    DFS(1,0);cin>>s>>t;s=" "+s;t=" "+t;
    for (int i=1;i<=n;i++) {
        if (s[i]=='?') {A[0]++;continue;}
        if (col[i]==1) {
            if (s[i]=='0') s[i]='1';
            else s[i]='0';
        }
        if (s[i]=='0') A[1]++;
        else A[2]++;
    }
    for (int i=1;i<=n;i++) {
        if (t[i]=='?') {B[0]++;continue;}
        if (col[i]==1) {
            if (t[i]=='0') t[i]='1';
            else t[i]='0';
        }
        if (t[i]=='0') B[1]++;
        else B[2]++;
    }
    X=-B[1]+A[1]-B[0];
    Y=-B[1]+A[1]+A[0];
    Z=A[0]+B[0];
    DFS2(1,0);
    cout<<res<<"\n";
}
signed main()
{
    ios::sync_with_stdio(false);
    int t=0;cin>>t;init();
    while (t--) solve();
    return 0;
}