#include <bits/stdc++.h>
using namespace std;
const int mod = 998244353;
int add(int x, int y) { return (x + y) % mod; }
int mul(int x, int y) { return (x * 1ll * y) % mod; }
int main() {
  ios_base::sync_with_stdio(false);
  cin.tie(nullptr);
  string s, t;
  cin >> s >> t;
  int n = s.size(), m = t.size();
  int dp[2][n][26], ind = 0;
  memset(dp, 0, sizeof(dp));
  for (int i = 0; i < n; i++) {
    if (i < m) {
      if (s[0] == t[i]) dp[ind][i][s[0] - 'a'] = 1;
    } else {
      dp[ind][i][s[0] - 'a'] = 1;
    }
  }
  int ans = 0;
  if (0 >= m - 1) ans = dp[ind][0][t[0] - 'a'];
  ind ^= 1;
  for (int i = 1; i < n; i++) {
    for (int j = 0; j < n; j++)
      for (int z = 0; z < 26; z++) dp[ind][j][z] = 0;
    for (int j = 0; j < n; j++) {
      for (int z = 0; z < 26; z++) {
        if (j - 1 >= 0) {
          if (j - 1 < m) {
            if (s[i] == t[j - 1])
              dp[ind][j - 1][z] = add(dp[ind ^ 1][j][z], dp[ind][j - 1][z]);
          } else {
            dp[ind][j - 1][z] = add(dp[ind ^ 1][j][z], dp[ind][j - 1][z]);
          }
        }
        if (j + i < m) {
          if (s[i] == t[j + i]) {
            dp[ind][j][s[i] - 'a'] =
                add(dp[ind ^ 1][j][z], dp[ind][j][s[i] - 'a']);
          }
        } else if (j + i < n) {
          dp[ind][j][s[i] - 'a'] =
              add(dp[ind ^ 1][j][z], dp[ind][j][s[i] - 'a']);
        }
      }
    }
    if (i >= m - 1)
      for (int z = 0; z < 26; z++) ans = add(ans, dp[ind][0][z]);
    ind ^= 1;
  }
  cout << mul(ans, 2) << "\n";
  return 0;
}
