#include <bits/stdc++.h>
using namespace std;
long long int Pow(long long int a, long long int b, long long int md,
                  long long int ans = 1) {
  for (; b; b >>= 1, a = a * a % md)
    if (b & 1) ans = ans * a % md;
  return ans % md;
}
const long long int MAXN = 5e3 + 10;
const long long int INF = 8e18;
const long long int MOD = 1e9 + 7;
int M[MAXN], dp[MAXN][30], fact[MAXN], inv[MAXN], n, ans = 0;
char C[MAXN];
void add(int &x, int y) {
  y = (y % MOD + MOD) % MOD;
  x = (x + y) % MOD;
}
long long int Ct(int x, int y) {
  if (x < y || y < 0 || x < 0) return 0;
  return 1LL * fact[x] * inv[y] % MOD * inv[x - y] % MOD;
}
int main() {
  scanf("%d%s", &n, C + 1);
  fact[0] = inv[0] = 1;
  for (int i = 1; i < MAXN; i++) {
    fact[i] = 1LL * i * fact[i - 1] % MOD;
    inv[i] = Pow(fact[i], MOD - 2, MOD);
  }
  M[0] = 1;
  for (int i = 1; i <= n; i++) {
    for (int j = i; j; j--) {
      add(M[j], -dp[j][C[i] - 'a' + 1]);
      add(dp[j][C[i] - 'a' + 1],
          M[j - 1] - dp[j - 1][C[i] - 'a' + 1] - dp[j][C[i] - 'a' + 1]);
      add(M[j], dp[j][C[i] - 'a' + 1]);
    }
  }
  for (int i = 1; i <= n; i++) {
    int res = 0;
    for (int j = 1; j <= 26; j++) add(res, dp[i][j]);
    add(ans, 1LL * res * Ct(n - 1, i - 1) % MOD);
  }
  printf("%d\n", ans);
  return 0;
}
