#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<ll, ll> pii;
const int mod = 1e9 + 7;
bool chmin(int &x, int y) {
	if (x > y) {
		x = y;
		return 1;
	}
	return 0;
}
bool chmax(int &x, int y) {
	if (x < y) {
		x = y;
		return 1;
	}
	return 0;
}


struct DSU {
    std::vector<int> f, siz;
    
    DSU() {}
    DSU(int n) {
        init(n);
    }
    
    void init(int n) {
        f.resize(n);
        std::iota(f.begin(), f.end(), 0);
        siz.assign(n, 1);
    }
    
    int find(int x) {
        while (x != f[x]) {
            x = f[x] = f[f[x]];
        }
        return x;
    }
    
    bool same(int x, int y) {
        return find(x) == find(y);
    }
    
    bool merge(int x, int y) {
        x = find(x);
        y = find(y);
        if (x == y) {
            return false;
        }
        siz[x] += siz[y];
        f[y] = x;
        return true;
    }
    
    int size(int x) {
        return siz[find(x)];
    }
};


void solve() {
	int n, m;
	cin >> n >> m;
	vector<string>g(n);
	for (int i = 0; i < n; i ++) {
		cin >> g[i];
	}
	vector<int> cnt(m);
	for (int i = 0; i < n; i ++) {
		for (int j = 0; j < m; j ++) {
			if (g[i][j] == '1') {
				cnt[j] ++;
				cnt[m - j - 1] ++;
			}
		}
	}
	ll ans = 1;

	for (int i = 0; i < m; i ++) {
		if (cnt[i] > 2) {
			cout << 0 << '\n';
			return;
		}
	}

	vector<int> pos(m, -1);
	DSU dsu(n);
	for (int i = 0; i < n; i ++) {
		for (int j = 0; j < m; j ++) {
			if (g[i][j] == '1') {
				if (pos[j] != -1) {
					dsu.merge(pos[j], i);
				} else {
					pos[m - j - 1] = i;
					pos[j] = i;
				}
			}
		}
	}
	map<int, bool> mp;
	for (int i = 0; i < n; i ++) {
		int p = dsu.find(i);
		//cout << p << " \n"[i == n - 1];
		if (!mp.count(p)) {
			mp[p] = 1;
			ans = ans * 2 % mod;
		}  
	}
	cout << ans << '\n';

}

int main() {
	ios::sync_with_stdio(false);
	cin.tie(nullptr); cout.tie(nullptr);
	int t = 1;
	cin >> t;
	while (t --) {
		solve();
	}

	return 0;
}