SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Finding splice sites: TopHat and SOAPsplice export empty junction files jgw RNA Sequencing 0 01-31-2012 12:04 PM
Empty junction.bed TopHat vinay052003 Bioinformatics 0 11-24-2011 05:53 PM
GATK IndelRealigner returns empty output file bjfan Bioinformatics 3 06-08-2011 11:43 AM
TopHat - empty junction database altodor RNA Sequencing 3 02-21-2011 05:22 PM
Tophat error -Junction database is empty repinementer Bioinformatics 1 08-02-2010 06:47 AM

Reply
 
Thread Tools
Old 02-14-2013, 04:43 AM   #1
sachitad
Junior Member
 
Location: Kathmandu, Nepal

Join Date: Oct 2012
Posts: 4
Default Tophat2 returns empty junction file: Warning: junction database is empty!

I am using Tophat2 V2.0.7 with this command:
Code:
tophat2 -p 1 --library-type fr-unstranded --segment-length 12 /data/rathi/Homo_sapiens/UCSC/hg19/Sequence/Bowtie2Index/genome Galaxy164-\[FASTQ_Groomer_on_data_152\].fastqsanger
#Running:
Code:
[2013-02-14 03:31:11] Beginning TopHat run (v2.0.7)
-----------------------------------------------
[2013-02-14 03:31:11] Checking for Bowtie
		  Bowtie version:	 2.0.6.0
[2013-02-14 03:31:11] Checking for Samtools
		Samtools version:	 0.1.18.0
[2013-02-14 03:31:11] Checking for Bowtie index files
[2013-02-14 03:31:11] Checking for reference FASTA file
	Warning: Could not find FASTA file /data/rathi/Homo_sapiens/UCSC/hg19/Sequence/Bowtie2Index/genome.fa
[2013-02-14 03:31:11] Reconstituting reference FASTA file from Bowtie index
  Executing: /usr/bin/bowtie2-inspect /data/rathi/Homo_sapiens/UCSC/hg19/Sequence/Bowtie2Index/genome > ./tophat_out/tmp/genome.fa
[2013-02-14 03:33:21] Generating SAM header for /data/rathi/Homo_sapiens/UCSC/hg19/Sequence/Bowtie2Index/genome
	format:		 fastq
	quality scale:	 phred33 (default)
[2013-02-14 03:33:54] Preparing reads
	 left reads: min. length=36, max. length=36, 2 kept reads (0 discarded)
[2013-02-14 03:33:54] Mapping left_kept_reads to genome genome with Bowtie2 
[2013-02-14 03:34:26] Mapping left_kept_reads_seg1 to genome genome with Bowtie2 (1/3)
[2013-02-14 03:34:59] Mapping left_kept_reads_seg2 to genome genome with Bowtie2 (2/3)
[2013-02-14 03:35:32] Mapping left_kept_reads_seg3 to genome genome with Bowtie2 (3/3)
[2013-02-14 03:36:05] Searching for junctions via segment mapping
Warning: junction database is empty!
[2013-02-14 03:37:27] Reporting output tracks
-----------------------------------------------
[2013-02-14 03:38:48] Run complete: 00:07:37 elapsed
The only thing I get ouput in junctions.bed is:
Code:
track name=junctions description="TopHat junctions"
What am I doing wrong?
sachitad is offline   Reply With Quote
Old 02-14-2013, 05:15 AM   #2
iris_aurelia
Member
 
Location: Netherlands

Join Date: Jul 2012
Posts: 22
Default

It seems that you only have 2 reads in your dataset (2 kept reads). I would try a bigger set of reads to start with.
iris_aurelia is offline   Reply With Quote
Old 02-14-2013, 09:21 PM   #3
sachitad
Junior Member
 
Location: Kathmandu, Nepal

Join Date: Oct 2012
Posts: 4
Default

I tried with all the possible data but am getting the same error about junction. Is it because of the reference genome?
sachitad is offline   Reply With Quote
Old 02-15-2013, 12:35 AM   #4
iris_aurelia
Member
 
Location: Netherlands

Join Date: Jul 2012
Posts: 22
Default

Have you checked if the created reference fasta file wasn't empty? Mayby you could try to recreate the reference file upon analysis using bowtie2-inspect?
iris_aurelia is offline   Reply With Quote
Old 02-15-2013, 01:40 AM   #5
sachitad
Junior Member
 
Location: Kathmandu, Nepal

Join Date: Oct 2012
Posts: 4
Default

I already have .fa file in Bowtie2Index directory and it's not empty:

These are files with reference gnome:

genome.1.bt2 genome.2.bt2 genome.4.bt2 genome.fa.fai genome.rev.2.bt2 reads_2.fq
genome.3.bt2 genome.fa genome.rev.1.bt2 reads_1.fq

Last edited by sachitad; 02-15-2013 at 01:45 AM.
sachitad is offline   Reply With Quote
Reply

Tags
tophat, 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 03:52 PM.


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