#pragma GCC Optimize("Ofast")
#include<bits/stdc++.h>
using namespace std;
#define pb push_back
#define int long long
#define lowbit(i) (i&(-i))
const int mod=1e9+7,maxn=5e6+5,M=1e5+5;
inline int read(){
char ch=getchar();int x=0;bool f=0;
for(;!isdigit(ch);ch=getchar())if(ch=='-')f=1;
for(;isdigit(ch);ch=getchar())x=(x<<1)+(x<<3)+(ch^48);
if(f==1)x=-x;return x;
}
int T,n,m,x;
int sz[maxn],st[maxn];
int top,ch[M][27];
vector<int>sum[27];
vector<int>a[maxn];
char s[maxn];
void print(){
puts("YES");
vector<int>tmp;
for(int j=0;j<26;j++){
for(auto i:sum[j])tmp.pb(i*26+j);
}
sort(tmp.begin(),tmp.end());
for(int i=sz[1]-1;i>=0;i--){
s[sz[1]-i-1]=tmp[i]%26+'a';
}
for(int i=1;i<=n;i++){
for(int j=0;j<sz[i];j++)
putchar(s[a[i][j]]);
putchar('a');
puts("");
}
}
int get(int x){
int l=0,r=top,res=0;
while(l<=r){
int mid=(l+r)/2;
if(sz[st[mid]]<=x)res=mid,l=mid+1;
else r=mid-1;
}return st[res];
}
void solve(){
n=read();top=0;
scanf("%s",s);
//if(T==53){cout<<n<<endl;cout<<s<<endl;}
sz[1]=strlen(s);a[1].resize(sz[1]);
for(int i=0;i<sz[1];i++){
ch[1][s[i]-'a']++;
sum[s[i]-'a'].pb(0);
a[1][i]=i;
}
st[++top]=1;bool flag=0;
//cout<<sz[1]<<endl;
for(int i=2;i<=n;i++){
scanf("%s",s);
//if(T==53){cout<<s<<endl;}
sz[i]=strlen(s);
a[i].resize(sz[i]);
vector<int>tmp;
tmp.clear();
tmp.resize(27);
for(int j=0;j<sz[i];j++){
a[i][j]=a[i-1][j%sz[i-1]];
ch[i][s[j]-'a']++;
tmp[s[j]-'a']++;
}
int t=sz[i],z=get(t);
while(z){
int p=t/sz[z];
t=t%sz[z];
for(int j=0;j<26;j++){
tmp[j]=tmp[j]-ch[z][j]*p;
if(tmp[j]<0){flag=1;break;}
}z=get(t);
}
for(int j=0;j<26;j++)if(tmp[j]>0){
if(tmp[j]>sum[j].size()){flag=1;break;}
for(int k=0;k<tmp[j];k++)sum[j][k]++;
}
while(top&&sz[st[top]]>sz[i])top--;
st[++top]=i;
}
if(flag){puts("NO");
else print();
for(int i=1;i<=n;i++){
for(int j=0;j<=26;j++)ch[i][j]=0;
a[i].clear();
}
for(int i=0;i<=26;i++)sum[i].clear();
}
signed main(){
T=read();
while(T--){
solve();
}
}