xiaochuanle / NECAT

Nanopore data assembler
135 stars 24 forks source link

运行 necat.pl correct Necat.txt 报错 #8

Open baizixiansheng opened 5 years ago

baizixiansheng commented 5 years ago

我的配置文件是: (base) lge@dell:~/Chen_/bio-soft$ cat Necat.txt PROJECT=S-necat ONT_READ_LIST=read_list.txt GENOME_SIZE=1200000000 THREADS=8 MIN_READ_LENGTH=5000 OVLP_FAST_OPTIONS="-n 500 -z 20 -b 2000 -e 0.5 -j 0 -u 1 -a 1000" OVLP_SENSITIVE_OPTIONS="-n 500 -z 10 -e 0.5 -j 0 -u 1 -a 1000" CNS_FAST_OPTIONS="-a 2000 -x 4 -y 12 -l 1000 -e 0.5 -p 0.8 -u 0" CNS_SENSITIVE_OPTIONS="-a 2000 -x 4 -y 12 -l 1000 -e 0.5 -p 0.8 -u 0" TRIM_OVLP_OPTIONS="-n 100 -z 10 -b 2000 -e 0.5 -j 1 -u 1 -a 400" ASM_OVLP_OPTIONS="-n 100 -z 10 -b 2000 -e 0.5 -j 1 -u 0 -a 400" NUM_ITER=2 CNS_OUTPUT_COVERAGE=100 CLEANUP=1 USE_GRID=false GRID_NODE=0 FSA_OL_FILTER_OPTIONS="--max_overhang=-1 --min_identity=-1 --coverage=40" FSA_ASSEMBLE_OPTIONS="" FSA_CTG_BRIDGE_OPTIONS="--dump --read2ctg_min_identity=80 --read2ctg_min_coverage=4 --read2ctg_m0 --ctg_min_length=1000 --read2ctg_min_aligned_length=5000 --select_branch=best" (base) lge@dell:~/Chen_/bio-soft$ cat read_list.txt /home/lge/Chen_/Data/N01.filt.fq 具体如下: (base) lge@dell:~/Chen_/bio-soft$ ./NECAT/Linux-amd64/bin/necat.pl correct Necat.txt Smartmatch is experimental at /home/lge/Chen_/bio-soft/NECAT/Linux-amd64/bin/Plgd/Project.pm lin 2019-07-30 15:14:46 [Info] Start correcting rawreads. 2019-07-30 15:14:46 [Info] Start filtering rawreads for consensus. 2019-07-30 15:14:46 [Info] Run script: /home/lge/Chen_/bio-soft/Sty-necat/scripts/cns_pprr.sh 2>necat/scripts/cns_pprr.sh.log [Tue Jul 30 15:14:46 2019] INFO: process fastq file /home/lge/Chen_/Data/N01.filt.fq (preprocess 2019-07-30 17:24:26 [Info] End filtering rawreads for consensus. 2019-07-30 17:24:26 [Info] Start correcting iteratively for cns 1.. 2019-07-30 17:24:26 [Info] Start aliging for cns 1. 2019-07-30 17:24:26 [Info] Start making vol for cns 1. 2019-07-30 17:24:26 [Info] Run script: /home/lge/Chen_/bio-soft/Sty-necat/scripts/cns_1_mk_vol.sSty-necat/scripts/cns_1_mk_vol.sh.log file_list: /home/lge/Chen_/bio-soft/Sty-necat/1-consensus/raw_reads/raw_read_list.txt [Tue Jul 30 17:39:50 2019] INFO: pack /home/lge/Chen_/bio-soft/Sty-necat/1-consensus/raw_reads/r: 2505406 reads, 69390818071 bps (makedb/main.c, 44) 2019-07-30 17:39:50 [Info] End making vol for cns 1. 2019-07-30 17:39:50 [Info] Start aligning vol for cns 1. 2019-07-30 17:39:50 [Info] Parallelly start aligning vol 0 for cns 1. 2019-07-30 17:39:50 [Info] Parallelly start aligning vol 1 for cns 1. 2019-07-30 17:39:50 [Info] Parallelly start aligning vol 2 for cns 1. 2019-07-30 17:39:50 [Info] Parallelly start aligning vol 3 for cns 1. 2019-07-30 17:39:50 [Info] Parallelly start aligning vol 4 for cns 1. 2019-07-30 17:39:50 [Info] Parallelly start aligning vol 5 for cns 1. 2019-07-30 17:39:50 [Info] Parallelly start aligning vol 6 for cns 1. 2019-07-30 17:39:50 [Info] Parallelly start aligning vol 7 for cns 1. 2019-07-30 17:39:50 [Info] Parallelly start aligning vol 8 for cns 1. 2019-07-30 17:39:50 [Info] Parallelly start aligning vol 9 for cns 1. 2019-07-30 17:39:50 [Info] Parallelly start aligning vol 10 for cns 1. 2019-07-30 17:39:50 [Info] Parallelly start aligning vol 11 for cns 1. 2019-07-30 17:39:50 [Info] Parallelly start aligning vol 12 for cns 1. 2019-07-30 17:39:50 [Info] Parallelly start aligning vol 13 for cns 1. 2019-07-30 17:39:50 [Info] Parallelly start aligning vol 14 for cns 1. 2019-07-30 17:39:50 [Info] Parallelly start aligning vol 15 for cns 1. 2019-07-30 17:39:50 [Info] Parallelly start aligning vol 16 for cns 1. 2019-07-30 17:39:50 [Info] Parallelly start aligning vol 17 for cns 1. 2019-07-30 17:39:50 [Info] Parallelly start aligning vol 18 for cns 1. 2019-07-30 17:39:50 [Info] Parallelly start aligning vol 19 for cns 1. 2019-07-30 17:39:50 [Info] Parallelly start aligning vol 20 for cns 1. 2019-07-30 17:39:50 [Info] Parallelly start aligning vol 21 for cns 1. 2019-07-30 17:39:50 [Info] Parallelly start aligning vol 22 for cns 1. 2019-07-30 17:39:50 [Info] Parallelly start aligning vol 23 for cns 1. 2019-07-30 17:39:50 [Info] Parallelly start aligning vol 24 for cns 1. 2019-07-30 17:39:50 [Info] Parallelly start aligning vol 25 for cns 1. 2019-07-30 17:39:50 [Info] Parallelly start aligning vol 26 for cns 1. 2019-07-30 17:39:50 [Info] Parallelly start aligning vol 27 for cns 1. 2019-07-30 17:39:50 [Info] Parallelly start aligning vol 28 for cns 1. 2019-07-30 17:39:50 [Info] Parallelly start aligning vol 29 for cns 1. 2019-07-30 17:39:50 [Info] Parallelly start aligning vol 30 for cns 1. 2019-07-30 17:39:50 [Info] Parallelly start aligning vol 31 for cns 1. 2019-07-30 17:39:50 [Info] Parallelly start aligning vol 32 for cns 1. 2019-07-30 17:39:50 [Info] Parallelly start aligning vol 33 for cns 1. 2019-07-30 17:39:50 [Info] Parallelly start aligning vol 34 for cns 1. 2019-07-30 17:39:50 [Info] Run script: /home/lge/Chen_/bio-soft/Sty-necat/scripts/cns_1_al_vol_0ome/lge/Chen_/bio-soft/Sty-necat/scripts/cns_1_al_vol_0.sh.log /home/lge/Chen_/bio-soft/Sty-necat/scripts/cns_1_al_vol_0.sh: line 11: 8142 Segmentation fault d) /home/lge/Chen_/bio-soft/NECAT/Linux-amd64/bin/oc2pmov -j 0 -e 0.5 -m 500 -d 0.25 -q 500 -a 1 -k 15 -z 10 -u 1 -t 1 -n 500 -i 1 -t 8 /home/lge/Chen_/bio-soft/Sty-necat/1-consensus/cns_iter1ome/lge/Chen_/bio-soft/Sty-necat/1-consensus/cns_iter1/PackedData/pm_result_0 2019-07-30 17:42:41 [Error] Failed to run script, 139, /home/lge/Chen_/bio-soft/Sty-necat/scripts/cns_1_al_vol_0.sh cns_1_al_vol0.sh如下: `#!/bin/bash export PATH=/home/lge/Chen/bio-soft/NECAT/Linux-amd64/bin:/home/lge/Chen/bio-soft/NECAT/Linux-amd64/bin:$PATH retVal=0 if [ $retVal -eq 0 ]; then /home/lge/Chen/bio-soft/NECAT/Linux-amd64/bin/oc2pmov -j 0 -e 0.5 -m 500 -d 0.25 -q 500 -a 1000 -s 3 -b 2000 -k 15 -z 10 -u 1 -t 1 -n 500 -i 1 -t 8 /home/lge/Chen_/bio-soft/Sty-necat/1-consensus/cnsiter1/PackedData 0 /home/lge/Chen/bio-soft/Sty-necat/1-consensus/cns_iter1/PackedData/pm_result_0 temp_result=$? if [ $retVal -eq 0 ]; then retVal=$tempresult fi fi echo $retVal > /home/lge/Chen/bio-soft/Sty-necat/scripts/cns_1_al_vol_0.sh.done` 请问肖老师,这个该怎么解决?

baizixiansheng commented 5 years ago

我试了一下添加到环境变量,然后继续跑了

baizixiansheng commented 5 years ago

但是还是不行,请问老师,这个可以怎么解决?