#include <bits/stdc++.h>
using namespace std;
const int N = 5e3;
int mod = 1e9 + 7;
int mul(int a, int b) { return ((a % mod) * (b % mod)) % mod; }
int add(int a, int b) { return ((a % mod) + (b % mod)) % mod; }
int power(int a, int b) {
  int res = 1;
  while (b > 0) {
    if (b & 1) {
      res = mul(res, a);
    }
    b = b / 2;
    a = mul(a, a);
  }
  return res;
}
int sub(int a, int b) { return ((a % mod) - (b % mod) + mod) % mod; }
int divi(int a, int b) { return mul(a, power(b, mod - 2)); }
int dp[N][N];
int main() {
  ios_base::sync_with_stdio(0);
  cin.tie(0);
  cout.tie(0);
  ;
  string s, t;
  cin >> s >> t;
  int n1 = s.size();
  int n2 = t.size();
  for (int i = 0; i < n1; i++) {
    for (int j = 0; j < n2; j++) {
      if (s[i] == t[j]) {
        if (i > 0 && j > 0)
          dp[i][j] = add(dp[i - 1][j - 1], 1);
        else
          dp[i][j] = 1;
      }
      if (j > 0) dp[i][j] = add(dp[i][j], dp[i][j - 1]);
    }
  }
  int ans = 0;
  for (int i = 0; i < n1; i++) ans = add(ans, dp[i][n2 - 1]);
  cout << ans;
}
