0

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

4

0 回答 0