#include <bits/stdc++.h>
using namespace std;
using LL = long long;
const int N = 1e6 + 10;
const LL mod = 1e9 + 7;

string S[N];
int n, m, check[N], fa[N];
int find(int x) {
    if(x != fa[x]) fa[x] = find(fa[x]);
    return fa[x];
}

void solve()
{
    
    cin >> n >> m;
    vector <int> ck(m);
    for(int i = 1; i <= n; i++) cin >> S[i];
    for(int i = 1; i <= n; i++) {
        int f = 1;
        for(int j = 0; j < m; j++) {
            if(S[i][j] == '1' && ck[j]) {
                f = 0;
                break;
            }
        }
        if(f) {
            for(int j = 0; j < m; j++) {
                if(S[i][j] == '1') ck[j] = i;
            }
        }
        else {
            for(int j = 0; j < m; j++) {
                if(S[i][j] == '1' && ck[m - j - 1]) {
                    cout << 0 << '\n';
                    return;
                }
            }
            for(int j = 0; j < m; j++) {
                if(S[i][j] == '1') ck[m - j - 1] = i;
            }
        }
    }

    vector <int> belo(m);
    for(int i = 1; i <= n; i++) fa[i] = i;

    for(int i = 1; i <= n; i++) {
        for(int j = 0; j < m; j++) {
            if(S[i][j] == '1') {
                if(belo[j]) {
                    int rj = find(belo[j]);
                    fa[rj] = i;
                }
                else belo[j] = i;
                if(belo[m - j - 1]) {
                    int rj = find(belo[m - j - 1]);
                    fa[rj] = i;
                }
                else belo[m - j - 1] = i;
            }
        }
    }

    LL ans = 1;
    for(int i = 1; i <= n; i++) {
        if(find(i) == i) ans *= 2, ans %= mod;
    }
    cout << ans << '\n';
}
int main()
{
    ios::sync_with_stdio(0);
    cin.tie(0);
    int t;
    cin >> t;
    while(t--) solve();
    return 0;
}