#include <bits/stdc++.h>
using namespace std;
const int mod = 1000000007;
long long v[1010];
char s[1010];
long long dp[1010][1010];
int cnt;
long long work(int n) {
  for (int i = 0; i <= n; i++)
    for (int j = 0; j <= i; j++) dp[i][j] = 0;
  dp[0][0] = 1;
  for (int i = 1; i <= n; i++)
    for (int j = 1; j <= i; j++)
      dp[i][j] = (dp[i - 1][j - 1] + dp[i - 1][j] * j) % mod;
  long long ret = 0;
  for (int i = 1; i <= n; i++) ret += dp[n][i];
  return ret % mod;
}
int main() {
  int m, n;
  long long ans = 1;
  scanf("%d%d", &m, &n);
  memset(dp, 0, sizeof(dp));
  memset(v, 0, sizeof(v));
  for (int i = 1; i <= n; i++) {
    scanf("%s", s + 1);
    for (int j = 1; j <= m; j++) v[j] = (v[j] << 1) | (s[j] - '0');
  }
  sort(v + 1, v + 1 + m);
  int cnt = 0;
  v[0] = v[1];
  v[m + 1] = -1;
  m++;
  for (int i = 1; i <= m; i++)
    if (v[i] == v[i - 1])
      cnt++;
    else {
      (ans *= work(cnt)) %= mod;
      cnt = 1;
    }
  printf("%I64d", ans);
  return 0;
}
