#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;
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]++;
    // cerr<<x<<" "<<sz1[x][0]<<" "<<sz2[x][0]<<" qoijrt\n";
    for (int v:e[x]) {
        if (v==fa) continue;
        DFS2(v,x);
        
        for (int A1=0;A1<=sz1[v][0];A1++) {
            for (int B1=0;B1<=sz2[v][0];B1++) {
                for (int A11=0;A11<=A[0]-sz1[v][0];A11++) {
                    for (int B11=0;B11<=B[0]-sz2[v][0];B11++) {
                        // 0 的个数是 A11+A1+A[1] =  B11+B1+B[1]
                        // cerr<<sz1[v][0]<<" "<<sz2[v][0]<<" \n";
                        if (A11+A1+A[1] ==  B11+B1+B[1]) {
                            res+=C(sz1[v][0],A1)*C(sz2[v][0],B1)%Mod
                                *C(A[0]-sz1[v][0],A11)%Mod*C(B[0]-sz2[v][0],B11)%Mod*abs(A1-B1+sz1[v][1]-sz2[v][1]);
                            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]++;
    }
    DFS2(1,0);
    cout<<res<<"\n";
}
signed main()
{
    ios::sync_with_stdio(false);
    int t=0;cin>>t;init();
    while (t--) solve();
    return 0;
}