#include <bits/stdc++.h>
using namespace std;
const long long mod = 1e9 + 7;
bool debug = true;
bool debug1 = true;
int A[1000000 + 1];
long long M = 1000000007;
long long rev(long long b) {
  long long t = 1;
  for (long long p = M - 2; p > 0; p >>= 1) {
    if (p & 1) (t *= b) %= M;
    (b *= b) %= M;
  }
  return t;
}
long long ANS[1000000 + 1];
int main() {
  ios_base::sync_with_stdio(false);
  cin.tie(0);
  string s;
  cin >> s;
  int sz = s.size();
  int a = 0;
  long long ans = 0;
  long long modInv2 = 1;
  modInv2 = rev(2);
  ANS[0] = 0;
  for (int i = 1; i < sz; i++) {
    ANS[i] = (2 * ANS[i - 1] + 1) % mod;
  }
  for (int i = 0; i < sz; i++) {
    if (s[i] == 'a')
      a++;
    else {
      ans = (ans + ANS[a]) % mod;
    }
  }
  cout << ans << endl;
  return 0;
}
