美文网首页
juicer挂载遇到的报错

juicer挂载遇到的报错

作者: 徒唤奈何_c5f0 | 来源:发表于2024-07-02 09:48 被阅读0次

conda配好环境之后,运行结果如下

bash ../01.juicer/juicer.sh -y 460m.3d_MboI.txt -z 460m.polished2.fasta -s MboI -p chr.size -D ../01.juicer -t 16 --assembly
Using 460m.3d_MboI.txt as site file
(-: Aligning files matching /mnt/e/bio_data/cy_genome_make/01.hic.gua.2/fastq/*_R*.fastq*
 to genome 460m.polished2.fasta with no fragment delimited maps.
[E::hts_open_format] Failed to open file "/mnt/e/bio_data/cy_genome_make/01.hic.gua.2/splits/*.bam" : No such file or directory
samtools merge: fail to open "/mnt/e/bio_data/cy_genome_make/01.hic.gua.2/splits/*.bam": No such file or directory
***! Some problems occurred somewhere in creating sorted align files.

经群友指点之后发现,在fastq文件夹中的hic测序文件格式,需要遵循_R1.fastq或者_R1.fastq.gz的命名格式。而且每次需要将splits文件夹内的内容清空才能运行。

bash ../01.juicer/juicer.sh -y 460m.3d_MboI.txt -z 460m.polished2.fasta -s MboI -p chr.size -D ../01.juicer -t 16 --assembly
Using 460m.3d_MboI.txt as site file
(-: Looking for fastq files...fastq files exist
(-: Created /mnt/e/bio_data/cy_genome_make/01.hic.gua.2/splits.
(-: Aligning files matching /mnt/e/bio_data/cy_genome_make/01.hic.gua.2/fastq/*_R*.fastq*
 to genome 460m.polished2.fasta with no fragment delimited maps.

修改之后成功运行。

相关文章

网友评论

      本文标题:juicer挂载遇到的报错

      本文链接:https://www.haomeiwen.com/subject/mjevcjtx.html