Go Back   SEQanswers > Applications Forums > RNA Sequencing

Similar Threads
Thread Thread Starter Forum Replies Last Post
Knowledge about SOAPdenovo-trans pbrand Bioinformatics 5 02-07-2017 09:25 PM
SOAPdenovo-trans without paired reads rexxi Bioinformatics 4 03-09-2015 11:14 AM
Questions SOAPdenovo trans JManuelMonroyKuhn De novo discovery 0 11-05-2014 08:37 AM
SOAPdenovo-Trans: Seg fault kenietz De novo discovery 5 05-31-2013 05:53 PM
Comparing Soapdenovo and Soapdenovo-Trans Kate.W Bioinformatics 0 08-06-2012 03:04 PM

Thread Tools
Old 03-01-2016, 12:04 PM   #1
Rahul shelke
Junior Member
Location: Guwahati,India

Join Date: Mar 2015
Posts: 8
Default Fastx and Soapdenovo trans


I have filtered my illumina pair-end reads (Forward lib-24 million reads, Reverse Lib 24 million) using FASTX_Quality_Filter by applying the Q20 score to 90 percent of bases. (75 bp reads, insert size 200 bp)

But after filtering, I am observing around 18 million reads in a forward library and 20 million reads in a reverse library. I can see here 2 million bases difference between two libraries. Can I use above libraries for making transcriptome assembly purpose given that the number of reads are unequal?

Regards Rahul
Rahul shelke is offline   Reply With Quote
Old 03-01-2016, 12:17 PM   #2
Senior Member
Location: uk

Join Date: Mar 2009
Posts: 667

You will have problems with downstream analysis because the Forward reads and Reverse reads in your files will be out of sync, unless the reads where one of the pair has been filtered out are in separate 'unpaired_reads' files.

Last edited by mastal; 03-01-2016 at 12:22 PM.
mastal is offline   Reply With Quote
Old 03-01-2016, 01:40 PM   #3
Senior Member
Location: East Coast USA

Join Date: Feb 2008
Posts: 7,057

@Rahul: You should have used a paired-end aware trimming program to begin with.

At this stage you can use the program from BBMap suite as follows:

$ in1=r1.fq in2=r2.fq out1=fixed1.fq out2=fixed2.fq outsingle=singletons.fq
This will re-sync the order of your trimmed files and separate the single reads that remain into a separate file.
GenoMax is offline   Reply With Quote

fastx toolkit, quality assessment, rnaseq alignment, soapdenovo-trans

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 07:49 AM.

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