SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
Help with trim_galore output using rrbs option Dipro Bioinformatics 0 05-19-2015 01:28 PM
trim_galore cutadapt error daveiles Bioinformatics 13 02-14-2014 06:02 AM
Problem running trim_galore and cutadapt jadon Bioinformatics 4 01-14-2014 08:19 AM
How come cutadapt and trim_galore don't work at all for me?? prs321 Bioinformatics 7 06-11-2013 05:46 PM
Trouble With Trim_Galore processing paired reads Dave-bo_Baggins Bioinformatics 4 05-09-2013 09:41 AM

Reply
 
Thread Tools
Old 11-01-2015, 06:28 AM   #1
rvann
Junior Member
 
Location: Toronto

Join Date: Sep 2015
Posts: 3
Default trim_galore - when to stop trimming?

I've just used trim_galore (arguments below) to trim the adapter sequences off of fastq files from Illumina hiseq 4000 run TruSeq prep. Cutadapt seemed to work well, running it in the default mode to auto-detect adapters and remove them, as well as remove any bases with phred score < 5, but my fastQC reports for some files show that Illumina Single End PCR primer or TruSeq Adapter, Index 7, remain in certain samples (0.15 % and 0.53 %, respectively).

Do I have to run cutadapt again and feed it these specific sequences to remove? I have many samples and searching through each report for specific adapters to remove in a second cutadapt run is not ideal.

Was I not stringent enough in trimming?

Do I need to get rid of the remaining contaminants to perform differential gene expression analysis?

trim_galore --paired -q 5 -o /output/path/ --fastqc_args "--outdir /fastqc/output/path/" sample_R1.fastq.gz sample_R2.fastq.gz
rvann is offline   Reply With Quote
Old 11-01-2015, 11:22 AM   #2
Brian Bushnell
Super Moderator
 
Location: Walnut Creek, CA

Join Date: Jan 2014
Posts: 2,707
Default

I suggest including all Illumina adapters in the first pass, rather than trying to autodetect and then running a second pass to clean the stuff that was not removed. All standard Illumina adapters are contained in the file "bbmap/resources/adapters.fa" in the BBMap suite.
Brian Bushnell is offline   Reply With Quote
Old 11-08-2015, 05:00 PM   #3
fkrueger
Senior Member
 
Location: Cambridge, UK

Join Date: Sep 2009
Posts: 625
Default

Hi rvann,

sorry for replying so late. For adapter trimming it is most important to capture and remove read through adapters which obviously worked fine. If FastQC still picks up certain indexed adapter you will most likely find that these are full-length adapter sequence (which are missing the A at the start) and are thus not removed during the adapter removal process. Full length adapter sequences won't align to your genome of interest anyway and are normally not a problem. Cheers, Felix
fkrueger is offline   Reply With Quote
Reply

Tags
cutadapt, gene expression profiling, rna seq, trim adapter, trim_galore

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 12:38 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO