#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 j=0; j<m; j++ ){
            if( S[ i ][ j ]=='1' ){
                if( A[ j+1 ]==1 ){
                    reverse( S[ i ].begin() , S[ i ].end() );
                    break;
                }
            }
        }

        for( int j=0; j<m; j++ ){
            if( S[ i ][ j ]=='1' ){
                if( A[ j+1 ]==1 ){
                    ok=true;
                    break;
                }
                else A[ j+1 ]=1;
            }
        }
    }

    if( ok ) cout<<0<<endl;
    else{
        for( int i=1; i<=n; i++ ){
            for( int j=0; j<m; j++ ){
                if( S[ i ][ j ]=='1' ){
                    int a1=j+1 , a2=m-j;

                    G[ a1 ].emplace_back( i );
                    G[ a2 ].emplace_back( i );

//                    if( nex[ a1 ]==0 ) nex[ a1 ]=i;
//                    else{
//                        if( find( nex[ a1 ] ) != find( i ) ) add( nex[ a1 ] , i );
//                        nex[ a1 ]=i;
//                    }
//
//                    if( nex[ a2 ]==0 ) nex[ a2 ]=i;
//                    else{
//                        if( find( nex[ a2 ] ) != find( i ) ) add( nex[ a2 ] , i );
//                        nex[ a2 ]=i;
//                    }
                }
            }
        }

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