玩转基因组浏览器之tdf文件

本文介绍了如何通过igvtools将bam文件转换为tdf文件,以减小内存占用并快速查看测序深度信息。tdf是IGV推荐的二进制格式,与bedgraph相似,但文件大小更小。通过调整窗口大小,可以在空间和分辨率之间做出权衡。文章还展示了tdf文件与wig、bam文件的比较,并提供了将tdf转为bedgraph的方法。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

欢迎关注”生信修炼手册”!

将bam文件导入IGV之后,可以直观的查看测序深度的分布情况, 但是直接导入bam文件会占用比较大的内存,如果只是想要查看测序深度信息,有很多其他的代替方案。

tdf是IGV官方推荐的一种二进制格式,类似bedgraph格式,用窗口的方式来记录测序深度信息。相比bam文件,tdf文件会小很多,导入和查看也更快速。可以通过igvtools来生成tdf文件,命令如下

igvtools count input.bam out.tdf hg19.chrom.sizes

需要三个参数,第一个参数为输入文件,支持bam, sam等格式,第二个参数为输出文件,支持tdf和wig两种格式,第三个参数为基因组的ID或者保存染色体大小的chrom.sizes文件。在igvtools的安装目录,有个lib/genomes文件夹,保存了很多物种的chrom.sizes文件

当你提供了基因组ID时,软件会自动在该目录下查找对应的文件,默认的ID是hg18, 当然你也可以自己准备这个文件,内容如下所示

\t分隔的两列,第一列为染色体名称,第二列为染色体长度。在我的测试中,bam文件大小为1.3G, 输出的tdf文件大小为17M, 而wig文件为59M,转换为二进制的bigwig文件大小为25M。可以看到,tdf的文件大小最小。输出的tdf文件可以直接导入IGV进行查看,示意如下

上图的wig,tdf来自同一个bam文件,从峰型可以看出,三者基本一致。只所以说基本一致,是因为wig和tdf都是以窗口的方式来统计测序深度的,而bam则以单个碱基为单位。在IGV中进一步放大,可以看到下图

wig和tdf不再是平滑的峰,变成了一个个很宽的柱子。相比之下,bam文件依然是以碱基为单位的柱子。通过下面的命令,可以将tdf文件转换为bedgraph这种纯文本格式

igvtools tdftobedgraph  input.tdf out.bedgraph

内容如下

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值