#include <bits/stdc++.h>
#define int long long
#define io ios::sync_with_stdio(0),cin.tie(0),cout.tie(0)
#define pb push_back
using namespace std;
const int N=2e6+5,mod1=1e9+7,mod2=998244353;
typedef pair<int,int> pi;
int st[N];
vector<int>g[N];
int dfs(int u)
{
    for(auto v:g[u])
    {
        if(st[v]==st[u])return 0;
        if(st[v])continue;
        st[v]=3^st[u];
        if(!dfs(v))return 0;
    }
    return 1;
}
void solve()
{
    int n,m;
    cin>>n>>m;
    for(int i=1;i<=2*n;i++)g[i].clear(),st[i]=0;
    vector<int>p[n+5];
    for(int i=1;i<=n;i++)
    {
        string s;
        cin>>s;
        s=' '+s;
        for(int j=1;j<=m;j++)
            if(s[j]=='1')p[j].pb(i);
    }
    if(m&1)
    {
        if(p[m/2+1].size()>1) 
        {
            cout<<0<<'\n';
            return ;
        }
    }
    for(int i=1;i<=n;i++)
    {
        g[i].pb(i+n);
        g[i+n].pb(i);
    }
    for(int j=1;j<=m/2;j++)
    {
        int x=j,y=m-j+1;
        if(p[x].size()+p[y].size()>2)
        {
            cout<<0<<'\n';
            return ;
        }
        if(p[x].size()+p[y].size()==2&&p[x].size()==1)
        {
            int u=p[y][0],v=p[x][0];
            if(u==v)continue;
            g[u].pb(v+n);
            g[v].pb(u+n);
            g[u+n].pb(v);
            g[v+n].pb(u);
        }
        else if(p[x].size()==2)
        {
            int u=p[x][0],v=p[x][1];
            g[u].pb(v);
            g[v].pb(u);
            g[u+n].pb(v+n);
            g[v+n].pb(u+n);
        }
        else if(p[y].size()==2)
        {
            int u=p[y][0],v=p[y][1];
            g[u].pb(v);
            g[v].pb(u);
            g[u+n].pb(v+n);
            g[v+n].pb(u+n);
        }
    }
    int ans=1;
    for(int i=1;i<=2*n;i++)
    {
        if(st[i])continue;
        st[i]=1;
        if(!dfs(i))ans=0;
        else ans=ans*2%mod1;
    }
    cout<<ans<<'\n';
}
signed main()
{
    io;
    int t=1;
    cin>>t;
    while(t--)
    {
        solve();
    }
    return 0;
}