生物软件4-Multiqc整合多个质控结果

Multiqc安装

安装conda

conda下载网址:https://docs.anaconda.com/anaconda/install/linux/

##网络下载conda
wget https://repo.anaconda.com/archive/Anaconda3-2021.05-Linux-x86_64.sh
bash Anaconda3-2021.05-Linux-x86_64.sh

安装好后重启虚拟机,使用conda --version查看


conda --version结果

使用conda安装Mulitiqc

##安装python2环境
conda create --name python2 python=2.7 -c https://mirrors.ustc.edu.cn/anaconda/cloud/bioconda/ -y
##激活python2
conda activate python2
##使用conda下载multiqc
conda install multiqc -c https://mirrors.tuna.tsinghua.edu.cn/anaconda/cloud/bioconda

Multiqc整合使用

数据准备

使用prefetch下载目标SRR文件,这里以SRR21608304、SRR19963056、SRR19963069为例

##prefetch批量下载数据
prefetch SRR21608304 SRR19963056 SRR19963069
##prefetch下载的数据在家目录下的ncbi文件夹,去ncbi文件夹查看结果
cd ncbi/public/sra/
##查看下载结果
ll
prefetch批量下载数据
查看下载结果
##使用fastq-dump解压sra文件
 fastq-dump --gzip --split-3 SRR19963069.sra SRR19963056.sra SRR21608304.sra 
解压后结果

使用fastqc进行数据质量评估

如下方命令,使用fastqc进行质量评估并使用-o参数将结果保存到家目录下MultiQC_test(要实现在家目录下建立该文件夹),当然也可以将6个压缩文件同时进行评估,为了方便观看,我将它们分开评估。

##使用fastqc进行质量评估
fastqc SRR19963056_1.fastq.gz SRR19963056_2.fastq.gz -o ~/MultiQC_test/
fastqc SRR19963069_1.fastq.gz SRR19963069_2.fastq.gz -o ~/MultiQC_test/
fastqc SRR21608304_1.fastq.gz SRR21608304_2.fastq.gz -o ~/MultiQC_test/
##切换到存放结果的目录
cd ~/MultiQC_test/
##查看结果
ll
此时生成了对应的html文件和zip文件

使用multiqc整合结果

直接在当前目录下执行multiqc .整合结果

multiqc .
multiqc .

整合结束生成了一个html文件和一个data文件


ll查看结果

查看整合结果

打开html文件或将html下载到本地并打开,即可见整合结果,在网页中可以看到分析的种种可视化结果,这里不一一展示了


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

推荐阅读更多精彩内容