Skip to content

生成bigwig

By Hongjiang, Feb 27, 2022

.bam.bedGraph等文件转换为 .bigwig文件,缩小文件体积,以便在浏览器中加载。

使用bamCoverage转换

Deeptools官方文档

安装Deeptools

一般使用conda安装deeptools。

conda activate $YOUR_ENV
conda install deeptools

在使用此程序处理.bam文件前,需要使用samtools index建立索引。

Samtools index

samtools index -b reads.bam

一个最简单的命令为:

bamCoverage -b reads.bam -o coverage.bw

必须参数

-b,也可写作--bam,后接输入的.bam文件。

若不指定输出参数,则为stdout,使用 > 来处理。

输出参数

-o,也可写作--outFileName,后接输出的.bigWig文件。

-of,也可以写作--outFileFormat,后接输出文件格式。可选bigwigbedgraph两种,默认为bigwig

使用bedGraphToBigWig转换

UCSC官方文档

二进制程序已编译好,可以直接使用。

此程序仅支持Linux使用,UCSC同时提供MacOS版本

wget http://hgdownload.soe.ucsc.edu/admin/exe/linux.x86_64/bedGraphToBigWig

赋予执行权限

chmod a+x bedGraphToBigWig

下载基因组大小文件

bedGraphToBigWig程序需要基因组大小文件作为参考,使用以下指令下载hg19参考文件。

wget https://genome.ucsc.edu/goldenPath/help/hg19.chrom.sizes

同理,使用以下指令下载hg38参考文件。

wget https://genome.ucsc.edu/goldenPath/help/hg38.chrom.sizes

运行程序

bedGraphToBigWig程序需要输入三个变量,第一个是.bedGraph文件地址,第二个是基因组参考文件.chrom.sizes地址,第三个是.bigwig(.bw)文件输出地址。

参考指令:

./bedGraphToBigWig in.bedGraph chrom.sizes myBigWig.bw

注:.bedGraph文件需要sort -k1,1 -k2,2n

注:UCSC也提供bigWigToBedGraph软件。

UCSC软件库