#include <bits/stdc++.h>
using namespace std;
const int MAX = 5000 + 10;
int dp[MAX][MAX];
int mod = 1e9 + 7;
int main() {
  string s, t;
  cin >> s >> t;
  for (int i = 0; i < s.size(); i++) {
    for (int j = 0; j < t.size(); j++) {
      dp[i + 1][j + 1] = dp[i + 1][j];
      if (s[i] == t[j]) dp[i + 1][j + 1] += dp[i][j] + 1;
      dp[i + 1][j + 1] %= mod;
    }
  }
  int ans = 0;
  for (int i = 0; i < s.size(); i++) {
    ans += dp[i + 1][t.size()];
    ans %= mod;
  }
  cout << ans << endl;
  return 0;
}
