tophat2 runs with error:: segment-based junction search failed with err =-6
0
0
Entering edit mode
20 months ago
pxd • 0

The command I use is :

tophat -p 8 -r 101 -G gencode.v41.annotation.gtf -o GM12878_1 pxd ENCFF001REJ.fastq ENCFF001REK.fastq 

[2022-08-18 17:30:19] Beginning TopHat run (v2.1.1)
-----------------------------------------------
[2022-08-18 17:30:19] Checking for Bowtie
          Bowtie version:    2.2.5.0
[2022-08-18 17:30:19] Checking for Bowtie index files (genome)..
[2022-08-18 17:30:19] Checking for reference FASTA file
[2022-08-18 17:30:19] Generating SAM header for pxd
[2022-08-18 17:30:21] Reading known junctions from GTF file
[2022-08-18 17:31:20] Preparing reads
     left reads: min. length=101, max. length=101, 97505371 kept reads (42681 discarded)
    right reads: min. length=101, max. length=101, 97521988 kept reads (26064 discarded)
[2022-08-18 18:43:47] Building transcriptome data files GM12878_1/tmp/gencode.v41.annotation
[2022-08-18 18:45:18] Building Bowtie index from gencode.v41.annotation.fa
[2022-08-18 19:15:59] Mapping left_kept_reads to transcriptome gencode.v41.annotation with Bowtie2 
[2022-08-18 22:19:22] Mapping right_kept_reads to transcriptome gencode.v41.annotation with Bowtie2 
[2022-08-19 01:24:03] Resuming TopHat pipeline with unmapped reads
[2022-08-19 01:24:04] Mapping left_kept_reads.m2g_um to genome pxd with Bowtie2 
[2022-08-19 02:39:17] Mapping left_kept_reads.m2g_um_seg1 to genome pxd with Bowtie2 (1/4)
[2022-08-19 02:41:14] Mapping left_kept_reads.m2g_um_seg2 to genome pxd with Bowtie2 (2/4)
[2022-08-19 02:43:19] Mapping left_kept_reads.m2g_um_seg3 to genome pxd with Bowtie2 (3/4)
[2022-08-19 02:45:21] Mapping left_kept_reads.m2g_um_seg4 to genome pxd with Bowtie2 (4/4)
[2022-08-19 02:48:02] Mapping right_kept_reads.m2g_um to genome pxd with Bowtie2 
[2022-08-19 04:09:07] Mapping right_kept_reads.m2g_um_seg1 to genome pxd with Bowtie2 (1/4)
[2022-08-19 04:11:20] Mapping right_kept_reads.m2g_um_seg2 to genome pxd with Bowtie2 (2/4)
[2022-08-19 04:14:00] Mapping right_kept_reads.m2g_um_seg3 to genome pxd with Bowtie2 (3/4)
[2022-08-19 04:16:34] Mapping right_kept_reads.m2g_um_seg4 to genome pxd with Bowtie2 (4/4)
[2022-08-19 04:19:50] Searching for junctions via segment mapping
    [FAILED]
Error: segment-based junction search failed with err =-6
  Reason: tried: '/usr/local/opt/boost/lib/libboost_thread-mt.dylib' (no such file), '/usr/local/lib/libboost_thread-mt.dylib' (no such file), '/usr/lib/libboost_thread-mt.dylib' (no such file)

I installed the boost using conda and I can not find the libboost_thread-my.dylib in the directories, and I want to know how to solve it? Thanks in advance.

tophat2 • 400 views
ADD COMMENT

Login before adding your answer.

Traffic: 2364 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