#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[2000005];
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;
    vector<string>q(n+1);
    for(int i=1;i<=2*n;i++)g[i].clear(),st[i]=0;
    for(int i=1;i<=n;i++)cin>>q[i];
    vector<int>b(m+5);
    vector<int>p[n+5];
    for(int i=1;i<=n;i++)
        for(int j=0;j<m;j++)
        {
            if(q[i][j]=='1')
            {
                b[j+1]++;
                if(b[j+1]>2)
                {
                    cout<<0<<'\n';
                    return ;
                }
                p[j+1].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(x==y)continue;
        if(p[x].size()==2)
        {
            int u=p[x][0],v=p[x][1];
            g[u].pb(v+n);
            g[v].pb(u+n);
            g[u+n].pb(v);
            g[v+n].pb(u);
        }
        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);
        }
        else if(p[x].size()+p[y].size()==2)
        {
            int u=p[y][0],v=p[x][0];
            g[u].pb(v+n);
            g[v].pb(u+n);
            g[u+n].pb(v);
            g[v+n].pb(u);
        }
    }
    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;
}