#include<bits/stdc++.h>
using namespace std;
const int N=5e6+9;
int Q,n,m,deg[N],dfn[N],tot,son[N],sz[N],top[N],fa[N],d[N],rt1,rt2;
int id1[N],id2[N],cnt,ls[N<<1],rs[N<<1],id=1;
vector<int>e[N],t[N];
void add(int x,int y){
	deg[y]++;
	e[x].push_back(y);
}
void prework(){
	memset(dfn,0,sizeof(dfn));
	memset(sz,0,sizeof(sz));
	memset(top,0,sizeof(top));
	memset(fa,0,sizeof(fa));
	memset(d,0,sizeof(d));
	memset(son,0,sizeof(son));
	memset(deg,0,sizeof(deg));
	for(int i=1;i<=n;i++){
	  	t[i].clear();
  	}
  	for(int i=1;i<=cnt;i++){
		e[i].clear(); 
  	}
	tot=cnt=0;
	for(int i=1;i<=id;i++){
		ls[i]=rs[i]=0;
	} 
	rt1=rt2=0;
	for(int i=1;i<=n;i++){
		id1[i]=id2[i]=0;
	} 
	id=1;
}
void builds(int &p,int l,int r){
    p=++id;
    if(l==r){
		id1[l]=p; 
		return ;
	} 
	int mid=(l+r)>>1;
    builds(ls[p],l,mid);
	builds(rs[p],mid+1,r);
    add(ls[p],p);
	add(rs[p],p);
}
void buildt(int &p,int l,int r){
    p=++id;
    if(l==r){
		id2[l]=p; 
		return ;
	} 
	int mid=(l+r)>>1;
    buildt(ls[p],l,mid);
	buildt(rs[p],mid+1,r);
    add(p,ls[p]);
	add(p,rs[p]);
}
void add1(int p,int l,int r,int L,int R,int c){
    assert(L);
    if(l==L&&r==R){	
		add(p,c); 
		return ;
	} 
	int mid=(l+r)>>1;
    if(R<=mid) add1(ls[p],l,mid,L,R,c);
    else if(L>mid) add1(rs[p],mid+1,r,L,R,c);
    else add1(ls[p],l,mid,L,mid,c),add1(rs[p],mid+1,r,mid+1,R,c);
}
void add2(int p,int l,int r,int L,int R,int c){
    if(l==L&&r==R){
		add(c,p); 
		return ;
	} 
	int mid=(l+r)>>1;
    if(R<=mid) add2(ls[p],l,mid,L,R,c);
    else if(L>mid) add2(rs[p],mid+1,r,L,R,c);
    else add2(ls[p],l,mid,L,mid,c),add2(rs[p],mid+1,r,mid+1,R,c);
}
void dfs1(int x){
  	d[x]=d[fa[x]]+1;
	sz[x]=1;
  	for(int y:t[x]){ 
	  	if(y==fa[x]) continue; 
	    fa[y]=x;
		dfs1(y); 
		sz[x]+=sz[y];
	    if(sz[y]>sz[son[x]]||!son[x]) son[x]=y;
  	}
}
void dfs2(int x,int tp) {
  	dfn[x]=++tot;
  	top[x]=tp; 
	if(!son[x]) return ;
	dfs2(son[x],tp);
  	for(int y:t[x]){
  		if(y!=fa[x]&&y!=son[x]) dfs2(y,y);
  	}
}
void work(int x,int y,int c){
  	if(d[x]<d[y]) swap(x,y);
  	if(dfn[y]<=dfn[x]&&dfn[x]<dfn[y]+sz[y]){
	    x=fa[x];
	    while(top[x]!=top[y]){
	      	add1(rt1,1,n,dfn[top[x]],dfn[x],c);
	      	add2(rt2,1,n,dfn[top[x]],dfn[x],c);
	      	x=fa[top[x]];
	    }
    	add1(rt1,1,n,dfn[y]+1,dfn[x],c);
    	add2(rt2,1,n,dfn[y]+1,dfn[x],c);
  	}else{
    	x=fa[x],y=fa[y];
    	while(top[x]!=top[y]){
	      	if(d[top[x]]<d[top[y]]) swap(x,y);
	      	add1(rt1,1,n,dfn[top[x]],dfn[x],c);
	      	add2(rt2,1,n,dfn[top[x]],dfn[x],c);
	      	x=fa[top[x]];
    	}
    	if(dfn[x]<dfn[y]) swap(x,y);
    	add1(rt1,1,n,dfn[y],dfn[x],c);
    	add2(rt2,1,n,dfn[y],dfn[x],c);
  	}
}
bool topo(){//拓扑排序 
  	queue<int>q; 
  	for(int i=1;i<=cnt;i++){
  		if(!deg[i]) q.push(i);
  	}
  	while(!q.empty()){
    	int x=q.front();q.pop();
    	for(int y:e[x]){
			if(!(--deg[y])) q.push(y);
		} 
  	}
  	for(int i=1;i<=cnt;i++){
  		if(deg[i]) return 0;
  	}
  	return 1;
}
int main(){
	ios::sync_with_stdio(false);
	cin.tie(0),cout.tie(0);
	cin>>Q;
	while(Q--){
		cin>>n;
		prework();
		for(int i=2;i<=n;i++){
			int x,y;
			cin>>x>>y;
		  	t[x].push_back(y);
			t[y].push_back(x);
		}
		dfs1(1);
		dfs2(1,1);
		builds(rt1,1,n);
		buildt(rt2,1,n);
		cnt=id;
		cin>>m;
		for(int i=1;i<=m;i++){
		  	int s,t;
		  	cin>>s>>t;
			cnt++;
		  	add(id2[dfn[t]],cnt);
		  	add(cnt,id1[dfn[s]]);
		  	add(cnt,id2[dfn[s]]);
		  	add(id1[dfn[t]],cnt);
		  	work(s,t,cnt);
		}
		if(topo())cout<<"Yes\n";
		else cout<<"No\n";
	}
	return 0;
}