执行ALLHiC分组不成功

执行下列命令行:

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,这是什么原因呢


请先 登录 后评论
  • 0 关注
  • 0 收藏,88 浏览
  • 池yi鱼 提出于 2025-03-17 11:25

相似问题