SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
tophat Error running running 'prep_reads' victoryhe Bioinformatics 2 10-17-2011 05:53 AM
Error when running TopHat clk RNA Sequencing 3 09-29-2011 01:10 PM
Tophat Running Error Wei-HD Bioinformatics 12 02-11-2010 03:42 PM
TopHat running error pfranchini Bioinformatics 2 08-10-2009 07:46 AM
TopHat running error ? baohua100 Bioinformatics 3 07-22-2009 04:58 PM

Reply
 
Thread Tools
Old 09-02-2009, 05:12 PM   #1
day_sleeperz
Junior Member
 
Location: Adelaide

Join Date: Aug 2009
Posts: 1
Exclamation Tophat running error 1?

Hi , I tried to run tophat to align my Solexa reads against the latest GrCh37 assembly. I encounter this error. Could someone help me? Thanks

tophat -i 69 --solexa1.3-quals --coverage-search --fill-gap --microexon-search --GFF Ensembl_h_sapiens.GFF3 h_sapiens_37_asm s_2_sequence_fastq.txt

[Wed Sep 2 16:43:33 2009] Beginning TopHat run (v1.0.11)
-----------------------------------------------
[Wed Sep 2 16:43:33 2009] Preparing output location ./tophat_out/
[Wed Sep 2 16:43:33 2009] Checking for Bowtie index files
[Wed Sep 2 16:43:33 2009] Checking for reference FASTA file
Warning: Could not find FASTA file /home/neuro/Documents/Sonny/bowtie-0.10.1/indexes/h_sapiens_37_asm.fa
[Wed Sep 2 16:43:33 2009] Reconstituting reference FASTA file from Bowtie index
[Wed Sep 2 17:01:26 2009] Checking for Bowtie
Bowtie version: 0.10.1.0
[Wed Sep 2 17:01:27 2009] Checking reads
seed length: 65bp
format: fastq
quality scale: --solexa1.3-quals
[Wed Sep 2 17:03:59 2009] Reading known junctions from GFF file
[Wed Sep 2 17:05:23 2009] Mapping reads against h_sapiens_37_asm with Bowtie
[Wed Sep 2 18:45:28 2009] Joining segment hits
Splitting reads into 2 segments
[Wed Sep 2 18:46:04 2009] Mapping reads against h_sapiens_37_asm with Bowtie
[Wed Sep 2 19:07:29 2009] Mapping reads against h_sapiens_37_asm with Bowtie
[Wed Sep 2 19:26:40 2009] Searching for junctions via segment mapping
[Wed Sep 2 19:28:49 2009] Retrieving sequences for splices
[Wed Sep 2 19:30:36 2009] Indexing splices
Warning: Empty input file
Reference file does not seem to be a FASTA file
[FAILED]
Error: Splice sequence indexing failed with err = 1
day_sleeperz is offline   Reply With Quote
Old 10-03-2011, 07:10 AM   #2
hrajasim
Member
 
Location: DC Metro

Join Date: Aug 2009
Posts: 27
Default Figured out yet?

I am running the latest version of tophat 1.3.2 but hit this same error=1 for the first time. I had run tophat 1.2 successfully on 50+ samples. ANy resolution to this issue yet?

I am getting this error only for the second iteration of tophat. I ran first iteration of tophat on each individual sample without the -j option. Then, merged all .juncs.bed files into a master juncs.bed file. Now, when I try tophat with the -j option on each samples (2nd iteration), I get this error.

Thanks for any pointers.
hrajasim is offline   Reply With Quote
Old 10-04-2011, 08:19 AM   #3
hrajasim
Member
 
Location: DC Metro

Join Date: Aug 2009
Posts: 27
Default

Okay this issue was resolved by ensuring a single junctions header at the top and all others removed. Bed needs to be converted into .juncs file using the script bed_to_juncs before feeding into the tophat -j option.
hrajasim is offline   Reply With Quote
Reply

Tags
error, 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 12:21 AM.


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