#pragma GCC optimize("Ofast")
#include <iostream>
#include <vector>
#include <algorithm>
#include <map>
#include <queue>
#include <cstdio>
#include <ctime>
#include <assert.h>
#include <chrono>
#include <random>
#include <numeric>
#include <set>
#include <deque>
#include <stack>
#include <sstream>
#include <utility>
#include <cstring>
#include <unordered_map>
#include <unordered_set>
#include <tuple>
#include <array>
using namespace std;
typedef long long int ll;
typedef unsigned long long ull;

mt19937_64 rng(chrono::steady_clock::now().time_since_epoch().count());
ll myRand(ll B) {
    return (ull)rng() % B;
}
inline double time() {
    return static_cast<double>(chrono::duration_cast<chrono::nanoseconds>(chrono::steady_clock::now().time_since_epoch()).count()) * 1e-9;
}

constexpr ll mod = 1e9+7;
string s[51];

int dp[51][51][22][27];
bool done[51][51][22][27];

void add(int &x, int y){
    x += y;
    if(x >= mod)x -= mod;
}
int mul(int x, int y){
    return 1LL*x*y%mod;
}

int main(){
    cin.tie(nullptr);
    ios::sync_with_stdio(false);
    int n; cin >> n;
    for(int i=0;i<n;i++){
        cin >> s[i];
    }
    auto solve=[&](auto self,int l,int r,int idx,int c)->int{
        if(l == r) return 1;
        // メモ化
        if(done[l][r][idx][c]) return dp[l][r][idx][c];
        done[l][r][idx][c] = true;
        // これ以上文字は使えないよ
        if(c >= 26){
            dp[l][r][idx][c] = 0;
            return 0;
        }
        // 文字列長のチェック
        for(int j=l;j<r;j++){
            if(l<j and (int)s[j].size() == idx){
                dp[l][r][idx][c] = 0; return 0;
            }
            if((int)s[j].size() < idx){
                dp[l][r][idx][c] = 0; return 0;
            }
        }
        // lを一つズラす
        // (上の1つ目のif文で i<j⇨s_{i}<s_{j} が違反されることはない)
        if((int)s[l].size() == idx){
            int res = self(self,l+1,r,idx,c);
            dp[l][r][idx][c] = res;
            return res; 
        }
        // s[l][idx]がc+1の時
        int res = self(self,l,r,idx,c+1);
        // s[l][idx]がcの時
        for(int j=l;j<r;j++){
            if(s[j][idx] != (char)('a'+c) and s[j][idx] != '?')break;
            add(res, mul(self(self,l,j+1,idx+1,0), self(self,j+1,r,idx,c+1)));
        }
        dp[l][r][idx][c] = res;
        return res;
    };
    cout << solve(solve,0,n,0,0) << endl;
}