Go Back   SEQanswers > Applications Forums > RNA Sequencing

Similar Threads
Thread Thread Starter Forum Replies Last Post
Tophat error running from bam2fastx? Triple_W Bioinformatics 5 01-12-2015 05:15 AM
Tophat2: Error running bowtie: JonB Bioinformatics 5 01-18-2014 03:09 AM
TopHat2 error running 'prep_reads' Kaskar Bioinformatics 5 07-20-2013 01:16 AM
TopHat2 Error running 'prep_reads' srmaruyama Bioinformatics 0 05-30-2013 03:34 PM
Tophat2 Error running 'long_spanning_reads': dvanic Bioinformatics 16 04-26-2013 12:56 PM

Thread Tools
Old 02-29-2016, 08:55 AM   #1
Junior Member
Location: michigan

Join Date: Feb 2016
Posts: 1
Default tophat2/2.1.0 error running from bam2fastx

I am running tophat2.1.0 (bowtie on SE RNAseq data and encountering the following error:

[2016-02-29 09:38:34] Beginning TopHat run (v2.1.0)
[2016-02-29 09:38:34] Checking for Bowtie
Bowtie version:
[2016-02-29 09:38:34] Checking for Bowtie index files (transcriptome)..
[2016-02-29 09:38:34] Checking for Bowtie index files (genome)..
[2016-02-29 09:38:34] Checking for reference FASTA file
[2016-02-29 09:38:34] Generating SAM header for PGSC_DM_v4.03_index
[2016-02-29 09:38:40] Reading known junctions from GTF file
[2016-02-29 09:38:44] Preparing reads
left reads: min. length=30, max. length=40, 47957373 kept reads (864 discarded)
[2016-02-29 09:48:09] Using pre-built transcriptome data..
[2016-02-29 09:48:11] Mapping left_kept_reads to transcriptome known with Bowtie2
Error running:
/opt/software/TopHat2/2.1.0--GCC-4.4.5/bin/bam2fastx --all tophat_Kalkaska/tmp/left_kept_reads.bam|/opt/software/bowtie2/2.2.6--GCC-4.4.5/bin/bowtie2 -k 60 -D 15 -R 2 -N 0 -L 20 -i S,1,1.25 --gbar 4 --mp 6,2 --np 1 --rdg 5,3 --rfg 5,3 --score-min C,-62,0 -p 1 --sam-no-hd -x transcriptome_data/known -|/opt/software/TopHat2/2.1.0--GCC-4.4.5/bin/fix_map_ordering --bowtie2-min-score 55 --read-mismatches 3 --read-gap-length 10 --read-edit-dist 10 --read-realign-edit-dist 11 --sam-header tophat_Kalkaska/tmp/known.bwt.samheader.sam - - tophat_Kalkaska/tmp/left_kept_reads.m2g_um.bam | /opt/software/TopHat2/2.1.0--GCC-4.4.5/bin/map2gtf --sam-header tophat_Kalkaska/tmp/PGSC_DM_v4.03_index_genome.bwt.samheader.sam transcriptome_data/known.fa.tlst - tophat_Kalkaska/tmp/left_kept_reads.m2g.bam > tophat_Kalkaska/logs/m2g_left_kept_reads.out

When I run the error output, I get:

Error at parsing .tlst line (invalid strand):
53551 PGSC0003DMT400030180 ST4.03ch12. 56298573-56298656
(ERR): bowtie2-align died with signal 13 (PIPE)

Looking at the GTF file (PGSC_DM_V403_genes.gff from SpudDB), there are a fair number of entries where the strand is recorded as .. Online forums report the same issue but no solutions, save for deleting these entries, where are very numerous. Are there any other solutions I could try?

Thank you,
natalie_1 is offline   Reply With Quote
Old 02-29-2016, 09:02 AM   #2
Senior Member
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,795

For reference cross-posted:
GenoMax is offline   Reply With Quote

error correction, gtf, rnaseq analysis, 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 05:34 PM.

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