#include <bits/stdc++.h>
#define SZ(x) (int) x.size() - 1
#define all(x) x.begin(), x.end()
#define ms(x, y) memset(x, y, sizeof x)
#define F(i, x, y) for (int i = (x); i <= (y); i++)
#define DF(i, x, y) for (int i = (x); i >= (y); i--)
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
template <typename T> void chkmax(T &x, T y) { x = max(x, y); }
template <typename T> void chkmin(T &x, T y) { x = min(x, y); }
template <typename T> void read(T &x) {
	x = 0; int f = 1; char c = getchar();
	for (; !isdigit(c); c = getchar()) if (c == '-') f = -f;
	for (; isdigit(c); c = getchar()) x = (x << 1) + (x << 3) + (c ^ 48);
	x *= f;
}
const int N = 2e5 + 10;
int n, tot, dfn[N], rdfn[N], s[N], ff[N], t[N];
vector <int> v[N];
set <int> vv[N];
void dfs(int x, int fa) {
	// cout << "# " << x << " " << v[x].size() << endl;
	ff[x] = fa;
	if (v[x].size() == 1) t[dfn[x] = ++tot] = x;
	else dfn[x] = tot + 1;
	for (int i: v[x])
		if (i != fa) dfs(i, x);
	rdfn[x] = tot;
}
void zhk() {
	tot = 0;
	read(n);
	F(i, 1, n) s[i] = 0, v[i].clear(), vv[i].clear();
	// cout << "N: " << n << endl;
	F(i, 1, n - 1) {
		int x, y; read(x), read(y);
		// cout << x << " " << y << endl;
		vv[x].insert(y);
		vv[y].insert(x);
	}
	if (n == 2) {
		puts("YES\n1 2\n1 2");
		return;
	}
	while (true) {
		bool flag = true;
		F(i, 1, n) {
			if (vv[i].size() == 2) {
				// cout << "!! " << i << endl;
				flag = false;
				int a = *vv[i].begin(), b = *vv[i].rbegin();
				vv[a].erase(i);
				vv[b].erase(i);
				vv[i].clear();
				vv[a].insert(b);
				vv[b].insert(a);
			}
		}
		if (flag) break;
	}
	F(i, 1, n)
		for (int j: vv[i])
			v[i].push_back(j);
	int rt = -1;
	F(i, 1, n)
		if (v[i].size() && v[i].size() != 1) {
			rt = i;
			break;
		}
	if (rt == -1) {
		F(i, 1, n) if (v[i].size()) cout << i << " "; cout << '\n';
		F(i, 1, n) if (v[i].size()) cout << i << " "; cout << '\n'; cout << '\n';
		return;
	}
	dfs(rt, 0);
	// cout << "! " << rt << " " << tot << endl;
	// F(i, 1, tot) cout << t[i] << " "; cout << endl;
	// F(i, 1, n)
	vector <pair <int, int>> ans;
	// puts("beg:");
	// F(i, 1, n) {
	// 	for (int j: v[i])
	// 		cout << i << " " << j << endl;
	// }
	// 	cout << "! " << rt << " " << v[rt].size() << " " << n << " " << tot << endl;
	if (tot % 2 == 0) {
		F(i, 1, n) {
			if (i == rt) continue;
			if (v[i].empty()) continue;
			if ((dfn[i] + rdfn[i] + 1) % 2) continue;
			int id = ((dfn[i] + rdfn[i]) / 2 - 1) % (tot / 2) + 1;
			// cout << "! " << i << " " << id << endl;
			// cout << "! " << (dfn[i] + rdfn[i]) / 2 << endl;
			// cout << "! " << i << " " << dfn[i] << " " << rdfn[i] << " " << id << endl;
			if (!s[id]) s[id] = i;
			else s[id] = -1;
		}
		F(i, 1, tot / 2) {
			if (s[i] == 0) {
				// cout << "@ " << i << endl;
				// cout << "! " << tot << endl;
				F(j, 1, tot / 2 + 1) {
					if (t[(i - j + tot) % tot + 1] == t[(i + j - 1) % tot + 1]) ans.emplace_back(t[(i - j + tot) % tot + 1], rt);
					else ans.emplace_back(t[(i - j + tot) % tot + 1], t[(i + j - 1) % tot + 1]);
				}
				break;
			}
			if (s[i] != -1) {
				// cout << "~ " << i << endl;
				// cout << "! " << i << " " << s[i] << endl;
				int ff;
				F(j, 1, tot / 2) {
					int A = t[(i - j + tot) % tot + 1], B = t[(i + j - 1) % tot + 1];
					if (dfn[s[i]] <= dfn[A] && rdfn[A] <= rdfn[s[i]]) {
						if (j == 1) ff = 1;
						else if (ff == 2) {
							ans.emplace_back(A, rt);
							ff = 3;
						}
					} else {
						if (j == 1) ff = 2;
						else if (ff == 1) {
							ans.emplace_back(s[i], rt);
							ff = 3;
						}
					}
					ans.emplace_back(A, B);
					// if (t[(i - j - a + tot) % tot + 1] == t[(i + j + b - 1) % tot + 1]) ans.emplace_back(t[(i - j - a + tot) % tot + 1], rt);
					// else ans.emplace_back(t[(i - j - a + tot) % tot + 1], t[(i + j + b - 1) % tot + 1]);
				}
				// puts("!");
				// cout << "! " << tot << endl;
				// int len = (rdfn[s[i]] - dfn[s[i]] + 1) / 2;
				// F(j, 1, len) ans.emplace_back(t[(i - j + tot) % tot + 1], t[(i + j - 1) % tot + 1]);
				// // cout << tot << " " << dfn[s[i]] << "  " << rdfn[s[i]] << " " << len << endl;
				// // ans.emplace_back(s[i], ff[s[i]]);
				// // if (m & 1)
				// ans.emplace_back(t[(i - len - 1 + tot) % tot + 1], s[i]);
				// F(j, len + 1, tot / 2) ans.emplace_back(t[(i - j - 1 + tot) % tot + 1], t[(i + j - 1) % tot + 1]);
				break;
			}
		}
	} else {
		F(i, 1, n) {
			if (v[i].empty()) continue;
			if (i == rt) continue;
			if ((dfn[i] + rdfn[i] + 1) % 2) {
				if (tot % 2) {
					int id = ((dfn[i] + rdfn[i]) / 2 + ((tot + 1) / 2) - 2) % tot + 1;
					// cout << "!! " << i << " " << dfn[i] << " " << id << endl;
					if (!s[id]) s[id] = i;
					else s[id] = -1;
				}
				continue;
			}
			int id = ((dfn[i] + rdfn[i]) / 2 - 1) % tot + 1;
			// cout << "! " << i << " " << id << endl;
			// cout << "! " << (dfn[i] + rdfn[i]) / 2 << endl;
			// cout << "! " << i << " " << dfn[i] << " " << rdfn[i] << " " << id << endl;
			if (!s[id]) s[id] = i;
			else s[id] = -1;
		}
		F(i, 1, tot) {
			// cout << "! " << i << " " << s[i] << endl;
			if (s[i] == 0) {
				// cout << "! " << tot << endl;
				F(j, 1, tot / 2 + 1) {
					if (t[(i - j + tot) % tot + 1] == t[(i + j - 1) % tot + 1]) ans.emplace_back(t[(i - j + tot) % tot + 1], rt);
					else ans.emplace_back(t[(i - j + tot) % tot + 1], t[(i + j - 1) % tot + 1]);
				}
				break;
			}
			if (s[i] != -1) {
				int a = 0;
				bool flag = true;
				F(j, 1, tot / 2) {
					int A = t[(i - j - a + tot) % tot + 1], B = t[(i + j - 1) % tot + 1];
					if (dfn[s[i]] <= dfn[A] && rdfn[A] <= rdfn[s[i]] && flag) {
						ans.emplace_back(A, rt);
						a++;
						A = t[(i - j - a + tot) % tot + 1];
						flag = false;
					}
					ans.emplace_back(A, B);
					// if (t[(i - j - a + tot) % tot + 1] == t[(i + j + b - 1) % tot + 1]) ans.emplace_back(t[(i - j - a + tot) % tot + 1], rt);
					// else ans.emplace_back(t[(i - j - a + tot) % tot + 1], t[(i + j + b - 1) % tot + 1]);
				}
				if (flag) {
					ans.emplace_back(s[i], rt);
				}
				// puts("!");
				// cout << "! " << tot << endl;
				// int len = (rdfn[s[i]] - dfn[s[i]] + 1) / 2;
				// F(j, 1, len) ans.emplace_back(t[(i - j + tot) % tot + 1], t[(i + j - 1) % tot + 1]);
				// // cout << tot << " " << dfn[s[i]] << "  " << rdfn[s[i]] << " " << len << endl;
				// // ans.emplace_back(s[i], ff[s[i]]);
				// // if (m & 1)
				// ans.emplace_back(t[(i - len - 1 + tot) % tot + 1], s[i]);
				// F(j, len + 1, tot / 2) ans.emplace_back(t[(i - j - 1 + tot) % tot + 1], t[(i + j - 1) % tot + 1]);
				break;
			}
		}
	}
	// cout << ans.size() << '\n';
	for (auto i: ans) cout << i.first << " " << i.second << '\n'; cout << '\n';
	assert(ans.size() == tot / 2 + 1);
}
signed main() {
	int _ = 1;
	cin >> _;
	F(i, 1, _) {
		// if (i == 331) {
		// 	int n; cin >> n;
		// 	cout << n << endl;
		// 	F(i, 1, n - 1) {
		// 		int x, y; cin >> x >> y;
		// 		cout << x << " " << y << endl;
		// 	}
		// 	return 0;
		// }
		zhk();
	}
	return 0;
}