怎么使用trim_galore對NGS數(shù)據(jù)進行質(zhì)量過濾，相信很多沒有經(jīng)驗的人對此束手無策，為此本文總結了問題出現(xiàn)的原因和解決方法，通過這篇文章希望你能解決這個問題。

創(chuàng)新互聯(lián)專注于浦口網(wǎng)站建設服務及定制，我們擁有豐富的企業(yè)做網(wǎng)站經(jīng)驗。 熱誠為您提供浦口營銷型網(wǎng)站建設，浦口網(wǎng)站制作、浦口網(wǎng)頁設計、浦口網(wǎng)站官網(wǎng)定制、微信小程序定制開發(fā)服務，打造浦口網(wǎng)絡公司原創(chuàng)品牌,更為您提供浦口網(wǎng)站排名全網(wǎng)營銷落地服務。

cutadapt軟件可以對NGS數(shù)據(jù)進行質(zhì)量過濾，F(xiàn)astQC軟件可以查看NGS數(shù)據(jù)的質(zhì)量分布，trim_galore將這兩個軟件封裝到一起，使用起來更加的方便。

該軟件會對數(shù)據(jù)進行以下4步處理

1. 去除reads 3’端的低質(zhì)量堿基

illumina平臺的測序數(shù)據(jù)，通常3’端質(zhì)量較差。trim_galore首先會過濾掉3’端的低質(zhì)量堿基，本質(zhì)上是調(diào)用了cutadapt的質(zhì)量過濾算法。下圖是過濾前后堿基質(zhì)量的分布圖

可以看到，過濾掉低質(zhì)量堿基后，序列的整體質(zhì)量顯著提高。

2. 去除adapter序列

過濾掉低質(zhì)量的堿基之后，trim_galore會調(diào)用cutadapt在reads的3’端查找adapter 序列并去除。通常情況下，我們需要指定對應的adapter序列，如果沒有指定的化，trim_galore會自動查找以下3種類型的adapter

Illumina:   AGATCGGAAGAGC
Small RNA:  TGGAATTCTCGG
Nextera:    CTGTCTCTTATA

默認讀取前一百萬條序列，通過這一百萬條序列判斷adapter屬于上述三種的哪一種，然后進行去除。如果你不希望軟件自動判斷，也可以通過--illumina, --nextera, --small_rna參數(shù)指定對應的adapter類型。

3. 去除長度太短的序列

經(jīng)過上述兩步處理之后，有可能剩余的序列長度很短，這部分短序列也會被去除。默認情況下，如果序列長度少于20bp, 這條序列會被丟掉。

4. 其它過濾

對于所有的輸入序列，以上3個步驟是肯定會執(zhí)行的。除此之，trim_galore還支持一些其他的過濾措施，以滿足個性化的需求。

hardtrim5參數(shù)用于從序列的3’端切除堿基，示意如下

before:         CCTAAGGAAACAAGTACACTCCACACATGCATA
--hardtrim5 20: CCTAAGGAAACAAGTACACT

通過hardtrim5參數(shù)可以將序列截取成固定長度。與之對應的，還有一個hardtrim3參數(shù)，從序列的5’端切除堿基，示意如下

before:         CAAATGTTATTTTTAAGAAAATGGAAAAT
--hardtrim3 20:          TTTTTAAGAAAATGGAAAAT

軟件的安裝也很方便，首先需要確保cutadapt和fastqc這兩個軟件已經(jīng)安裝，并且可執(zhí)行文件位于PAH環(huán)境變量定義的路徑種。然后下載trim_galore的源代碼包，解壓即可，代碼如下

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

在軟件的安裝目錄有一個名為trim_galore的可執(zhí)行文件。

對于單端測序數(shù)據(jù)，基本用法如下

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

對于雙端測序數(shù)據(jù)，基本用法如下

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

看完上述內(nèi)容，你們掌握怎么使用trim_galore對NGS數(shù)據(jù)進行質(zhì)量過濾的方法了嗎？如果還想學到更多技能或想了解更多相關內(nèi)容，歡迎關注創(chuàng)新互聯(lián)行業(yè)資訊頻道，感謝各位的閱讀！

分享題目：怎么使用trim_galore對NGS數(shù)據(jù)進行質(zhì)量過濾
標題鏈接：http://aaarwkj.com/article12/jpoidc.html

成都網(wǎng)站建設公司_創(chuàng)新互聯(lián)，為您提供品牌網(wǎng)站設計、做網(wǎng)站、微信公眾號、網(wǎng)站營銷、App開發(fā)、標簽優(yōu)化

廣告

聲明：本網(wǎng)站發(fā)布的內(nèi)容（圖片、視頻和文字）以用戶投稿、用戶轉(zhuǎn)載內(nèi)容為主，如果涉及侵權請盡快告知，我們將會在第一時間刪除。文章觀點不代表本網(wǎng)站立場，如需處理請聯(lián)系客服。電話：028-86922220；郵箱：631063699@qq.com。內(nèi)容未經(jīng)允許不得轉(zhuǎn)載，或轉(zhuǎn)載時需注明來源： 創(chuàng)新互聯(lián)