s = input()
mod = 10**9+7

n = len(s)
rem = dict()
ans = n* pow(2,n-1,mod)
ans %= mod
for i,si in enumerate(s):
    if(not si in rem):
        rem[si] = 0
    ans -= rem[si] * pow(2,n-i-1,mod)
    ans %= mod
    rem[si] += pow(2,i,mod)
    rem[si] %= mod
print(ans)