#include <bits/stdc++.h>
using namespace std;
char S[5003];
char T[5003];
int s, t;
int dp[5003][5003];
int mod = 1000000007;
int main() {
  scanf("%s", S);
  scanf("%s", T);
  s = strlen(S);
  t = strlen(T);
  for (int i = 0; i < s; i++)
    if (S[i] == T[0]) dp[i][0] = 1;
  for (int j = 0; j < t; j++)
    dp[0][j] = dp[0][j - 1] + ((S[0] == T[j]) ? 1 : 0);
  long long ans = 0;
  for (int i = 1; i < s; i++) {
    for (int j = 1; j < t; j++) {
      long long x = dp[i][j - 1];
      if (S[i] == T[j]) {
        x += (1 + dp[i - 1][j - 1]);
      }
      x %= mod;
      dp[i][j] = (int)x;
    }
    ans += dp[i][t - 1];
    ans %= mod;
  }
  ans += dp[0][t - 1];
  ans %= mod;
  cout << ans << endl;
  return 0;
}
