#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <string>
#include <map>
#include <queue>
#include <bitset>
#include <vector>
#include <stack>
#include <cmath>
#include <set>
#include <cassert>
#include <random>
#include <chrono>
#include <ctime>
#include <cstdlib>
#define pii pair<int,int> 
#define mp make_pair
#define pb emplace_back
#define debug() cout<<"qwq"<<endl
typedef long long ll;
typedef double db;
typedef long double ldb;
typedef unsigned long long ull;
const int DMAX = 1000000 + 5;
const double INF = 1e9;
const int MOD = 1000000000 + 7;
using namespace std;
mt19937 rnd(chrono::steady_clock::now().time_since_epoch().count());
template<class T> inline void read(T &x){
    T f=1; x=0; char ch=getchar();
    while(ch<'0' || ch>'9'){
        if(ch=='-') f=-1; 
        ch=getchar();
    }
    while(ch<='9' && ch>='0'){ x=x*10+ch-'0';ch=getchar(); }
    x*=f;
}
int T;
char s[DMAX];
int n, m;
vector<vector<int> > a;
int sum[DMAX];
bool vst[DMAX];
void clear(){
    for(int i = 0; i < m; i++) vst[i] = 0;
}
void solve(){
    read(n), read(m);
    a.clear();
    a.resize(n);
    for(int i = 0; i < m; i++) vst[i] = 0;
    for(int i = 0; i < n; i++) a[i].resize(m);
    for(int i = 0; i < m; i++) sum[i] = 0;
    for(int i = 0; i < n; i++){
        scanf("%s", s);
        for(int j = 0; j < m; j++){
            a[i][j] = s[j] - '0';
            if(a[i][j]){
                sum[j] ++;
            }
        }
    }  
    int qwe = 0;
    for(int j = 0; j < m; j++){
        qwe += sum[j];
        if(sum[j] >= 3){
            puts("0");
            clear();
            return ;
        }
    }
    if(qwe > m){
        puts("0");
        clear();
        return ;
    }
    for(int j = 0; j < m; j++){
        int x = 0;
        if(j != m - 1 - j){
            x = sum[j] + sum[m - 1 - j];
            if(x >= 3){
                puts("0");
                clear();
                return ;
            }
        }
        else{
            x = sum[j];
            if(x >= 2){
                puts("0");
                clear();
                return ;
            }
        }
    }
    int ans = 1;
    for(int i = 0 ; i < n; i++){
        bool flag = 1;
        for(int j = 0; j < m; j++){
            if(a[i][j]){
                if(vst[j]) flag = 0;
                // vst[j] =  vst[m - 1 - j] = 1;
            }
        }
        for(int j = 0; j < m; j++){
            if(a[i][j]){
                vst[j] =  vst[m - 1 - j]  = 1;
            }
        }
        if(flag){
            ans = 2ll * ans % MOD;
        }
    }
    printf("%d\n", ans);
    clear();
    return ;
}
int main(){
    read(T);
    while(T--){
        solve();
    }
    return 0;
}