#include <algorithm>
#include <functional>
#include <iostream>
#include <iomanip>
#include <vector>
#include <map>
#include <set>
#include <unordered_set>
#include <queue>
#include <numeric>
#include <bitset>
#include <cmath>

using namespace std;
using ll = long long;
using ull = unsigned long long;
// const ll mod =  998244353 ;
const ll mod =  1000000007 ;
const ll inf = 1e18;
const int dx[] = {1, -1, 0, 0};
const int dy[] = {0, 0, 1, -1};


#define REP(i, n) for (ll i = 0; i < (n); ++i)
#define REP_FROM(i, j, n) for (ll i = (j); i < (n); ++i)
#define REP_REV(i, n) for (ll i = n-1; i >= 0; --i)
#define all(x) (x).begin(),(x).end()
#define sz(x) int(x.size())

template<typename T> void print_vec(vector<T> a) {int n = a.size();REP(i, n) {cout >> a[i]; if(i!=n-1) cout << " "; else cout << endl;}}
ll power(ll base, ll exponent){ if(exponent % 2){ return power(base, exponent - 1) * base % mod; }else if(exponent){ ll root_ans = power(base, exponent / 2); return root_ans * root_ans % mod; }else{ return 1; }}
ll inverse(ll x){ return power(x, mod - 2);}
ll gcd(ll a, ll b) {if(a < b) gcd(b, a);ll r;while(r=a%b) {a=b;b=r;}return b;}
template<typename T> ll sum(T begin, T end) {return accumulate(begin, end, 0ll);}
struct combination{ vector<ll> fact, inv; combination(int sz) : fact(sz + 1), inv(sz + 1){ fact[0] = 1; for(int i = 1; i <= sz; i++){ fact[i] = fact[i - 1] * i % mod; } inv[sz] = power(fact[sz], mod - 2); for(int i = sz - 1; i >= 0; i--){ inv[i] = inv[i + 1] * (i + 1) % mod; } } ll C(int p, int q) const{ if(q < 0 || p < q) return 0; return (fact[p] * inv[q] % mod * inv[p - q] % mod); } };

using Pair = pair<ll, ll>;

template <ll Modulus> struct ModInt {
    ll a;
    constexpr ModInt(const ll x = 0) noexcept  : a(x % mod){}
    constexpr ll &value() noexcept {return a;}
    constexpr const ll &value() const noexcept {return a;}
    constexpr ModInt operator+(const ModInt x) const noexcept {
        return ModInt(*this) += x;
    }
    constexpr ModInt operator-(const ModInt x) const noexcept {
        return ModInt(*this) -= x;
    }
    constexpr ModInt operator*(const ModInt x) const noexcept {
        return ModInt(*this) *= x;
    }
    constexpr ModInt operator/(const ModInt x) const noexcept {
        return ModInt(*this) /= x;
    }
    constexpr ModInt &operator+=(const ModInt x) noexcept {
        a += x.a;
        if (a >= Modulus) {
            a -= Modulus;
        }
        return *this;
    }
    constexpr ModInt &operator-=(const ModInt x) noexcept {
        if (a < x.a) {
            a += Modulus;
        }
        a -= x.a;
        return *this;
    }
    constexpr ModInt &operator*=(const ModInt x) noexcept {
        a = a * x.a % Modulus;
        return *this;
    }
    constexpr ModInt &operator/=(ModInt x) noexcept {
        ll exp = Modulus - 2;
        while (exp) {
            if (exp % 2) {
                *this *= x;
            }
            x *= x;
            exp /= 2;
        }
        return *this;
    }

};

using mint = ModInt<mod>;

vector<vector<ll>> dp;
vector<ll> a, b;

ll dfs(ll i, ll j) {
    if(dp[i][j] > -inf) {
        return dp[i][j];
    }
    ll s = (i + j) % 2 == 0? 1: 0;
    if(j == sz(b)) {
        return dp[i][j] = dfs(i+1, j) + a[i] * s;
    }
    if(i == sz(a)) {
        return dp[i][j] = dfs(i, j+1) + b[j] * s;
    }
    if(s){
        return dp[i][j] = max(dfs(i+1, j) + a[i] * s, dfs(i, j+1) + b[j] * s);
    }
    return dp[i][j] = min(dfs(i+1, j) + a[i] * s, dfs(i, j+1) + b[j] * s);
}

signed main() {
    ios::sync_with_stdio(false);
    cout << fixed << setprecision(15);
    string s;
    cin >> s;
    int n = sz(s);
    vector<mint> a(n+1), b(n+1), c(n+1), x(n+1);
    x[0] = 1;
    REP(i, n) {
        a[i+1] = a[i];
        b[i+1] = b[i];
        c[i+1] = c[i];
        x[i+1] = x[i];
        if(s[i] == 'A') {
            a[i+1] += x[i];
        }
        else if(s[i] == 'B') {
            b[i+1] += a[i];
        }
        else if(s[i] == 'C') {
            c[i+1] += b[i];
        }
        else {
            a[i+1] *= 3;
            b[i+1] *= 3;
            c[i+1] *= 3;
            x[i+1] *= 3;
            a[i+1] += x[i];
            b[i+1] += a[i];
            c[i+1] += b[i];
        }
    }
    cout << c[n].value() << endl;
    return 0;
}


Measured memory_bytes: 3664