SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Initial QC and grooming for Illumina HiSeq2000 paired end RNAseq on Galaxy lindseykelly RNA Sequencing 5 07-30-2014 01:09 PM
Paired end Illumina prep problems? James Sample Prep / Library Generation 4 09-27-2010 11:10 AM
Problems with Modified Paired End Sample Prep Protocol vas72985 Sample Prep / Library Generation 3 07-14-2010 07:22 PM
454 Paired End orientation problems pr0t3us Bioinformatics 5 06-03-2010 12:25 AM
Problems with bwa on Q2 trimmed paired end reads curious_mapper Bioinformatics 2 05-06-2010 02:44 PM

Reply
 
Thread Tools
Old 12-13-2012, 03:16 AM   #1
quique_vzquez
Junior Member
 
Location: Salamanca, Spain

Join Date: Dec 2012
Posts: 4
Default Problems with TopHat and Paired End RNASeq

Hi, I'm very new at RNASeq Analysis and I think that I'm doing something wrong.
I've tried to align that data: http://trace.ncbi.nlm.nih.gov/Traces...mple=SRS167019.

First I execute sratoolkit to extract the fastq files with this command:

Quote:
fastq-dump --split-3 ./1_SRR097897.sra
Then I have two files: 1_SRR097897_1.fastq and 1_SRR097897_2.fastq.

The next step that I've done is to execute TopHat to generate the bam file with the reads aligned. The command was:

Quote:
../tophat2 -N 0 Pombe 1_SRR097897_1.fastq 1_SRR097897_2.fastq
Then I use the IGB to visualize the results and than shows the + and the - strands are the same



If I only use one of the _n files to generate the output I have two strands with polarity, but depending on the file to select the polarity is reversed.



I have also tried running TopHat with different parameters enabling --library-type or a GFF annotation file with identical results.
I could stay only with the results obtained through _2.fastq file but I wonder to know what I'm doing wrong, for other occasions.

Thanks for your time!!
quique_vzquez is offline   Reply With Quote
Old 12-13-2012, 03:05 PM   #2
rflrob
Member
 
Location: Berkeley, CA

Join Date: May 2010
Posts: 50
Default

I think it makes sense that the different paired ends will be on opposite strands (since you basically read the second end forwards, but on the other strand). At the level of zoom you're at in the included figures, I'm not sure that you could count on seeing the different spacing of the read1s vs the read2s (in my experiments, something like 100-300 bp). It doesn't seem like the best solution, but you could perhaps go through the BAM file and change all of the "SEQ being reverse complemented" bits in the FLAG field to get both ends pointing the same way.
rflrob is offline   Reply With Quote
Old 12-14-2012, 08:30 AM   #3
quique_vzquez
Junior Member
 
Location: Salamanca, Spain

Join Date: Dec 2012
Posts: 4
Default

Thanks for your answer.
Finaly I change the reads in one of the FASTQ files to their reverse complementary form and I can align them with polarity but I haven't understood yet why it works :S
quique_vzquez 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 08:22 AM.


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