跟着Nature Ecology&Evolution学作图:R语言ggmsa包展示多序列比对结果
2023-02-19 12:27:40 时间
论文
https://www.nature.com/articles/s41559-022-01771-6#code-availability
论文没有权限下载
但是查看数据代码链接的时候发现github主页上提供了论文的下载链接
论文中的图做的都非常好看,而且提供数据和代码,我们可以找来学习
数据代码链接
https://github.com/sebepedroslab/chromatin-evolution-analysis
今天的推文我们学习一下论文中提供的画多序列比对的代码。我没有在论文中找到对应的图,只是github的链接里有数据和代码
论文中国提供的代码
library(stringr)
library(seqinr)
library(msa)
library(Biostrings)
list_ali = list(
"data/20220620/archaeal_Ntails/tails_archaea-sub-MXKK.g.fasta",
"data/20220620/archaeal_Ntails/tails_archaea-sub-mcl1.g.fasta",
"data/20220620/archaeal_Ntails/tails_archaea-sub-mcl2.g.fasta",
"data/20220620/archaeal_Ntails/tails_archaea-sub-mcl4.g.fasta"
)
for (i in 1:length(list_ali)) {
fai = list_ali[[i]]
ali = Biostrings::readAAMultipleAlignment(fai, format = "fasta")
# plot alignment
msa::msaPrettyPrint(
ali, askForOverwrite=FALSE,
shadingMode = "functional",
shadingModeArg = "chemical",
showConsensus = "none",
showLogo = "none",
shadingColors = "blues",
logoColors = "chemical",
psFonts = TRUE,
paperWidth = 32, paperHeight = 10,
output = "pdf",
file = sprintf("%s.colour.pdf",basename(fai)))
}
这里作图用到的是msa
这个R包
读取多序列比对的数据使用的是Biostrings
这个R包
上面的代码是写了一个简单的循环,做了四个数据的图,我试着做其中一个图,但是遇到了报错
Error in texi2dvi(texfile, quiet = !verbose, pdf = identical(output, "pdf"), :
unable to run pdflatex on 'tails_archaea-sub-MXKK.g.fasta.colour.tex'
LaTeX errors:
! Paragraph ended before \inf@@get was complete.
<to be read again>
\par
l.26 ...cal/Temp/RtmpagEAWD/seq397c6ecb2093.fasta}
! Misplaced alignment tab character &.
\msfline ->\par &
& & & @
l.26 ...cal/Temp/RtmpagEAWD/seq397c6ecb2093.fasta}
! Misplaced alignment tab character &.
\msfline ->\par & &
& & @
l.26 ...cal/Temp/RtmpagEAWD/seq397c6ecb2093.fasta}
! Misplaced alignment tab character &.
\msfline ->\par & & &
& @
l.26 ...cal/Temp/RtmpagEAWD/seq397c6ecb2093.fasta}
搞不懂是什么原因
R语言里做多序列比对的图还有更好的选择使用 ggmsa 这个R包
安装
devtools::install_github("YuLab-SMU/ggmsa")
这里R语言必须是4.1以上的
输入数据是比对好的fasta文件的路径
作图代码
library(ggmsa)
fai<-"data/20220620/archaeal_Ntails/tails_archaea-sub-MXKK.g.fasta"
pdf(file = "aligned_fasta.pdf",width = 9.4,height = 4)
ggmsa(fai)+
geom_seqlogo() + geom_msaBar()
dev.off()
image.png
这个我认为比论文中提供的代码出图要好看的多
这个是论文中的代码出的图
相关文章
- PS简洁水彩手绘+颗粒感及AI绘画教程!PS全版本软件下载地址
- CAD 2020软件免费下载及安装教程CAD软件全版本下载
- AutoCAD 2022软件免费下载及安装教程 pdf编辑器全版本下载
- 基于ERNIELayout&pdfplumber-UIE的多方案学术论文信息抽取
- 零售、跨境外贸中的6大主要EDI报文
- 【BSP视频教程】BSP视频教程第25期:CAN/CANFD/CANopen专题,CAN知识点干货分享, 收发执行过程和错误帧处理(2023-01-03)
- 红中群内每日分享题目解析——第一天
- 桌面客户端性能提升,优化使用资源消耗
- EazyDraw for Mac(矢量图绘制软件)v10.10.2中文版
- 绿盟 2019物联网安全年报
- 技术汇总:第十六章:关于登录与退出的token
- 【玩转腾讯云】必须保障169.254网段的路由网关正确
- 读 NebulaGraph源码 | 查询语句 LOOKUP 的一生
- 连接格式优化,支持自定义
- 实用工具分享~IDM下载器idm多个版本(电脑、手机、浏览器插件都有)
- 紧急推送:高速下载,破除下载权限,包括百毒文库,网盘,视频...idm多个版本
- 【Mac】AutoCAD 2019软件免费下载及安装教程 CAD软件全版本下载
- 下载神器IDM_6.37绿色特别版分享idm多个版本(电脑、手机、浏览器插件都有)
- 【生信技能树培训笔记】R语言基础(20230112更新)
- 破解网盘限速 !idm多个版本(电脑、手机、浏览器插件都有)