SEQanswers

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 06:15 AM
Tophat2: Error running bowtie: JonB Bioinformatics 5 01-18-2014 04:09 AM
TopHat2 error running 'prep_reads' Kaskar Bioinformatics 5 07-20-2013 02:16 AM
TopHat2 Error running 'prep_reads' srmaruyama Bioinformatics 0 05-30-2013 04:34 PM
Tophat2 Error running 'long_spanning_reads': dvanic Bioinformatics 16 04-26-2013 01:56 PM

Reply
 
Thread Tools
Old 02-29-2016, 09:55 AM   #1
natalie_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 2.2.6.0) 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: 2.2.6.0
[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
[FAILED]
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
natalie_1 is offline   Reply With Quote
Old 02-29-2016, 10:02 AM   #2
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,548
Default

For reference cross-posted: https://www.biostars.org/p/179178/
GenoMax is offline   Reply With Quote
Reply

Tags
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 01:16 PM.


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