#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(x!=p[x]) return p[x]=find(p[x]);
	      return 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];
		//	cout<<mp[i][j]<<endl;
		}
	}
//	for(int i=1;i<=m;i++){
//		for(int j=1;j<=n;j++){
//			if(mp[j][i]=='1') cout<<mp[j][i]<<' ';
//		//	cout<<mp[i][j]<<endl;
//		}
//		cout<<endl;
//	}
	map<int,int>cnt;
//	for(int i=1;i<=m;i++) cout<<cnt[i]<<' ';
	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]++;
		   	 //  cnt[i]+=1;
		   	//  cout<<i<<' '<<cnt[i]<<'\n';
		   	  id[i].push_back(j);
		   }
		}
	}
////	for(int i=1;i<=m;i++) cout<<cnt[i]<<' '<<(int)id[i][0]<<endl;
	if(m%2){
		int id=m/2+1;
		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];
	//	int rowl1=id[i][],row2=id[ot];
		if(sum>2){
			cout<<0<<'\n';
			return;
		}
		if(sum<2) continue;
	//	cout<<i<<' '<<ot<<' '<<cnt[i]<<' '<<cnt[ot]<<endl;
		int r1=-n,r2=-n;
		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);
		//	cout<<"???"<<cur<<' '<<(int)id[cur].size()<<endl;
			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(p[i]==i) u++;
	cout<<qsumi(2,u/2)<<'\n';
}
signed main()
{
	ios::sync_with_stdio(false);
	cin.tie(0),cout.tie(0);
	int t;
	cin>>t;
	while(t--) solve();
	return 0;
}