三代转录组测序数据,如果要进行多个数据合并(比如加测一次),以合并*subreads.bam为例,
1. 需要采用SMRT中的pbmerge, 不能采用samtools 的merge 。
2. 需要采用pbindex 对bam 文件进行index, 也不能采用samtools index 构建index
构建的index 文件(pbi后缀)非常的重要,因为后面的分析中,有些步骤虽然采用的是subreads.bam 文件,但是如果没有index 的pbi 文件,命令执行就会失败,而且没有提示(大坑)。
如果觉得我的文章对您有用,请随意打赏。你的支持将鼓励我继续创作!