/**
*    author:  boutarou
*    created: 14.08.2020 21:56:06
**/

#include <bits/stdc++.h>
using namespace std;
#define rep(i,n) for(int i = 0; i < int(n); i++)
using ll = long long;
using P = pair<int, int>;

const ll MOD = 1e9 + 7;

long long mod_pow(long long n, long long p, long long m) {
    if (p == 0) return 1;
    if (p % 2 == 0) {
        long long t = mod_pow(n, p / 2, m);
        return (t * t) % m;
    }
    return n * mod_pow(n, p - 1, m) % m;
}

int main() {
    string s;
    cin >> s;
    vector<int>a(26, 0);
    ll ans = 0;
    rep(i, s.size()) {
        ans += mod_pow(2, s.size() - 1 - i, MOD) * (mod_pow(2, i, MOD) - a[s[i] - 'a']);
        // cout << ans << endl;
        ans %= MOD;
        a[s[i] - 'a'] += mod_pow(2, i, MOD);
        a[s[i] - 'a'] %= MOD;
    }
    ans %= MOD;
    if (ans < 0) ans += MOD;
    cout << ans << endl;
    return 0;
}