#include<bits/stdc++.h>
using namespace std;
#define fast ios::sync_with_stdio(false),cin.tie(0),cout.tie(0)
#define int long long
#define PII pair<int,int>
int r,c,x,cnt,ii,jj,iii,jjj;
const int N=1e6+7;
const int mod=1e9+7;
int f[N],f1[N],f2[N],vis[N],a[N];
string s;
vector<int> v[N];
void init(){
    for(int i=1;i<=r;i++)
        f[i]=f1[i]=f2[i]=i,v[i].clear(),vis[i]=0;
}
int qpow(int a,int m){
    int res=1,x=a;
    while(m)
    {
        if(m&1) res=res*x%mod;//
        x=x*x%mod;m>>=1;
    }
    return (res+mod)%mod;
}
void solve(){
    cin>>r>>c;
    init();
    for(int i=1;i<=r;i++)
    {
        cin>>s;
        for(int j=0;j<s.length();j++)
        {
            x=s[j]-'0';
            v[i].push_back(x);
        }        
    }

    //cout<<r<<"!"<<c<<endl;
    /*for(int i=1;i<=r;i++){
        for(int j=1;j<=c;j++)
        {
            cout<<v[i][j-1];
        }
        cout<<endl;
    }*/
    for(int j=1;j<=c/2;j++)
    {
        cnt=0;
        for(int i=1;i<=r;i++)
        {
            if(v[i][j-1]==1) cnt++;//,cout<<i<<"!"<<j<<"!!"
            if(v[i][c-j]==1) cnt++;//,cout<<i<<"!"<<c-i+1<<"!!"
        }
        if(cnt>=3) 
        {
            cout<<0<<endl;return ;
        }
    }
    cnt=0;
    if(c&1)
    {
        for(int i=1;i<=r;i++)
        {
            if(v[i][c/2]==1) cnt++;
        }
        if(cnt>=2)
        {
            cout<<0<<endl;return ;
        }
    }
    for(int j=1;j<=c/2;j++)
    {
        cnt=0;ii=jj=iii=jjj=0;
        for(int i=1;i<=r;i++)
        {
            if(v[i][j-1]==1)
            {
                cnt++;
                if(cnt==1)
                {
                    ii=i,jj=j;
                }
                else 
                {
                    iii=i,jjj=j;
                }
            }
            if(v[i][c-j]==1)
            {
                cnt++;
                if(cnt==1)
                {
                    ii=i,jj=c-j+1;
                }
                else 
                {
                    iii=i,jjj=c-j+1;
                }                
            }
            if(cnt==2) break;
        }
        if(cnt<2) continue;
        else{            
            if(ii>iii) swap(ii,iii);
            if(jj>jjj) swap(jj,jjj);

            if(ii==iii) continue;
            else{//if(j==1)cout<<ii<<"!"<<jj<<"!!"; 
                if(jj==jjj)
                {
                    if(f2[iii]==ii) 
                    {
                        cout<<0<<endl;return;
                    }
                    f1[iii]=ii; //
                }
                else{
                    if(f1[iii]==ii) 
                    {
                        cout<<0<<endl;return;
                    }
                    f2[iii]=ii;                  
                }
            }
        }
    }
    //int tp=0;
    for(int i=1;i<=r;i++)
    {      

        int tp=f2[i],flg=0,tmp=f1[i];
        while(f2[i+1]==tp||i<=r-1) 
        {
            i++;
            if(f1[i]==tmp) 
            {
                cout<<0<<endl;return;
            }
            tmp=f1[i];
        }
    }
    for(int i=1;i<=r;i++)
    {
        f[i]=min(f1[f2[i]],f2[f1[i]]);
        vis[f[i]]=1;
    }
    cnt=0;
    for(int i=1;i<=r;i++)
    {
        if(vis[i]) cnt++;
    }
    //cout<<f[2]<<"!"<<endl;
    cout<<qpow(2ll,cnt)<<endl;
}
signed main(){
    fast;
    int t=1;
    cin>>t;
    while(t--){
        solve();
    }
}