前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >可视化基因组结构重排的工具~plotsr

可视化基因组结构重排的工具~plotsr

作者头像
用户7010445
发布2022-05-23 15:48:40
8300
发布2022-05-23 15:48:40
举报

论文

https://academic.oup.com/bioinformatics/advance-article/doi/10.1093/bioinformatics/btac196/6569079?login=false

plotsr: visualizing structural similarities and rearrangements between multiple genomes

github主页

https://github.com/schneebergerlab/plotsr

安装

可以直接使用conda安装

代码语言:javascript
复制
conda install -c bioconda plotsr

用自己的数据试试

首先是做基因组比对

代码语言:javascript
复制
seqkit faidx tunisia_genomic.fna NC_045127.1 > tunisiaChr01.fa
seqkit faidx ../contig2scaffold/ragtag_output/ragtag.scaffold.fasta NC_045127.1_RagTag > ysChr01.fa

minimap2 -ax asm5 -t 8 --eqx tunisiaChr01.fa ysChr01.fa | samtools sort -O BAM - > A_B.bam

结构变异鉴定

代码语言:javascript
复制
syri -c A_B.bam -r tunisiaChr01.fa -q ysChr01.fa -F B --prefix A_B

画图

代码语言:javascript
复制
plotsr --sr A_Bsyri.out --genomes genome.txt -o output_plot.png

这里我的genome.txt文件如下

image.png

这里会报错

image.png

暂时搞不懂是什么原因

运行下github上example里的数据试试

代码语言:javascript
复制
plotsr --sr col_lersyri.filtered.out --sr ler_cvisyri.filtered.out --sr cvi_erisyri.filtered.out --genomes genomes.txt -o output_plot.pdf

运行这个命令能够得到结果

image.png

他的genomes.txt文件是

image.png

我把两个不同物种的染色体的名字改成一样的再试试

代码语言:javascript
复制
minimap2 -ax asm5 -t 8 --eqx tunisiaChr01.fa ysChr01.fa | samtools sort -O BAM - > A_B.bam

samtools index A_B.bam
syri -c A_B.bam -r tunisiaChr01.fa -q ysChr01.fa -F B
plotsr --sr syri.out --genomes genome.txt -o output_plot.pdf

还是没有搞定,有时间再来看看吧

本文参与 腾讯云自媒体同步曝光计划,分享自微信公众号。
原始发表:2022-04-27,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 小明的数据分析笔记本 微信公众号,前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体同步曝光计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
目录
  • 论文
  • github主页
  • 安装
  • 用自己的数据试试
  • 首先是做基因组比对
  • 结构变异鉴定
  • 画图
  • 运行下github上example里的数据试试
  • 我把两个不同物种的染色体的名字改成一样的再试试
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档