This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Error running bowtie: .. (When I use tophat2)

Hi

I was install bowtie2 (apt-get install bowtie, bowtie2, samtools)

and version is

/usr/bin/bowtie2-align version 2.1.0

I ran tophat2 for RNA-seq mapping, but this error message popped up.

[2015-12-10 14:56:25] Beginning TopHat run (v2.0.9)

-----------------------------------------------
[2015-12-10 14:56:25] Checking for Bowtie
          Bowtie version:     2.1.0.0
[2015-12-10 14:56:25] Checking for Samtools
        Samtools version:     0.1.19.0
[2015-12-10 14:56:26] Checking for Bowtie index files (transcriptome)..
[2015-12-10 14:56:26] Checking for Bowtie index files (genome)..
[2015-12-10 14:56:26] Checking for reference FASTA file
[2015-12-10 14:56:26] Generating SAM header for /home/genome/Homo_sapiens/hg19
    format:         fastq
    quality scale:     phred33 (default)
[2015-12-10 14:56:36] Reading known junctions from GTF file
[2015-12-10 14:56:38] Preparing reads
     left reads: min. length=20, max. length=101, 17432242 kept reads (15275 discarded)
    right reads: min. length=20, max. length=101, 17439959 kept reads (7558 discarded)
[2015-12-10 15:07:27] Using pre-built transcriptome data..
[2015-12-10 15:07:32] Mapping left_kept_reads to transcriptome hg19_transcripts with Bowtie2
    [FAILED]
Error running bowtie:
Error reading _ebwt[] array: 50459520, 50652672
Error: Encountered internal Bowtie 2 exception (#1)
Command: /usr/bin/bowtie2-align -q -k 60 -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 C,-14,0 -p 1 --sam-no-hd -x /home/genome/Homo_sapiens/hg19_transcripts -

My command is

tophat2 -o output -G  GTFfile --transcriptome-index /home/genome/Homo_sapiens/hg19_transcripts --library-type fr-firststrand /home/genome/Homo_sapiens/hg19 sample1.fastq sample2.fastq

Before running tophat, bowtie2-build also I did.

How can I solve this Error?

rna-seq software-error alignment

0 answers

No answers yet.

Log in to answer this question.