SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
how to merge paired ends provided in separate files during (bfast) alignement david.tamborero Bioinformatics 3 01-03-2012 04:31 AM
Bfast alignement with paired end reads in separate files david.tamborero Bioinformatics 2 11-29-2011 07:49 AM
Mira assembler: Medium sized genomes;How to use 2 separate files for paired-end reads ndeshpan Bioinformatics 3 05-23-2011 05:59 PM
Why are Illumina paired-end SRA datasets made up of 3 FASTQ files? Bio.X2Y Illumina/Solexa 9 12-21-2010 11:36 AM

Reply
 
Thread Tools
Old 10-13-2009, 12:35 PM   #1
jyli
Member
 
Location: North Carolina

Join Date: Nov 2008
Posts: 21
Default tophat on paired files

Dear all,

I am using tophat and ran into a simple error. When I put in more pairs in the command, it failed complaining not being able to open the files.

Any comment is appreciated...

What was said in the documentation:
******************
<reads1_1[,...,readsN_1]> A comma-separated list of files containing reads in FASTQ or FASTA format. When running TopHat with paired-end reads, this should be the *_1 ("left") set of files.
<[reads1_2,...readsN_2]> A comma-separated list of files containing reads in FASTA or FASTA format. Only used when running TopHat with paired end reads, and contains the *_2 ("right") set of files. The *_2 files *MUST* appear in the same order as the *_1 files.
******************

My command:

********************************************
tophat --solexa1.3-quals -m 1 -r 130 -G /datastore/nextgenproc/projects/exonDB/custom_exonDB.04152009.hg18.exon.gff /home/jyli/Bowtie/indexes/hg18_genome
../raw_data/MCF7_CP7_00300_003_s_4_1.txt, ../raw_data/MCF7_CP7_00300_008_s_7_1.txt ../raw_data/MCF7_CP7_00300_003_s_4_2.txt, ../raw_data/MCF7_CP7_00300_008_s_7_2.txt
****************************************

Error here...
#########################################
[Tue Oct 13 16:10:37 2009] Beginning TopHat run (v1.0.10)
-----------------------------------------------
[Tue Oct 13 16:10:37 2009] Preparing output location ./tophat_out/
[Tue Oct 13 16:10:37 2009] Checking for Bowtie index files
[Tue Oct 13 16:10:37 2009] Checking for reference FASTA file
[Tue Oct 13 16:10:37 2009] Checking for Bowtie
Bowtie version: 0.10.0.2
[Tue Oct 13 16:10:37 2009] Checking reads
Error: could not open file
#############################################
jyli is offline   Reply With Quote
Old 10-13-2009, 12:59 PM   #2
Cole Trapnell
Senior Member
 
Location: Boston, MA

Join Date: Nov 2008
Posts: 212
Default

You appear to have the reads separated by spaces and commas. The left reads should be a single comma separated list, and same for the right. Like this:

tophat --solexa1.3-quals -m 1 -r 130 -G /datastore/nextgenproc/projects/exonDB/custom_exonDB.04152009.hg18.exon.gff /home/jyli/Bowtie/indexes/hg18_genome
../raw_data/MCF7_CP7_00300_003_s_4_1.txt,../raw_data/MCF7_CP7_00300_008_s_7_1.txt ../raw_data/MCF7_CP7_00300_003_s_4_2.txt,../raw_data/MCF7_CP7_00300_008_s_7_2.txt
Cole Trapnell is offline   Reply With Quote
Old 10-14-2009, 06:01 AM   #3
jyli
Member
 
Location: North Carolina

Join Date: Nov 2008
Posts: 21
Default

Thank you so much, it worked out properly!
jyli 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 06:31 PM.


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