//#define dxx
#ifdef dxx
#define dbg(...) fprintf(stderr,__VA_ARGS__)
#define dex(a) dbg(#a"=%lld onL%d infun %s\n",a,__LINE__,__FUNCTION__)
#include<cstdlib>
#define pause sys##tem("pause")
#define _GLIBCXX_DEBUG
#endif

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#include<set>
#include<map>

using ll=long long;
using std::max;
using std::min;
using std::abs;
template<class T> void cmax(T&a,T b){a=max(a,b);}
template<class T> void cmin(T&a,T b){a=min(a,b);}
template<class T> T sqr(T a){return a*a;}

namespace gph{
    struct edge{
        int t,id;
    };
    std::vector<edge> G[4000006];
    int st[4000006],ind[4000006],oud[4000006];
    void ade(int s,int t,int id){
        G[s].push_back({t,id});
        ++ind[t];
        ++oud[s];
        //printf("%d %d %d\n",s,t,id);
    }
    std::vector<int> ans;
    bool evis[4000006];
    void dfs(int x,int eid){
        for(int e=st[x];e<(int)G[x].size();e=max(e+1,st[x])){
            int t=G[x][e].t;
            if(!evis[G[x][e].id]){
                evis[G[x][e].id]=1;
                st[x]=e+1;
                dfs(t,G[x][e].id);
            }
        }
        if(eid) ans.push_back(eid);
    }bool geteuler(int N,int ecnt){
        for(int i=1;i<=N;++i) if(ind[i]!=oud[i]) return 0;
        ans.clear();
        dfs(1,0);
        if((int)ans.size()==ecnt) return 1;
        return 0;
    }void clear(int N,int ecnt){
        for(int i=1;i<=N;++i) G[i].clear();
        memset(evis+1,0,sizeof*evis*ecnt);
        memset(st+1,0,sizeof*st*N);
        memset(ind+1,0,sizeof*ind*N);
        memset(oud+1,0,sizeof*oud*N);
    }
}

namespace xm{
    const ll mod=(ll)1e17+13,mod2=(ll)1e16+13;
    const int bas=2347,bas2=127;

    std::map<std::pair<ll,ll>,int> mp;
    std::set<std::pair<std::pair<ll,ll>,int> > hshmst;
    std::vector<ll> ha[1000006],hb[1000006];
    std::vector<ll> ha2[1000006],hb2[1000006];
    ll bpw[1000006],bpw2[1000006];
    int oneans[1000006];
    char tsz[1000006];

    void getbpw(){
        bpw[0]=1;
        for(int i=1;i<1000006;++i) bpw[i]=(__int128)bpw[i-1]*bas%mod;
        bpw2[0]=1;
        for(int i=1;i<1000006;++i) bpw2[i]=(__int128)bpw2[i-1]*bas2%mod;
    }void gethash(std::vector<ll>&hsh,std::vector<ll>&hsh2){
        hsh.clear();
        for(int i=1;tsz[i];++i)
            hsh.push_back(i==1?tsz[i]-'a'+1:
                    ((__int128)hsh.back()*bas%mod+tsz[i]-'a'+1)%mod);
        hsh2.clear();
        for(int i=1;tsz[i];++i)
            hsh2.push_back(i==1?tsz[i]-'a'+1:
                    ((__int128)hsh2.back()*bas2%mod2+tsz[i]-'a'+1)%mod2);
        //for(ll x:hsh) printf("%lld ",x);
        //putchar(10);
    }std::pair<ll,ll> getrange(std::vector<ll>&hsh,std::vector<ll>&hsh2,int l,int r){
        --l;
        --r;
        //printf("[%d %d]\n",l,r);
        return std::make_pair(ll(l?(hsh[r]-(__int128)hsh[l-1]*bpw[r-l+1]%mod+mod)%mod:
            hsh[r]),ll(l?(hsh2[r]-(__int128)hsh2[l-1]*bpw2[r-l+1]%mod2+mod2)%mod2:
            hsh2[r]));
    }void _(){
        int N,M;

        scanf("%d%d",&N,&M);
        for(int i=1;i<=N;++i){
            scanf("%s",tsz+1);
            gethash(ha[i],ha2[i]);
        }
        for(int i=1;i<=N;++i){
            scanf("%s",tsz+1);
            gethash(hb[i],hb2[i]);
        }

        for(int d=1;d<M;++d){
            int vcnt=0;

            mp.clear();
            for(int i=1;i<=N;++i){
                auto pref=getrange(ha[i],ha2[i],1,d),suff=getrange(ha[i],ha2[i],d+1,M);
                if(!mp[pref]) mp[pref]=++vcnt;
                if(!mp[suff]) mp[suff]=++vcnt;
                gph::ade(mp[pref],mp[suff],i);
                //printf("%d %d\n",pref,suff);
            }
            for(int i=1;i<=N;++i){
                auto pref=getrange(hb[i],hb2[i],1,M-d),suff=getrange(hb[i],hb2[i],M-d+1,M);
                if(!mp[pref]) mp[pref]=++vcnt;
                if(!mp[suff]) mp[suff]=++vcnt;
                gph::ade(mp[pref],mp[suff],i+N);
                //printf("%d %d\n",pref,suff);
            }
            if(gph::geteuler(vcnt,N*2)){
                std::vector<int> ans1,ans2;
                for(int x:gph::ans)
                    if(x<=N) ans1.push_back(x);
                    else ans2.push_back(x-N);
                for(auto it=ans1.rbegin();it!=ans1.rend();++it)
                    printf("%d ",*it);
                putchar(10);
                for(auto it=ans2.rbegin();it!=ans2.rend();++it)
                    printf("%d ",*it);
                putchar(10);
                gph::clear(vcnt,N*2);
                return;
            }
            gph::clear(vcnt,N*2);
        }

        hshmst.clear();
        for(int i=1;i<=N;++i)
            hshmst.emplace(std::make_pair(ha[i].back(),ha2[i].back()),i);
        for(int i=1;i<=N;++i){
            auto it=hshmst.lower_bound(std::make_pair(std::make_pair(hb[i].back(),hb2[i].back()),0));
            if(it!=hshmst.end()&&it->first==std::make_pair(hb[i].back(),hb2[i].back())){
                oneans[i]=it->second;
                hshmst.erase(it);
            }else{
                puts("-1");
                return;
            }
        }
        for(int i=1;i<=N;++i) printf("%d ",i);
        putchar(10);
        for(int i=1;i<=N;++i) printf("%d ",oneans[i]);
        putchar(10);
    }
}

int main(){
    int t;
    xm::getbpw();
    scanf("%d",&t);
    while(t--) xm::_();
    return 0;
}