#define TLE ios::sync_with_stdio(0); cin.tie(0); cout.tie(0)
#include <bits/stdc++.h>
using namespace std;
#ifndef ONLINE_JUDGE
#include "../template/debug.h"
#else 
#define dbg(...) ;
#endif
using ll =long long;
#define int ll
#define inf 0x3f3f3f3f
#define LINF 0x3f3f3f3f3f3f3f3f
#define rep(i,a,n) for (int i=a;i<=n;i++)
#define per(i,n,a) for (int i=n;i>=a;i--)
#define all(x) (x).begin(),(x).end()
#define fi first
#define se second
#define SZ(x) ((int)(x).size())
typedef vector<int> VI;
#define pb push_back
typedef pair<int,int> PII;
const int mod=1e9+7;
// const int mod=998244353;
typedef double db;
mt19937 rnd(time(0));
ll powmod(ll a,ll b) {ll res=1;a%=mod; assert(b>=0); for(;b;b>>=1){if(b&1)res=res*a%mod;a=a*a%mod;}return res;}
ll gcd(ll a,ll b) { return b?gcd(b,a%b):a;}
const int N = 1e6+5;
vector<int >v1[N];
int n,m;
int fa[N],sz1[N],sz2[N];
int vis[N];
int find(int x){
	if(x==fa[x]) return x;
	return fa[x]=find(fa[x]);
}
void solve(){
	cin>>n>>m;
	int flag=1;
	for(int i=1;i<=n;i++){
		string s;
		cin>>s;
		// cout<<s<<'\n';
		for(int j=0;j<m;j++){
			int p=s[j]-'0';
			// cout<<p<<'\n';
			if(p){
				v1[j+1].push_back(i);
			}
			if(v1[j+1].size() + v1[m-j].size() > 2) flag=0;
		}
	}
	if(!flag) {
		cout<<0<<'\n';
		return ;
	}
	for(int i=1;i<=n;i++){
		fa[i]=i;
		sz1[i]=1;
		sz2[i]=0;
	}
	for(int i=1;i<=m;i++){
		if(v1[i].empty()) continue;
		int x,y;
		int tmp;
		if(v1[i].size()==1){
			x=v1[i][0];
			tmp=0;
			if(v1[m-i+1].size()==1){
				y=v1[m-i+1][0];
			}
			else{
				continue;
			}
		}
		else{
			tmp=1;
			x=v1[i][0];
			y=v1[i][1];
		}
		// cout<<x<<' '<<y<<'\n';
		int u=find(x);
		int v=find(y);
		if(u!=v){
			fa[u]=v;
			sz1[v]+=sz1[u];
			sz1[u]=0;
			sz2[v]+=sz2[u];
			sz2[u]=0;
		}
		else{
			vis[v]=1;
		}
		sz2[v]+=tmp;
	}
	int ans=1;
	for(int i=1;i<=n;i++) {
		if(i==find(i)) {
			// cout<<i<<" "<<sz1[i]<<' '<<sz2[i]<<' '<<sz3[i]<<'\n';
			ans*=2;
			ans%=mod;
			if(sz1[i]==2 && sz2[i]==2) {
				ans=0;
			}
			if(sz1[i]==2 && sz2[i]==1 && vis[i]) {
				ans=0;
			}
			if(sz2[i]==sz1[i] && sz1[i]&1) {
				ans=0;
			} 
		}
	}
	cout<<ans<<'\n';
}
void clear(){
	for(int i=1;i<=m;i++) v1[i].clear();
	for(int i=1;i<=n;i++) sz1[i]=sz2[i]=vis[i]=0;
}
signed main(){
	TLE;
	int t;
	cin>>t;
	while(t--){
		solve();
		clear();
	}
	return 0;
}