#include <bits/stdc++.h>
#define fi first
#define se second
#define pii pair<int,int>
#define mp make_pair
#define pb push_back
#define space putchar(' ')
#define enter putchar('\n')
#define eps 1e-10
#define MAXN 500005
//#define ivorysi
using namespace std;
typedef long long int64;
typedef unsigned int u32;
typedef double db;
template<class T>
void read(T &res) {
    res = 0;T f = 1;char c = getchar();
    while(c < '0' || c > '9') {
    	if(c == '-') f = -1;
    	c = getchar();
    }
    while(c >= '0' && c <= '9') {
    	res = res * 10 +c - '0';
    	c = getchar();
    }
    res *= f;
}
template<class T>
void out(T x) {
    if(x < 0) {x = -x;putchar('-');}
    if(x >= 10) {
	out(x / 10);
    }
    putchar('0' + x % 10);
}
const int MOD = 1000000007;
int N,M;
int a[205],b[205],fac[205],invfac[205];
int col[(1 << 20) + 5][21],num[(1 << 20) + 5],c[21];
int to[21],e[(1 << 20) + 5],cnt[(1 << 20) + 5];
int dp[(1 << 19) + 5],sum[(1 << 19) + 5],len[(1 << 19) + 5];
int fa[21];
int inc(int a,int b) {
    return a + b >= MOD ? a + b - MOD : a + b;
}
int mul(int a,int b) {
    return 1LL * a * b % MOD;
}
void update(int &x,int y) {
    x = inc(x,y);
}
int lowbit(int x) {
    return x & (-x);
}
int getfa(int x) {
    return fa[x] == x ? x : fa[x] = getfa(fa[x]);
}
int fpow(int x,int c) {
    int res = 1,t = x;
    while(c) {
        if(c & 1) res = mul(res,t);
        t = mul(t,t);
        c >>= 1;
    }
    return res;
}
void Solve() {
    read(N);read(M);
    for(int i = 1 ; i <= M ; ++i) {
        read(a[i]);read(b[i]);
        to[a[i]] |= (1 << b[i] - 1);
        to[b[i]] |= (1 << a[i] - 1);
    }
    fac[0] = 1;
    for(int i = 1 ; i <= 200 ; ++i) fac[i] = mul(fac[i - 1],i);
    invfac[200] = fpow(fac[200],MOD - 2);
    for(int i = 199 ; i >= 0 ; --i) invfac[i] = mul(invfac[i + 1],i + 1);
    for(int S = 0 ; S < (1 << N) ; ++S) {
        if(S != 0) cnt[S] = cnt[S - lowbit(S)] + 1;
        for(int j = 1 ; j <= N ; ++j) {
            if(S >> (j - 1) & 1) {
                int T = to[j] & S;
                e[S] = e[S ^ (1 << j - 1)] + cnt[T];break;
            }
        }
    }
    for(int S = 0 ; S < (1 << N - 1) ; ++S) {
        for(int i = 1 ; i <= N ; ++i) fa[i] = i;
        for(int i = 1 ; i < N ; ++i) {
            if(S >> (i - 1) & 1) {
                fa[getfa(a[i])] = getfa(b[i]);
            }
        }
        int tot = 0;
        for(int i = 1 ; i <= N ; ++i) {
            if(fa[i] == i) {col[S][i] = ++tot;c[tot] = 0;}
        }
        for(int i = 1 ; i <= N ; ++i) {
            col[S][i] = col[S][getfa(i)];
            int t = col[S][i];
            c[t] |= (1 << i - 1);
        }
        for(int i = 1 ; i <= tot ; ++i) {
            num[S] += e[c[i]] - (cnt[c[i]] - 1);
        }
    }
    dp[0] = 1;sum[0] = 0;
    int ALL = (1 << N - 1) - 1;
    for(int S = 0 ; S < (1 << N - 1) ; ++S) {
        for(int j = 1 ; j < N ; ++j) {
            if(!(S >> (j - 1) & 1)) {
                int a = num[ALL ^ S] - num[ALL ^ S ^ (1 << j - 1)];
                int T = S ^ (1 << j - 1);
                int n = len[S];
                len[T] = n + a + 1;
                int c = mul(dp[S],mul(fac[n + a],invfac[n]));
                int s = mul(sum[S],mul(fac[n + a + 1],invfac[n + 1]));
                update(s,mul(c,cnt[S] + 1));
                update(dp[T],c);update(sum[T],s);
            }
        }
    }
    out(sum[(1 << N - 1) - 1]);enter;
}
int main() {
#ifdef ivorysi
    freopen("f1.in","r",stdin);
#endif
    Solve();
    return 0;
}


Measured memory_bytes: 63744