在NGS基础:测序原始数据下载一文中提到可以使用SRA-toolkit
中的命令fastq-dump
从NCBI下载原始测序数据,命令如下。
nohup fastq-dump -v --split-3 --gzip SRR5908360 &
nohup fastq-dump -v --split-3 --gzip SRR5908361 &
这个代码,给我们4个提示:
fastq-dump
不只可以转换下载好的sra
文件为fastq
文件,还可以顺带下载sra
文件。
只需提供SRR
号,就可以获得FASTQ
序列。
不需要先调用prefetch
下载,然后再转换。
其它参数解释见引用文章。
每一行命令后面&
号表示把命令放入后台运行,当前终端可以继续输入其它命令;
此处也相当于实现了一个手动并行下载多样本,配合for
可以自动并行下载。
nohup
表示让程序在终端因人为原因或网络原因断开后不挂断,适用于运行时间比较长的命令,一般与&
连用,形式如nohup 你的命令 &
(注意空格的存在)。
如果程序运行输出错误信息,则会写入当前目录下nohup.out
文件里面,供后续查看和调试。
经常会有一些培训班“拿来主义”比较严重,以上推文和生信宝典的其它推文都被发现过直接用于某些培训班的教材,但从未申请过授权,也未引用过出处。
更有甚者,盗版易生信早期培训教案和视频,用于自己的课程或在全网发布,希望大家多多举报。
言归正传,通常我们运行程序前,会有个预判,如前面那个例子,运行时间比较长,会使用nohup 我的命令 &
的形式进行运行,从而保证程序不受网络或终端异常退出的影响。
但有时也会有误判,如没想到某个程序运行了半个小时还没结束,或数据传输时网太慢,需要传输很久,这时怎么办?中止程序,然后加上nohup
再从头运行?还是有更好的办法?
下面看这个例子:马上要去吃午饭了,把文件同步到另一个服务器,饭后回来继续操作:
ysx@ehbio:~/test/Bigwig$ rsync -av * [email protected]:/tmp
[email protected]'s password:
sending incremental file list
test1Y_DK10.bw
输入密码后,发现同步速度太慢了,1
分钟只同步了1
个文件,后面还有99
个文件,待会离开后,如果网断了,终端退出,程序终止怎么办?同步不能完成,饭后怎么愉快的工作?
还好我们有下面的方案,一步步跟着操作,补救一下。
第一步,按ctrl+z
把程序挂起,操作后屏幕会出现如下提示([1]
中的1
表示命令的作业号,后面会用到):
^Z
[1]+ 已停止 rsync -av * [email protected]:/tmp
第二步(可选),用jobs
命令查看下任务状态,跟刚才的屏幕提示一致,程序被暂时终止,作业号还是1
:
ysx@ehbio:~/test/Bigwig$ jobs
[1]+ 已停止 rsync -av * [email protected]:/tmp
第三步,使用bg %1
命令把作业号为1
的任务放入后台,并从停止状态变为运行状态,相当于加了&
后接着运行。再用jobs
查看,任务状态变成了运行中
,这一步很关键。如果没有运行bg %1
则程序处于停止状态,一直不会运行,吃几顿饭都不会运行。
ysx@ehbio:~/test/Bigwig$ bg %1
[1]+ rsync -av * [email protected]:/tmp &
ysx@ehbio:~/test/Bigwig$ jobs
[1]+ 运行中 rsync -av * [email protected]:/tmp &
第四步,运行disown -h %1
,表示在终端关闭时不对作业号为1
的程序发送终止信号,外部因素将不影响程序的运行。通过ps
命令查看下任务进程 (可选)。
ysx@ehbio:~/test/Bigwig$ disown -h %1
ysx@ehbio:~/test/Bigwig$ ps -auwx | grep 'rsync'
ysx 18214 0.0 0.0 117844 1720 ? S 09:43 0:01 rsync -av *.bw [email protected]:/tmp
ysx 18215 0.1 0.0 182376 8360 ? S 09:43 0:04 ssh -l ysx 46.93.19.14 rsync --server -vlogDtpre.iLsfxC . /tmp
ysx 18340 0.0 0.0 112724 984 pts/1 S+ 10:17 0:00 grep --color=auto rsync
通过以上4步就完成了对这次操作的事后补救。以后遇到同类问题,试一试这个新方案吧!
同时还有5点提示:
例子中使用的是rsync
同步,从节省时间来看,不是一个很好的例子。
因为把命令停掉再运行一次时,已经同步完整的数据不会再同步,时间损失不会太大。
这也是使用同步命令rsync
相比于scp
的一个好处。
更多同步方式见(Linux服务器数据定期同步和备份方式。
例子中的rsync
或其它涉及两个服务器交互的命令,都需要我们人为输入登录密码,因此直接加nohup &
运行是行不通的,无法接受密码的输入。
因此通过上面这个操作先在前台启动运行、输入密码,再放入后台不挂断运行。
从这个角度看,是一个不错的例子。
当然解决这个问题也有其它方式,具体见ssh免密码登录远程服务器。
如果程序运行时,已加了&
号,放入后台了,则只需运行jobs
获得作业号,再运行disown
不挂断即可。
程序作业号不一定都是1
,如果之前就有程序在后台运行,作业号相应的会自加。
后面用到作业号时也需要相应修改,不要刻板总用1
。
nohup
和disown
都可以使程序不挂断,可以获得一样的效果,但原理不太一致。
nohup
可以使程序忽略挂断信号(SIGHUP
)或者使程序脱离终端的控制,从而终端不能再对其发送挂断信号(SIGHUP
);
disown
则是内生于shell
,告诉shell
在终止时不对对应程序发送挂断信号(SIGHUP
)。
10000+:菌群分析 宝宝与猫狗 梅毒狂想曲 提DNA发Nature Cell专刊 肠道指挥大脑
系列教程:微生物组入门 Biostar 微生物组 宏基因组
专业技能:学术图表 高分文章 生信宝典 不可或缺的人
一文读懂:宏基因组 寄生虫益处 进化树
必备技能:提问 搜索 Endnote
文献阅读 热心肠 SemanticScholar Geenmedical
扩增子分析:图表解读 分析流程 统计绘图
16S功能预测 PICRUSt FAPROTAX Bugbase Tax4Fun
在线工具:16S预测培养基 生信绘图
科研经验:云笔记 云协作 公众号
编程模板: Shell R Perl
生物科普: 肠道细菌 人体上的生命 生命大跃进 细胞暗战 人体奥秘
为鼓励读者交流、快速解决科研困难,我们建立了“宏基因组”专业讨论群,目前己有国内外5000+ 一线科研人员加入。参与讨论,获得专业解答,欢迎分享此文至朋友圈,并扫码加主编好友带你入群,务必备注“姓名-单位-研究方向-职称/年级”。PI请明示身份,另有海内外微生物相关PI群供大佬合作交流。技术问题寻求帮助,首先阅读《如何优雅的提问》学习解决问题思路,仍未解决群内讨论,问题不私聊,帮助同行。
学习16S扩增子、宏基因组科研思路和分析实战,关注“宏基因组”
点击阅读原文,跳转最新文章目录阅读