#include <bits/stdc++.h>
using namespace std;
const int MAX = 5010;
const int MOD = 1000000007;
const int INF = 1e8;
const long long LLINF = 0x3f3f3f3f3f3f3f3f;
const long double EPS = 1e-7;
int n, m;
string s, t;
long long tab[MAX][MAX];
long long dp(int i, int j) {
  int res = 0;
  if (i >= n or j >= m) return 0;
  if (tab[i][j] != -1) return tab[i][j];
  if (s[i] == t[j]) {
    res = 1 + dp(i + 1, j + 1) + dp(i, j + 1);
  } else {
    res = dp(i, j + 1);
  }
  return tab[i][j] = res % MOD;
}
int main() {
  ios_base::sync_with_stdio(false);
  cin.tie(NULL);
  cout.tie(NULL);
  ;
  cin >> s >> t;
  memset(tab, -1, sizeof(tab));
  long long res = 0;
  n = s.size();
  m = t.size();
  for (int i = 0; i < (int)n; i++) {
    res = (res + dp(i, 0)) % MOD;
    ;
  }
  cout << res << '\n';
  return 0;
}
