#include <bits/stdc++.h>
using namespace std;
const int mod = 1e9 + 7;
int ncr[1002][1002], bell[1002], sz[1002], par[1002];
int root(int x) { return (par[x] == x) ? x : par[x] = root(par[x]); }
int main() {
  int n, m, ans = 1;
  char s[55][1002];
  scanf("%d%d", &m, &n);
  for (int i = 0; i <= m; i++) {
    ncr[i][0] = 1;
    for (int j = 1; j <= i; j++) {
      ncr[i][j] = (ncr[i - 1][j - 1] + ncr[i - 1][j]) % mod;
    }
  }
  bell[0] = 1;
  for (int i = 1; i <= m; i++) {
    par[i] = i;
    for (int j = 0; j < i; j++) {
      bell[i] += (1LL * ncr[i - 1][j] * bell[j]) % mod;
      bell[i] %= mod;
    }
  }
  for (int i = 1; i <= n; i++) {
    cin >> s[i] + 1;
  }
  for (int i = 1; i <= m; i++) {
    for (int j = 1; j < i; j++) {
      int cnt = 0;
      for (int k = 1; k <= n; k++) {
        if (s[k][i] == s[k][j]) {
          cnt++;
        }
      }
      if (cnt != n) continue;
      par[root(i)] = par[root(j)];
    }
  }
  for (int i = 1; i <= m; i++) sz[root(i)]++;
  for (int i = 1; i <= m; i++) {
    ans = (1LL * ans * bell[sz[i]]) % mod;
  }
  printf("%d\n", ans);
  return 0;
}
