美文网首页
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