Question: Getting: Error reading _plen[] array: 1776466139, 4294967292 with TopHat2/bowtie2
gravatar for rubic
5.3 years ago by
United States
rubic190 wrote:


I think seem to be doing something very basically wrong in running TopHat2 so it's probably an easy answer.


Here's my command:

tophat --read-mismatches 1 --output-dir /n/user/rubi/data/tophat.alns/s1/s1 --num-threads 8 --transcriptome-index /n/user/rubi/refdata/tophat.2.gencode.vM2.annotation.2wayconspseudos.mm10/mm10 /n/user/rubi/refdata/bowtie.2.gencode.vM2.annotation.2wayconspseudos.mm10/mm10 /n/user/rubi/data/fq/s1.fq



Here's the log file:


[2015-03-20 00:07:34] Beginning TopHat run (v2.0.4)
[2015-03-20 00:07:34] Checking for Bowtie
Bowtie version:
[2015-03-20 00:07:34] Checking for Samtools
Samtools version:  0.1.12.a
[2015-03-20 00:07:34] Checking for Bowtie index files
[2015-03-20 00:07:34] Checking for Bowtie index files
[2015-03-20 00:07:34] Checking for reference FASTA file
[2015-03-20 00:07:34] Generating SAM header for /n/user/rubi/refdata/tophat.2.gencode.vM2.annotation.2wayconspseudos.mm10/mm10
format:  fastq
quality scale:  phred33 (default)
[2015-03-20 00:07:35] Reading known junctions from GTF file
[2015-03-20 00:08:03] Preparing reads
 left reads: min. length=20, max. length=59, 128844351 kept reads (192369 discarded)
[2015-03-20 00:49:52] Using pre-built transcriptome index..
[2015-03-20 00:50:20] Mapping left_kept_reads to transcriptome mm10 with Bowtie2 
[2015-03-20 03:33:13] Resuming TopHat pipeline with unmapped reads
[2015-03-20 03:33:15] Mapping left_kept_reads.m2g_um to genome mm10 with Bowtie2 
Error running bowtie:
Error reading _plen[] array: 1776466139, 4294967292
Error: Encountered internal Bowtie 2 exception (#1)
Command: /n/sw/bowtie2-2.0.0-beta6/bowtie2-align -q -k 20 -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 L,-0.6,-0.6 -p 8 --sam-no-hd -x /n/user/rubi/refdata/bowtie.2.gencode.vM2.annotation.2wayconspseudos.mm10/mm10 - 

As far as I can tell BOWTIE_INDEXES is defined as an environment variable:

rubi@rclogin09 10:25 ~ >echo $BOWTIE_INDEXES



And the transcriptome and genome indix folders seem to be found - according to the log. So what's wrong?


bowtie rna-seq tophat • 3.2k views
ADD COMMENTlink modified 5.1 years ago by Biostar ♦♦ 20 • written 5.3 years ago by rubic190

This error has something to do with reading the index into memory. My guess would be that either you're running out of memory or that there's something wrong with one of the indices. In either case, try to delete and remake then genome index. BTW, you should really upgrade your versions of tophat2 and bowtie2.

ADD REPLYlink written 5.3 years ago by Devon Ryan95k
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: 1492 users visited in the last hour