DNA序列重组可视化方法

这是高老师文章里的图，每一个小图都是通过比较两条比对序列差异而做出来的，从图A和B中可以看出，CF_YL21有可能是Mont和Oz两条序列重组而来。高老师自己做出这图后，深知手工统计序列差异、画图、拼图之困难。问我有没有办法读两条序列，生成一个图。

很多东西都没有现成的，需要自己造轮子，于是我写了个R包来做这个事情：

1. library(seqcombo)

2. 
   

3. fas = list.files(system.file("examples","GVariation", package="seqcombo"),

4. 
   

5. pattern="fas", full.names=TRUE)

6. 
   

7. x = seqdiff(fas[1])

8. 
   

9. plot(x)

一行代码读文件，然后plot画出来，就是这么简单。

这里有三个文件，我们可以三个文件一起读，三个图一起画。

1. x = lapply(fas, seqdiff)

2. 
   

3. plts = lapply(x, plot)

然后可以很容易把三张图拼在一起

1. plot_grid(plotlist=plts, ncol=1, labels=LETTERS[1:3])

简直完美，跟文章中的图一毛一样。

这个包我取名 seqcombo，90后大概不知道何为combo，如果你知道，没错，说的就是你，你已经暴露年龄了！因为这个图是可以辅助看序列重组的，所以取这个名字，这个功能写于2016年，当然这个包目前可不止这个可视化方法，起码得有三把刷子，才对得起combo这名字，以后会陆续介绍其它的功能。

往期精彩

• 爱要高声喊出来

• （翻）云（覆）雨图

• 邪恶的我画出了善良的饼图

• 好饿，连刚需都解决不了的医学院

• reverse time scale

• 文章发表：两种方法在进化树上可视化数据

• 公众号「biobabble」是什么意思？

• learn R wisdom with fortunes

点击阅读原文，可揭开序列宝宝的内幕。