#include<bits/stdc++.h>
using namespace std;
#define FOR(i,s,t) for(int i=(s),_t=(t); i<=_t; ++i)
#define DOR(i,s,t) for(int i=(s),_t=(t); i>=_t; --i)
typedef long long ll;
const int N=2005;
const int Mod=1e9+7;
int a[N];
ll f[N][2],g[N][2],h[N];
ll F[N];
int main() {
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    string s;
    cin >> s;
    int n=s.size();
    FOR(i,1,n) {
        a[i]=(s[i-1]-'A'+n-i)%4;
    }
    vector<int> b;
    int m=0;
    FOR(i,1,n) {
        b.push_back(a[i]);
        while(b.size()>=4 && b[b.size()-2]==b.back() && b[b.size()-3]==b.back() && b[b.size()-4]==b.back()) {
            b.pop_back();
            b.pop_back();
            b.pop_back();
            b.pop_back();
            ++m;
        }
    }
    f[0][0]=f[0][1]=g[0][0]=g[0][1]=1;
    h[0]=1;
    FOR(i,1,m) {
        FOR(k,0,i-1) {
            g[i][0]=(g[i][0]+4*h[k]*f[i-1-k][1])%Mod;
            g[i][1]=(g[i][1]+3*h[k]*f[i-1-k][1])%Mod;
        }
        FOR(k,0,i-1) {
            f[i][0]=(f[i][0]+f[k][0]*g[i-k][0])%Mod;
            f[i][1]=(f[i][1]+f[k][1]*g[i-k][1])%Mod;
        }
        FOR(k,0,i) {
            h[i]=(h[i]+f[i-k][1]*f[k][1])%Mod;
        }
    }
    // cerr << f[2] << endl;
    F[0]=1;
    FOR(i,0,b.size()) {
        DOR(j,m,0) {
            FOR(k,1,j) {
                F[j]=(F[j]+f[k][i>0]*F[j-k])%Mod;
            }
        }
    }
    cout << F[m] << '\n';
    return 0;
}