#include<bits/stdc++.h>
using namespace std;
const int N=1000005,P=1e9+7;
int r,c,col[N];
vector<int> G[N];
int dfs(int u){
    for(auto v:G[u]){
        if(!col[v]) col[v]=3-col[u],dfs(v);
        else if(col[v]+col[u]!=3) return 1;
    }
    return 0;
}
void solve(){
    scanf("%d%d",&r,&c);
    string b[r];
    for(int i=0;i<r;i++) cin>>b[i];
    int ct=0;
    if(c%2==1){
        for(int i=0;i<r;i++) ct+=(b[i][(c+1)/2-1]=='1');
        if(ct>1) return puts("0"),void();
    }
    for(int i=0;i<r;i++) G[i].clear();
    for(int i=0;i<(c/2);i++){
        vector<int> v;
        for(int j=0;j<r;j++){
            if(b[j][i]=='1') v.push_back(j);
            if(b[j][c-1-i]=='1') v.push_back(j);
            if(v.size()>2) return puts("0"),void();
        }
        if(v.size()==2&&v[0]!=v[1]) G[v[0]].push_back(v[1]),G[v[1]].push_back(v[0]); 
    }
    for(int i=0;i<r;i++) col[i]=0;
    int as=1;
    for(int i=0;i<r;i++) if(!col[i]){
        col[i]=1,as=2ll*as%P;
        if(dfs(i)) return puts("0"),void();
    }
    printf("%d\n",as);
}
int main(){
    int T;
    scanf("%d",&T);
    while(T--) solve();
}