SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > Illumina/Solexa



Similar Threads
Thread Thread Starter Forum Replies Last Post
Introducing the Trimmomatic tonybolger Bioinformatics 189 08-16-2018 11:22 AM
Install Trimmomatic? Palgrave Bioinformatics 10 11-02-2014 05:47 AM
Problem with trimmomatic amango Bioinformatics 9 12-29-2013 09:43 AM
Trimmomatic explanation tahamasoodi Bioinformatics 11 03-20-2013 09:26 PM
question about Trimmomatic dejavu2010 Bioinformatics 4 02-27-2012 10:27 AM

Reply
 
Thread Tools
Old 11-30-2012, 10:25 PM   #1
janaahan13
Junior Member
 
Location: CA,USA

Join Date: Sep 2012
Posts: 6
Default Question about trimmomatic

Hi,

I'm a first time user to trimmomatic and I have paired-end RNA seq data. My question why do I get paired and unpaired output for each read at the output?

java -classpath <path to trimmomatic jar> org.usadellab.trimmomatic.TrimmomaticPE [-threads <threads>] [-phred33 | -phred64] [-trimlog <logFile>] <input 1> <input 2> <paired output 1> <unpaired output 1> <paired output 2> <unpaired output 2> <step 1>

I have used CLC genomics for trimming and using that I'll get two output files namely read1_trimmed.fastq and read2_trimmed.fastq. Why do I get 4 output files in trimmmomatic?
janaahan13 is offline   Reply With Quote
Old 12-01-2012, 08:45 AM   #2
janaahan13
Junior Member
 
Location: CA,USA

Join Date: Sep 2012
Posts: 6
Default

I got the answer from Federico.

the answer is: because Trimmomatic doesn't discard unpaired reads, as CLC does.
Trimmomatic trims both reads in a pair. If one of the reads is of high quality, and the second one is of low quality, only the first is kept. This may happen in cases when the two sequencing cycle groups have generally different noise levels (recently, it happened in our institute that the laser died during the second cycle, and trimming removed most of the _2 reads, allowing us to keep the _1).

Most of the other non-CLC trimming tools do the same as Trimmomatic. If you really don't want them, you can ignore _U1 and _U2. Trimmomatic just give you an extra chance to use your unpaired sequences.

Federico
janaahan13 is offline   Reply With Quote
Old 03-02-2016, 01:07 AM   #3
Rahul shelke
Junior Member
 
Location: Guwahati,India

Join Date: Mar 2015
Posts: 8
Default

For assembly which reads should be used? paired one or unpaired...
Rahul shelke is offline   Reply With Quote
Old 03-02-2016, 08:13 AM   #4
westerman
Rick Westerman
 
Location: Purdue University, Indiana, USA

Join Date: Jun 2008
Posts: 1,104
Default

Quote:
Originally Posted by Rahul shelke View Post
For assembly which reads should be used? paired one or unpaired...
You can use both. Depends on the assembler.
westerman is offline   Reply With Quote
Reply

Tags
bioinfomatics, rna-seq, trimmomatic

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 01:19 PM.


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