fastqc 批处理文件

当我们获取到许多的测序数据的fastq文件,我们为了方便,通过shell编程写一个批处理脚本来对许多文件进行质控。
1 首先在创建一个文件夹存放fastq文件或者fastq.gz文件,将fastq文件和fastq.gz文件放进去

mkdir ~/fastq
mv *fastq *fastq.gz ~/fastq
图一

2 创建一个results目录存放质控结果

mkdir ~/results

3 通过编辑器(vim/gedit/nano)创建fastqc.sh,其内容如下

#!/bin/bash
for fastq_file in ~/Seqs/fastq/*
do
    fastqc $fastq_file
       mv ~/Seqs/fastq/*fastqc.html ~/Seqs/fastq/*fastqc.zip ~/results
done
echo "Ok!"

(注:我的fastq文和fastq.gz文件存放在~/Seqs/fastq下)
4 写完之后退出保存,更改fastqc.sh文件的权限为可执行

chmod 777 fastqc.sh

5 运行脚本

sh fastqc.sh

或者

./fastqc.sh

运行及结果:


运行过程

质控分析结果
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容