Entering edit mode
5.6 years ago
Sun Jiji
•
0
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
illumina iGenomes, UCSC hg38
illumina iGenomes, NCBI GRCh38
illumina iGenomes, Ensembl GRCh37
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: 1.2.2.0
[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
[FAILED]
Error: Couldn't build bowtie index with err = -11
TopHat has been retired, see note at TopHat site:
Or this tweet from one of the original TopHat authors:
I use tophat2 and bowtie2, and use the illumina NCBI GRCh38 index, it really works.Thank you.