怎么使用trim_galore对NGS数据进行质量过滤

怎么使用trim_galore对NGS数据进行质量过滤，相信很多没有经验的人对此束手无策，为此本文总结了问题出现的原因和解决方法，通过这篇文章希望你能解决这个问题。

创新互联专注于浦口网站建设服务及定制，我们拥有丰富的企业做网站经验。 热诚为您提供浦口营销型网站建设，浦口网站制作、浦口网页设计、浦口网站官网定制、微信小程序定制开发服务，打造浦口网络公司原创品牌,更为您提供浦口网站排名全网营销落地服务。

cutadapt软件可以对NGS数据进行质量过滤，FastQC软件可以查看NGS数据的质量分布，trim_galore将这两个软件封装到一起，使用起来更加的方便。

该软件会对数据进行以下4步处理

1. 去除reads 3’端的低质量碱基

illumina平台的测序数据，通常3’端质量较差。trim_galore首先会过滤掉3’端的低质量碱基，本质上是调用了cutadapt的质量过滤算法。下图是过滤前后碱基质量的分布图

可以看到，过滤掉低质量碱基后，序列的整体质量显著提高。

2. 去除adapter序列

过滤掉低质量的碱基之后，trim_galore会调用cutadapt在reads的3’端查找adapter 序列并去除。通常情况下，我们需要指定对应的adapter序列，如果没有指定的化，trim_galore会自动查找以下3种类型的adapter

Illumina:   AGATCGGAAGAGC
Small RNA:  TGGAATTCTCGG
Nextera:    CTGTCTCTTATA

默认读取前一百万条序列，通过这一百万条序列判断adapter属于上述三种的哪一种，然后进行去除。如果你不希望软件自动判断，也可以通过--illumina, --nextera, --small_rna参数指定对应的adapter类型。

3. 去除长度太短的序列

经过上述两步处理之后，有可能剩余的序列长度很短，这部分短序列也会被去除。默认情况下，如果序列长度少于20bp, 这条序列会被丢掉。

4. 其它过滤

对于所有的输入序列，以上3个步骤是肯定会执行的。除此之，trim_galore还支持一些其他的过滤措施，以满足个性化的需求。

hardtrim5参数用于从序列的3’端切除碱基，示意如下

before:         CCTAAGGAAACAAGTACACTCCACACATGCATA
--hardtrim5 20: CCTAAGGAAACAAGTACACT

通过hardtrim5参数可以将序列截取成固定长度。与之对应的，还有一个hardtrim3参数，从序列的5’端切除碱基，示意如下

before:         CAAATGTTATTTTTAAGAAAATGGAAAAT
--hardtrim3 20:          TTTTTAAGAAAATGGAAAAT

软件的安装也很方便，首先需要确保cutadapt和fastqc这两个软件已经安装，并且可执行文件位于PAH环境变量定义的路径种。然后下载trim_galore的源代码包，解压即可，代码如下

wget https://github.com/FelixKrueger/TrimGalore/archive/0.5.0.tar.gz
tar xzvf 0.5.0.tar.gz

在软件的安装目录有一个名为trim_galore的可执行文件。

对于单端测序数据，基本用法如下

trim_galore  --quality 20 -a AGATCGGAAGAGC --length 20  -o out_dir  input.fq

对于双端测序数据，基本用法如下

trim_galore --paired --quality 20 -a AGATCGGAAGAGC   -a2 AGATCGGAAGAGC --length 20  -o out_dir  R1.fq.gz R2.fq.gz

看完上述内容，你们掌握怎么使用trim_galore对NGS数据进行质量过滤的方法了吗？如果还想学到更多技能或想了解更多相关内容，欢迎关注创新互联行业资讯频道，感谢各位的阅读！