#include <bits/stdc++.h>
#ifdef LOCAL
#include <debug.hpp>
#else
#define debug(...) void(0)
#endif

struct LowLink {
    int n, m = 0, t = 0, b = 0;
    std::vector<std::vector<int>> G;
    std::vector<std::pair<int, int>> es;
    std::vector<int> ord, low, tecc, bcc, tmp;
    std::vector<bool> is_articulation, is_bridge;

    LowLink(int n) : n(n), G(n), ord(n, -1), low(n), tecc(n, -1), is_articulation(n, false) {}

    void add_edge(int u, int v) {
        assert(0 <= u and u < n);
        assert(0 <= v and v < n);
        G[u].emplace_back(m);
        G[v].emplace_back(m);
        es.emplace_back(u, v);
        is_bridge.emplace_back(false);
        bcc.emplace_back();
        m++;
    }

    void build() {
        for (int i = 0; i < n; i++) {
            if (ord[i] != -1) continue;
            dfs1(i, 0, -1);
            dfs2(i, -1);
        }
    }

    std::pair<int, int> operator[](int k) const { return es[k]; }

private:
    void dfs1(int v, int k, int pre) {
        ord[v] = k++, low[v] = ord[v];
        int cnt = 0;
        for (int& e : G[v]) {
            if (e == pre) continue;
            int u = es[e].first ^ es[e].second ^ v;
            if (ord[u] == -1) {
                cnt++;
                dfs1(u, k, e);
                low[v] = std::min(low[v], low[u]);
                if (pre != -1 and ord[v] <= low[u]) is_articulation[v] = true;
                if (ord[v] < low[u]) is_bridge[e] = true;
            } else
                low[v] = std::min(low[v], ord[u]);
        }
        if (pre == -1 and cnt > 1) is_articulation[v] = true;
    }

    void dfs2(int v, int pre) {
        if (pre == -1) tecc[v] = t++;
        for (int& e : G[v]) {
            if (e == pre) continue;
            int u = es[e].first ^ es[e].second ^ v;
            if (tecc[u] == -1 or ord[u] < ord[v]) tmp.emplace_back(e);
            if (tecc[u] >= 0) continue;
            if (ord[v] >= low[u])
                tecc[u] = tecc[v];
            else
                tecc[u] = t++;
            dfs2(u, e);
            if (ord[v] <= low[u]) {
                while (true) {
                    int ne = tmp.back();
                    tmp.pop_back();
                    bcc[ne] = b;
                    if (ne == e) break;
                }
                b++;
            }
        }
    }
};

using namespace std;

typedef long long ll;
#define all(x) begin(x), end(x)
constexpr int INF = (1 << 30) - 1;
constexpr long long IINF = (1LL << 60) - 1;
constexpr int dx[4] = {1, 0, -1, 0}, dy[4] = {0, 1, 0, -1};

template <class T> istream& operator>>(istream& is, vector<T>& v) {
    for (auto& x : v) is >> x;
    return is;
}

template <class T> ostream& operator<<(ostream& os, const vector<T>& v) {
    auto sep = "";
    for (const auto& x : v) os << exchange(sep, " ") << x;
    return os;
}

template <class T, class U = T> bool chmin(T& x, U&& y) { return y < x and (x = forward<U>(y), true); }

template <class T, class U = T> bool chmax(T& x, U&& y) { return x < y and (x = forward<U>(y), true); }

template <class T> void mkuni(vector<T>& v) {
    sort(begin(v), end(v));
    v.erase(unique(begin(v), end(v)), end(v));
}

template <class T> int lwb(const vector<T>& v, const T& x) { return lower_bound(begin(v), end(v), x) - begin(v); }

void solve() {
    int n, m;
    cin >> n >> m;
    LowLink LL(n);
    for (int i = 0; i < m; i++) {
        int u, v;
        cin >> u >> v;
        LL.add_edge(--u, --v);
    }
    vector<int> a(n);
    cin >> a;

    sort(all(a));
    if (a.front() == a.back()) {
        cout << "Yes\n";
        cout << a << '\n';
        return;
    }
    LL.build();
    int N = LL.t;
    auto& tecc = LL.tecc;
    auto& is_bridge = LL.is_bridge;
    vector<vector<int>> G(N), group(N);
    for (int i = 0; i < m; i++) {
        if (not is_bridge[i]) continue;
        auto [u, v] = LL[i];
        u = tecc[u], v = tecc[v];
        G[u].emplace_back(v);
        G[v].emplace_back(u);
    }
    for (int i = 0; i < n; i++) group[tecc[i]].emplace_back(i);

    vector<int> sub(N);
    auto dfs_sz = [&](auto self, int v, int p) -> int {
        sub[v] = group[v].size();
        for (int& u : G[v]) {
            if (u == p) continue;
            sub[v] += self(self, u, v);
        }
        return sub[v];
    };
    auto dfs_search = [&](auto self, int v, int p, int mid) -> int {
        for (int& u : G[v]) {
            if (u == p) continue;
            if (sub[u] > mid) return self(self, u, v, mid);
        }
        return v;
    };
    int c = dfs_search(dfs_search, 0, -1, dfs_sz(dfs_sz, 0, -1) / 2);
    vector<int> cnt(n, 0);
    for (int& x : a) cnt[x - 1]++;

    int A = -1, B = -1;
    for (int i = 0, sum = 0; i < n; i++) {
        sum += cnt[i];
        if (A == -1 and sum >= (n + 1) / 2) A = i;
        if (B == -1 and (~n & 1) and sum >= (n + 2) / 2) B = i;
    }

    auto f = [&](int val) -> vector<int> {
        if (val == -1) return {};
        vector<vector<pair<int, int>>> route(2);
        for (int i = 0, sum = 0; i < val; i++) {
            if (cnt[i] == 0) continue;
            sum += cnt[i];
            route[0].emplace_back(sum, i);
        }
        for (int i = n - 1, sum = 0; i > val; i--) {
            if (cnt[i] == 0) continue;
            sum += cnt[i];
            route[1].emplace_back(sum, i);
        }
        vector nxt(N, vector<int>(2, -1));
        vector maxi(N, vector<int>(2, 0));
        int nc = -1;
        vector<int> par(N, -1);
        auto dfs = [&](auto self, int v, int p) -> void {
            par[v] = p;
            sub[v] = group[v].size();
            vector<vector<int>> OK(2);
            for (int& u : G[v]) {
                if (u == p) continue;
                self(self, u, v);
                for (int i = 0; i < 2; i++) {
                    if (nxt[v][i] == -1) nxt[v][i] = u;
                    if (chmax(maxi[v][i], maxi[u][i])) nxt[v][i] = u;
                    if (maxi[u][i] == int(route[i].size())) OK[i].emplace_back(u);
                }
                sub[v] += sub[u];
            }
            for (int& e : OK[0]) {
                if (nc != -1) break;
                for (int& f : OK[1]) {
                    if (e != f) {
                        nc = v;
                        break;
                    }
                }
            }
            for (int i = 0; i < 2; i++) {
                if (maxi[v][i] < int(route[i].size()) and sub[v] == route[i][maxi[v][i]].first) maxi[v][i]++;
                if (maxi[v][i] < int(route[i].size()) and sub[v] > route[i][maxi[v][i]].first) maxi[v][i] = -1;
            }
        };
        vector<vector<int>> ok(2);
        for (int& ch : G[c]) {
            dfs(dfs, ch, c);
            for (int i = 0; i < 2; i++) {
                if (maxi[ch][i] == int(route[i].size())) ok[i].emplace_back(ch);
            }
        }
        if (nc != -1) {
            c = nc;
            for (int i = 0; i < 2; i++) ok[i].clear();
            for (int& ch : G[c]) {
                if (ch == par[c]) continue;
                for (int i = 0; i < 2; i++) {
                    if (maxi[ch][i] == int(route[i].size())) ok[i].emplace_back(ch);
                }
            }
        }
        // debug(c, val);
        // debug(route[0], route[1]);
        // debug(ok[0], ok[1]);
        int p = -1, q = -1;
        if (route[0].empty()) {
            if (not ok[1].empty()) q = ok[1][0];
        } else if (route[1].empty()) {
            if (not ok[0].empty()) p = ok[0][0];
        } else {
            for (int& u : ok[0]) {
                for (int& v : ok[1]) {
                    if (u != v) {
                        p = u, q = v;
                        break;
                    }
                }
                if (p != -1) break;
            }
        }
        if (p == -1 and q == -1) return {};
        // debug(val, p, q);
        vector<int> ans(n, -1);
        auto dfs2 = [&](auto self, int v, int p, int idx, bool on, int& cur,
                        const vector<pair<int, int>>& path) -> vector<int> {
            if (v == -1) return {};
            vector<int> tmp;
            for (int& vv : group[v]) tmp.emplace_back(vv);
            for (int& u : G[v]) {
                if (u == p) continue;
                auto ch = self(self, u, v, idx, on & (u == nxt[v][idx]), cur, path);
                if (tmp.size() < ch.size()) swap(tmp, ch);
                for (int& vv : ch) tmp.emplace_back(vv);
            }
            if (on and cur < int(path.size()) and sub[v] == path[cur].first) {
                assert(int(tmp.size()) == cnt[path[cur].second]);
                for (int& vv : tmp) ans[vv] = path[cur].second;
                tmp.clear();
                cur++;
            }
            return tmp;
        };
        int init = 0;
        dfs2(dfs2, p, c, 0, true, init, route[0]);
        assert(init == int(route[0].size()));
        init = 0;
        dfs2(dfs2, q, c, 1, true, init, route[1]);
        assert(init == int(route[1].size()));
        for (int i = 0; i < n; i++) {
            if (ans[i] == -1) {
                ans[i] = val;
            }
        }
        return ans;
    };

    for (int val : {A, B}) {
        auto ans = f(val);
        if (ans.empty()) continue;
        cout << "Yes\n";
        for (int i = 0; i < n; i++) cout << ans[i] + 1 << (i + 1 == n ? '\n' : ' ');
        return;
    }
    cout << "No\n";
}

int main() {
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    int T;
    cin >> T;
    for (; T--;) solve();
    return 0;
}