s = input()
ngs = len(s)
num = 26
dp = [[0]*2 for i in range(num)]
#dp[0][0] = 1

mod = 10**9+7

for i in range(ngs):
    idx = ord(s[i])-ord("a")
    dp_n = [[0]*2 for i in range(num)]
    dp_n[idx][0] += 1
    dp_n[idx][1] += 1
    #print(dp_n)
    for j in range(num):
        if j == idx:
            dp_n[idx][0] += dp[j][0]*2
            dp_n[idx][1] += dp[idx][1]*2
            dp_n[idx][0] %= mod
            dp_n[idx][1] %= mod
        elif dp[j][0] != 0:
            dp_n[idx][0] += dp[j][0]+dp[j][1]
            dp_n[idx][0] %= mod
            dp_n[idx][1] += dp[j][1]
            dp_n[idx][1] %= mod
            dp_n[j][0] = dp[j][0]
            dp_n[j][0] %= mod
            dp_n[j][1] = dp[j][1]
            dp_n[j][1] %= mod
        #else:
        #    dp_n[j] = dp[i][j]
    #print(dp_n)
    dp = dp_n[:]

ans = 0
#print(dp)
for i in range(num):
    ans += dp[i][0]
    ans %= mod

print(ans)