Question: TopHat Error: Couldn't build bowtie index with err = -11
gravatar for Sun Jiji
6 months ago by
Sun Jiji0
Sun Jiji0 wrote:

When I was running tophat, I met this problem:

Error: Couldn't build bowtie index with err = -11

It seemed that it was because of the differences between the annotation and index.

But I tried different annotation files, all return the same error -11

Here are annotation and index I tried

  1. illumina iGenomes, UCSC hg38

  2. illumina iGenomes, NCBI GRCh38

  3. illumina iGenomes, Ensembl GRCh37

  4. Gencode Release 28 (GRCh38.p12) (annotation) and bowtie pre-built index NCBI GRCh38

Here is the call

root@DESKTOP-2SILU7S:/mnt/c/ubuntu# tophat --bowtie1 -p 8 -o ./tophat_out_Control_1 -G ./bowtie_indexes/GRCh38/Annotation/Archives/archive-2015-08-11-09-31-31/Genes.gencode/genes.gtf ./bowtie_indexes/GRCh38/Sequence/BowtieIndex/genome ./Control_1_1_filtered.fastq

Here is the error

[2018-08-16 14:41:23] Beginning TopHat run (v2.1.1)
[2018-08-16 14:41:23] Checking for Bowtie
                  Bowtie version:
[2018-08-16 14:41:23] Checking for Bowtie index files (genome)..
[2018-08-16 14:41:23] Checking for reference FASTA file
[2018-08-16 14:41:23] Generating SAM header for ./bowtie_indexes/GRCh38/Sequence/BowtieIndex/genome
[2018-08-16 14:41:40] Reading known junctions from GTF file
[2018-08-16 14:42:09] Preparing reads
         left reads: min. length=18, max. length=101, 56278536 kept reads (13726 discarded)
Warning: short reads (<20bp) will make TopHat quite slow and take large amount of memory because they are likely to be mapped in too many places
[2018-08-16 15:01:09] Building transcriptome data files ./tophat_out_Control_1/tmp/genes
[2018-08-16 15:02:25] Building Bowtie index from genes.fa
Error: Couldn't build bowtie index with err = -11
bowtie rna-seq tophat • 309 views
ADD COMMENTlink modified 6 months ago • written 6 months ago by Sun Jiji0

TopHat has been retired, see note at TopHat site:

Please note that TopHat has entered a low maintenance, low support stage as it is now largely superseded by HISAT2 which provides the same core functionality (i.e. spliced alignment of RNA-Seq reads), in a more accurate and much more efficient way.

Or this tweet from one of the original TopHat authors:

ADD REPLYlink modified 6 months ago • written 6 months ago by h.mon23k

I use tophat2 and bowtie2, and use the illumina NCBI GRCh38 index, it really works.Thank you.

ADD REPLYlink written 6 months ago by Sun Jiji0
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2342 users visited in the last hour