現在、genbankのpttファイルを読み込んで、genoplotRを使用してRでゲノムをプロットするために使用しました
plot_gene_map(dna_segs=list(mo),xlims=xlims,annotations=annotMED,annotation_height=5,main="Region",gene_type="side_blocks",dna_seg_scale=TRUE, scale=FALSE)
また、対応するソートされたbamファイルを読み込み、rbamtoolsを使用してカバレッジプロットを作成しました
start<-130000
end<-140000
coords<-as.integer(c(0,start,end))
range<-bamRange(reader,coords)
bamClose(reader)
ad<-alignDepth(range)
add<-data.frame(ad@pos, ad@depth)
count<-bamCountAll(reader,verbose=TRUE)
ggplot(add, aes(x=ad.pos, y=ad.depth)) + geom_histogram(stat="identity", position="dodge") + theme_bw()
これらの 2 つの図を 1 つのグラフに重ねて、R で基本的なゲノム アラインメント ビューアーを作成したいと思います。
どんな助けでも大歓迎です!
ありがとう