Tophat Usage difficulty?
2
0
Entering edit mode
2.1 years ago
jaqx008 ▴ 110

Hello All, I am having difficulty understanding and executing the tophat usage after looking up the help. its just not as simple as running other aligners like bowtie and I am needing help. below is the command I put together from what I understand. ( my reads are single end btw). can someone tell me what I am doing wrong and how the command is usually written so it execute without complain?

tophat --GTF Genome.Bw2.index.fa --library-type fr-unstranded --num-threads 24 Genome.Bw2.index Male.fastq --output-dir

Error message

[2019-10-02 10:46:19] Building transcriptome data files ./tophat_out/tmp/Genome.Bw2.index [2019-10-02 10:46:27] Building Bowtie index from Genome.Bw2.index.fa [FAILED] Error: Couldn't build bowtie index with err = 1

alignment tophat mapping bowtie2 • 626 views
ADD COMMENT
1
Entering edit mode

Consider other alternative tools such as HISAT2 instead of Tophat.

<script async="" src="https://platform.twitter.com/widgets.js" charset="utf-8"></script>

ADD REPLY
0
Entering edit mode

Thanks for your responses and suggestions. I figured the command works if the last imputes on the command include a reference genome index folled by the fastq RNAseq file. Again thank you

ADD REPLY
2
Entering edit mode
2.1 years ago
27atcggcta27 ▴ 170

here is a link to the manual for tophat TOPHAT2 MANUAL

I highly suggest you use hisat2 over tophat though. It is a much better option.

ADD COMMENT
0
Entering edit mode
2.1 years ago
Mensur Dlakic ★ 14k

See if making the output directory (and tmp directory within it) changes your luck. Do this in a directory where your .fasta file is.

mkdir -p tophat_out/tmp

ADD COMMENT

Login before adding your answer.

Traffic: 1720 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6