#include<bits/stdc++.h>
#define int long long
using namespace std;
const int mod=1e9 + 7;
int qsumi(int a,int b){
	int res=1;
	while(b){
		if(b&1) res=res*a%mod;
		b>>=1;
		a=a*a%mod;
	}
	return res;
}
void solve()
{
	int n,m;
	cin>>n>>m;
	vector<int>p(2*n+10); 
	std::function<int(int)>find=[&](int x)->int{
	      if (p[x]==x) {
		return x;
	}
	return p[x]=find(p[x]);
	};
	for(int i=1;i<=2*n+1;i++) p[i]=i;
	vector<vector<char>>mp(n+10,vector<char>(m+10));
	for(int i=1;i<=n;i++){
		for(int j=1;j<=m;j++){
			cin>>mp[i][j];
		}
	}
	vector<int>cnt(m+1,0);
	map<int,vector<int>>id;
	for(int i=1;i<=m;i++){
		for(int j=1;j<=n;j++){
		   if(mp[j][i]=='1'){
		   	cnt[i]++;
		   	  id[i].push_back(j);
		   }
		}
	}
	if(m%2){
		int id=(m+1)/2;
		if(cnt[id]>1){
			cout<<0<<'\n';
			return;
		}
	}
	for(int i=1;i<=m/2;i++){
		int ot=m-i+1;
		int sum=cnt[i]+cnt[ot];
		if(sum>2){
			cout<<0<<'\n';
			return;
		}
		if(sum<2) continue;
		int r1,r2;
		if(cnt[i]==cnt[ot]) r1=id[i][0],r2=id[ot][0];
		else{
			int cur=i;
			if(cnt[cur]<cnt[ot]) swap(cur,ot);
			r1=id[cur][0],r2=id[cur][1];
		}
		if(r1==r2){
			continue;
		}
		int r11=r1+n,r22=r2+n;
		int fa1=find(r1),fa2=find(r2);
		int fa11=find(r11),fa22=find(r22);
		p[fa1]=fa22,p[fa2]=fa11; 
	}
	for(int i=1;i<=n;i++){
		int fa1=find(i),fa2=find(i+n);
		if(fa1==fa2){
			cout<<0<<'\n';
			return;
		}
	}
	int u=0;
	for(int i=1;i<=2*n;i++)if(find(p[i])==i) u++;
	cout<<qsumi(2,u/2)%mod<<'\n';
}
signed main()
{
	ios::sync_with_stdio(false);
	cin.tie(0),cout.tie(0);
	int t;
	cin>>t;
	while(t--) solve();
	return 0;
}