前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >批量对多个测序文件进行fastqc

批量对多个测序文件进行fastqc

作者头像
Y大宽
发布2019-05-29 09:02:24
3.6K0
发布2019-05-29 09:02:24
举报
文章被收录于专栏:Y大宽Y大宽

现在一共是728*2=1456个测序文件,需要全部进行质控。 fastqc的命令很简单,直接跟文件即可,参数里面主要用-o(输出路径)和-t(线程,一般用2或4)

代码语言:javascript
复制
$ ls SRR85182*.gz|while read id;do  nohup fastqc -t 4 $id ;done& 

上面这个代码是错的,&放错了位置,后面会说原因,但是竟然还是可以运行。 这个代码对测序文件一个个进行,全部完成是无法忍受的时间。所以需要改进。

当我google后,发现了这个命令

代码语言:javascript
复制
ls SRR85182*.gz|while read id;do  nohup fastqc  $id &;done

但奇怪的是我一直报错报错....

代码语言:javascript
复制
-bash: syntax error near unexpected token `;'

于是又google,总体是产生这种错题提示的原因上千种,没有统一解决方案。后来发现这个描述觉得有点用

The first and foremost reason why you might experience this error message is that of bad syntax in your code or you not following the exact format of the commands. Each command has a predefined format which you can see in its documentation. Several parameters are optional which others are mandatory. Furthermore, extra care should be taken for extra space, use of double quotes, and the mandatory parameters required. If any of them are missing or have been declared incorrectly, you will not be able to execute your code.

所以我就去查了nohup用法,但没发现错误。 最后我直接把&放到done的后面,“令人惊喜”的可以运行了,但又沮丧的发现仍然是一个个任务进行。

折腾了几个小时,最后求助Jimmy,被告知

&放错了位置。

我说放前面就报错,并复制了错误信息-bash: syntax error near unexpected token;'`。 然后Jimmy说,错误信息的原因是

nohup&之间需要用()括起来。 同时告知不要用ls SRR85182*这种方式选择样本不建议,建议用循环。

暂时我现在这么运行。

代码语言:javascript
复制
$ ls SRR85182*.gz | while read id;do (nohup fastqc $id &);done

应该seq是可以的。

我又去重新选择关键词去查询nohup和&的用法,没有 找到用()的相关内容。

所以也有疑惑,为什么sra到fq的时候不用()就可以呢? 什么是可以省略?

本文参与 腾讯云自媒体分享计划,分享自作者个人站点/博客。
原始发表:2019.05.28 ,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 作者个人站点/博客 前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体分享计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档