重测序之快速比对(速度比谁都快哦~)

原理

讲一下

这个比较重要,就不瞎比比了,screen建立作业窗口,再分别运行比对个体

都不喜欢讲想法,那就看原理。

python 复制代码
import os
##比对之前如果不建立索引,那你猜我会怎么骂你呢?
def create_screen_window(screen_name):
    os.system(f"screen -dmS {screen_name}")
def execute_code_in_screen(screen_name, code):
    os.system(f"screen -S {screen_name} -X stuff '{code}\n'")
    os.system(f"screen -S {screen_name} -X stuff 'exit\n'")
for i in range(52,82):
    screen_name = "SRR131783"+str(i)
    create_screen_window(screen_name)
    code_to_execute = "bwa-mem2 mem sidai.fa SRR131783"+str(i)+"_1.fastq.gz SRR131783"+str(i)+"_2.fastq.gz > SRR131783"+str(i)+".sam"
    execute_code_in_screen(screen_name, code_to_execute)

基本上每一个参数都是重点没有可浪费的,每一步都是重点。充分理解screen就明白了

相关推荐
CodeCraft Studio31 分钟前
CAD文件处理控件Aspose.CAD教程:使用 Python 将绘图转换为 Photoshop
python·photoshop·cad·aspose·aspose.cad
Python×CATIA工业智造2 小时前
Frida RPC高级应用:动态模拟执行Android so文件实战指南
开发语言·python·pycharm
onceco3 小时前
领域LLM九讲——第5讲 为什么选择OpenManus而不是QwenAgent(附LLM免费api邀请码)
人工智能·python·深度学习·语言模型·自然语言处理·自动化
狐凄4 小时前
Python实例题:基于 Python 的简单聊天机器人
开发语言·python
悦悦子a啊5 小时前
Python之--基本知识
开发语言·前端·python
2401_826097625 小时前
JavaEE-Linux环境部署
java·linux·java-ee
(:满天星:)6 小时前
第31篇:块设备与字符设备管理深度解析(基于OpenEuler 24.03)
linux·运维·服务器·网络·centos
爱莉希雅&&&6 小时前
shell编程之awk命令详解
linux·服务器·git
笑稀了的野生俊6 小时前
在服务器中下载 HuggingFace 模型:终极指南
linux·服务器·python·bash·gpu算力
Naiva6 小时前
【小技巧】Python+PyCharm IDE 配置解释器出错,环境配置不完整或不兼容。(小智AI、MCP、聚合数据、实时新闻查询、NBA赛事查询)
ide·python·pycharm