fastqc

Posted 2021-02-04 emanlee

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