#include <bits/stdc++.h>
const int mod=1e9+7;
using namespace std;
int n,m,l[5010],r[5010],R[5010],C[5010],f[5010][5010],lk[5010],rk[5010],lj[5010],rj[5010];
void qmo(int &x){
	x+=(x>>31)&mod;
} 
void add(int x,int y,int w){
	qmo(w);
	qmo(f[x][y]+=w-mod);
	qmo(R[x]+=w-mod);
	qmo(C[y]+=w-mod);
}
void clr(int j,int lk,int rk){
	while(l[j]<=r[j]&&l[j]<lk) add(j,l[j],mod-f[j][l[j]]),l[j]++;
	while(l[j]<=r[j]&&r[j]>rk) add(j,r[j],mod-f[j][r[j]]),r[j]--;
}
void sol(){
	scanf("%d%d",&n,&m);
	for(int i=0;i<=n;i++) R[i]=C[i]=0,lk[i]=lj[i]=l[i]=0,rk[i]=rj[i]=r[i]=n;
	for(int i=0;i<=n;i++) for(int j=0;j<=n;j++) f[i][j]=0;
	add(0,0,1);
	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==3) lk[r]=max(lk[r],l);
		if(x==2) rk[r]=min(rk[r],l-1),lj[r]=max(lj[r],l);
	}
	for(int i=2;i<=n;i++){
		for(int k=0;k<i-1;k++) add(i-1,k,C[k]); 
		for(int j=1;j<i-1;j++) add(i-1,j,R[j]);
		for(int j=0;j<i;j++){
			if(j<lj[i]||rj[i]<j) clr(j,1,0);
			else clr(j,lk[i],rk[i]);
		}
//		for(int j=0;j<i;j++) for(int k=0;k<=j;k++) if(f[j][k]) printf("f[%d][%d][%d]=%d\n",i,j,k,f[j][k]);
	}
	int ans=0;
	for(int j=1;j<=n;j++) ans=(ans+6ll*R[j])%mod;
	ans=(ans+3ll*R[0])%mod;
	printf("%d\n",ans); 
}
int main(){
	int _;
	scanf("%d",&_);
	while(_--) sol();
}