SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Tophat/Cufflinks treating paired end reads as single Helical Bioinformatics 5 05-24-2016 09:31 AM
Convert paired-end, longer reads to single-end, shorter reads mlodato Bioinformatics 7 01-30-2016 06:47 AM
Creating psuedo paired-end sequencing reads from single-end reads cburke04 Bioinformatics 6 01-14-2015 07:10 PM
Uniquely mapped reads and difference for single end and paired end reads gene_x Bioinformatics 2 01-13-2015 01:55 PM
TopHat -paired end vs single end reads adarshjose RNA Sequencing 10 06-12-2012 07:15 PM

Reply
 
Thread Tools
Old 12-10-2015, 02:57 AM   #1
JonB
Member
 
Location: Norway

Join Date: Jan 2010
Posts: 83
Default Mixing paired-end and single-end reads in Tophat: Do I have to reverse the SE-reads f

This is a cross-post from the Tuxedo Google group: https://groups.google.com/forum/#!to...rs/BtI3_fp_n-o

And a follow-up from a previous post


Hope my question makes sense. I have paired-end reads, but after trimming some reads are left as single-end. I know that tophat can take a mix of paired-end and single-end reads, but do I have to reverse the single-end reads from the rightmost pair (usually called the R2-file or something similar) if they have no mate?

Thanks,

Jon
JonB is offline   Reply With Quote
Old 12-10-2015, 05:08 AM   #2
SylvainL
Senior Member
 
Location: Geneva

Join Date: Feb 2012
Posts: 174
Default

Not sure I understood, but all the reads are from 5' to 3' so you should not have to reverse the SE...
SylvainL is offline   Reply With Quote
Old 12-11-2015, 12:28 AM   #3
JonB
Member
 
Location: Norway

Join Date: Jan 2010
Posts: 83
Default

The paired reads are sequenced like this
R1------> <------R2
and I guess tophat takes this into account when it is fed paired reads. But when it only gets the R2 reads as single end, without knowing that they were part of a pair, will it still be ok?
JonB is offline   Reply With Quote
Old 12-11-2015, 01:37 AM   #4
SylvainL
Senior Member
 
Location: Geneva

Join Date: Feb 2012
Posts: 174
Default

If your R1 read is mapping to the forward strand, then the R2 will map to the reverse one...

Actually, the only interest to reverseComplement (not only reverse) would be if you did stranded RNAseq... but even in this case, if you map separately only the R2 reads, you can simply change the strand afterwards...
SylvainL is offline   Reply With Quote
Old 12-11-2015, 01:53 AM   #5
JonB
Member
 
Location: Norway

Join Date: Jan 2010
Posts: 83
Default

I should have mentioned that I do have stranded sequencing. So I guess I will map the R2 reads separately.
Thanks!
JonB is offline   Reply With Quote
Old 12-11-2015, 02:08 AM   #6
SylvainL
Senior Member
 
Location: Geneva

Join Date: Feb 2012
Posts: 174
Default

Ok, then it makes sense to reverseComplement (or easiest, as mentioned, is to change the strand after the mapping for R2 alone...)...
Actually, maybe to confuse you a bit more, the R2 (without reverseComplementing them, or witghout changing the strand after mapping), will be on the same strand than your gene (the coding sequence). R1 will be in the opposite one
SylvainL is offline   Reply With Quote
Reply

Tags
mapping, paired-end, single-end reads, tophat2

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


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