#include <cstdio>
#include <vector>
#include <stack>
#include <cstring>
using namespace std;
const int MOD=998244353;
int clo,dfn[10007],low[10007],vis[10007],tot,sz[10007],cnt[10007];
long long jc[25000007],dp[5007][5007],f[10007];
stack<int> st;
vector<int> g[10007],e[10007];
void tarjan(int u,int fa){
	dfn[u]=low[u]=++clo;
	st.push(u);
	for(int i=0;i<g[u].size();i++){
		int v=g[u][i];
		if(v==fa) continue;
		if(!dfn[v]){
			tarjan(v,u);
			low[u]=min(low[u],low[v]);
		}else low[u]=min(low[u],dfn[v]);
	}
	if(dfn[u]==low[u]){
		tot++;
		while(st.top()!=u){
			int x=st.top();
			st.pop();
			sz[tot]++;
			vis[x]=tot;
		}
		int x=st.top();
		st.pop();
		sz[tot]++;
		vis[x]=tot;
	}
}
void dfs(int u,int fa){
	dp[u][sz[u]]=jc[sz[u]*(sz[u]-1)/2-cnt[u]];
	for(int i=0;i<e[u].size();i++){
		int v=e[u][i];
		if(v==fa) continue;
		dfs(v,u);
		memset(f,0,sizeof(f));
		for(int j=1;j<=sz[u];j++){
			for(int k=1;k<=sz[v];k++){
				f[j]=(f[j]+MOD-(dp[u][j]*dp[v][k])%MOD)%MOD;
				(f[j+k]+=((dp[u][j]*dp[v][k]%MOD)*jc[j*k-1])%MOD)%=MOD;
			}
			dp[u][j]=f[j];
		}
		for(int j=sz[u];j<=sz[u]+sz[v];j++) dp[u][j]=f[j];
		sz[u]+=sz[v];
	}
}
int x[10007],y[10007];
int main(){
	jc[0]=1;
	for(int i=1;i<=25000000;i++) jc[i]=jc[i-1]*2%MOD;
	int n,m;
	scanf("%d%d",&n,&m);
	for(int i=1;i<=m;i++){
		scanf("%d%d",&x[i],&y[i]);
		g[x[i]].push_back(y[i]);
		g[y[i]].push_back(x[i]);
	}
	tarjan(1,0);
	for(int i=1;i<=m;i++){
		int u=x[i],v=y[i];
		if(vis[u]==vis[v]) cnt[vis[u]]++;
		else{
			e[vis[u]].push_back(vis[v]);
			e[vis[v]].push_back(vis[u]);
		}
	}
	dfs(1,0);
	long long ans=0;
	for(int i=1;i<=sz[1];i++) (ans+=dp[1][i])%=MOD;
	printf("%lld",ans);
	return 0;
}