Entering edit mode
2.3 years 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.