SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Need help in mapping Single end and paired-end reads Lizex RNA Sequencing 0 02-17-2013 10:48 PM
Aligning paired end read files that have different numbers of reads after trimming d1antho RNA Sequencing 2 01-21-2013 09:41 AM
Can paired-end mapping produce more reads than single-end ? warrenemmett Bioinformatics 13 03-20-2012 11:10 PM
should I reverse paired end reads before mapping? supermario Bioinformatics 2 03-06-2012 08:59 PM
BFAST mapping paired end reads. tanghz Bioinformatics 10 04-29-2011 06:29 AM

Reply
 
Thread Tools
Old 12-09-2015, 09:33 AM   #1
JonB
Member
 
Location: Norway

Join Date: Jan 2010
Posts: 83
Default Mapping paired-end reads, include also the SE-reads after trimming?

Hi,

I have trimmed paired-end reads with Trimmomatic, and some reads are therefore left as only single-end (one pair filtered out entirely). I am using tophat2 to map against a reference genome and my goal is differential gene expression analysis.

My question is if I should include also the single-end reads in the mapping? I am thinking that more reads are better, but are there any drawbacks of including both PE and SE reads?

Thanks!

Jon
JonB is offline   Reply With Quote
Old 12-09-2015, 09:44 AM   #2
Brian Bushnell
Super Moderator
 
Location: Walnut Creek, CA

Join Date: Jan 2014
Posts: 2,707
Default

Quote:
Originally Posted by JonB View Post
I am thinking that more reads are better, but are there any drawbacks of including both PE and SE reads?
More reads are better! The only drawback is that it makes the workflow slightly more complicated; depending on the aligner, you will generally need to map the paired reads and the singletons in two different passes, then merge the sam files. But I think the BWA can handle both in one pass by examining the read names.
Brian Bushnell is offline   Reply With Quote
Old 12-10-2015, 02:03 AM   #3
JonB
Member
 
Location: Norway

Join Date: Jan 2010
Posts: 83
Default

Thanks! I think also Tophat can take a mix of paire-end and single-end reads.
JonB is offline   Reply With Quote
Reply

Tags
mapping, tophat2, trimming

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:40 PM.


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