#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 dis[ N*2 ];

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

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 );
    if( b>a ) swap( a , b );
    dis[ a ]=b;
}

vector<int> G[ N ];

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

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

    for( int i=1; i<=n; i++ ){
        cin>>S[ i ];
        for( int j=0; j<m; j++ ){
            if( S[ i ][ j ]=='1' ){
//                G[ m-j ].emplace_back( i );
                G[ j+1 ].emplace_back( i );
            }
        }
    }

    bool ok=false;
    int cnt=0;
    for( int i=1; i<=m; i++ ){
        for( int e=1; e<G[ i ].size(); e++ ){
            int a=G[ i ][ 0 ] , b=G[ i ][ e ];
            //if( a==b ) continue;
            //cout<<a<<" "<<b<<endl;
            add( a+n , b ) , add( a , b+n );
            if( find( a+n )==find( b+n ) || find( a )==find( b ) ) ok=true;
        }

        int o=0;
        for( auto x:G[ i ] ){
            if( find( x )==x ) o++;
        }
        if( o!=G[ i ].size() ) ok=true;

    }

    for( int i=1; i<=n; i++ ){
        reverse( S[ i ].begin() , S[ i ].end() );
        for( int e=0; e<m; e++ ){
            if( S[ i ][ e ]=='1' ){
                for( auto son:G[ e+1 ] ){
                    if( son==i ) continue;
                    //cout<<i<<" "<<son<<endl;
                    add( son , i ); add( son+n , i+n );
                    if( find( son+n )==find( i ) || find( i+n )==find( son ) ) ok=true;
                    //if( find( son+n )==find( i+n ) || find( i )==find( son ) ) continue;
                }

                int o=0;
                for( auto x:G[ e+1 ] ){
                    if( find( x )==x ) o++;
                }
                if( o!=G[ e+1 ].size() ) ok=true;

            }
        }
        reverse( S[ i ].begin() , S[ i ].end() );
    }

    if( ok ) cout<<0<<endl;
    else {
        int ans = 1;
        init( n );

        for( int i=1; i<=n; i++ ){
            for( int j=0; j<m; j++ ){
                if( S[ i ][ j ]=='1' ){
                    G[ m-j ].emplace_back( i );
                    //G[ j+1 ].emplace_back( i );
                }
            }
        }

        for( int i=1; i<=m; i++ ){
            for( int j=1; j<G[ i ].size(); j++ ){
                add( G[ i ][ j ] , G[ i ][ 0 ] );
            }
        }

        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;
}