#include <algorithm>
#include <cstdint>
#include <iostream>
#include <random>
#include <unordered_map>
#include <vector>

using namespace std;

constexpr uint32_t N = 1e6, mod = -37u / 2;

vector<int> todo, q;
uint64_t head[N], tail[N];
unordered_map<uint64_t, vector<int>> byh;

uint32_t pw[N], remap[26][2], h[2][N + 1][2];

int nx[N];
unordered_map<uint64_t, vector<uint64_t>> to;

uint64_t hget(int x, int l, int n) { return (h[x][l][0] * (uint64_t)pw[n] + mod - h[x][l + n][0]) % mod << 32 |
                                            (h[x][l][1] * (uint64_t)pw[n] + mod - h[x][l + n][1]) % mod; }

bool extend(int v) {
  int skipped = 0;
  for (auto it = to[tail[v]].end(); it != to[tail[v]].begin(); ) {
    auto x = *--it;
    while (byh[x].size()) {
      auto y = byh[x].back(); byh[x].pop_back();
      to[tail[v]].erase(it);
      todo.push_back(y);
      return extend(nx[v] = y);
    }
    if (skipped++) return 0;
  }
  return 1;
}

void solve() {
  int n, m; cin >> n >> m;
  string s[2];
  for (int i = 0; i < 2 * n; ++i) { string c; cin >> c, s[i >= n] += c; }
  for (int x: {0, 1}) for (int y: {0, 1})
  for (int i = 0, p = 0; i < n * m; ++i) h[x][i + 1][y] = p = (p * 2ull + remap[s[x][i] - 'a'][y]) % mod;
  for (int d = 0; d < m; ++d) {
    to = {};
    byh = {};
    for (int i = 0; i < n; ++i) head[i] = hget(0, i * m, d), tail[i] = hget(0, i * m + d, m - d), byh[head[i]].push_back(i), nx[i] = -1;
    for (int i = 0; i < n; ++i) to[hget(1, i * m, m - d)].push_back(hget(1, i * m + m - d, d));
    byh[head[0]].erase(find(byh[head[0]].begin(), byh[head[0]].end(), 0));
    todo = {{}};
    if (!extend(0)) todo = {};
    else if (int x = todo.back(); to[tail[x]].empty() || to[tail[x]][0] != head[0]) todo = {};//, cout << ":( " << x << '\n';
    else to[tail[x]] = {}, nx[x] = 0;
    for (int i = 0; i < todo.size(); ++i) {
      int x = todo[i], y = nx[x];
      if (!extend(x)) { todo = {}; break; }
      if (nx[x] == y) continue;
      --i;
      x = todo.back();
      if (to[tail[x]].empty() || to[tail[x]][0] != head[y]) { todo = {}; break; }
      to[tail[x]] = {}, nx[x] = y;
    }
    if (todo.size() != n) continue;
    byh = {};
    for (int i = 0; i < n; ++i) byh[hget(1, i * m, m - d) * 2 + hget(1, i * m + m - d, d)].push_back(i);
    todo = q = {};
    for (int i = 0; todo.push_back(i), i = nx[i]; ) ;
    for (int i = 0; i < n; ++i) {
      auto h = tail[todo[i]] * 2 + head[todo[i + (i != n - 1 ?: 1 - n)]];
      q.push_back(byh[h].back()), byh[h].pop_back();
    }
    for (auto x: todo) cout << x + 1 << ' '; cout << '\n';
    for (auto x: q) cout << x + 1 << ' '; cout << '\n';
    return;
  }
  cout << "-1\n";
}

int main() {
  cin.tie(0)->sync_with_stdio(0);
  generate(*remap, *end(remap), mt19937_64{163467});
  for (int i = 0, p = 1; i < N; ++i) pw[i] = p, (p += p) %= mod;
  for (int tc = (cin >> tc, tc); tc--; ) solve();
}