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.
修改之后成功运行。
网友评论