#include<bits/stdc++.h>
using namespace std;
constexpr int maxn=1e6+5;
string s[maxn];
int fa[maxn<<1];
inline int getfa(int x){return fa[x]==x?x:fa[x]=getfa(fa[x]);}
inline void con(int u,int v){
    // cout<<u<<' '<<v<<'\n';
    u=getfa(u),v=getfa(v);
    if(u==v) return;
    fa[v]=u;
}
inline void solve(){
    int n,m;
    cin>>n>>m;
    for(int i=1;i<=n;i++) cin>>s[i],s[i]=' '+s[i];
    for(int i=1;i<=2*n;i++) fa[i]=i;
    if(m&1){
        int cnt=0;
        for(int i=1;i<=n;i++) cnt+=(s[i][(m+1)/2]-'0');
        if(cnt>1){
            cout<<"0\n";
            return;
        }
    }
    for(int j=1;j<m-j+1;j++){
        vector<pair<int,int>> v;
        for(int i=1;i<=n;i++){
            if(s[i][j]=='1') v.emplace_back(i,0);
            if(s[i][m-j+1]=='1') v.emplace_back(i,1);
        }
        if(v.size()>2) return cout<<"0\n",void();
        if(v.size()<2||v[0].first==v[1].first) continue;
        auto tr=[&](pair<int,int> p,int rev){
            if(rev) p.second^=1;
            return p.first*2-p.second;
        };
        con(tr(v[0],0),tr(v[1],1));
        con(tr(v[0],1),tr(v[1],0));
    }
    int cnt=0;
    for(int i=1;i<=2*n;i++) if(fa[i]==i) cnt++;
    cnt/=2;
    int ans=1;
    while(cnt--) ans=ans*2%1000000007;
    cout<<ans<<'\n';
}
int main(){
    std::ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);
    int t;cin>>t;
    while(t--) solve();
}