#include<bits/stdc++.h>
using namespace std;using ll=long long;using ld=long double;using pii=pair<int,int>;
const int M=998244353,N=1e6+10;const ll inf=1e17;const double eps=1e-10;
ll b[2][N];
struct Hsr
{
	ll res[2];
	bool operator==(Hsr A){return A.res[0]==res[0]&&A.res[1]==res[1];}
};
bool operator<(Hsr A,Hsr B)
{
	if(A.res[0]^B.res[0])return A.res[0]<B.res[0];
	return A.res[1]<B.res[1];
};
struct Hash{
	string s;int n;vector<vector<ll>>l;int mod[2]={998244353,int(1e9+9)},fac[2]={9527,3960};
	Hsr GhL(int L,int R){if(L>R)return Hsr{{0,0}};Hsr C;for(int e=0;e<2;e++)C.res[e]=
		(l[e][R]-l[e][L-1]*b[e][R-L+1]%mod[e]+mod[e])%mod[e];return C;}
	int tran(char c){return c-'a'+1;}
	void ini(){cin>>s;n=s.length();l.resize(2,vector<ll>(n+2));
		for(int e=0;e<2;e++)l[e][0]=0;
		for(int e=0;e<2;e++)for(int i=1;i<=n;i++){l[e][i]=(l[e][i-1]*fac[e]+tran(s[i-1]))%mod[e];}}
};
struct S2
{	
	int n,m,siz=0;
	vector<int>gd,cov,l,r;
	vector<vector<int>>es;vector<int>ec;
	bool ini(int _n,vector<pii>vr)
	{
		n=_n;m=vr.size();
		es.resize(n+1);ec.resize(n+1,0);
		gd.resize(m+1,0);cov=l=r=gd;
		for(int i=1,u,v;i<=m;i++)
		{
			u=vr[i-1].first,v=vr[i-1].second;
			es[u].push_back(i);ec[u]++;gd[i]=u^v;
		}	
		return solve();
	}
	int tar,thd;
	void spr2(int p,int x)
	{
		cov[x]=1;siz++;
		if(thd==p){r[x]=tar,l[tar]=x;tar=thd=0; }
		while(ec[p])
		{
			int k=es[p][--ec[p]];
			if(k==x&&gd[k])continue;
			if(!cov[k])
			{
				int y=r[x];r[x]=k,l[k]=x;
				if(y)l[y]=0,tar=y,thd=p;
				spr2(gd[k]^p,k);
			}			
		}
	}
	vector<int>R;
	bool solve()
	{	
		int x=es[1][0];tar=x;thd=1;
		spr2(gd[x]^1,x);if(siz!=m)return 0;		
		R.push_back(1);int p=1,c=siz;
		while(c--){p=gd[x]^p;R.push_back(p);x=r[x];}
		return 1;
	}
};
struct S
{	
	int n,m;vector<Hash>A,B;
	void ini()
	{
		cin>>n>>m;
		A.resize(n+1);B.resize(n+1);
		for(int i=1;i<=n;i++)A[i].ini();
		for(int i=1;i<=n;i++)B[i].ini();
	}
	bool Sol(int p)
	{
		S2 C;vector<pii>es;vector<int>cnt(n*6+1,0);int cc=n*2;
		{
			map<Hsr,int>mp;
			for(int i=1;i<=n;i++)
			{
				auto a=A[i].GhL(1,p);int id=mp[a];
				if(!id)mp[a]=id=++cc;es.push_back({i,id});cnt[id]++;
			}
			for(int i=1;i<=n;i++)
			{
				auto a=B[i].GhL(m-p+1,m);int id=mp[a];
				if(!id)mp[a]=id=++cc;es.push_back({id,i+n});cnt[id]--;
			}
		}
		{
			map<Hsr,int>mp;
			for(int i=1;i<=n;i++)
			{
				auto a=A[i].GhL(p+1,m);int id=mp[a];
				if(!id)mp[a]=id=++cc;es.push_back({id,i});cnt[id]++;
			}
			for(int i=1;i<=n;i++)
			{
				auto a=B[i].GhL(1,m-p);int id=mp[a];
				if(!id)mp[a]=id=++cc;es.push_back({i+n,id});cnt[id]--;
			}			
		}
		for(int i=1;i<=cc;i++)if(cnt[i])return 0;
		if(!C.ini(cc,es))return 0;
		auto&R=C.R;int siz=R.size()-1;vector<int>ra,rb;
		for(int i=siz;i;i--)
		{
			int id=R[i];if(id>n*2)continue;
			if(id<=n)ra.push_back(id);else rb.push_back(id);
		}
		for(int i=0;i<n;i++)cout<<ra[i]<<" \n"[i==n-1];
		for(int i=0;i<n;i++)cout<<rb[i]-n<<" \n"[i==n-1];		
		return 1;
	}
	void solve()
	{
		for(int i=0;i<m;i++)if(Sol(i))return;cout<<"-1\n";
	}
};
void precal()
{
	int mod[2]={998244353,int(1e9+9)},fac[2]={9527,3960},P=N-5;
	for(int e=0;e<2;e++)b[e][0]=1;
	for(int e=0;e<2;e++)for(int i=1;i<=P;i++)b[e][i]=b[e][i-1]*fac[e]%mod[e];	
}
int main()
{
	//freopen("1.in","r",stdin);
	//cout<<fixed<<setprecision(12);
	ios::sync_with_stdio(0);cin.tie(0);precal();
	int t=1;cin>>t;
	clock_t a=clock();
	while(t--){S SS;SS.ini();SS.solve();}
}