#include<bits/stdc++.h>
using namespace std;
#define mid ((l+r)>>1)
const int N=5e5+5;
vector<int>e[N],E[N];
struct scan{
	int s,t;
}a[N];
int top[N],du[N],dfn[N],son[N],size[N],tot,lson[N],rson[N],line[N],fa[N],dep[N];
void add(int u,int v) {
	E[u].push_back(v);
	du[v]++;
}//线段树建DAG
void build1(int &rt,int l,int r) {
	rt=++tot;
	if(l==r) return;
	build1(lson[rt],l,mid);
	build1(rson[rt],mid+1,r);
	add(lson[rt],rt);
	add(rson[rt],rt);
}//入树 
void build2(int &rt,int l,int r) {
	rt=++tot;
	if(l==r) return;
	build2(lson[rt],l,mid);
	build2(rson[rt],mid+1,r);
	add(rt,lson[rt]);
	add(rt,rson[rt]);
}//出树 
void update1(int rt,int l,int r,int x,int y){
	if(l==r) return add(y,rt);
	if(x<=mid) update1(lson[rt],l,mid,x,y);
	else update1(rson[rt],mid+1,r,x,y);
}
void update2(int rt,int l,int r,int x,int y){
	if(l==r) return add(rt,y);
	if(x<=mid) update2(lson[rt],l,mid,x,y);
	else update2(rson[rt],mid+1,r,x,y);
}
//单点 
void query1(int rt,int l,int r,int x,int y,int z){
	if(l>y||r<x) return;
	if(x<=l && r<=y) return add(rt,z);
	query1(lson[rt],l,mid,x,y,z);query1(rson[rt],mid+1,r,x,y,z);
}
void query2(int rt,int l,int r,int x,int y,int z){
	if(l>y||r<x) return;
	if(x<=l&&r<=y) return add(z,rt);
	query2(lson[rt],l,mid,x,y,z);query2(rson[rt],mid+1,r,x,y,z);
}
//区间修改 
void dfs1(int u) {
	son[u]=0;
	size[u]=1;
	for(auto v:e[u]) {
		if(v==fa[u]) continue;
		fa[v]=u;
		dep[v]=dep[u]+1;
		dfs1(v);
		size[u]+=size[v];
		if(size[v]>size[son[u]]) son[u]=v;
	}
}
void dfs2(int u) {
	dfn[u]=++tot;
	if(u==son[fa[u]]) top[u]=top[fa[u]];
	else top[u]=u;
	if(son[u]) dfs2(son[u]);
	for(auto v:e[u]) {
		if(v==fa[u]||v==son[u])continue;
		dfs2(v);
	}
}//树链剖分
void solve() {
	int n,m,rt1=0,rt2=0;
	scanf("%d",&n);
	for(int i=1,u,v;i<n;i++) {
		scanf("%d%d",&u,&v);
		e[u].push_back(v);
		e[v].push_back(u);
	}
	dfs1(1);
	dfs2(1);
	scanf("%d",&m);
	tot=m;
	build1(rt1,1,n);
	build2(rt2,1,n);
	for(int i=1;i<=m;i++) scanf("%d%d",&a[i].s,&a[i].t);
	for(int i=1;i<=m;i++) update1(rt1,1,n,dfn[a[i].s],i),update2(rt2,1,n,dfn[a[i].t],i);
	
	for(int u,v,i=1;i<=m;i++){
		u=a[i].s;v=a[i].t;
		while(top[u]!=top[v])
		{
			if(dep[top[u]]<dep[top[v]]) swap(u,v);
			query1(rt1,1,n,dfn[top[u]]+(top[u]==a[i].s),dfn[u]-(u==a[i].s),i);
			query2(rt2,1,n,dfn[top[u]]+(top[u]==a[i].t),dfn[u]-(u==a[i].t),i);
			u=fa[top[u]];
		}
		if(dep[u]>dep[v]) swap(u,v);
		query1(rt1,1,n,dfn[u]+(u==a[i].s),dfn[v]-(v==a[i].s),i);
		query2(rt2,1,n,dfn[u]+(u==a[i].t),dfn[v]-(v==a[i].t),i);
	}
	int head=0,last=0;
	for(int i=1;i<=tot;i++)
		if(!du[i]) line[last++]=i;
	while(head<last){
		int u=line[head++];
		for(auto v:E[u])
			if(!(--du[v])) line[last++]=v;
	}
	if(head==tot) printf("Yes\n");
	else printf("No\n"); 
	
	for(int i=1;i<=tot;i++) du[i]=lson[i]=rson[i]=0,E[i].clear();
	for(int i=1;i<=n;i++)e[i].clear();
	tot=0;
}
int main(){
	int T; 
	scanf("%d",&T);
	while(T--) solve();
	return 0;
}