执行下列命令行:
k=5 #设置染色体数目,染色体基数为5
ALLHiC_partition -b sample.clean.bam -r draft.asm.fasta -e DPNII -k $k
## Extract CLM file and counts of restriction sites #(HindIII: AAGCTT; MboI: GATC)
allhic extract sample.clean.bam draft.asm.fasta --RE GATC
## 排序和方向优化
for K in `seq 1 $k`;
do
allhic optimize sample.clean.counts_GATC.${k}g${K}.txt sample.clean.clm
done
报错## 排序和方向优化
[root@90227036c754 20:53:04 /work/my_LQ28473/30.ALLHiC1]# for K in `seq 1 $k`;
> do
> allhic optimize sample.clean.counts_GATC.${k}g${K}.txt sample.clean.clm
> done
20:53:16 mustOpen | CRITIC open sample.clean.counts_GATC.5g1.txt: no such file or directory
20:53:16 mustOpen | CRITIC open sample.clean.counts_GATC.5g2.txt: no such file or directory
20:53:16 mustOpen | CRITIC open sample.clean.counts_GATC.5g3.txt: no such file or directory
20:53:16 mustOpen | CRITIC open sample.clean.counts_GATC.5g4.txt: no such file or directory
20:53:16 mustOpen | CRITIC open sample.clean.counts_GATC.5g5.txt: no such file or directory,这是什么原因呢