nf-core/hlatyping

'remap_to_hla' index declaration

kkyung412 opened this issue · 2 comments

Hi, I have a problem.
How can I solve the problem?

[0;35m[nf-core/hlatyping] Pipeline completed with errors
Error executing process > 'remap_to_hla (1)'

Caused by:
  Process `remap_to_hla (1)` terminated with an error exit status (1)

Command executed:

  samtools view -@ 1 -h -f 0x40 example_pe.bam > output_1.bam
  samtools view -@ 1 -h -f 0x80 example_pe.bam > output_2.bam
  samtools bam2fq output_1.bam > output_1.fastq
  samtools bam2fq output_2.bam > output_2.fastq
  yara_mapper -e 3 -t 1 -f bam /root/.nextflow/assets/nf-core/hlatyping/data/indices/yara/hla_reference_dna output_1.fastq output_2.fastq > output.bam
  samtools view -@ 1 -h -F 4 -f 0x40 -b1 output.bam > mapped_1.bam
  samtools view -@ 1 -h -F 4 -f 0x80 -b1 output.bam > mapped_2.bam

Command exit status:
  1

Command output:
  (empty)

Command error:
  [M::bam2fq_mainloop] discarded 0 singletons
  [M::bam2fq_mainloop] processed 44315 reads
  [M::bam2fq_mainloop] discarded 0 singletons
  [M::bam2fq_mainloop] processed 44315 reads
  Open failed on file /root/.nextflow/assets/nf-core/hlatyping/data/indices/yara/hla_reference_dna.txt.size: "No such file or directory"
  yara_mapper: Error while opening reference file.

Work dir:
  /data1/k2/ramu_optitype_out/results/work/99/143a4235e86bfbb45a9027550bc866

Tip: you can replicate the issue by changing to the process work dir and entering the command bash .command.run

Can you post your command line used and the .nextflow.log file?

We had a revisit to this issue in slack, and @pditommaso pointed out the issue is that the full index should be declared in the input of this process properly:

https://github.com/nf-core/hlatyping/blob/master/main.nf#L279