一、MultiQC介紹
NGS技術(shù)的進(jìn)步催生了新的實(shí)驗(yàn)設(shè)計秧均、分析類型和極高通量測序數(shù)據(jù)的生成锉罐。對于這些數(shù)據(jù)的質(zhì)量評估，每一步分析結(jié)果的評估是后續(xù)結(jié)果可信度的衡量和保障赏枚。不少生信工具都可以給樣品生成一個評估結(jié)果亡驰，如FastQC、Qualimap 和RSeQC等 (39個轉(zhuǎn)錄組分析工具饿幅，120種組合評估)凡辱。但是這時又出現(xiàn)了一個難題，那就是幾乎所有的質(zhì)控工具都是針對單個樣本生成一個報告栗恩，這就要求用戶自己去逐一查找各個QC結(jié)果透乾，這無疑是個十分耗時、重復(fù)又復(fù)雜的事，而且還不能快速看出所有樣本的異同续徽。
那能否把所有質(zhì)控結(jié)果整合在一起呢蚓曼？可以自己寫程序造輪子（我們之前就是這么做的）。但現(xiàn)在有了MultiQC钦扭，基于Python的小工具很好地解決了這個繁瑣的事情纫版，其強(qiáng)大的功能主要體現(xiàn)在以下三個方面：
1）能將測序數(shù)據(jù)的多個QC結(jié)果整合成一個HTLM網(wǎng)頁交互式報告，同時也能導(dǎo)出pdf文件客情；
2）支持多種分析類型的質(zhì)控結(jié)果查看其弊，如：RNAseq、Whole-Genome Seq膀斋、Bisulfite Seq梭伐、Hi-C和MultiQC_NGI；
3）支持整合68種軟件分析的結(jié)果仰担，而且支持的軟件還在持續(xù)增加糊识，也可以自己寫作一個插件，具體見下圖摔蓝。
二赂苗、安裝MultiQC
依賴python2.7+, 3.4+ 或者 3.5+
# pip安裝
pip install git+https://github.com/ewels/MultiQC.git #Installation with pip
# conda安裝
conda install -c bioconda multiqc # Installing with conda
三、運(yùn)行MultiQC
直接指定MultiQC要分析的文件路徑即可贮尉，若數(shù)據(jù)在當(dāng)前目錄下輸入multiqc .即可拌滋。
multiqc .
multiqc data/
multiqc data/ ../proj_one/analysis/ /tmp/results
multiqc data/*_fastqc.zip
multiqc data/sample_1*
使用--ignore忽略掉某些文件
multiqc . --ignore *_R2*
multiqc . --ignore run_two/
multiqc . --ignore */run_three/*/fastqc/*_R2.zip
四、MultiQC報告解讀（以RNA-Seq數(shù)據(jù)為例）
1. General Statistics
每一個樣本reads數(shù)量猜谚、比對層面的質(zhì)量評估整合統(tǒng)計表败砂，點(diǎn)擊Configure Columns可以選擇顯示或不顯示某些項。點(diǎn)擊Plot可以繪圖魏铅。
點(diǎn)擊Configure Columns選擇展示哪些項
點(diǎn)擊Plot可以對任意兩種屬性的評估結(jié)果做交互式二維圖昌犹，若各樣本均一性好，散點(diǎn)會比較集中览芳，反之會出現(xiàn)某些離散的點(diǎn)祭隔，這樣方便查看某些指標(biāo)異常的離群樣本。
2. featureCounts
利用featureCounts工具計算每個基因外顯子的reads數(shù)的結(jié)果展示路操。featureCounts不僅可以支持gene的定量疾渴，也支持exon, gene bodies, genomic bins, chromsomal locations的定量。功能類似的軟件是HTSeq屯仗。
軟件官網(wǎng)：http://bioinf.wehi.edu.au/featureCounts/
3. STAR
基于STAR比對工具的分析結(jié)果搞坝，STAR會將沒有paired mapping的reads都剔除，避免single reads比對到基因組上魁袜；并且STAR對lower-quality（采用more soft-clipped和錯配堿基）比對有較高的容忍度桩撮。
軟件官網(wǎng)：https://github.com/alexdobin/STAR
更多分析工具比較見：轉(zhuǎn)錄組分析工具大比拼
4. Cutadapt
用cutadapt軟件來對雙端測序數(shù)據(jù)進(jìn)行預(yù)處理敦第，去除接頭和低質(zhì)量堿基。
對測序數(shù)據(jù)進(jìn)行過濾時cutadapt對測序數(shù)據(jù)進(jìn)行識別店量、剪切并去除adapters, primers , poly_A等序列芜果，移除被adapter污染的reads部分（指由于插入片段長度不夠，測序儀讀到的測序引物等序列）融师。具體見NGS基礎(chǔ) - 高通量測序原理右钾。
軟件官網(wǎng)：https://cutadapt.readthedocs.io/en/stable/
5. FastQC
MultiQC將fastqc工具分析得到的10個結(jié)果分別整合成一個模塊，集中查看旱爆。
軟件官網(wǎng)：http://www.bioinformatics.babraham.ac.uk/projects/fastqc/
具體的關(guān)于FastQC報告解讀可以見歷史推文：NGS基礎(chǔ) - FASTQ格式解釋和質(zhì)量評估
MultiQC的可定制性也比較強(qiáng)舀射，更多功能值得進(jìn)一步探索。
