#include <bits/stdc++.h>
using namespace std;
inline bool EQ(double a, double b) { return fabs(a - b) < 1e-9; }
inline bool isLeapYear(long long int year) {
  return (year % 400 == 0) | (year % 4 == 0 && year % 100 != 0);
}
inline void normal(long long int &a) {
  a %= 1000000007;
  (a < 0) && (a += 1000000007);
}
inline long long int modMul(long long int a, long long int b) {
  a %= 1000000007, b %= 1000000007;
  normal(a), normal(b);
  return (a * b) % 1000000007;
}
inline long long int modAdd(long long int a, long long int b) {
  a %= 1000000007, b %= 1000000007;
  normal(a), normal(b);
  return (a + b) % 1000000007;
}
inline long long int modSub(long long int a, long long int b) {
  a %= 1000000007, b %= 1000000007;
  normal(a), normal(b);
  a -= b;
  normal(a);
  return a;
}
inline long long int modPow(long long int b, long long int p) {
  long long int r = 1;
  while (p) {
    if (p & 1) r = modMul(r, b);
    b = modMul(b, b);
    p >>= 1;
  }
  return r;
}
inline long long int modInverse(long long int a) {
  return modPow(a, 1000000007 - 2);
}
inline long long int modDiv(long long int a, long long int b) {
  return modMul(a, modInverse(b));
}
inline bool isInside(pair<int, int> p, long long int n, long long int m) {
  return (p.first >= 0 && p.first < n && p.second >= 0 && p.second < m);
}
inline bool isInside(pair<int, int> p, long long int n) {
  return (p.first >= 0 && p.first < n && p.second >= 0 && p.second < n);
}
inline bool isSquare(long long int x) {
  long long int s = sqrt(x);
  return (s * s == x);
}
inline bool isFib(long long int x) {
  return isSquare(5 * x * x + 4) || isSquare(5 * x * x - 4);
}
inline bool isPowerOfTwo(long long int x) {
  return ((1LL << (long long int)log2(x)) == x);
}
template <typename first, typename second>
ostream &operator<<(ostream &os, const pair<first, second> &p) {
  return os << "(" << p.first << ", " << p.second << ")";
}
template <typename T>
ostream &operator<<(ostream &os, const vector<T> &v) {
  os << "{";
  for (auto it = v.begin(); it != v.end(); ++it) {
    if (it != v.begin()) os << ", ";
    os << *it;
  }
  return os << "}";
}
template <typename T>
ostream &operator<<(ostream &os, const set<T> &v) {
  os << "[";
  for (auto it = v.begin(); it != v.end(); ++it) {
    if (it != v.begin()) os << ", ";
    os << *it;
  }
  return os << "]";
}
template <typename first, typename second>
ostream &operator<<(ostream &os, const map<first, second> &v) {
  os << "[";
  for (auto it = v.begin(); it != v.end(); ++it) {
    if (it != v.begin()) os << ", ";
    os << it->first << " = " << it->second;
  }
  return os << "]";
}
clock_t tStart = clock();
void faltu() { cerr << endl; }
template <typename T>
void faltu(T a[], int n) {
  for (int i = 0; i < n; ++i) cerr << a[i] << ' ';
  cerr << endl;
}
template <typename T, typename... hello>
void faltu(T arg, const hello &...rest) {
  cerr << arg << ' ';
  faltu(rest...);
}
struct func {
  bool operator()(pair<int, int> const &a, pair<int, int> const &b) {
    if (a.first == b.first) return (a.second < b.second);
    return (a.first < b.first);
  }
};
const long long int INF = 0x3f3f3f3f3f3f3f3f;
const long double EPS = 1e-9;
const int inf = 0x3f3f3f3f;
const int mx = (int)1e5 + 9;
struct Point {
  long long int s, p, x, y;
  string name;
  bool operator<(const Point &ob) const {
    return s == ob.s ? p < ob.p : s > ob.s;
  }
};
double distance(Point a, Point b) {
  return sqrt((a.x - b.x) * (a.x - b.x) + (a.y - b.y) * (a.y - b.y));
}
string s, s1, ss;
long long int n, m, a, b, i, j, d, cs = 0, counT = 0, k, ans = 0, l = 0,
                                   sum1 = 0, sum = 0, Max, Min, num;
vector<long long int> vc[30];
map<long long int, long long int> mp;
int frequency[30];
int mat[30][30];
bool vis[30];
void dfs(int t) {
  ss += (char)(t + 'a');
  vis[t] = true;
  for (auto x : vc[t]) {
    if (!vis[x]) dfs(x);
  }
}
int main() {
  int tt;
  scanf("%d", &tt);
  for (int t = 1; t <= tt; t++) {
    cin >> s;
    memset(frequency, 0, sizeof(frequency));
    memset(vis, 0, sizeof(vis));
    memset(mat, 0, sizeof(mat));
    for (int i = 0; i < (26); i++) vc[i].clear();
    for (i = 1; i < s.size(); i++) {
      a = s[i] - 'a';
      b = s[i - 1] - 'a';
      if (mat[a][b]) continue;
      mat[a][b] = mat[b][a] = 1;
      frequency[a]++;
      frequency[b]++;
      vc[a].push_back(b);
      vc[b].push_back(a);
    }
    bool flag = true;
    for (int i = 0; i < 26; i++)
      if (frequency[i] > 2) flag = false;
    ss = "";
    if (flag) {
      for (int i = 0; i < 26; i++) {
        if (!vis[i] and frequency[i] < 2) dfs(i);
      }
      if (ss.size() == 26) {
        cout << "YES\n";
        cout << ss << '\n';
      } else {
        cout << "NO\n";
      }
    } else {
      cout << "NO\n";
    }
  }
}
