#include<bits/stdc++.h>
using namespace std;
using ll=long long;
using i64=long long;
using u64=unsigned long long;
template<class T>
struct DSU{
    T n;
    vector<T> fa;
    vector<T> siz;
    DSU(int x):n(x),fa(x+1),siz(x+1){
        for(int i=0;i<=x;i++){
            fa[i]=i;
            siz[i]=1;
        }
    }
    T find(T x){
        if(x==fa[x]) return x;
        else{
            fa[x]=find(fa[x]);
            return fa[x];
        }
    }
    bool merge(T x,T y){
        x=find(x);
        y=find(y);
        if(x!=y){
            fa[y]=x;
            siz[x]+=siz[y];
            return 1;
        }
        return 0;
    }
    T ask_size(T x){
        x=find(x);
        return siz[x];
    }
};
const ll mod = 1e9+7;
i64 power(i64 a, i64 b, i64 m)
{
    i64 res = 1;
    a %= m;
    while (b > 0)
    {
        if (b & 1)
            res = (res % m * a % m) % m;
        a = (a % m * a % m) % m;
        b >>= 1;
    }
    return res;
}
void solve(){
    int n,m;
    cin>>n>>m;
    vector<string> s(n+1);
    for(int i=1;i<=n;i++) cin>>s[i],s[i]=' '+s[i];
    vector<pair<int,int>> ans;
    for(int i=1;i<=m/2;i++){
        int sum=0;
        int l=0;
        int r=0;
        for(int j=1;j<=n;j++){
            if(s[j][i]=='1'||s[j][m-i+1]=='1'){
                if(s[j][i]=='1'){
                    sum++;
                    if(!l) l=j;
                    else r=j;
                }
                if(s[j][m-i+1]=='1'){
                    sum++;
                    if(!l) l=j;
                    else r=j;
                }
            }
        }
        if(l&&r) ans.push_back({l,r});
        else if(l) ans.push_back({l,l});
        else if(r) ans.push_back({r,r});
        if(sum>2){
            cout<<"0"<<'\n';
            return;
        }
    }
    if(m%2){
        int sum=0;
        for(int i=1;i<=n;i++){
            if(s[i][(m+1)/2]=='1') sum++;
        }
        if(sum>1){
            cout<<"0"<<'\n';
            return;
        }
    }
    DSU<int> dsu(2*n+1);
    map<int,int> mp;
    for(auto [x,y]:ans){
        if(x==y){
            dsu.merge(x,y);
            dsu.merge(x+n,y+n);
        }
        else{
            dsu.merge(x+n,y);
            dsu.merge(x,y+n);
        }
    }
    for(int i=1;i<=n;i++){
        if(dsu.find(i)==dsu.find(i+n)){
            cout<<"0"<<"\n";
            return;
        }
    }
    ll sum=0;
    for(int i=1;i<=2*n;i++){
        if(dsu.find(i)==i) sum++;
    }
    sum/=2;
    cout<<power(2,sum,mod)<<"\n";
    
}
int main(){
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    int t;
    cin>>t;
    while(t--) solve();
    return 0;
}