#include<bits/stdc++.h>
using namespace std;
const int maxn=5050;
const int mod=998244353;
int n,m,u[maxn<<1],v[maxn<<1],pre[maxn<<1],now[maxn],son[maxn<<1],tot=1;
int dfn[maxn],low[maxn],ok[maxn<<1],scc[maxn];
int cnt[maxn],siz[maxn],pw2[maxn],tmp[maxn];
int f[maxn][maxn];
vector<int> G[maxn];
int power(int x,long long y)
{
    y%=(mod-1);
    int ans=1;
    while(y)
    {
        if(y&1)
            ans=1ll*ans*x%mod;
        y>>=1;
        x=1ll*x*x%mod;
    }
    return ans;
}
void put(int x,int y)
{
    pre[++tot]=now[x];
    now[x]=tot;
    son[tot]=y;
}
void tarjan(int x,int fa)
{
    dfn[x]=low[x]=++dfn[0];
    for(int p=now[x];p;p=pre[p])
    {
        int t=son[p];
        if(t==fa)
            continue;
        if(!dfn[t])
        {
            tarjan(t,x);
            low[x]=min(low[x],low[t]);
            if(low[t]>dfn[x])
                ok[p]=ok[p^1]=1;
        }
        else
            low[x]=min(low[x],dfn[t]);
    }
}
void dfs1(int x,int id)
{
    scc[x]=id;cnt[id]++;
    for(int p=now[x];p;p=pre[p])
    {
        int t=son[p];
        if(!scc[t]&&!ok[p])
            dfs1(t,id);
    }
}
void dfs2(int x,int fa)
{
    siz[x]=cnt[x];
    f[x][cnt[x]]=1;
    for(int t:G[x])
        if(t!=fa)
        {
            dfs2(t,x);
            for(int i=0;i<=siz[x]+siz[t];i++)
                tmp[i]=0;
            for(int i=1;i<=siz[x];i++)
                for(int j=1;j<=siz[t];j++)
                    tmp[i+j]=(tmp[i+j]+1ll*f[x][i]*f[t][j])%mod,tmp[i]=(tmp[i]-1ll*pw2[j]*f[x][i]%mod*f[t][j]%mod+mod)%mod;
            siz[x]+=siz[t];
            for(int i=1;i<=siz[x];i++)
                f[x][i]=tmp[i];
        }
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
        pw2[i]=power(2,1ll*i*(i-1)/2-(i-1));
    for(int i=1;i<=m;i++)
        scanf("%d%d",&u[i],&v[i]),put(u[i],v[i]),put(v[i],u[i]);
    tarjan(1,0);
    for(int i=1;i<=n;i++)
        if(!scc[i])
            dfs1(i,++scc[0]);
    for(int i=1;i<=m;i++)
        if(scc[u[i]]!=scc[v[i]])
            G[scc[u[i]]].push_back(scc[v[i]]),G[scc[v[i]]].push_back(scc[u[i]]);
    dfs2(1,0);
    int ans=0;
    for(int i=1;i<=n;i++)
        ans=(ans+1ll*f[1][i]*pw2[i])%mod;
    printf("%lld\n",1ll*ans*power(2,1ll*(mod-2)*(m-n+1))%mod);

    return 0;
}