FastQC及MultiQC整合使用

FastQC

FastQC 是一个基于Java写的测序数据质量评估软件。因为是用跨平台的语言Java写的，自然而然FastQC应是可以在不同系统运行的了。不过也许大多时候我们还是在Linux服务器上用的多吧。

安装

安装软件，方便的还是通过conda了，一行命令：

$ conda install -c bioconda fastqc -y

当然这需要你已经安装了anaconda的前提下。

若没有的安装anaconda的话，可以选择安装anaconda（这更方便点）或者用使用官方安装方式。
官方软件包下载：

https://www.bioinformatics.babraham.ac.uk/projects/fastqc/

官方安装说明：

https://raw.githubusercontent.com/s-andrews/FastQC/master/INSTALL.txt

这里需要提醒的是一般Linux系统是已经安装了Java的了。而若在Window下，则需要先安装好Java才可以使用FastQC。
若在有桌面环境的电脑下，FastQC其实提供了交互式的GUI操作界面，具体使用参考上面的官方安装说明

Win10 FastQC

使用

FastQC GUI操作界面看上图就发现了很简单。不过多数我们还是通过Linux命令行使用，这更灵活一点。下面以CentOS7 下的 FastQC v0.11.9为例进行简单讲解。

fastqc [-o output dir] [-t threads] [-f fastq|bam|sam] seqfile1 .. seqfileN

参数	说明
-o	指定输出目录
-t	线程数
-f	输入文件格式，默认是fastq的
seqfile	位置参数，输入文件，可输入多个文件或使用通配符匹配多个文件

例如下面例子，采用10个线程，输出地址是atac, 输入文件包裹通过*匹配多个文件作为输入

$ fastqc ATAC.1*day.r*.fq.gz -o atac/ -t 10

输出报告是html网页文件，需要传到win下方便查看。有多少个输入文件，就有多少个html报告。

若需要详细的参数说明可使用参数-h

$ fastqc -h

以及查看官方文档（也提供了一些质量报告例子供参考）：

https://www.bioinformatics.babraham.ac.uk/projects/fastqc/

fastqc

MultiQC

FastQC的输出是每一个输入文件对应一个输出报告，当有多个输入文件，会产生多个输出报告。这时一个一个查看，以及它们之间需要对比的也不是很方便。MultiQC的产生解决了这个问题，它可以将FastQC产生的多个输出报告，整合为一个，方便查看。

安装

MultiQC 是个Python包，可以通过pip下载

$ pip install multiqc

或

$ conda install -c bioconda multiqc

使用

MultiQC 有很多参数的，下面例子是个简单示例。

$ multiqc  atac/ -o all -n test

参数	说明
dir	atac/, 分析整合目录
-o	整合后的输出目录
-n	输入文件名字，默认multiqc_report

参考

https://multiqc.info/
http://www.bioinformatics.babraham.ac.uk/projects/fastqc/