Error running bowtie: .. (When I use tophat2)
0
0
Entering edit mode
9.0 years ago
leenaehyeon ▴ 10

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 • 3.3k views
ADD COMMENT

Login before adding your answer.

Traffic: 1150 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6