重点：这个报错主要和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软件的人提供一些帮助
网友评论