#include<bits/stdc++.h>
using namespace std;
#define fo(i,a,b) for(int i=a;i<=b;++i)
#define fod(i,a,b) for(int i=a;i>=b;--i)
#define p 1000000007
#define ll long long
int read(){
    int u=0;
    char c=getchar();
    while(c>'9'||c<'0')c=getchar();
    while(c>='0'&&c<='9'){
        u=u*10+c-48;
        c=getchar();
    }
    return u;
}
int n;
int s[100011][26],len[100011];
struct Sub{
    int A[26];
    int l,r;
    friend bool operator <(Sub x,Sub y){
        return (x.l<y.l);
    }
    void build(int x){
        l=0,r=len[x]-1;
        fo(j,0,25)A[j]=s[x][j];
    }
    bool Cut(Sub y){
        fo(j,0,25){
            if(y.A[j]>A[j])return 0;
            A[j]-=y.A[j];
        }
        l=y.r+1;
        return 1;
    }
};
Sub fir;
char s2[5000011];
set<Sub>F;
string ans[100011];
bool Add(Sub X){
    auto it=F.begin();
    while(it!=F.end()){
        if((*it).r<=X.r){
            int v=X.Cut(*it);
            if(!v){
                puts("NO");
                return 0;
            }
        }
        else{
            if(X.l<=X.r){
                Sub v=(*it);
                int u=v.Cut(X);
                if(!u){
                    puts("NO");
                    return 0;
                }
                F.erase(it);
                F.insert(X);
                F.insert(v);
            }
            return 1;
        }
        ++it;
    }
    return 1;
}
void qin(Sub x){
    int u=0;
    fo(i,x.l,x.r){
        while(x.A[u]==0)++u;
        if(ans[1][i]=='#'){
            ans[1][i]='a'+u;
        }
        --x.A[u];
    }
}
vector<int>zha;
void solve(){
    n=read();
    fo(i,1,n){
        scanf("%s",s2);
        len[i]=strlen(s2);
        memset(s[i],0,sizeof(s[i]));
        fo(j,0,len[i]-1){
            ++s[i][s2[j]-'a'];
        }
        ans[i].clear();
        fo(j,1,len[i])ans[i].push_back('#');
    }
    F.clear();
    Sub Init;
    Init.build(1);
    F.insert(Init);
    zha.clear();
    zha.push_back(1);
    fo(i,2,n){
        while(zha.size()&&len[i]<=len[zha.back()])zha.pop_back();
        if(!zha.size()){
            Sub X;
            X.build(i);
            int v=Add(X);
            if(!v)return ;
            auto it=F.end();
            --it;
            while((*it).l>=len[i]){
                qin((*it));
                --it;
            }
            ++it;
            while(it!=F.end()){
                auto it2=it;
                ++it2;
                F.erase(it);
                it=it2;
            }
        }
        else{
            int t=zha.back(),t2=zha.size()-1;
            int h=len[i]%len[t];
            Sub X;
            X.build(i);
            int w=len[i]/len[t];
            while(h){
                X.l=0,X.r=h-1;
                fo(j,0,25){
                    X.A[j]=X.A[j]-s[t][j]*w;
                    if(X.A[j]<0){
                        puts("NO");
                        return ;
                    }
                }
                int v=Add(X);
                if(!v)return ;
                if(h<len[zha[0]])break;
                while(len[zha[t2]]>h)--t2;
                t=zha[t2];
                int h2=h%len[t];
                w=h/len[t];
                h=h2;
            }
        }
        zha.push_back(i);
    }
    while(F.size()){
        qin(*F.begin());
        F.erase(F.begin());
    }
    fo(i,2,n){
            fo(j,0,len[i-1]-1){
                int u=j;
                while(u<len[i]){
                    ans[i][u]=ans[i-1][j];
                    u+=len[i-1];
                }
            }
    }
    fo(i,1,n){
        int h[26]={0};
        fo(j,0,len[i]-1)h[ans[i][j]-'a']++;
        fo(j,0,25){
            if(s[i][j]!=h[j]){
                puts("NO");
                return ;
            }
        }
    }
    puts("YES");
    fo(i,1,n){
        cout<<ans[i]<<"\n";
    }
}
int main(){
    int T=read();
    while(T){
        --T;
        solve();
    }
    return 0;
}