#include<bits/stdc++.h>
using namespace std;
#define LL long long
struct edge{
	LL to,nt;
}a[10000005];
LL n,i,j,k,m,t,cnt=0,ans=0;
LL sum[1000005],stat[1000005],nxt[1000005];
string s[1000005];
map<LL,bool> mp;
bool vis[1000005];
bool flag;
const LL mod=1e9+7;
void add(LL x,LL y){
	//printf("%lld %lld\n",x,y);
	a[++cnt].to=y;a[cnt].nt=nxt[x];nxt[x]=cnt;
}
void dfs(LL x){
	vis[x]=true;
	for(LL i=nxt[x];i;i=a[i].nt){
	    if(vis[a[i].to]==false) stat[a[i].to]=stat[x]^1,dfs(a[i].to);
	    else if(stat[a[i].to]==stat[x]) flag=true;
	}
}
int main(){
	cin>>t;
	while(t--){
		for(i=1;i<=n;i++)
		  nxt[i]=0,stat[i]=0,vis[i]=false;
		for(i=1;i<=m;i++)
		  sum[i]=0;
		for(i=1;i<=cnt;i++)
		  a[i].to=0,a[i].nt=0;
		mp.clear();
		cnt=0;
		cin>>n>>m;
		for(i=1;i<=n;i++){
			s[i]="";
			cin>>s[i];s[i]=" "+s[i];
		}
		for(i=1;i<=n;i++){
			for(j=1;j<=m;j++)
			  if(s[i][j]=='1') sum[j]++;
		}
		flag=false;
		for(i=1;i<=(m)/2;i++){
			if(sum[i]+sum[m-i+1]<=2){
		  		if(sum[i]+sum[m-i+1]==2){
		  			LL x=0,y=0;
		  			for(j=1;j<=n;j++)
		  		  	  if(s[j][i]=='1' || s[j][m-i+1]=='1'){
		  		  		if(x==0) x=j;
		  		  		else y=j;
				  	  }
				    if(y>0){
				    	add(x,y),add(y,x);
				    	if(x>y) swap(x,y);
				    	if(mp[x*n+y]==true) flag=true; 
				    	mp[x*n+y]=true;
					}
				}
		  	}
		  	else{
		  		flag=true;break;
			  }
		}
		if(m%2==1){
		    if(sum[(m+1)>>1]>1){
				flag=true;
			}
		}
		if(flag==true){
			printf("0\n");
			continue;
		}
		LL ans=1;
		for(i=1;i<=n;i++)
		  if(vis[i]==false) dfs(i),ans*=2,ans%=mod;
		if(flag==true){
			printf("0\n");
			continue;
		}
		printf("%lld\n",ans);
	}
	return 0;
}