下拉菜单实现树状结构 您所在的位置:网站首页 motif分析图怎么看 下拉菜单实现树状结构

下拉菜单实现树状结构

#下拉菜单实现树状结构| 来源: 网络整理| 查看: 265

我们在日常分析中,有时候会比较不同物种间motif序列的保守性。今天小编教大家使用R包“motifStack ”绘制美观的motif序列结构图!

5d33ba71851e4642bb3e17c364dcd1fc.png ## 安装R包 if (!requireNamespace("BiocManager", quietly = TRUE)) install.packages("BiocManager") BiocManager::install("motifStack")

安装好R包后,我们就可以准备输入文件啦~

比如我想绘制拟南芥中HSF结合的motif结构。

在JASPAR(JASPAR 2020: An open-access database of transcription factor binding profiles)中搜索“HSF”。

de2958b816f3a51dbdeac672213d983c.png

找到自己感兴趣的motif,点击ID号。

f44f3ef858ef2d23c372ccb400095a9b.png

57357abfed108011160e61fcaa1fc956.png

下载JASPAR 格式的矩阵。

506e2b51a660536d9640760bf56e72a6.png

我们需要将它处理成以下格式(空格分隔),并命名为"motif_ID.pcm"。

a16d046398ab73141e20552100f50f9f.png

如果你想展示你自己鉴定的motif结构的话,可以使用小编写的脚本(https://github.com/biozhp/motifStack_input)。

使用脚本前需要准备两个输入文件:

输入文件一:mofit序列(第一列为ID,第二列为序列,Tab分隔)。

d55b7ba4830058cc07d033745e18f4fa.png

输入文件二:mofit ID。

cb2dadcfe8f6afaa30222201b57fe59d.png ## 输入mofit序列, motif ID及输出文件位置 sh ./run.sh motif.seq motif.id out_path

执行脚本完毕后即可获得矩阵文件,开始绘制motif序列结构图。

## 加载R包 library("motifStack") ## 导入motif文件 pcm


【本文地址】

公司简介

联系我们

今日新闻

    推荐新闻

    专题文章
      CopyRight 2018-2019 实验室设备网 版权所有