#include <bits/stdc++.h>
using namespace std;
const int mod = 1e9 + 7;
int add(int a, int b) { return a + b < mod ? a + b : a + b - mod; }
int mul(int a, int b) { return int(1LL * a * b % mod); }
const int N = 5005;
int C[N][N];
int dp[N][N];
int n;
char s[N];
int nxt[N][26];
int main() {
  for (int i = 0; i < N; i++) {
    C[i][0] = C[i][i] = 1;
    for (int j = 1; j < i; j++) {
      C[i][j] = add(C[i - 1][j], C[i - 1][j - 1]);
    }
  }
  scanf("%d %s", &n, s);
  for (int i = 0; i < 26; i++) nxt[n][i] = n;
  for (int i = n - 1; i >= 0; i--) {
    for (int j = 0; j < 26; j++) nxt[i][j] = nxt[i + 1][j];
    nxt[i][int(s[i] - 'a')] = i;
  }
  for (int j = n; j >= 1; j--) {
    int tot = 0;
    for (int i = n - 1; i >= 0; i--) {
      int &ans = dp[i][j];
      ans = 0;
      int me = int(s[i] - 'a');
      if (nxt[i + 1][me] != n) {
        tot = add(tot, mod - dp[nxt[i + 1][me]][j + 1]);
      }
      ans = tot;
      ans = add(ans, C[n - 1][j - 1]);
      tot = add(tot, dp[i][j + 1]);
    }
  }
  int ans = 0;
  for (int c = 0; c < 26; c++) {
    int it = nxt[0][c];
    if (it == n) continue;
    ans = add(ans, dp[it][1]);
  }
  printf("%d\n", ans);
}
