#include<bits/stdc++.h>
#define LB lower_bound
using namespace std;using ll=long long;
const int N=1e6+5,mod=1e9+7;
int n,m;
string A[N];
vector<int> S[N];
int fa[N*2];
int GF(int x){return fa[x]^x?fa[x]=GF(fa[x]):x;}
void merge(int x,int y){
	fa[GF(x)]=GF(y);
}
void Solve(){
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++) cin>>A[i];
	iota(fa+1,fa+2*n+1,1);
	for(int i=0;i<m;i++) S[i].clear();
	for(int i=1;i<=n;i++){
		for(int j=0;j<m&&j<=m-j-1;j++)if(A[i][j]=='1'||A[i][m-j-1]=='1') S[j].push_back(i);
	}
	for(int i=0;i<m&&i<=m-i-1;i++){
		if(S[i].size()>2){puts("0");return;}
		if(S[i].size()<=1) continue;
		int x=S[i][0],y=S[i][1];
		for(int j:{x,y}){
			if(A[j][i]=='1'&&A[j][m-i-1]=='1'){
				puts("0");return;
			}
		}
		if(A[x][i]==A[y][i]) merge(x,y+n),merge(x,y+n);
		else merge(x,y),merge(x+n,y+n);
	}
	for(int i=1;i<=n;i++) if(GF(i)==GF(i+n)){puts("0");return;}
	int cnt=0;
	for(int i=1;i<=2*n;i++) if(fa[i]==i) cnt++;
	cnt>>=1;
	ll ans=1;
	for(int i=1;i<=cnt;i++) ans=ans*2%mod;
	printf("%lld\n",ans);
}
int main(){
	int t=1;
	scanf("%d",&t);
	while(t--) Solve();
	return 0;
}