我有多个 .bam 文件(与 bwa 对齐),我想生成一个子集 .bam 文件,其中仅包含 50bp 以上的读取(读取长度,而不是对齐长度)。有没有办法从 bam 文件中做到这一点?还是我应该回去对 fastq 进行子集化,然后重新对齐?
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)