SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Samtools flagstat on Tophat generated .bam file rg_gis Bioinformatics 6 02-27-2015 09:35 PM
tophat problem: no accepted_hits.bam generated RNAer Bioinformatics 9 08-30-2013 09:41 AM
TOPHAT EMPTY accepted_hits.bam ISSUE waterboy Bioinformatics 1 11-16-2010 08:48 AM
Split accepted_hits.bam file after Tophat run? hong_sunwoo Bioinformatics 6 10-18-2010 12:06 AM
Tophat: Is accepted_hits.sam normalized? jling Bioinformatics 2 11-03-2009 04:22 PM

Reply
 
Thread Tools
Old 07-19-2011, 12:18 PM   #1
RNAer
Junior Member
 
Location: NY

Join Date: Mar 2011
Posts: 5
Default tophat problem: no accepted_hits.bam generated

I am using tophat to map the SOLiD single-end 50mer reads to
S.cerevisae genome. The commands I used are:

$ bowtie-build -C reads/sgd_genome.fa reads/sgd_genome
$ qsub -V -o WT.1_tophat.O -e WT.1_tophat.E -N WT.1_tophat -l
nodes=1pn1,walltime=10:00:00 <<< " tophat --color --output-dir
cufflinks/WT.1 --GTF reads/gene_features.gff --quals reads/genome
reads/WT.1.csfasta reads/WT.1.qual"

The job ran without any error but there is no accepted_hits.bam
generated and all the BED files are empty.
$ ls -l WT.1
total 16K
-rw-rw-r-- 1 zxu7 zxu7 0 07/19 15:15:24 deletions.bed
-rw-rw-r-- 1 zxu7 zxu7 0 07/19 15:15:24 insertions.bed
-rw-rw-r-- 1 zxu7 zxu7 0 07/19 15:15:24 junctions.bed
-rw-rw-r-- 1 zxu7 zxu7 70 07/19 12:30:42 left_kept_reads.info
drwxrwxr-x 2 zxu7 zxu7 8.0K 07/18 14:06:23 logs

$ cat WT.1_tophat.E
[Mon Jul 18 10:36:07 2011] Beginning TopHat run (v1.3.1)
-----------------------------------------------
[Mon Jul 18 10:36:07 2011] Preparing output location cufflinks/WT.1/
[Mon Jul 18 10:36:07 2011] Checking for Bowtie index files
[Mon Jul 18 10:36:07 2011] Checking for reference FASTA file
Warning: Could not find FASTA file cufflinks/genome.fa
[Mon Jul 18 10:36:07 2011] Reconstituting reference FASTA file from Bowtie index
Executing: /home/zxu7/Softwares/bowtie-0.12.7/bowtie-inspect
cufflinks/genome > cufflinks/WT.1/tmp/genome.fa
[Mon Jul 18 10:36:08 2011] Checking for Bowtie
Bowtie version: 0.12.7.0
[Mon Jul 18 10:36:08 2011] Checking for Samtools
Samtools Version: 0.1.14
[Mon Jul 18 10:36:08 2011] Generating SAM header for cufflinks/genome
[Mon Jul 18 10:36:11 2011] Preparing reads
format: fasta
[Mon Jul 18 10:36:11 2011] Reading known junctions from GTF file
Left reads: min. length=50, count=32186882
[Mon Jul 18 10:46:11 2011] Mapping left_kept_reads against genome with Bowtie
[Mon Jul 18 11:22:32 2011] Processing bowtie hits
[Mon Jul 18 11:56:23 2011] Mapping left_kept_reads_seg1 against genome
with Bowtie (1/2)
[Mon Jul 18 12:26:07 2011] Mapping left_kept_reads_seg2 against genome
with Bowtie (2/2)
[Mon Jul 18 13:01:47 2011] Searching for junctions via segment mapping
[Mon Jul 18 13:13:27 2011] Retrieving sequences for splices
[Mon Jul 18 13:13:28 2011] Indexing splices
[Mon Jul 18 13:13:29 2011] Mapping left_kept_reads_seg1 against
segment_juncs with Bowtie (1/2)
[Mon Jul 18 13:38:03 2011] Mapping left_kept_reads_seg2 against
segment_juncs with Bowtie (2/2)
[Mon Jul 18 13:48:21 2011] Joining segment hits
[Mon Jul 18 14:06:53 2011] Reporting output tracks
-----------------------------------------------
Run complete [03:32:48 elapsed]

By the way, the gff files is in GFF version 3 format. Does anyone know
what's wrong with my tophat command?
RNAer 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 02:02 AM.


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