#include <bits/stdc++.h>
#include <boost/multiprecision/cpp_int.hpp>
using namespace std;
namespace mp = boost::multiprecision;
int main(){
    int T;
    cin>>T;
    for(int j=0;j<T;++j){
        int A,B,C,D;
        cin>>A>>B>>C>>D;
        if(A!=B){
            cout<<0<<endl;
            continue;
        }
        int mod=1e9+7;
        long long a=1,b=A;
        for(int i=1;i<A;++i){
            a=(a*(D+i))%mod;
            b=(b*i)%mod;
            a=(a*(A+D+1+i))%mod;
            b=(b*i)%mod;
        }
        for(int i=1;i<=C;++i){
            a=(a*(A+B+D+i))%mod;
            b=(b*i)%mod;
        }
        int counter=mod-2;
        int ans=1;
        while(counter){
            if(counter%2){
                ans=(ans*b)%mod;
            }
            b=(b*b)%mod;
            counter/=2;
        }
        ans=(ans*a)%mod;
        cout<<ans<<endl;
    }
}