记录一个BioNano Solve软件运行时的一个报错

2019-02-24  本文已影响34人  xuzhougeng

重点:这个报错主要和perl的环境变量 PERL5LIB有关,不同版本的模块管理让人头疼。

在运行hybridScaffold.pl时出现下面的报错

运行代码为

perl /opt/biosoft/Solve3.3_10252018/HybridScaffold/10252018/hybridScaffold.pl\
    -n ref.fa \
    -b EXP_REFINEFINAL1.cmap \
    -c hybridScaffold_config.xml \
    -r /opt/biosoft/Solve3.3_10252018/RefAligner/7915.7989rel/RefAligner \
    -o hybrid_assembly_output \
    -B 2 -N 2 \
    -f

报错信息为

Beginning MergeNGS_BN...
*Using conflict-cut BioNano and conflict-cut sequence CMAP*
Running command: /usr/bin/perl ./scripts/MergeNGS_BN.pl -outputDir /data/xzg_data/1800_assembly/assembly/bionano/hybrid_assembly_output/mergeNGS_BN -refaligner /opt/biosoft/Solve3.3_10252018/RefAligner/7915.7989rel/RefAligner -logFile mergeNGS_BN_script.log -ngs_cmap_fn /data/xzg_data/1800_assembly/assembly/bionano/hybrid_assembly_output/assignAlignType/cut_conflicts/ref_BSPQI_0kb_0labels_cut.cmap -bng_cmap_fn /data/xzg_data/1800_assembly/assembly/bionano/hybrid_assembly_output/assignAlignType/cut_conflicts/EXP_REFINEFINAL1_bppAdjust_cut.cmap -id_shift 100000 -readparameters /data/xzg_data/1800_assembly/assembly/bionano/hybrid_assembly_output/align1/align1.errbin -xmlFile /data/xzg_data/1800_assembly/assembly/bionano/hybridScaffold_config.xml
MergeNGS_BN complete in 4 seconds.

Beginning extraction of used and not used NGS contigs in hybrid scaffold pairmerge step...
Running command: /usr/bin/perl scripts/find_used_not_used_ngs.pl /data/xzg_data/1800_assembly/assembly/bionano/hybrid_assembly_output/mergeNGS_BN 100000 /data/xzg_data/1800_assembly/assembly/bionano/hybrid_assembly_output/assignAlignType/cut_conflicts/ref_BSPQI_0kb_0labels_cut.cmap /opt/biosoft/Solve3.3_10252018/RefAligner/7915.7989rel/RefAligner
ERROR: getUsedIds: cannot open /data/xzg_data/1800_assembly/assembly/bionano/hybrid_assembly_output/mergeNGS_BN/step1.merge.pairs.txt: No such file or directory

ERROR: ERROR: Extraction of NGS contigs cannot be completed.

该步骤原本可以运行,为了避免是数据集的原因,因此在另一台服务器上安装相同的软件处理相同的数据集,结果没有任何问题。

在多方排查无果之后,闹钟突然想到会不会是PERL5LIB这个环境变量的原因。于是我检查了原本服务器上该环境变量

echo $PERL5LIB
/opt/sysoft/perl-5.28.0/lib/5.28.0:

由于上述流程的perl是v5.16.3,但是加载了5.28.0的模块,所以导致运行异常。

希望能为哪些使用BioNano软件的人提供一些帮助

上一篇下一篇

猜你喜欢

热点阅读