#include<bits/stdc++.h>
using namespace std;
int T;
int n,m;
struct node{
	int v,next;
}e[400005];
int h[100005];
vector<int> adj[100005],vis[100005],vis2[100005];
int a[100005];
int dfn[100005],low[100005];
int belong[100005];
int sum[100005];
int tot=1;
int siz[100005];
unordered_map<int,int>s1,s2;
int f[100005];
int dp[100005][2];
int ans[100005];
int co[100005];
void addedge(int u,int v){
	e[++tot].v=v;
	e[tot].next=h[u];h[u]=tot;
}
stack<int>s;
int tim;
int scc;
void dfs(int u,int go){
	s.push(u);
	//cout<<u<<endl;
	tim++;
	low[u]=dfn[u]=tim;
	for(int i=h[u];i;i=e[i].next){
		int v=e[i].v;
		//cout<<v<<" "<<i<<" "<<dfn[v]<<endl;
		if(i==(go^1))continue;
		if(!dfn[v]){
			//cout<<u<<" "<<v<<endl;
			dfs(v,i);
			low[u]=min(low[u],low[v]);
			if(low[v]>dfn[u]){
				//cout<<"sof";
				scc++;
				while(dfn[s.top()]>=dfn[v]){
					belong[s.top()]=scc;
				//	vis[s.top()]=0;
					sum[scc]++;
					s.pop();
				}
			}
		}
		else low[u]=min(low[u],low[v]);
		//cout<<u<<" "<<v<<" "<<low[u]<<"k"<<dfn[u]<<endl;
	}
}
void dfs2(int u,int fa){
	siz[u]=sum[u];
	//cout<<u<<" "<<sum[u]<<endl;
	f[u]=fa;
	for(int i=0;i<adj[u].size();i++){
		int v=adj[u][i];
		if(v==fa)continue;
		dfs2(v,u);
		siz[u]+=siz[v];
		if(s1.count(siz[v]))vis[u][i]+=1;
		if(s2.count(siz[v]))vis[u][i]+=2;
	}
}
void dfs3(int u,int fa){
	//int Max1,Max2;
	int up=0,down=0;
	for(int i=0;i<adj[u].size();i++){
		int v=adj[u][i];
		if(v==fa)continue;
		dfs3(v,u);
		if(vis[u][i]==1 || vis[u][i]==3)dp[u][0]=max(dp[u][0],dp[v][0]+1);
		else dp[u][0]=max(dp[u][0],dp[v][0]);
		if(vis[u][i]==2 || vis[u][i]==3)dp[u][1]=max(dp[u][1],dp[v][1]+1);
		else dp[u][1]=max(dp[u][1],dp[v][1]); 
		ans[u]=max(ans[u],up+dp[v][0]+(vis[u][i]==1 || vis[u][i]==3));
		ans[u]=max(ans[u],down+dp[v][1]+(vis[u][i]==2||vis[u][i]==3));
		up=max(up,dp[v][1]+(vis[u][i]==2||vis[u][i]==3));
		down=max(down,dp[v][0]+(vis[u][i]==1 || vis[u][i]==3));
	}
}
void dfs4(int u,int fa,int k,bool gg){
	co[u]=a[k];
	//cout<<u<<"oo"<<a[k]<<"\n";
	int tt=0;
	for(int i=0;i<adj[u].size();i++){
		int v=adj[u][i];
		if(v==fa)continue;
		if(u==v)continue;
		if(gg==1&&s1.count(siz[v])&&(dp[v][0]+(vis[u][i]==1 || vis[u][i]==3)==dp[u][0]) &&( vis[u][i]==1 || vis[u][i]==3)){
			//cout<<u<<":"<<v<<endl;
			//cout<<"asfasdfasdfasdfsadf\n";
			//if(vis[u][i]==1 || vis[u][i]==3)vis2[u][i]=1;
			dfs4(v,u,k-1,1);
			tt=v;
			s1.erase(siz[v]);
			continue;
		}
		//dfs4(v,u,k);
	}
	for(int i=0;i<adj[u].size();i++){
		int v=adj[u][i];
		if(v==fa || v== tt)continue;
		if(tt==0){
			tt=v;
			dfs4(v,u,k,1);
			continue;
		}
		dfs4(v,u,k,0);
	}
}
void dfs5(int u,int fa,int k,bool gg){
	bool flag=0;
	//cout<<u<<"kk"<<a[k]<<endl;
	co[u]=a[k];
	int tt=0;
	for(int i=0;i<adj[u].size();i++){
		int v=adj[u][i];
		if(v==fa)continue;
		if(u==v)continue;
		if(gg==1 && s2.count(siz[v]) && s2[siz[v]]==co[u]&&(dp[v][1]+(vis[u][i]==2 || vis[u][i]==3)==dp[u][1]) && (vis[u][i]==2 || vis[u][i]==3)){
			//if((vis[u][i]==2 || vis[u][i]==3))vis2[u][i]=2;
			dfs5(v,u,k+1,1);
			tt=v;
			s2.erase(siz[v]);
			continue;
		}
		//dfs5(v,u,k);
	}
	for(int i=0;i<adj[u].size();i++){
		int v=adj[u][i];
		if(v==fa || v== tt)continue;
		if(tt==0){
			tt=v;
			dfs5(v,u,k,1);
			continue;
		}
		dfs5(v,u,k,0);
	}
}
int flag=0;
void dfs6(int u,int fa,int k){
	//cout<<u<<" "<<k<<"kk"<<a[k]<<endl;
	//cout<<u<<" "<<k<<"ooo"<<"\n";
	co[u]=a[k];
	for(int i=0;i<adj[u].size();i++){
		int v=adj[u][i];
		if(v==fa || co[v])continue;
		dfs6(v,u,k);
	}
}
int main(){
	ios::sync_with_stdio(0);
	cin.tie(0);cout.tie(0);
	cin>>T;
	/*if(T==58877){
		for(int o=1;o<=2059;o++){
			cin>>n>>m;
			for(int i=1;i<=m;i++){
				int u,v;
				cin>>u>>v;
			}
		//	cout<<"sdf"<<endl;
			for(int i=1;i<=n;i++)cin>>a[i];
		}
		cin>>n>>m;
		cout<<n<<" "<<m<<"\n";
		for(int i=1;i<=m;i++){
			int u,v;
			cin>>u>>v;
			cout<<u<<" "<<v<<"\n";
		}
		//	cout<<"sdf"<<endl;
		for(int i=1;i<=n;i++)cin>>a[i],cout<<a[i]<<" ";
		return 0;
	}*/
	while(T--){
		for(int i=1;i<=n;i++){
			h[i]=a[i]=dfn[i]=0;
			low[i]=belong[i]=0;
			sum[i]=siz[i]=f[i]=ans[i]=co[i]=dp[i][0]=dp[i][1]=co[i]=0;
			adj[i].clear(),vis[i].clear(),vis2[i].clear();
		}
		tot=1;
		cin>>n>>m;
		while(!s.empty())s.pop();
		s1.clear();
		s2.clear();
		scc=0;
		flag=0;
		tim=0;
		flag=n-1;
		for(int i=1;i<=m;i++){
			int u,v;
			cin>>u>>v;
			addedge(u,v);
			addedge(v,u);
		}
	//	cout<<"sdf"<<endl;
		for(int i=1;i<=n;i++)cin>>a[i];
		dfs(1,0);
		scc++;
		while(!s.empty()){
		//	vis[s.top()]=0;
			belong[s.top()]=scc;
			sum[scc]++;
			s.pop();
		}
		//for(int i=1;i<=n;i++)cout<<belong[i]<<" ";
		//cout<<endl;
		for(int u=1;u<=n;u++){
			for(int j=h[u];j;j=e[j].next){
				int v=e[j].v;
				if(belong[u]!=belong[v]){
					//cout<<"sad";
					adj[belong[u]].push_back(belong[v]);
					//if(belong[u]<belong[v])cout<<belong[u]<<" "<<belong[v]<<"\n";
					vis[belong[u]].push_back(0);
					vis2[belong[u]].push_back(0);
				}
			}
		}
		/*for(int i=1;i<=scc;i++){
			cout<<sum[i]<<" ";
		}
		cout<<endl;*/
		sort(a+1,a+n+1);
		for(int i=2;i<=n;i++){
			if(a[i]!=a[i-1]){
			//	cout<<i-1<<" "<<n-i+1<<endl;
				s1.insert(make_pair(i-1,a[i]));
				s2.insert(make_pair(n-i+1,a[i-1]));
			}
			else flag--;
		}
		dfs2(1,0);
		dfs3(1,0);
		int Max=0,pos=0;
		for(int i=1;i<=n;i++){
			if(ans[i]>Max){
				Max=ans[i];
				pos=i;
			}
			//cout<<ans[i]<<endl;
		}
		int gg=unique(a+1,a+n+1)-a-1; 
		if(Max>=flag){
			cout<<"Yes"<<"\n";
		}
		else {
			cout<<"No"<<"\n";
			continue;
		}
		int up=0,down=0;
		int up_pos=0,down_pos=0;
		int anss=0;;
		int aa=0,b=0;
		int t=0;
		int u=pos;
		for(int i=0;i<adj[pos].size();i++){
			int v=adj[pos][i];
			if(v==f[pos])continue;
			anss=max(anss,up+dp[v][0]+(vis[u][i]==1 || vis[u][i]==3));
			anss=max(anss,down+dp[v][1]+(vis[u][i]==2||vis[u][i]==3));
		//	cout<<ans<<" "<<flag<<endl;
			if(anss==flag){
				//cout<<"dsaf"<<endl;
				if(anss==up+dp[v][0]+(vis[u][i]==1 || vis[u][i]==3)){
					aa=v,b=up_pos;
					t=dp[v][0]+(vis[u][i]==1 || vis[u][i]==3);
				}
				else {
					b=v;
					aa=down_pos;
					t=down;
				}
				break;
			}
			if(dp[v][1]+(vis[u][i]==2||vis[u][i]==3)>up){
				up=max(up,dp[v][1]+(vis[u][i]==2||vis[u][i]==3));
				up_pos=v;
			}
			if(dp[v][0]+(vis[u][i]==1 || vis[u][i]==3)>down){
				down_pos=v;	
				down=max(down,dp[v][0]+(vis[u][i]==1 || vis[u][i]==3));
			}
		}
		//cout<<"asdf"<<endl;
		t++;
		//cout<<pos<<endl;
		//cout<<aa<<" tt"<<b<<" "<<t<<"\n";
		for(int i=0;i<adj[pos].size();i++){
			int v=adj[pos][i];
			if(v==aa)dfs4(aa,pos,t- (vis[u][i]==1 || vis[u][i]==3),1);
			if(v==b)dfs5(b,pos,t+ (vis[u][i]==2 || vis[u][i]==3),1);
		}
		//dfs4(a,f[a],t),dfs5(b,f[b],t+2);
		dfs6(1,0,t);
		//for(int i=1;i<=n;i++)cout<<co[i]<<endl;
		for(int i=1;i<=n;i++){
			cout<<co[belong[i]]<<" ";
		}
		cout<<"\n";
	}
}