/*
_/_/_/_/    _/_/_/_/_/  _/_/_/
_/      _/      _/    _/      _/
_/      _/      _/    _/      _/
_/      _/      _/    _/      _/
_/      _/      _/    _/  _/  _/
_/      _/  _/  _/    _/    _/_/
_/_/_/_/      _/_/     _/_/_/_/_/

_/_/_/_/    _/    _/  _/      _/
_/      _/   _/  _/   _/_/  _/_/
_/      _/    _/_/    _/ _/_/ _/
_/      _/     _/     _/  _/  _/
_/      _/    _/_/    _/      _/
_/      _/   _/  _/   _/      _/
_/_/_/_/    _/    _/  _/      _/

_/_/_/_/_/ _/_/_/_/_/ _/_/_/_/_/
    _/         _/     _/
    _/         _/     _/
    _/         _/     _/_/_/_/
    _/         _/     _/
    _/         _/     _/
    _/     _/_/_/_/_/ _/_/_/_/_/

_/_/_/_/_/ _/_/_/_/_/ _/_/_/_/_/
    _/         _/     _/
    _/         _/     _/
    _/         _/     _/_/_/_/
    _/         _/     _/
    _/         _/     _/
    _/     _/_/_/_/_/ _/_/_/_/_/
*/
#include <cassert>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <ctime>
#include <cctype>

#include <algorithm>
#include <random>
#include <bitset>
#include <queue>
#include <functional>
#include <set>
#include <map>
#include <vector>
#include <chrono>
#include <iostream>
#include <iomanip>
#include <limits>
#include <numeric>

#define LOG(FMT...) fprintf(stderr, FMT)

using namespace std;

typedef long long ll;
typedef unsigned long long ull;

// mt19937 rng(chrono::steady_clock::now().time_since_epoch().count());

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

template<class T>
ostream &operator<<(ostream &os, const vector<T> &v) {
  if (!v.empty()) {
    os << v.front();
    for (int i = 1; i < v.size(); ++i)
      os << ' ' << v[i];
  }
  return os;
}

const int P = 1000000007;

int norm(int x) { return x >= P ? (x - P) : x; }

void add(int &x, int y) { if ((x += y) >= P) x -= P; }

void sub(int &x, int y) { if ((x -= y) < 0) x += P; }

void exGcd(int a, int b, int &x, int &y) {
  if (!b) {
    x = 1;
    y = 0;
    return;
  }
  exGcd(b, a % b, y, x);
  y -= a / b * x;
}

int inv(int a) {
  int x, y;
  exGcd(a, P, x, y);
  return norm(x + P);
}

const int N = 150;

int fac[N];
char g[N][N];
char s[N];
int in[N], out[N];
int mat[N][N];

int det(int n) {
  int ret = 1;
  for (int i = 1; i <= n; ++i) {
    int p = i;
    while (p <= n && !mat[p][i])
      ++p;
    if (p > n)
      return 0;
    if (p != i) {
      for (int j = i; j <= n; ++j)
        swap(mat[p][j], mat[i][j]);
      ret = norm(P - ret);
    }
    ret = ret * (ll) mat[i][i] % P;
    int nv = inv(mat[i][i]);
    for (int j = i; j <= n; ++j)
      mat[i][j] = mat[i][j] * (ll) nv % P;
    for (int j = p + 1; j <= n; ++j)
      for (int k = n; k >= i; --k)
        mat[j][k] = norm(P + mat[j][k] - mat[j][i] * (ll) mat[i][k] % P);
  }
  return ret;
}

int main() {
#ifdef ELEGIA
  freopen("test.in", "r", stdin);
  int nol_cl = clock();
#endif
  ios::sync_with_stdio(false);
  cin.tie(nullptr);

  fac[0] = 1;
  for (int i = 1; i != N; ++i) fac[i] = fac[i - 1] * (ull) i % P;
  int T;
  cin >> T;
  while (T--) {
    int n, k;
    cin >> n >> k;
    cin >> (s + 1);
    for (int i = 1; i <= n; ++i) cin >> (g[i] + 1);
    memset(in, 0, sizeof(in));
    memset(out, 0, sizeof(out));
    memset(mat, 0, sizeof(mat));
    auto adde = [&](int i, int j) {
      ++in[j];
      ++out[i];
      ++mat[i][i];
      sub(mat[i][j], 1);
    };
    for (int i = 1; i <= n; ++i)
      if (s[i] == 'S')
        adde(0, i);
      else if (s[i] == 'R') adde(i, 0);
    for (int i = 1; i <= n; ++i)
      for (int j = 1; j <= n; ++j)
        if (g[i][j] == '1') adde(i, j);
    for (int i = 1; i <= n; ++i)
      if (in[i] == 0) mat[i][i] = 1;
    int ans = det(n);
    for (int i = 1; i <= n; ++i) if (in[i]) ans = ans * (ull) fac[in[i] - 1] % P;
    cout << ans << '\n';
  }

#ifdef ELEGIA
  LOG("Time: %dms\n", int ((clock()
          -nol_cl) / (double)CLOCKS_PER_SEC * 1000));
#endif
  return 0;
}