#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll mod=1000000007;
vector<int>v[1000050];
int sumc[2000050]={0};
int root[2000050]={0};
int vis[2000050]={0};
int find(int x){
	if(x!=root[x]){
		root[x]=find(root[x]);
	}
	return root[x];
}
void merge(int a,int b){
	int ra=find(a);
	int rb=find(b);
	root[ra]=root[rb];
	return ;
}
ll qm(ll a,ll x){
	ll ans=1;
	while(x){
		if(x%2){
			ans=(ans%mod*a%mod)%mod;
		}
		x/=2;
		a=(a%mod)*(a%mod)%mod;
	}
	return ans%mod;
}
void solve(){
	int r,c;
	scanf("%d %d",&r,&c);
	for(int i=1;i<=c;i++){sumc[i]=0;}
	for(int i=1;i<=max(r,c);i++){
		v[i].clear();
	}
	for(int i=1;i<=r;i++){
		string s;
		cin>>s;
		for(int j=1;j<=c;j++){
			int tem=(int)s[j-1]-48;
			sumc[j]+=tem;
			if(tem==1){
				v[j].push_back(i);//列的1的行
			}
		}
	}
	for(int i=1;i<=c;i++){
		if(sumc[i]+sumc[c-i+1]>2){
			printf("0");
			return ;
		}
	}
	for(int i=1;i<=2*max(r,c);i++){
		root[i]=i;
	}
	for(int i=1;i<=c;i++){
		if(sumc[i]+sumc[c-i+1]==2){
			if(sumc[i]==0) continue;
			if(sumc[i]==1){//相同
				merge(v[c-i+1][0],v[i][0]);
				merge(v[c-i+1][0]+r,v[i][0]+r);
			}else{//相反
				merge(v[i][1],v[i][0]+r);
				merge(v[i][1]+r,v[i][0]);
			}
		}
	}
	for(int i=1;i<=2*r;i++){
		vis[i]=0;
	}
	for(int i=1;i<=r;i++){
		if(find(i)==find(i+r)){
			printf("0");return;
		}
		vis[find(i)]=1;
	}
	for(int i=2*r;i>r;i--){
		vis[find(i)]=1;
	}
	ll siz=0;
	for(int i=1;i<=2*r;i++){
		siz+=(vis[i]==0);
	}
	ll ans=1;
	ans=qm(2,siz);
	printf("%lld",ans);
	
}
int main(){ 
	int t;
	scanf("%d",&t);
	while(t--){
		solve();
		if(t) printf("\n");
	}	
	
	
	return 0;
}