# include <bits/stdc++.h>
using namespace std;
using ll=long long;
constexpr int cof[]={1,3,6,6},mod=1e9+7;
int L1[5010],L2[5010],R1[5010],R2[5010],f[5010][5010],s1[5010],s2[5010];
template<typename T>
inline void upd(int &x,const T &y){x=(x+y)%mod;}
inline void chkmin(int &x,const int &y){x=(x<y?x:y);}
inline void chkmax(int &x,const int &y){x=(x>y?x:y);}
int solve()
{
	int n,m,l,r,x,x1=0,x2=0,y1=0,y2=0,ans=0;
	scanf("%d%d",&n,&m);
	fill(L1+1,L1+n+1,0);
	fill(L2+1,L2+n+1,0);
	iota(R1+1,R1+n+1,0);
	iota(R2+1,R2+n+1,0);
	for(int i=0;i<m;i++)
	{
		scanf("%d%d%d",&l,&r,&x);
		if(x==1) chkmin(R1[r],l-1),chkmin(R2[r],l-1);
		else if(x==2) chkmax(L1[r],l),chkmin(R2[r],l-1);
		else chkmax(L1[r],l+1),chkmax(L2[r],l);
	}
	for(int i=2;i<=n;i++) chkmax(L1[i],L1[i-1]),chkmax(L2[i],L2[i-1]);
	for(int i=n-1;i;i--) chkmin(R1[i],R1[i+1]),chkmin(R2[i],R2[i+1]);
	for(int i=1;i<=n;i++)
		if(L1[i]>R1[i] || L2[i]>R2[i])
			return 0;
	s1[0]=s2[0]=f[0][0]=1;
	for(int i=0;i<n;i++)
	{
		for(int j=max(x1,1);j<=x2;j++)
		{
			int v=s1[j];
			upd(f[i][j],v);
			if(i>=y1 && i<=y2) upd(s1[i],v),upd(s2[j],v);
		}
		if(i) for(int j=y1;j<=y2;j++)
		{
			int v=s2[j];
			upd(f[i][j],v);
			if(i>=x1 && i<=x2) upd(s1[i],v),upd(s2[j],v);
		}
		for(;x1<L1[i+1];x1++)
			for(int j=y1;j<=y2;j++)
				upd(s1[x1],mod-f[x1][j]),upd(s2[j],mod-f[x1][j]);
		for(;x2<R1[i+1];x2++)
			for(int j=y1;j<=y2;j++)
				upd(s1[x2+1],f[x2+1][j]),upd(s2[j],f[x2+1][j]);
		for(;y1<L2[i+1];y1++)
			for(int j=x1;j<=x2;j++)
				upd(s1[j],mod-f[j][y1]),upd(s2[y1],mod-f[j][y1]);
		for(;y2<R2[i+1];y2++)
			for(int j=x1;j<=x2;j++)
				upd(s1[j],f[j][y2+1]),upd(s2[y2+1],f[j][y2+1]);
	}
	for(int i=x1;i<=x2;i++)
		for(int j=y1;j<=y2;j++)
			upd(ans,(ll)f[i][j]*cof[(i>0)+(j>0)+1]);
	for(int i=0;i<=n;i++) memset(f[i],0,sizeof(int)*(n+1));
	memset(s1,0,sizeof(int)*(n+1));
	memset(s2,0,sizeof(int)*(n+1));
	return ans;
}
int main()
{
	int T;
	cin>>T;
	while(T--) printf("%d\n",solve());
	return 0;
}