#include "bits/stdc++.h"
using namespace std;
const int N=1e5+20,MAX=26,MOD=(int)(1e9)+7;
string a;
long long n,q,l,r,dp[N][MAX],f[MAX],p[N],ie[N],i,j,s,odd,ans;

int power(int a,int b,int MOD)
{
	if(b==0) return 1;
	int x=power(a,b>>1,MOD);
	x=(1LL*x*x)%MOD;
	return (1LL*x*(b&1?a:1))%MOD;
}

int InverseEuler(int n,int MOD)
{
	return power(n,MOD-2,MOD);
}

int main()
{
	p[0]=1, ie[0]=InverseEuler(p[0],MOD);
	for(i=1;i<N;i++)
	{
		p[i]=(p[i-1]*i)%MOD;
		ie[i]=InverseEuler(p[i],MOD);
	}
	
	cin>>a>>q;
	
	n=a.size();
	a="$"+a;
	for(i=1;i<=n;i++)
	{
		for(j=0;j<MAX;j++) dp[i][j]=dp[i-1][j];
		dp[i][a[i]-'a']++;
	}
	
	while(q--)
	{
		cin>>l>>r;
		
		for(odd=s=j=0;j<MAX;j++)
		{
			f[j]=dp[r][j]-dp[l-1][j];
			if(f[j]&1) f[j]--,odd++;
			f[j]/=2;
			s+=f[j];
		}
		
		ans=p[s];
		for(j=0;j<MAX;j++) ans=(ans*ie[f[j]])%MOD;
		if(odd) ans=(ans*odd)%MOD;
		cout<<ans<<"\n";
	}
}