#include <bits/stdc++.h>

using namespace std;

#define int long long
#define deg( x ) cout<<""#x"="<<x<<endl
#define endl '\n'
#define pll pair<int,int>
#define fr frist
#define se second

const int N=1e6+10;
const int mod=1e9+7;

string S[ N ];

int A[ N ];

int dis[ N ];
int cnt[ N ];

void init( int n ){
    for( int i=1; i<=n; i++ ) dis[ i ]=i,cnt[ i ]=1;
}

int find( int x ){
    if( x==dis[ x ] ) return x;
    dis[ x ]=find( dis[ x ] );
    return find( dis[ x ] );
}

void add( int i , int j ){
    int a=find( i ) , b=find( j );
    dis[ a ]=b;
   // cnt[ b ]+=cnt[ a ];
}

vector<int> G[ N ];

void solve(  ){
    int n,m;
    cin>>n>>m;

    for( int i=1; i<=m; i++ ) G[ i ].clear();
    for( int i=1; i<=m; i++ ) A[ i ]=0;
    init( n );

    bool ok=false;
    for( int i=1; i<=n; i++ ){
        cin>>S[ i ];
        for( int e=0; e<m; e++ ){
            if( S[ i ][ e ]=='1' ){
                int a1=e+1 , a2=m-e;

                G[ a1 ].emplace_back( i );
                if( a1!=a2 ) G[ a2 ].emplace_back( i );
            }
        }
    }

    if( m&1 ){
        if( G[ (m+1)/2 ].size()>1 ) ok=true;
    }

    for( int i=1; i<=m; i++ ){
        if( G[ i ].size() >= 3 ){
            ok=true;
            break;
        }
    }

    if( ok ) cout<<0<<endl;
    else{
        for( int i=1; i<=m; i++ ){
            for( int e=1; e<G[ i ].size(); e++ ){
                add( G[ i ][ e ] , G[ i ][ 0 ] );
            }
        }

        int ans=1;
        for( int i=1; i<=n; i++ ){
            if( dis[ i ]==i ){
                //if( cnt[ find( i ) ]==1 )
                ans*=2, ans%=mod;
                //else ans*=cnt[ find( i ) ], ans%=mod;
            }
        }

        cout<<ans<<endl;
    }

}

signed main() {
    ios::sync_with_stdio( 0 );
    cin.tie( 0 ); cout.tie( 0 );

    int T=1;
    cin>>T;

    while( T-- ){
        solve(  );
    }

    return 0;
}