#include <bits/stdc++.h>
using namespace std;
#pragma G++ optimize(2)
#define debug(x) cout << "[debug] " #x << " = " << x << '\n';
#define double long double
#define int long long
#define xiaowen ac
typedef unsigned long long ull;
typedef long long ll;
typedef pair<int, int> PII;
const int N = 2e6 + 7;
const int mod =1e9+7;
int a[N];
int fa[N];
string s[N];
int find(int x) {
	if (fa[x]==x) {
		return x;
	}
	return fa[x]=find(fa[x]);
}
void unionn(int x,int y) {
	x=find(x);
	y=find(y);
	if (x==y) {
		return;
	}
	fa[x]=y;
}
int quick_pow(int a,int b) {
	int ans=1;
	while (b) {
		if (b&1) {
			ans=ans*a%mod;
		}
		a=a*a%mod;
		b>>=1;
	}
	return ans%mod;
}
void solved() {
	int n,m;
	cin>>n>>m;
	for (int i=1; i<=2*n; i++) {
		fa[i]=i;
	}
	for (int i=1; i<=n; i++) {
		cin>>s[i];
		s[i]=' '+s[i];
	}
	vector<int>v[2*m];
	for (int i=1; i<=n; i++) {
		for (int j=1; j<=m; j++) {
			if (s[i][j]=='1') {
				v[j].push_back(i);
			}
		}
	}
	if (m&1) {
		if (v[(m+1)/2].size()>=2) {
			cout<<"0\n";
			return;
		}
	}
	for (int i=1; i<=m/2; i++) {
		int j=m-i+1;
		int ans1=v[i].size();
		int ans2=v[j].size();
		if (ans1+ans2>2) {
			cout<<"0\n";
			return;
		} else if (ans1+ans2<=1) {
			continue;
		} else {
			if (ans1&&ans2) {
				if (v[i][0]==v[j][0]) {
					continue;
				}
				unionn(v[i][0],v[j][0]);
				unionn(v[i][0]+n,v[j][0]+n);
			} else {
				if (ans1) {
//					if (v[i][0]==v[i][1]) {
//						continue;
//					}
					unionn(v[i][0],v[i][1]+n);
					unionn(v[i][0]+n,v[i][1]);
				} else {
//					if (v[j][0]==v[j][1]) {
//						continue;
//					}
					unionn(v[j][0],v[j][1]+n);
					unionn(v[j][0]+n,v[j][1]);
				}
			}
		}
	}
	int ans=0;
	for (int i=1; i<=n; i++) {
		if (find(i)==find(i+n)) {
//			cout<<i<<"****"<<"\n";
			cout<<"0\n";
			return;
		}

	}
	for (int i=1; i<=2*n; i++) {
		find(i);
		if (fa[i]==i) {
			ans+=1;
		}
	}
	cout<<quick_pow(2,ans/2)%mod<<"\n";
}
signed main() {
	ios::sync_with_stdio(false);
	cin.tie(nullptr), cout.tie(nullptr);
	int T;
	// T = 1;
	cin >> T;
	while (T--) {
		solved();
	}
	return 0;
}