#include <iostream>
using namespace std;
const int mod=1e9+7;
const int N = 1e6+10;
int p[N];
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 solve()
{
	int r,c,ans=0;
	cin>>r>>c;
	int a[r+10],b[c+10];
	string s;
	for(int i=0;i<=r;i++) p[i]=i;
	for(int i=0;i<=c;i++) a[i]=0,b[i]=0;
	for(int i=1;i<=r;i++) {
		cin>>s;
		s=" "+s;
		for(int j=1;j<=c;j++) {
			if(s[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 {
					p[i]=p[b[j]];
				}
			}
		}
	}
	for(int i=1;i<=r;i++) {
		if(p[i]==i) ans++;
		p[i]=0;
	}
	cout<<power(2,ans)<<endl;
	return ;
}
int main()
{
    int t = 1;
    cin >> t;
    while (t--)
        solve();
    return 0;
}