#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*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 );
	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();
	for( int i=1; i<=m; i++ ) A[ i ]=0;
	init( n*2 );
	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' ){
				A[min(j+1,m-j)]++;
				if(A[min(j+1,m-j)]>=3) ok=true;
//                G[ m-j ].emplace_back( i );
				G[ j+1 ].emplace_back( i );
			}
		}
	}
	
	
	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;
			
		}
	}
	
	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;
					if( (find( son+n )==find( i ) )||(find( i+n )==find( son ) )) ok=true;
					add( son+n , i+n ) , add( son , i );
					if( find( son )==find( son+n ) || find( i )==find( i+n ) ) ok=true;
					
					
				}
			}
		}
		reverse( S[ i ].begin() , S[ i ].end() );
	}
	
	for( int i=1; i<=n; i++ ){
		if( find( i )==find( i+n ) ) ok=true;
	}
	
	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;
}