首页 > 解决方案 > 对 bam 文件进行二次采样以仅包含长度超过 X 的读取

问题描述

我有多个 .bam 文件(与 bwa 对齐),我想生成一个子集 .bam 文件,其中仅包含 50bp 以上的读取(读取长度,而不是对齐长度)。有没有办法从 bam 文件中做到这一点?或者我应该回去对fastq进行子集化,然后重新对齐?

标签: alignmentsubsetbam

解决方案


推荐阅读