SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
cut 5' and 3' ends of paired-ends reads sebl Bioinformatics 4 12-08-2015 01:48 AM
Trinity denovo assemply error dinesh21.ksd Bioinformatics 0 07-15-2015 09:42 PM
DeNovo ChIPseq and Trinity blcUT General 3 11-13-2014 09:18 PM
Results of denovo assembly from trinity and velvet. vishwesh Bioinformatics 1 01-05-2014 10:03 PM

Reply
 
Thread Tools
Old 02-02-2016, 10:22 AM   #1
nimrod337
Member
 
Location: Wisconsin

Join Date: Jun 2013
Posts: 22
Default Trinity denovo transcript multiple reps paired ends

Sorry if this is in the wrong place. There are quite a few spots this could have gone. I've done some poking around, and haven't quite found an answer to my question, which is:

I am attempting to assemble a denovo transcript using paired end illumina 3000 reads from 4 different treatments, with 3 biological reps per treatment.

To do an assembly de novo, is it kosher (or better yet, is it required) to join all the reps from all the treatments into one "super-reads" file for R1 and another for the R2? If so, do I have to cat them in the same order to allow the correct R1 read to pair with its correct R2 read? Is there anyway to avoid this, and just give trinity all of the read files individually?
nimrod337 is offline   Reply With Quote
Old 02-02-2016, 11:04 AM   #2
nimrod337
Member
 
Location: Wisconsin

Join Date: Jun 2013
Posts: 22
Default

I figured it out. I can add multiple files, so long as they are separated by a comma, not a space. So my new question is: can I load paired end reads using --left, --right, and also include my singleton reads in the same run?

How would that effect my --SS_lib_type <string> argument?
nimrod337 is offline   Reply With Quote
Old 02-03-2016, 10:00 AM   #3
westerman
Rick Westerman
 
Location: Purdue University, Indiana, USA

Join Date: Jun 2008
Posts: 1,104
Default

I haven't done this but the answer as per the authors (via the Trinity mailing list a year or so ago) is to add the singleton reads to the end of the left reads.
westerman is offline   Reply With Quote
Old 02-04-2016, 10:09 AM   #4
nimrod337
Member
 
Location: Wisconsin

Join Date: Jun 2013
Posts: 22
Default

So I can just add the singleton files to the end of my --left file1.fastq,file2.fastq,etc.....,singleton1.fastq,singleton2.fastq ? This will result in 12 fastq files after --right, and 24 files after --left. I'll give it a shot, and see what happens.

Thanks for the reply!
nimrod337 is offline   Reply With Quote
Reply

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 02:31 PM.


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