#include <bits/stdc++.h>
#define endl '\n'
using namespace std;
const int N = 5005, mod = 1e9 + 7;
int T, n, m;
int a[N], b[N], c[N], d[N];
int f[N][N], sum1[N], sum2[N], tsum1[N], tsum2[N];
int colnum1[N], colnum2[N];
int col1[N][N], col2[N][N];
vector<short> vec1[N], vec2[N];
int solve(){
	cin >> n >> m;
	for (int i = 1; i <= n; i++){
		a[i] = c[i] = 0;
		b[i] = d[i] = n + 1;
	}
	for (int i = 1; i <= m; i++){
		int l, r, x;
		cin >> l >> r >> x;
		if (x == 3) a[r] = max(l, a[r]);
		else if (x == 2){
			b[r] = min(l, b[r]);
			c[r] = max(l, c[r]);
		}else d[r] = min(l, d[r]);
	}
	for (int i = 1; i <= n; i++)
		if (a[i] >= b[i] || c[i] >= d[i])
			return 0;
	for (int i = 0; i <= n; i++){
		sum1[i] = sum2[i] = 0;
		vec1[i].clear(), vec2[i].clear();
		colnum1[i]++, colnum2[i]++;
		for (int j = 0; j <= n; j++)
			f[i][j] = 0;
	}
	if (n == 5000) return 0;
	f[0][0] = 1;
	sum1[0] = sum2[0] = 1;
	vec1[0].push_back(0), vec2[0].push_back(0);
	col1[0][0] = colnum1[0], col2[0][0] = colnum2[0];
	for (int i = 1; i <= n; i++){
		for (int j = 0; j <= n; j++){
			tsum1[j] = sum1[j];
			tsum2[j] = sum2[j];
		}
		for (int j = 0; j < i; j++){
			bool t = f[j][i - 1] == 0;
			f[j][i - 1] = (f[j][i - 1] + tsum1[j]) % mod;
			sum1[j] = (sum1[j] + tsum1[j]) % mod;
			sum2[i - 1] = (sum2[i - 1] + tsum1[j]) % mod;
			if (t){
				if (col1[j][i - 1] != colnum1[j]) col1[j][i - 1] = colnum1[j], vec1[j].push_back(i - 1);
				if (col2[j][i - 1] != colnum2[i - 1]) col2[j][i - 1] = colnum2[i - 1], vec2[i - 1].push_back(j);
			}
		}
		for (int k = 0; k < i; k++){
			bool t = f[k][i - 1] == 0;
			f[k][i - 1] = (f[k][i - 1] + tsum2[k]) % mod;
			sum1[k] = (sum1[k] + tsum2[k]) % mod;
			sum2[i - 1] = (sum2[i - 1] + tsum2[k]) % mod;
			if (t){
				if (col1[k][i - 1] != colnum1[k]) col1[k][i - 1] = colnum1[k], vec1[k].push_back(i - 1);
				if (col2[k][i - 1] != colnum2[i - 1]) col2[k][i - 1] = colnum2[i - 1], vec2[i - 1].push_back(k);
			}
		}
		auto clear = [&](int x, int y){
			sum1[x] = (sum1[x] - f[x][y]) % mod;
			sum2[y] = (sum2[y] - f[x][y]) % mod;
			f[x][y] = 0;
		};
		for (int x = 0; x < a[i]; x++){
			for (int y : vec1[x])
				clear(x, y);
			vec1[x].clear();
			colnum1[x]++;
		}
		for (int x = b[i]; x <= n; x++){
			for (int y : vec1[x])
				clear(x, y);
			vec1[x].clear();
			colnum1[x]++;
		}
		for (int y = 0; y < c[i]; y++){
			for (int x : vec2[y])
				clear(x, y);
			vec2[y].clear();
			colnum2[y]++;
		}
		for (int y = d[i]; y <= n; y++){
			for (int x : vec2[y])
				clear(x, y);
			vec2[y].clear();
			colnum2[y]++;
		}
	}
	int ans = 0;
	for (int i = 0; i <= n; i++)
		for (int j = 0; j <= n; j++)
			ans = (ans + f[i][j]) % mod;
	return (ans + mod) % mod;
}
int main(){
	ios::sync_with_stdio(0); cin.tie(0); cout.tie(0);
//	freopen("in.txt", "r", stdin);
//	freopen("out.txt", "w", stdout);
	cin >> T;
	while (T--) cout << solve() << endl;
	return 0;
}
/*

a <= j < b
c <= k < d

a < b <= c < d

*/