#include<bits/stdc++.h>
using namespace std;
const int N=5e3+5,mod=1e9+7;
int t,n,m,Lj[N],Rj[N],Lk[N],Rk[N],sx[N],sy[N],f[N][N],ans;
deque<int>q[N];
void upd(int x,int y,int v){
	q[y].push_back(x);
	(sx[x]+=v)%=mod,(sy[y]+=v)%=mod,(f[x][y]+=v)%=mod;
}
void del(int x,int y){
	(sx[x]-=f[x][y])%=mod,(sy[y]-=f[x][y])%=mod,f[x][y]=0;
}
void work(int i){
	for(int k=0;k<i;k++){
		int lj=Lj[i],rj=Rj[i];
		if(!(Lk[i]<=k&&k<=Rk[i])) lj=n;
		while(q[k].size()&&q[k].front()<lj)
			del(q[k].front(),k),q[k].pop_front();
		while(q[k].size()&&q[k].back()>rj)
			del(q[k].back(),k),q[k].pop_back();
	}
}
signed main(){
	scanf("%d",&t);
	while(t--){
		scanf("%d%d",&n,&m),ans=0;
		for(int i=0;i<=n;i++)
			sx[i]=sy[i]=Lj[i]=Lk[i]=0,Rj[i]=Rk[i]=i-1,
			fill(f[i],f[i]+1+n,0),q[i].clear();
		for(int i=1,l,r,x;i<=m;i++){
			scanf("%d%d%d",&l,&r,&x);
			if(x==1) Rj[r]=min(Rj[r],l-1);
			if(x==2) Lj[r]=max(Lj[r],l),Rk[r]=min(Rk[r],l-1);
			if(x==3) Lk[r]=max(Lk[r],l);
		}
		upd(0,0,3),work(1);
		for(int i=1;i<n;i++){
			for(int j=0;j<i;j++) upd(i,j,(sx[j]+sy[j])%mod);
			work(i+1);
		}
		for(int i=0;i<n;i++)
			for(int j=0;j<n;j++) (ans+=f[i][j])%=mod;
		printf("%d\n",ans);
	}
	return 0;
}