#include <bits/stdc++.h>
using namespace std;
#define int long long
mt19937_64 rng(chrono::steady_clock::now().time_since_epoch().count());
const int MOD=1000000007;
void Solve() {
    int r,c;
    cin>>r>>c;
    string mp[r];
    for(int i=0;i<r;i++)cin>>mp[i];
    vector<vector<int>> adj(r+1);
    vector<bool> vis(r+1);
    for(int i=0;i<=(c-1)/2;i++) {
        vector<int> v;
        for(int j=0;j<r;j++) {
            if(mp[j][i]=='1') {
                v.push_back(j);
            }
            if(mp[j][c-i-1]=='1') {
                v.push_back(j);
            }
            if(v.size()>2) {cout<<0<<endl;return;}
        }
        if(v.size()==2&&v[0]!=v[1]) {
            adj[v[0]].emplace_back(v[1]);
            adj[v[1]].emplace_back(v[0]);
        }
    }
    function<void(int u)> dfs=[&](int u){
        vis[u]=1;
        for(auto v:adj[u]) {
            if(!vis[v]) dfs(v);
        }
    };
    int ans=1;
    for(int i=0;i<r;i++) {
        if(!vis[i]) {
            dfs(i);
            ans=ans*2%MOD;
        }
    }
    cout<<ans<<endl;
}

signed main() {
    int Task;
    cin >> Task;
    while (Task--) Solve();
}