import sys
input = sys.stdin.readline
from collections import *

S = input()[:-1]
N = len(S)
sub_acc = [[0] for _ in range(26)]
MOD = 10**9+7

for i in range(N):
    for j in range(26):
        sub_acc[j].append((sub_acc[j][-1]+(pow(2, i, MOD) if ord(S[i])-ord('a')==j else 0))%MOD)
        
ans = 0

for i in range(N):
    ans += (pow(2, i, MOD)-sub_acc[ord(S[i])-ord('a')][i])*pow(2, N-1-i, MOD)
    ans %= MOD

print(ans)