我正在尝试编写一个循环遍历目录中的 SAM 文件并用于samtools view
将这些 SAM 文件转换为 BAM 文件的 python 脚本。这就是我所拥有的,但我正在努力解决如何将 SAM 文件 (i) 输入到最后一行的 bash 命令中。
import os
import glob
for i in glob.iglob('/file/path/to/files/*.sam'):
os.system("samtools view -S -b %i > %i.bam" %i)
感谢您的帮助!