0

現在、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 で基本的なゲノム アラインメント ビューアーを作成したいと思います。

どんな助けでも大歓迎です!

ありがとう

4

1 に答える 1