欢迎关注”生信修炼手册”!
IGV支持动态查看gwas分析结果,对于gwas结果而言,要求至少要包含以下几列
CHR
BP
SNP
P
对列的顺序没有要求,IGV通过文件名后缀来识别文件格式,gwas结果对应的后缀可以是以下几种
linear
logistic
assoc
qassoc
gwas
从GWAS catalog数据库中随机下载了一个gwas的结果,内容示意如下
导入IGV之后,首先以曼哈顿图的形式展示所有位点的关联分析结果
曼哈顿图的含义可以查看以前的文章
还可以根据需要调整基因组区域,达到locuszoom的效果,示意如下
locuszoom的含义可以查看以前的文章
通过基因名称来检索,可以方便的查看基因上突变位点的关联结果。遗憾的一点是,如果只是单独导入gwas结果,无法直接以rs号进行检索,但是IGV支持导入多个track, 此时可以再导入一个rs号对应的bed文件,实现rs号的快速定位,这样就可以方便的以基因组位置,rs号,基因名称来查看gwas结果。
·end·
—如果喜欢,快分享给你的朋友们吧—
原创不易,欢迎收藏,点赞,转发!生信知识浩瀚如海,在生信学习的道路上,让我们一起并肩作战!
本公众号深耕耘生信领域多年,具有丰富的数据分析经验,致力于提供真正有价值的数据分析服务,擅长个性化分析,欢迎有需要的老师和同学前来咨询。
领取专属 10元无门槛券
私享最新 技术干货