#include "bits/stdc++.h"
using namespace std;
const int Len = 1e6 + 5 , mod = 1e9 + 7;
int ans,fa[Len << 1],flg[Len],flw[Len],N,n,m;
inline void makeSet(int x){for(int i = 1 ; i <= x ; i ++) flg[i] = flw[i] = 0 , fa[i] = i;}
int findSet(int x){return fa[x] == x ? fa[x] : fa[x] = findSet(fa[x]);}
inline int swp(int x){if(x > n) return x - n;return x + n;}
inline void unionSet(int x,int y)
{
	int u = findSet(x) , v = findSet(y);
	if(u != v) fa[u] = v;
	u = findSet(swp(x)) , v = findSet(swp(y));
	if(u != v) fa[u] = v;
}
vector<int> vc[Len];
char s[Len];
signed main()
{
	int T;scanf("%d",&T);
	while(T --)
	{
		scanf("%d %d",&n,&m);ans = 1;
		for(int i = 1 ; i <= m ; i ++)
		{
			vector<int> clr;
			swap(clr , vc[i]);
		}
		for(int i = 1 ; i <= n ; i ++)
		{
			scanf("%s",s + 1);
			for(int j = 1 ; j <= m ; j ++) if(s[j] == '1') vc[j].push_back(i);
		}
		N = (n << 1);
		makeSet(N);
		for(int i = 1 ; i <= m / 2 + (m % 2) ; i ++)
		{
			int SZ = (int)vc[i].size() + (i != m - i + 1) * (int)vc[m - i + 1].size();
			if(SZ >= 3 || (i == m - i + 1 && SZ >= 2)) ans = 0;
			else if(SZ == 2)
			{
				if((int)vc[i].size() == 2) unionSet(vc[i][0] , vc[i][1]);
				else if((int)vc[m - i + 1].size() == 2) unionSet(vc[m - i + 1][0] , vc[m - i + 1][1]);
				else 
				{
					if(vc[m - i + 1][0] != vc[i][0]) unionSet(swp(vc[m - i + 1][0]) , vc[i][0]);
				}
			}
		}
		for(int i = 1 ; i <= n ; i ++) if(findSet(i) == findSet(i + n)) ans = 0;
		int ret = 0;
		for(int i = 1 ; i <= N ; i ++)
		{
			int f = findSet(i);
			if(flw[f]) continue;
			flw[f] = 1;
			ret ++;
		}
		for(int i = 1 ; i <= (ret >> 1) ; i ++) ans = 1ll * ans * 2 % mod;
		printf("%d\n",ans);
	}
	return 0;
}