生成bigwig
By Hongjiang, Feb 27, 2022
将.bam
,.bedGraph
等文件转换为 .bigwig
文件,缩小文件体积,以便在浏览器中加载。
使用bamCoverage转换
安装Deeptools
一般使用conda安装deeptools。
在使用此程序处理.bam
文件前,需要使用samtools index
建立索引。
一个最简单的命令为:
必须参数
-b
,也可写作--bam
,后接输入的.bam
文件。
若不指定输出参数,则为stdout,使用 >
来处理。
输出参数
-o
,也可写作--outFileName
,后接输出的.bigWig
文件。
-of
,也可以写作--outFileFormat
,后接输出文件格式。可选bigwig
、bedgraph
两种,默认为bigwig
。
使用bedGraphToBigWig转换
二进制程序已编译好,可以直接使用。
此程序仅支持Linux使用,UCSC同时提供MacOS版本。
赋予执行权限
下载基因组大小文件
bedGraphToBigWig
程序需要基因组大小文件作为参考,使用以下指令下载hg19参考文件。
同理,使用以下指令下载hg38参考文件。
运行程序
bedGraphToBigWig
程序需要输入三个变量,第一个是.bedGraph
文件地址,第二个是基因组参考文件.chrom.sizes
地址,第三个是.bigwig(.bw)
文件输出地址。
参考指令:
注:.bedGraph
文件需要sort -k1,1 -k2,2n
。
注:UCSC也提供bigWigToBedGraph
软件。