#include<bits/stdc++.h>
#define ll long long
#define MOD 1000000007
using namespace std;
int sakuzyo,n,m,mnlmt[5003][4],mxlmt[5003][4],k1,k2,k3,k4,k5,k6,k7,k8,k9;
ll ans,vv[5003]; 
struct dpv{
	int nxth;
	int nxtw;
	ll val;
}dp[5003][5003];
ll sumh[5003],sumw[5003];
int headh[5003],headw[5003];
void setval(int X,int Y,ll v){
	dp[X][Y].val=v;
	sumh[X]=(sumh[X]+v)%MOD;
	sumw[Y]=(sumw[Y]+v)%MOD;
	dp[X][Y].nxth=headh[X];
	dp[X][Y].nxtw=headw[Y];
	headh[X]=Y;
	headw[Y]=X;
	return;
}
void clnh(int X){
	for(int i=headh[X];i!=-1;i=dp[X][i].nxth){
		sumw[i]=(sumw[i]-dp[X][i].val)%MOD;
		dp[X][i].val=0;
	}
	headh[X]=-1;
	sumh[X]=0;
	return;
}
void clnw(int X){
	for(int i=headw[X];i!=-1;i=dp[i][X].nxtw){
		sumh[i]=(sumh[i]-dp[i][X].val)%MOD;
		dp[i][X].val=0;
	}
	headw[X]=-1;
	sumw[X]=0;
	return;
}
void getlmt(int X){
	for(int i=X-1;i>=mnlmt[X][1];i--)clnh(i);
	for(int i=0;i<mxlmt[X][2];i++)clnh(i);
	for(int i=X-1;i>=mnlmt[X][2];i--)clnw(i);
	for(int i=0;i<mxlmt[X][3];i++)clnw(i);
	return;
}
void sol(){
	cin>>n>>m;
	for(int i=0;i<=n+1;i++){
		for(int j=0;j<4;j++){
			mnlmt[i][j]=n+1;
			mxlmt[i][j]=0;
		}
	}
	for(int i=1;i<=m;i++){
		cin>>k1>>k2>>k3;
		mnlmt[k2][k3]=min(mnlmt[k2][k3],k1);
		mxlmt[k2][k3]=max(mxlmt[k2][k3],k1);
	}
	for(int i=0;i<=n;i++){
		sumh[i]=sumw[i]=0;
		headh[i]=headw[i]=-1;
		for(int j=0;j<=n;j++){
			dp[i][j].nxth=dp[i][j].nxtw=-1;
			dp[i][j].val=0;
		}
	}
	setval(0,0,3);
	getlmt(1);
	for(int i=2;i<=n;i++){
		for(int j=0;j<i-1;j++)vv[j]=(sumh[j]+sumw[j])%MOD;
		for(int j=0;j<i-1;j++)setval(i-1,j,vv[j]);
		getlmt(i);
	}
	ans=0;
	for(int i=0;i<=n;i++)ans=(ans+sumh[i])%MOD;
	ans%=MOD;
	ans+=MOD;
	ans%=MOD;
	cout<<ans<<'\n';
	return;
}
signed main(){
	ios::sync_with_stdio(false);
	cin>>sakuzyo;
	while(sakuzyo--)sol();
	return 0;
}