这是一个小笔记,如何把bed直接转成bw(bigwig)文件。思路有很多,讲其中一种思路:
1:首先把bed文件转成bedgraph文件
bedtools genomecov -i a.bed -g /asnas/xuchh_group/liulinli/Database/hg38/hg38.chrom.sizes.bed -bg > CTCF2_q10.pair.bam.chr_start_end.bedgraph
(这里a.log就是需要转换的bed文件,-g 是输入基因组染色体大小的文件)
这里得到的bedgraph的格式:
chr start end value
然后得到的bedgraph文件,可以接着用bedgraphtobigwig命令转换:
bigwig其实是bedgraph的二进制文件。然后就可以做后续的分析啦:
bedGraphToBigWig CTCF2_q10.pair.bam.chr_start_end.bedgraph /asnas/xuchh_group/liulinli/Database/hg38/hg38.chrom.sizes.bed > CTCF2_q10.pair.bw
Ref:
1:https://bedtools.readthedocs.io/en/latest/content/tools/genomecov.html
2:https://cloud.tencent.com/developer/article/1455755
网友评论