#include <bits/stdc++.h>
using namespace std;
const long long mod = 1e9 + 7;
const long long mod1 = 998244353;
const double pie = 3.1415926535;
long long power(long long x, long long y) {
  if (y == 0)
    return 1;
  else {
    if (y % 2 == 0)
      return power(x * x, y / 2);
    else
      return x * power(x * x, (y - 1) / 2);
  }
}
long long mod_power(long long x, long long y, long long m) {
  long long r = 1;
  x = x % m;
  if (x == 0) return 0;
  while (y > 0) {
    if (y & 1) r = (r * x) % m;
    y = y / 2;
    x = (x * x) % m;
  }
  return r;
}
long long gcd(long long x, long long y) {
  if (y == 0) return x;
  return gcd(y, x % y);
}
vector<long long> extended_Euclid(long long a, long long b) {
  vector<long long> v(3);
  if (b == 0) {
    v[0] = a;
    v[1] = 1;
    v[2] = 0;
    return v;
  } else {
    vector<long long> v1 = extended_Euclid(b, a % b);
    v[0] = v1[0];
    v[1] = v1[2];
    v[2] = v1[1] - (a / b) * v1[2];
    return v;
  }
}
int main() {
  ios_base::sync_with_stdio(false);
  cin.tie(NULL);
  cout.tie(NULL);
  int TESTS = 1;
  cin >> TESTS;
  while (TESTS--) {
    int x;
    cin >> x;
    string s;
    cin >> s;
    long long ans = (int)s.size();
    for (int i = 1; i <= x; i++) {
      long long a = ans - i;
      long long b = i;
      ans = i + ((int)s[i - 1] - 48) * (a);
      ans = (ans + mod) % mod;
      if (s.size() < x) {
        string p = s.substr(i);
        for (int j = 1; j < (int)s[i - 1] - 48; j++) {
          s += p;
        }
      }
    }
    cout << ans << endl;
  }
}
