bam格式转bigWig

生信编程日常

发布于 2020-07-22 00:24:47

7.4K0

文章被收录于专栏：生物信息学、python、R、linux生物信息学、python、R、linux

当我们需要在UCSC browser中可视化的时候，将bam文件转化为bigwig文件会更加方便。而deeptools中的bamCoverage可以方便的实现这个功能。

1. 安装

conda或者pip安装：

conda install -c bioconda deeptools
# 或者
pip install deeptools

github源代码安装：

 git clone https://github.com/deeptools/deepTools.git
wget https://github.com/deeptools/deepTools/archive/1.5.12.tar.gz
tar -xzvf
python setup.py install --prefix /User/Tools/deepTools2.0

一定要注意一下的依赖包必须满足要求：

如果不满足要求的话，bam2Coverage可能会引发如下错误： "The XXX file does not have BAM or CRAM format"； "ModuleNotFoundError: No module named 'numpy.core._multiarray_umath'"等等。

2. bam2Coverage

必需参数：

需指定输入文件，输出文件名。输出文件格式默认为bigwig，如果想输出bedgraph文件需指定。

非必需参数：

此外还有Read coverage normalization options和Read coverage normalization options，有具体需要可以查阅。

ChIP-seq的例子：

bamCoverage --bam a.bam -o a.SeqDepthNorm.bw \
    --binSize 10
    --normalizeUsing RPGC
    --effectiveGenomeSize 2150570000
    --ignoreForNormalization chrX
    --extendReads

RNA-seq的例子：

bamCoverage -b a.bam -o a.bw

当我们需要分离正负链时（单端）：

bamCoverage -b a.bam -o forward.bw --filterRNAstrand forward
bamCoverage -b a.bam -o reverse.bw --filterRNAstrand reverse

# 如果在2.2版本之前：
bamCoverage -b a.bam -o a.fwd.bw --samFlagExclude 16
bamCoverage -b a.bam -o a.rev.bw --samFlagInclude 16

相当于我们做了如下操作：

samtools view -@ 8 a.bam | \
awk -F '\t' '$2==0' | \
bamCoverage -b a.bam -o a.fwd.bw

本文参与腾讯云自媒体同步曝光计划，分享自作者个人站点/博客。

如有侵权请联系 cloudcommunity@tencent.com 删除

编程算法

本文分享自作者个人站点/博客前往查看

如有侵权，请联系 cloudcommunity@tencent.com 删除。

本文参与腾讯云自媒体同步曝光计划，欢迎热爱写作的你一起参与！

编程算法

登录后参与评论

0 条评论

热度