S = input()
N = len(S)
P = 10 ** 9 + 7
ans = 0
_sum = 0
Sum = [0] * 26
dp = [0] * N
_tum = 0
Tum = [0] * 26
dpc = [0] * N
dp[0] = 1
_sum = 1
c = ord(S[0]) - ord('a')
Sum[c] = 1
Tum[c] = 1
_tum = 1
dpc[0] = 1
for i in range(N-1):
    c = ord(S[i+1]) - ord('a')
    tmp = _sum - Sum[c] + _tum - Tum[c] + 1
    tmp %= P
    tmp += Sum[c]
    tmp %= P
    ans = (ans + tmp) % P
    _sum = (_sum + tmp) % P
    Sum[c] = (Sum[c] + tmp) % P
    Tum[c] = (Tum[c] + _tum + 1) % P
    _tum = (_tum + 1) * 2 % P - 1
print((ans+1) % P)