#include <bits/stdc++.h>
using namespace std;
long long gcd(long long a, long long b) {
  if (a == 0) return b;
  return gcd(b % a, a);
}
long long vpow(long long a, long long b, long long p) {
  if (b == 0) return 1;
  if (b == 1) return a;
  long long temp = vpow(a, b / 2, p);
  temp *= temp;
  temp %= p;
  if (b % 2 == 0) return temp;
  temp = temp * a;
  temp %= p;
  return temp;
}
long long inverseMod(long long a, long long b, long long p) {
  long long temp = vpow(b, p - 2, p);
  long long ans = temp * a;
  ans %= p;
  return ans;
}
long long minl(long long a, long long b) {
  if (a > b) return b;
  return a;
}
long long maxl(long long a, long long b) {
  if (a > b) return a;
  return b;
}
void solve() {
  long long x;
  cin >> x;
  string s;
  cin >> s;
  long long l = 0;
  while (l < x && s.size() < x) {
    string temp = s.substr(l + 1, s.size() - l - 1);
    for (__typeof(((long long)(s[l] - '0'))) i = (1);
         i < ((long long)(s[l] - '0')); i++) {
      s += temp;
    }
    l++;
  }
  if (l == x) {
    cout << s.size() << endl;
    return;
  } else {
    long long ans = s.size();
    while (l < x) {
      ans += ((long long)(s[l] - '0') - 1) *
             ((ans + 1000000007 - l - 1) % 1000000007);
      ans %= 1000000007;
      l++;
    }
    cout << ans << endl;
    return;
  }
}
void prep() {}
int main() {
  ios_base::sync_with_stdio(false);
  cin.tie(0);
  cout.tie(0);
  long long t = 1;
  cin >> t;
  prep();
  cout << fixed << setprecision(12);
  while (t--) solve();
  return 0;
}
