#include <bits/stdc++.h>
using namespace std;
const int Maxn = 1E3 + 5;
const int Mod = 1E9 + 7;
int N, M;
char a[55][Maxn];
int T[Maxn];
int siz[Maxn];
int C[Maxn][Maxn];
int f[Maxn];
int Ans;
int Find(int X) { return (~T[X] ? T[X] = Find(T[X]) : X); }
void Union(int X, int Y) {
  int Fx = Find(X), Fy = Find(Y);
  if (Fx != Fy) T[Fx] = Fy, siz[Fy] += siz[Fx];
}
void PrePare() {
  for (int i = 0; i <= N; ++i) {
    C[i][0] = 1;
    for (int j = 1; j <= i; ++j)
      C[i][j] = (C[i - 1][j - 1] + C[i - 1][j]) % Mod;
  }
  f[0] = 1;
  for (int i = 1; i <= N; ++i)
    for (int j = 1; j <= i; ++j)
      f[i] = (f[i] + (long long)C[i - 1][j - 1] * f[i - j]) % Mod;
}
int main() {
  scanf("%d %d", &N, &M);
  for (int i = 1; i <= M; ++i) scanf("%s", a[i] + 1);
  PrePare();
  memset(T, -1, sizeof(T));
  for (int i = 1; i <= N; ++i) siz[i] = 1;
  for (int i = 1; i <= N; ++i)
    for (int j = i + 1, k; j <= N; ++j) {
      for (k = 1; k <= M; ++k)
        if (a[k][i] != a[k][j]) break;
      if (k > M) Union(i, j);
    }
  Ans = 1;
  for (int i = 1; i <= N; ++i)
    if (T[i] == -1) Ans = (long long)Ans * f[siz[i]] % Mod;
  printf("%d", Ans);
  return 0;
}
