#include<bits/stdc++.h>
using namespace std;
#define lli long long int
#define pb push_back
const lli mod=1e9+7;
const lli mod1=998244353;
#define fi first
#define se second
//  #include <boost/multiprecision/cpp_int.hpp> 
//  using namespace boost::multiprecision; 
//  #include <ext/pb_ds/assoc_container.hpp> 
//  #include <ext/pb_ds/tree_policy.hpp> 
//  using namespace __gnu_pbds; 
//  #define ordered_set tree<lli, null_type,less<lli>, rb_tree_tag,tree_order_statistics_node_update> 
/*
lli power(lli a, lli b) {
    lli res = 1;
    while (b > 0) {
        if (b & 1)
            res = res * a;
        a = a * a;
        b >>= 1;
    }
    return res;
}
lli powermod(lli a, lli b) 
{
    lli res = 1;
    while (b > 0) {
        if (b & 1)
            res = ((res%mod)*(a%mod))%mod;
        a = (a * a)%mod;
        b >>= 1;
    }
    return res;
}*/
const lli N=100001;
bool vis[N]={false};
lli cnt[N]={0},cnt1[N]={0};
vector<lli> v[N];
lli po[N]={0};
int main()
{
    //ios_base::sync_with_stdio(false);
    //cin.tie(NULL);
    lli T;T=1;
    cin>>T;
    while(T--)
    {
        lli n,m,i,j;
        cin>>n>>m;
        lli a[n],b[n],c[m];
        for(i=0;i<n;i++)
        {
            cin>>a[i];
        }
        for(i=0;i<n;i++)
        {
            cin>>b[i];
            vis[b[i]]=true;
        }
        for(i=0;i<m;i++)
        {
            cin>>c[i];
            cnt[c[i]]--;
        }
        if(vis[c[m-1]]==false)
        {
            cout<<"NO\n";
            for(i=0;i<n;i++)
                vis[b[i]]=false;
            for(i=0;i<=n;i++)
                cnt[i]=cnt1[i]=po[i]=0;
            for(i=0;i<n;i++)
            {
                v[b[i]].clear();
            }
            continue;
        }
        for(i=0;i<n;i++)
        {
            if(a[i]!=b[i])
            {   
                cnt[b[i]]++;
                v[b[i]].pb(i+1);
            }
        }
        lli flag=0;
        for(i=1;i<=n;i++)
        {
            if(cnt[i]>0)
                flag=1;
        }
        if(flag)
        {
            cout<<"NO\n";
            for(i=0;i<n;i++)
                vis[b[i]]=false;
            for(i=0;i<=n;i++)
                cnt[i]=cnt1[i]=po[i]=0;
            for(i=0;i<n;i++)
            {
                v[b[i]].clear();
            }
            continue;
        }
        lli ans=-1;
        cout<<"YES\n";
        for(i=0;i<n;i++)
        {
            if(b[i]==c[m-1])
                {ans=i+1;break;}
        }
        for(i=0;i<m-1;i++)
        {
            if(v[c[i]].empty())
            {
                cout<<ans<<" ";
            }
            else
            {
                cout<<v[c[i]].back()<<" ";
                v[c[i]].pop_back();
            }
        }
        cout<<ans<<"\n";
        for(i=0;i<n;i++)
            vis[b[i]]=false;
        for(i=0;i<=n;i++)
            cnt[i]=cnt1[i]=po[i]=0;
        for(i=0;i<n;i++)
        {
            v[b[i]].clear();
        }
    }
    return 0;
}