#include<bits/stdc++.h>
using namespace std;
const int N=15,M=(1<<13)+5,mod=998244353,inv2=499122177;
int a[N];
int n,m,ans;
int g[M],f[M];
int col[N],b[N];
int pt[M];
int A[N][N];
int ksm(int a,int b){
  int res=1;
  while(b){
    if(b&1)res=1ll*res*a%mod;
    a=1ll*a*a%mod,b>>=1;
  }
  return res;
}
int det(int cnt){
  int res=1;
  for(int i=1;i<=cnt;i++){
    int p=0;
    for(int j=i;j<=cnt;j++)
      if(A[j][i]){p=j;break;}
    if(p==0)return 0;
    if(p!=i){
      res=mod-res;
      for(int j=i;j<=cnt;j++)swap(A[p][j],A[i][j]);
    }
    res=1ll*res*A[i][i]%mod;
    int inv=ksm(A[i][i],mod-2);
    for(int j=1;j<=cnt;j++){
      if(j==i)continue;
      int xs=1ll*A[j][i]*inv%mod;
      for(int k=i;k<=cnt;k++)
        A[j][k]=(A[j][k]+mod-1ll*A[i][k]*xs%mod)%mod;
    }
  }
  return res;
}
void work(int cnt){
  int res=1;
  for(int i=1;i<=cnt;i++)res=1ll*res*g[col[i]]%mod;
  if(!res)return;
  for(int i=1;i<cnt;i++){
    int sd=0;
    for(int j=1;j<=n;j++)if(col[i]>>(j-1)&1)sd+=pt[b[j]^(b[j]&col[i])];
    A[i][i]=sd;
    for(int j=i+1;j<cnt;j++){
      sd=0;
      for(int k=1;k<=n;k++)if(col[i]>>(k-1)&1)sd+=pt[(b[k]^(b[k]&col[i]))&col[j]];
      A[i][j]=A[j][i]=(mod-sd)%mod;
    }
  }
  ans=(ans+1ll*res*det(cnt-1)%mod)%mod;
}
void dfs(int x,int cl,int st,int cnt){
  if(x>n){
    if(!g[col[cl]])return;
    work(cl);
    return;
  }
  if(st>n){
    if(g[col[cl]])dfs(x,cl+1,1,0);
    return;
  }
  if(!cnt){
    for(;a[st];st++);
    a[st]=cl;col[cl]|=(1<<st-1);
    dfs(x+1,cl,st+1,1);
    a[st]=0;col[cl]^=(1<<st-1);
    return;
  }
  if(a[st]){
    dfs(x,cl,st+1,cnt);
    return;
  }
  dfs(x,cl,st+1,cnt);
  a[st]=cl;col[cl]|=(1<<st-1);
  dfs(x+1,cl,st+1,cnt+1);
  a[st]=0;col[cl]^=(1<<st-1);
}
int dp[M][N];
int rec[N],tot;
int calc(int s){
  int res=0;tot=0;
  for(int i=1;i<=n;i++)if(s>>(i-1)&1)rec[++tot]=i;
  if(tot==1)return 1;
  if(tot==2)return 0;
  for(int i=0;i<(1<<tot);i++)
    for(int j=1;j<=tot;j++)
      dp[i][j]=0;
  dp[1][1]=1;
  for(int i=1;i<(1<<tot);i++){
    for(int j=1;j<=tot;j++){
      if((!(i>>(j-1)&1))||(!dp[i][j]))continue;
      for(int k=1;k<=tot;k++){
        if((i>>(k-1)&1)||(!(b[rec[j]]&(1<<rec[k]-1))))continue;
        if((i|(1<<k-1))==(1<<tot)-1){
          if(b[rec[1]]&(1<<rec[k]-1))res=(res+dp[i][j])%mod;
          continue;
        }
        dp[i|(1<<k-1)][k]=(dp[i|(1<<k-1)][k]+dp[i][j])%mod;
      }
    }
  }
  return 1ll*res*inv2%mod;
}
signed main(){
  cin>>n>>m;
  for(int i=1,x,y;i<=m;i++){
    cin>>x>>y;
    b[x]|=(1<<y-1);
    b[y]|=(1<<x-1);
  }
  for(int i=1;i<(1<<n);i++)g[i]=f[i]=calc(i),pt[i]=__builtin_popcount(i);
  for(int i=1;i<(1<<n);i++){
    int x=((1<<n)-1)^i;
    if(pt[i]==1)continue;
    for(int k=1;k<=n;k++){
      if(!(i>>(k-1)&1))continue;
      for(int j=x;j>i;j=x&(j-1))g[i|j]=(g[i|j]+1ll*g[i]*f[j|(1<<(k-1))]%mod)%mod;
    }
  }
  dfs(1,1,1,0);
  cout<<ans<<'\n';
  return 0;
}