前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >用脚本同时对多个fq文件进行数据质量检测

用脚本同时对多个fq文件进行数据质量检测

作者头像
戈贝尔光和热
发布2018-12-27 14:54:04
1.3K0
发布2018-12-27 14:54:04
举报
文章被收录于专栏:HUBU生信HUBU生信

当有多个fq文件要进行数据质量检测时,我们可通过建立一个脚本执行文件,执行该脚本,可同时批量对fq文件进行检测。 1.首先我们要建立一个脚本文件(后缀名为.sh)

代码语言:javascript
复制
lihan@lihan-virtual-machine:~/seqs/trim_out$ nano fastqc.sh

2.然后输入以下脚本:

代码语言:javascript
复制
#!/bin/sh 
for i in output* 
do 
echo $i 
fastqc $i 
done

Ctrl+X保存退出。这里output*是指以output开头的文件名。 3.执行脚本:

代码语言:javascript
复制
lihan@lihan-virtual-machine:~/seqs/trim_out$ sh fastqc.sh

就会自动检测所有以output开头的fq文件,运行结果如下:

代码语言:javascript
复制
output_forward_paired.fq.gz
Started analysis of output_forward_paired.fq.gz
Approx 5% complete for output_forward_paired.fq.gz
Approx 10% complete for output_forward_paired.fq.gz
Approx 15% complete for output_forward_paired.fq.gz
Approx 20% complete for output_forward_paired.fq.gz
Approx 25% complete for output_forward_paired.fq.gz
Approx 30% complete for output_forward_paired.fq.gz
Approx 35% complete for output_forward_paired.fq.gz
Approx 40% complete for output_forward_paired.fq.gz
Approx 45% complete for output_forward_paired.fq.gz
Approx 50% complete for output_forward_paired.fq.gz
Approx 55% complete for output_forward_paired.fq.gz
Approx 60% complete for output_forward_paired.fq.gz
Approx 65% complete for output_forward_paired.fq.gz
Approx 70% complete for output_forward_paired.fq.gz
Approx 75% complete for output_forward_paired.fq.gz
Approx 80% complete for output_forward_paired.fq.gz
Approx 85% complete for output_forward_paired.fq.gz
Approx 90% complete for output_forward_paired.fq.gz
Approx 95% complete for output_forward_paired.fq.gz
Analysis complete for output_forward_paired.fq.gz
output_forward_unpaired.fq.gz
Started analysis of output_forward_unpaired.fq.gz
Approx 20% complete for output_forward_unpaired.fq.gz
Approx 40% complete for output_forward_unpaired.fq.gz
Approx 60% complete for output_forward_unpaired.fq.gz
Approx 80% complete for output_forward_unpaired.fq.gz
Analysis complete for output_forward_unpaired.fq.gz
output_reverse_paired.fq.gz
Started analysis of output_reverse_paired.fq.gz
Approx 5% complete for output_reverse_paired.fq.gz
Approx 10% complete for output_reverse_paired.fq.gz
Approx 15% complete for output_reverse_paired.fq.gz
Approx 20% complete for output_reverse_paired.fq.gz
Approx 25% complete for output_reverse_paired.fq.gz
Approx 30% complete for output_reverse_paired.fq.gz
Approx 35% complete for output_reverse_paired.fq.gz
Approx 40% complete for output_reverse_paired.fq.gz
Approx 45% complete for output_reverse_paired.fq.gz
Approx 50% complete for output_reverse_paired.fq.gz
Approx 55% complete for output_reverse_paired.fq.gz
Approx 60% complete for output_reverse_paired.fq.gz
Approx 65% complete for output_reverse_paired.fq.gz
Approx 70% complete for output_reverse_paired.fq.gz
Approx 75% complete for output_reverse_paired.fq.gz
Approx 80% complete for output_reverse_paired.fq.gz
Approx 85% complete for output_reverse_paired.fq.gz
Approx 90% complete for output_reverse_paired.fq.gz
Approx 95% complete for output_reverse_paired.fq.gz
Analysis complete for output_reverse_paired.fq.gz
output_reverse_unpaired.fq.gz
Started analysis of output_reverse_unpaired.fq.gz
Analysis complete for output_reverse_unpaired.fq.gz

这时就会在目录下生成.html文件:

代码语言:javascript
复制
lihan@lihan-virtual-machine:~/seqs/trim_out$ ls
fastqc.sh                            output_reverse_paired_fastqc.html
output_forward_paired_fastqc.html    output_reverse_paired_fastqc.zip
output_forward_paired_fastqc.zip     output_reverse_paired.fq.gz
output_forward_paired.fq.gz          output_reverse_unpaired_fastqc.html
output_forward_unpaired_fastqc.html  output_reverse_unpaired_fastqc.zip
output_forward_unpaired_fastqc.zip   output_reverse_unpaired.fq.gz
output_forward_unpaired.fq.gz

4.之后可通过WinSCP或MobalXterm,下载到 windows,打开fastqc.html进行浏览

全文结束,欢迎在评论区讨论~

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

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

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

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

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