#pragma GCC optimize(2)
#include <bits/stdc++.h>
typedef long long ll;
using namespace std;
int head[1000001];
int vis[1000001];
int tot=0;
ll mod=1e9+7;
vector<int> vec[1000001];
ll fastpow(ll base,ll po){
	ll ans=1;
	while(po){
		if(po&1){
			ans*=base;
			ans%=mod;
		}
		base*=base;
		base%=mod;
		po>>=1;
	}
	return ans;
}
struct chainforwardstar{
	int to,next;
}e[2000001];
int r,c;
void addedge(int from,int to){
	e[++tot].next=head[from];
	head[from]=tot;
	e[tot].to=to;
}
int cnt[1000001];
void dfs(int x){
	vis[x]=1;
	for(int i=head[x];i;i=e[i].next){
		int to=e[i].to;
		if(!vis[to]){
			dfs(to);
		}
	}
}
void solve(){
	cin>>r>>c;
	string s;
	for(int p=1;p<=r;p++){
		cin>>s;
		for(int i=0;i<c;i++){
			if(s[i]=='1'){
				cnt[i+1]++;
				cnt[c-i]++;
				int idx=min(i+1,c-i);
				vec[idx].push_back(p);
			}
		}
	}
	for(int i=1;i<=c;i++){
		if(cnt[i]>2){
			cout<<0<<'\n';
			return;
		}
	}
	for(int i=1;i<=c;i++){
		if(vec[i].size()>=2){
			addedge(vec[i][0],vec[i][1]);
			addedge(vec[i][1],vec[i][0]);
		}
	}
	ll ct=0;
	for(int i=1;i<=r;i++){
		if(!vis[i]){
			dfs(i);
			ct++;
		}
	}
	cout<<fastpow(2LL,ct)<<'\n';
}
int main(){
	std::ios::sync_with_stdio(false);
	cin.tie(0);
	cout.tie(0);
	int t;
	cin>>t;
	while(t--){
		solve();
		for(int i=1;i<=r;i++){
			head[i]=0;
			vis[i]=0;
		}
		for(int i=1;i<=c;i++){
			cnt[i]=0;
			vec[i].clear();
		}
		tot=0;
	}
	return 0;
}