#include <bits/stdc++.h>
#define IOS ios::sync_with_stdio(false), cin.tie(0), cout.tie(0)
#define rep(i, l, r) for (int i = l; i <= r; i++)
#define frep(i, l, r) for (int i = l; i >= r; i--)
#define eb emplace_back
#define pb push_back
#define endl "\n"
#define int long long
using namespace std;
using ll = long long;
using PII = pair<int, int>;
const int N = 1e6 + 3;
const int mod = 1e9 + 7;
const int inf = 0x3f3f3f3f;
int n, m;
int vis[N], w[N];
bool tag;
vector<vector<PII>> ve;
void dfs(int x, int p)
{
    if (!vis[x])
    {
        vis[x] = 1;
        w[x] = p;
        for (auto k : ve[x])
        {
            if (k.second == 0)
                dfs(k.first, p);
            if (k.second == 1)
                dfs(k.first, p ^ 1);
        }
    }
    else
    {
        if (w[x] != p)
        {
            tag = false;
            return;
        }
    }
}
void solve()
{
    cin >> n >> m;
    ve.assign(n + 10, vector<PII>());
    string s[n + 10];
    rep(i, 1, n)
    {
        cin >> s[i];
        s[i] = '!' + s[i];
        vis[i] = 0;
        w[i] = -1;
    }
    vector<int> a(m + 10, 0);
    rep(i, 1, n)
    {
        rep(j, 1, (m + 1) / 2)
        {
            if (s[i][j] == '1')
                a[j]++;
            if (s[i][m - j + 1] == '1' && j != m - j + 1)
                a[j]++;
            if (a[j] > 2)
                return cout << "0" << endl, void();
        }
    }
    if (m % 2 == 1 && a[(m + 1) / 2] >= 2)
        return cout << "0" << endl, void();
    vector<PII> b(m + 10, {0, 0});
    vector<int> vis1(m + 10, 0);
    rep(i, 1, n)
    {
        rep(j, 1, (m + 1) / 2)
        {
            if (s[i][j] == '1')
            {
                if (!vis1[j])
                {
                    b[j] = {i, 0};
                    vis1[j] = 1;
                }
                else
                {
                    ve[b[j].first].push_back({i, b[j].second});
                    ve[i].push_back({b[j].first, b[j].second});
                    //cerr << b[j].first << "-" << i << "-" << (b[j].second) << endl;
                }
            }
            if (s[i][m - j + 1] == '1')
            {
                if (s[i][j] == '1')
                    continue;
                if (!vis1[j])
                {
                    b[j] = {i, 1};
                    vis1[j] = 1;
                }
                else
                {
                    ve[b[j].first].push_back({i, b[j].second ^ 1});
                    ve[i].push_back({b[j].first, b[j].second ^ 1});
                    //cerr << b[j].first << " " << i << " " << (b[j].second ^ 1) << endl;
                }
            }
        }
    }
    //cerr << "--" << endl;
    tag = true;
    int ans = 1;
    rep(i, 1, n)
    {
        if (!vis[i])
        {
            dfs(i, 0);
            ans = ans * 2 % mod;
        }
    }
    if (!tag)
        return cout << "0" << endl, void();
    cout << ans << endl;
}

signed main()
{
    IOS;
    int T = 1;
    cin >> T;
    while (T--)
        solve();
}