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

原理

讲一下

这个比较重要,就不瞎比比了,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就明白了

相关推荐
Turnsole_y44 分钟前
pycharm自动化测试初始化
python·selenium
郭式云源生法则1 小时前
归档及压缩、重定向与管道操作和综合使用,find精确查找、find处理查找结果、vim高级使用、vimdiff多文件使用
linux·运维·服务器
一张假钞1 小时前
Ubuntu 24.04 安装 Jenkins
linux·ci/cd·jenkins
tuokuac2 小时前
查看你电脑上某个端口正在被哪个进程占用
linux
weixin-a153003083162 小时前
[数据抓取-1]beautifulsoup
开发语言·python·beautifulsoup
AI量化投资实验室2 小时前
15年122倍,年化43.58%,回撤才20%,Optuna机器学习多目标调参backtrader,附python代码
人工智能·python·机器学习
倔强青铜三2 小时前
苦练Python第67天:光速读取任意行,linecache模块解锁文件处理新姿势
人工智能·python·面试
MANONGMN2 小时前
Linux 通配符与正则表达式(含实战案例+避坑指南)
linux·运维·正则表达式
带土12 小时前
18 .shell编程-正则表达式
linux·正则表达式
爱吃甜品的糯米团子2 小时前
Linux 学习笔记之进程管理、网络基础与常用软件安装
linux·网络·学习