重测序之快速比对(速度与服务器的配置挂钩。)

讲一下

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

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

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就明白了

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值