#include <bits/stdc++.h>
using namespace std;

typedef pair<int, int> PII;
const int mod=1e9+7;
const int N = 1000010;
vector<int> v[N];
bool st[N+10];
long long power(int a, int b){
	long long res = 1;
	while (b){
		if (b & 1) res = res * a % mod;
		b >>= 1;
		a = a * a % mod;
		
	}
	return res;
}
void dfs(int x) {
	for(auto f:v[x]) {
		if(st[f]==1) continue;
		st[f]=1;
		dfs(f);
	}
}
void solve()
{
	int r,c,ans=0,a[N]={0},b[N]={0};
	cin>>r>>c;
	string s[r+10];
	memset(st,0,sizeof st);
	for(int i=1;i<=r;i++) {
		cin>>s[i];
		s[i]=" "+s[i];
		for(int j=1;j<=c;j++) {
			if(s[i][j]=='1') {
				a[j]++;a[c-j+1]++;
				if(a[j]>2||a[c-j+1]>2) {
					cout<<0<<endl;
					return ;
				}
				if(b[j]==0&&b[c-j+1]==0) {
					b[j]=i;
					b[c-j+1]=i;
				}else {
					v[i].push_back(b[j]);
					v[b[j]].push_back(i);
				}
			}
		}
	}
	for(int i=1;i<=r;i++) {
		if(st[i]==0) {
			dfs(i);ans++;
		}
	}
	for(int i=1;i<=r;i++) {
		v[i].clear();
	}
	cout<<power(2,ans)<<endl;
}
int main()
{
    int t = 1;
    cin >> t;
    
    while (t--)
        solve();
    
    return 0;
}