fastqc

时间:2021-07-31 12:00:58

fastqc用于查看测序数据的质量。

1.下载:

http://www.bioinformatics.babraham.ac.uk/projects/download.html#fastqc

wget http://www.bioinformatics.babraham.ac.uk/projects/fastqc/fastqc_v0.11.8.zip

2.解压配置:

unzip fastqc_v0.11.8.zip

配置:

cd /data/software/FastQC

chmod +x fastqc  #由于fastqc是主程序,修改fastqc的使用权限

ln -s /data/software/FastQC/fastqc /usr/local/bin/fastqc

3.运行:

fastqc -h  #查看参数

#FastQC参数说明

# fastqc [-o output dir] [--(no)extract] [-f fastq|bam|sam] [-c contaminant file] seqfile1 .. seqfileN

# 主要是包括前面的各种选项和最后面的可以加入N个文件

# -o --outdir FastQC生成的报告文件的储存路径,生成的报告的文件名是根据输入来定的

# --extract 生成的报告默认会打包成1个压缩文件,使用这个参数是让程序不打包

# -t --threads 选择程序运行的线程数,每个线程会占用250MB内存,越多越快咯

# -c --contaminants 污染物选项,输入的是一个文件,格式是Name [Tab] Sequence,里面是可能的污染序列,如果有这个选项,FastQC会在计算时候评估污染的情况,并在统计的时候进行分析,一般用不到

# -a --adapters 也是输入一个文件,文件的格式Name [Tab] Sequence,储存的是测序的adpater序列信息,如果不输入,目前版本的FastQC就按照通用引物来评估序列时候有adapter的残留

# -q --quiet 安静运行模式,一般不选这个选项的时候,程序会实时报告运行的状况。

#运行

mkdir /data/hisat2/quality/

fastqc -o  /data/hisat2/quality/ -t 6 /data/hisat2/chrX_data/samples/ERR188044_chrX_1.fastq.gz