Biostar Beta. Not for public use.
Getting: Error reading _plen[] array: 1776466139, 4294967292 with TopHat2/bowtie2
Entering edit mode
15 months ago
rubic • 180
United States


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?

Entering edit mode

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.


Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.3