Go Back   SEQanswers > Applications Forums > RNA Sequencing

Similar Threads
Thread Thread Starter Forum Replies Last Post
what is a read group? wenhuang Bioinformatics 4 04-03-2014 06:46 AM
Add read groups in .bam file from 454 sequencing Maulik23 Bioinformatics 1 12-26-2012 09:03 AM
picard add read groups HGENETIC Bioinformatics 0 01-25-2012 04:19 AM
Add Readgroups to Tophat multiple sample input? Pamela Mukhopadhyay Bioinformatics 0 12-14-2011 09:55 PM
Adding read group and Platform information sbaheti Bioinformatics 2 09-25-2010 11:47 AM

Thread Tools
Old 04-17-2013, 10:13 PM   #1
Clare S
Junior Member
Location: Melbourne, Australia

Join Date: Jan 2010
Posts: 5
Default How to add in the right read group information when running Tophat?

Hi all,

I'm wondering about the best way to align reads from multiple lanes/runs when using Tophat.

As I understand it, if I have multiple fastq files for each sample, I should call tophat once per sample like

tophat ref sampleA_run1_R1.fastq,sampleA_run2_R1.fastq sampleA_run1_R2.fastq,sampleA_run2_R2.fastq

The problem is: I want each read in the resulting bam file to have read group information telling me which run/lane it came from, so that later differential expression analyses can correct for batch effects. But, if I run tophat using the above, the reads from run1 and run2 are merged together into one bam file with no distinguishing read group information.

Another alternative is to use tophat on each experimental condition separately and then merge the resulting bam files, something roughly like

tophat ref sampleA_run1_R1.fastq sampleA_run1_R2.fastq
tophat ref sampleA_run2_R1.fastq,sampleA_run2_R1.fastq
java -jar MergeSamFiles.jar INPUT=sampleA_run1.bam INPUT=sampleA_run2.bam OUTPUT=sampleA.bam

If I do this I could insert the correct read group information into each bam file before merging, so that the reads in the merged bam are distinguishable. But, if I do this, tophat is mapping each run separately, and doesn't get to use all available evidence to come up with exon-exon junctions in its iterative remapping strategy, so the mapping might not be ideal.

I have considered mapping everything once to get the exon-exon junction file, then feeding that back into tophat and mapping everything AGAIN in order to insert read group info - but this seems like it really shouldn't be necessary. Is there any way to get tophat to assign different read group information to different input fastq files? Or is there a better way?

Thanks for any help!
Clare S is offline   Reply With Quote

rna-seq, tophat

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 06:05 AM.

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