零基础小白笔记8 | 将bam文件转换为bw文件并进行可视化

❀前言:(书接上回)

零基础小白笔记1 | ChIP-seq原理、操作流程、分析流程

零基础小白笔记2 | 数据与样本信息处理

零基础小白笔记3 | 数据处理与质控

零基础小白笔记4 | 数据质控报告全解读

零基础小白笔记5 | 数据清洗与再次质控

零基础小白笔记6 | 使用bowtie2进行数据比对

零基础小白笔记7 | 对bam文件进行去重和过滤

一、将bam文件转换为bw文件:

1.bw文件,即bigwig文件,是一种常用的用于存储基因组测序数据的文件格式,它允许高效存储大量基因组范围的信号数据,并支持快速的数据检索与可视化。

2.使用工具:deeptools是一个用于深度分析测序数据的工具集,其中包含多个子命令,其中bamCoverage可以将BAM 文件转换成 BigWig 格式的输出文件。

conda install deptools

3.使用bamCoverage进行格式转换:

bamCoverage \
-b /your/input/file_sort_picard_filter.bam \
-o /your/output/file_sort_picard_filter.bigwig \
--binSize 10 --normalizeUsing RPKM \
--scaleFactor 1 --numberOfProcessors 6

(3.1)使用用法:

  • -b:用于指定输入的 BAM 文件;
  • -o:用于指定输出的Bigwig文件目录;
  • --binSize:用于设置每个输出 bin 的大小,即在生成 BigWig 文件时将测序覆盖度等信号数据按照每 10 个基因组单位进行划分和计算;
  • --normalizeUsing:进行归一化,确保在不同样本之间比较时考虑到测序深度的影响;
  • --scaleFactor :设置缩放因子为 1,表示不进行额外的信号缩放,保持原始信号的数值;
  • --numberOfProcessors:使用的线程数;

二、使用IGV进行可视化:

下载链接:Download IGV - IGV Desktop Application

1.田田目前使用的是windows系统 + mobaxterm远程终端连接linux服务器 + windows版本IGV;

所以bw文件在转换完成后需要先下载到本地才能在IGV上进行可视化分析;

2.使用步骤:

(2.1)从本地导入bigwig文件:

(2.2)设置参考基因组为mm10(本次实验动物模型),选择染色体及区域,最左侧轴设置显示区域的大小;

(2.3)设置纵轴大小一致:

3.详细操作步骤大家也可以自行查找相关视频或文章。

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

田不甜啊

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值