#include <bits/stdc++.h>
#define lb(x) (x&-x)
#define L(i,j,k) for(int i=(j);i<=(k);++i)
#define R(i,j,k) for(int i=(j);i>=(k);--i)
#define swap(a,b) (a^=b^=a^=b)

using namespace std;
using i64 = long long;

typedef pair<int, int> pii;
typedef long long ll;
typedef unsigned long long ull;
void chmin(int &x, int c) {
    x = min(x, c);
}
void chmax(int &x, int c) {
    x = max(x, c);
}

namespace fast_IO{
#define IOSIZE 100000
int precision=3,POW[10]={1,10,100,1000,10000,100000,1000000,10000000,100000000,1000000000};char ibuf[IOSIZE],obuf[IOSIZE],*p1=ibuf,*p2=ibuf,*p3=obuf;
#ifdef ONLINE_JUDGE
#define getchar()((p1==p2)and(p2=(p1=ibuf)+fread(ibuf,1,IOSIZE,stdin),p1==p2)?(EOF):(*p1++))
#define putchar(x)((p3==obuf+IOSIZE)&&(fwrite(obuf,p3-obuf,1,stdout),p3=obuf),*p3++=x)
#define isdigit(ch)(ch>47&&ch<58)
#define isspace(ch)(ch<33)
#endif
template<typename T>inline T read(){T s=0;int w=1;char ch;while(ch=getchar(),!isdigit(ch)&&(ch!=EOF))if(ch==45)w=-1;if(ch==EOF)return 0;while(isdigit(ch))s=s*10+ch-48,ch=getchar();return s*w;}template<typename T>inline bool read(T&s){s=0;int w=1;char ch;while(ch=getchar(),!isdigit(ch)&&(ch!=EOF))if(ch==45)w=-1;if(ch==EOF)return 0;while(isdigit(ch))s=s*10+ch-48,ch=getchar();return s*=w,1;}inline bool read(char&s){while(s=getchar(),isspace(s));return 1;}inline bool read(char*s){char ch;while(ch=getchar(),isspace(ch));if(ch==EOF)return 0;while(!isspace(ch))*s++=ch,ch=getchar();*s='\000';return 1;}template<typename T>inline void print(T x){if(x<0)putchar(45),x=-x;if(x>9)print(x/10);putchar(x%10+48);}inline void print(char x){putchar(x);}inline void print(char*x){while(*x)putchar(*x++);}inline void print(const char*x){for(int i=0;x[i];i++)putchar(x[i]);}inline bool read(std::string&s){s="";char ch;while(ch=getchar(),isspace(ch));if(ch==EOF)return 0;while(!isspace(ch))s+=ch,ch=getchar();return 1;}inline void print(std::string x){for(int i=0,n=x.size();i<n;i++)putchar(x[i]);}inline bool read(bool&b){char ch;while(ch=getchar(),isspace(ch));b=ch^48;return 1;}inline void print(bool b){putchar(b+48);}inline bool read(double&x){int a=0,b=0;char ch=getchar();bool f=0;while(ch<48||ch>57){if(ch==45)f=1;ch=getchar();}while(47<ch&&ch<58){a=(a<<1)+(a<<3)+(ch^48);ch=getchar();}if(ch!=46){x=f?-a:a;return 1;}ch=getchar();while(47<ch&&ch<58){b=(b<<1)+(b<<3)+(ch^48),ch=getchar();}x=b;while(x>1)x/=10;x=f?-a-x:a+x;return 1;}inline void print(double x){if(x<0){putchar(45),x=-x;}x=round((long double)x*POW[precision])/POW[precision],print((long long)x),putchar(46),x-=(long long)(x);for(int i=1;i<=precision;i++)x*=10,putchar(x+48),x-=(int)x;}template<typename T,typename...T1>inline int read(T&a,T1&...other){return read(a)+read(other...);}template<typename T,typename...T1>inline void print(T a,T1...other){print(a),print(other...);}struct Fast_IO{~Fast_IO(){fwrite(obuf,p3-obuf,1,stdout);}}io;template<typename T>Fast_IO&operator>>(Fast_IO&io,T&b){return read(b),io;}template<typename T>Fast_IO&operator<<(Fast_IO&io,T b){return print(b),io;}
#define cout io
#define cin io
#define endl '\n'
}using namespace fast_IO;
const int maxn = 1e6 + 10, mod = 1e9 + 7;
int N, M, p[maxn], f[maxn << 1];
int sz[maxn], ct[maxn];
string S[maxn];
vector<pii> cnt[maxn];
int gf (int x) {
	return p[x] == x ? x : p[x] = gf(p[x]);
}
void mg (int x, int y) {
	int px = gf (x);
	int py = gf (y);
	if (px != py) {
		p[px] = py, sz[py] += sz[px];
	}
}
int qmi (int a, int b) {
	int res = 1;
	while (b) {
		if (b & 1) {
			res = 1ll * res * a % mod;
		}
		a = 1ll * a * a % mod;
		b >>= 1;
	}
	return res;
}

void solve() {
  cin >> N >> M;
  L (i, 1, N * 2) p[i] = i, sz[i] = 1, f[i] = 0;
  L (i, 1, M) cnt[i].clear();
  int ctt = 0;
  L (i, 1, N) cin >> S[i], S[i] = ' ' + S[i];
  L (i, 1, N) L (j, 1, N) ctt += (S[i][j] == '1');
  L (i, 1, N) {
    L (j, 1, M) {
      if (S[i][j] == '0') continue;
      cnt[min(j, M - j + 1)].push_back(make_pair(j, i));
      ct[j] ++ ;
    }
  }
  L (i, 1, M) {
  	if (cnt[i].size() > 2) {
  		cout << "0\n";
  		return ;
		}
	} 
  L (i, 1, M) {
    if (cnt[i].size() < 2) continue;
    auto x = cnt[i][0];
    auto y = cnt[i][1];
//    cout << x.first << ' ' << x.second << '\n';
//    cout << y.first << ' ' << y.second << '\n';
    if (x.first != y.first) {
    	mg (x.second, y.second);
    	mg (x.second + N, y.second + N);
		} else {
//			cout << x.second << "->" << y.second << '\n';
			if (x.first == (M + 1) / 2 && (M & 1)) {
				cout << 0 << '\n';
				return ;
			}
			mg (x.second + N, y.second);
			mg (x.second, y.second + N);
		} 
  }
  L (i, 1, N) {
  	if (gf(i) == gf(i + N)) {
  		cout << "0\n";
  		return ;
		}
	}
  int c = 0;
  L (i, 1, N * 2) {
  	c += gf(i) == i;
	}
  cout << qmi(2, c / 2) << '\n';
}

signed main() {
  int T = 1;
  cin >> T;
  while (T--)solve();
  return 0;
}