#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 60,M = 10,mod = 1e9+7;
int addv(int x,int y){return (x+y>=mod)?(x+y-mod):(x+y);}
int subv(int x,int y){return (x>=y)?(x-y):(x-y+mod);}
void add(int& x,int y){x = addv(x,y);}
void sub(int& x,int y){x = subv(x,y);}
int qpow(int a,int n){
    int res = 1; a%=mod;
    while(n){
        if(n&1)res = 1ll*res*a%mod;
        n>>=1; a = 1ll*a*a%mod;
    }
    return res;
}
int qinv(int a){return qpow(a,mod-2);}
//
int T,n,buc[M],lim[M],ans;
string s; int a[N];

int fac[N],rfac[N];

int f[N],g[N];

int calc(int N,int x,int c){
    int flag = 1;
    for(int y=0;y<M;y++){
        lim[y] = c-buc[y]-(y>x);
        flag &= lim[y] >= 0;
    }
    if(!flag)return 0;
    memset(f,0,sizeof f); 
    f[c-buc[x]] = rfac[c-buc[x]];
    for(int y=0;y<M;y++)if(x != y){
        memset(g,0,sizeof g);
        for(int a=0;a<=N;a++)if(f[a]){
            for(int b=0;a+b<=N && b<=lim[y];b++){
                add(g[a+b],1ll*f[a]*rfac[b]%mod);
            }
        }
        memcpy(f,g,sizeof f);
    }
    return 1ll*fac[N]*f[N]%mod;
}

void solve(){
    cin>>s; n=s.length(); s=" "+s; ans = 0; for(int i=1;i<=n;i++)a[i] = s[i]-'0';
    memset(buc,0,sizeof buc);  
    
    for(int i=1;i<=n;i++)buc[a[i]]++;
    for(int i=1;i<M;i++)if(buc[ans] <= buc[i])ans = i;
    
    for(int i=1;i<=n;i++){
        int s = i==1;
        for(int j=s;j<a[i];j++){
            memset(buc,0,sizeof buc);
            for(int k=1;k<i;k++)buc[a[k]]++; buc[j]++;
            for(int x=1;x<M;x++)for(int k=buc[x];k<=n;k++){ 
                int w = calc(n-i,x,k);
                add(ans,1ll*w*x%mod);
            }
        }
    }
    for(int len=1;len<n;len++){
        for(int j=1;j<M;j++){
            memset(buc,0,sizeof buc);
            buc[j] = 1;
            for(int x=1;x<M;x++)for(int k=buc[x];k<=len;k++){
                int w = calc(len-1,x,k);
                add(ans,1ll*w*x%mod);
            }
        }
    }
    cout<<ans<<"\n";
}

int main(){
    fac[0] = 1; for(int i=1;i<=50;i++)fac[i] = 1ll*fac[i-1]*i%mod;
    for(int i=0;i<=50;i++)rfac[i] = qinv(fac[i]);
 
    cin>>T;
    while(T--)solve();

    return 0;
}