R语言如何在图中任何地方加字符-画序列比对结果图

biostat

在R群里看到大家在讨论问题，拷贝过来学些一下
如何在图片中画出序列比对的图形

1. a<-'ATCGGACTAA'
   
2. split_a<-unlist(strsplit(a,''))
   
3. b<-'TAGCCTGATT'
   
4. split_b<-unlist(strsplit(b,''))
   
5. l_a<-nchar(a)
   
6. plot(1:l_a,type='n')
   
7. text(x=seq(3,by=.3,length=l_a),y=3,split_a)
   
8. text(x=seq(3,by=.3,length=l_a),y=2.5,split_b)
   
9. segments(seq(3,by=.3,length=l_a),2.8,seq(3,by=.3,length=l_a),2.7)
   
10. text(x=4,y=5,cex=4,"I love you")
    
11. text(x=6,y=8,cex=5,col="red","I love you")

复制代码

运行结果：