#include<bits/stdc++.h>
using namespace std;
#define int long long
#define lowbit(x)  ((x) & - (x))
using PII = pair<int, int>;
using PIII = tuple<int, int, int>;
const int MOD = 1e9 + 7;
inline int mypow(int n, int k, int p = MOD) {
	int r = 1;
	for (; k; k >>= 1, n = n * n % p) {
		if (k & 1) r = r * n % p;
	}
	return r;
}
struct cmp {  //priority_queue
	bool operator()(const PIII& a, const PIII& b) {
		return std::get<0>(a) < std::get<0>(b); //down
	}
};

bool cmpsort(const PII& a, const PII& b) {  //sort
	return a.first < b.first; //up
}
int num[1000010]={0};
int fd(int k){ //find()
	//路径压缩
    if(num[k]==k)return k;
    return num[k]=fd(num[k]);
    //相当于num[k]=fd(num[k]); return num[k];
}
int n,m;
vector <string> mp;
string str;
inline void solve() {
    cin>>n>>m;
    for(int i=0;i<2*n;i++){
        num[i]=i;
    }
    for(int i=0;i<n;i++){
        cin>>str;
        mp.push_back(str);
    }
    int F=0;
    if(m%2==1){
        int cnt=0;
        for(int i=0;i<n;i++){
            if(mp[i][m/2]=='1'){
                cnt++;
            }
        }
        if(cnt>=2){
            cout<<0<<'\n';
            return ;
        }
    }
    for(int j=0;j<m/2;j++){
        int cnt=0;
        int f1=-1;
        int f2=-1;
        for(int i=0;i<n;i++){
            if(mp[i][j]=='1'||mp[i][m-1-j]=='1'){
                cnt++;
                if(f1==-1)f1=i;
                else f2=i;
            }
        }
        if(cnt==2){
            if(mp[f1][j]=='1'&&mp[f2][j]=='1'||mp[f1][m-j-1]=='1'&&mp[f2][m-j-1]=='1') F=1;
        }
        if(cnt>2){
            cout<<0<<'\n';
            return ;
        }
        else if(cnt==1) continue;
        else{
            if(F==1){
                if(fd(f1)!=fd(f2+n)){
                    num[fd(f1)]=fd(f2+n);
                }
                if(fd(f1+n)!=fd(f2)){
                    num[fd(f1+n)]=fd(f2);
                }
            }
            else{
                if(fd(f1)!=fd(f2)){
                    num[fd(f1)]=fd(f2);
                }
                if(fd(f1+n)!=fd(f2+n)){
                    num[fd(f1+n)]=fd(f2+n);
                }
            }
        }
    }
    for(int i=0;i<n;i++){
        if(fd(i)==fd(i+n)){
            cout<<0<<'\n';
            return ;
        }
    }   
    map<int,int> e;
    for(int i=0;i<2*n;i++){
        e[num[i]]++;
    }
    int ans=e.size()/2;
    cout<<mypow(2,ans)<<'\n';
}
signed main() {
	ios::sync_with_stdio(0), cin.tie(0); cout.tie(0);
	//cout << fixed << setprecision(12);
	int t = 1;
	cin >> t;
	while (t--) {
		solve();
        mp.clear();
        for(int i=0;i<=2*n;i++){
            num[i]=0;
        }
	}
}