Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

TRIM_READS error #210

Open
aperreault opened this issue Nov 21, 2024 · 0 comments
Open

TRIM_READS error #210

aperreault opened this issue Nov 21, 2024 · 0 comments
Labels
bug Something isn't working

Comments

@aperreault
Copy link

Description of the bug

I get the following error when I run this command: NXF_VER=23.10.1 nextflow run nf-core/hic -profile docker --input /media/aunderwood6/BCSP-extra/MCF7/HIC_MCF7_708.csv --genome GRCh38 --digestion mboi --outdir /media/aunderwood6/BCSP-extra/MCF7/HiC-analysis/HIC_MCF7_708

[nf-core/hic] Pipeline completed with errors- 
ERROR ~ Error executing process > 'NFCORE_HIC:HIC:HICPRO:HICPRO_MAPPING:TRIM_READS (HIC_MCF7_708)' 

Caused by:
  Process `NFCORE_HIC:HIC:HICPRO:HICPRO_MAPPING:TRIM_READS (HIC_MCF7_708)` terminated with an error exit status (1) 

Command executed: 
  zcat HIC_MCF7_708_0_R2.unmapped.fastq.gz > tmp.fastq 
  cutsite_trimming --fastq tmp.fastq \ 
      --cutsite AAGCTAGCTT \ 
      --out HIC_MCF7_708_0_R2_trimmed.fastq 
  gzip HIC_MCF7_708_0_R2_trimmed.fastq 
  /bin/rm -f tmp.fastq 

  cat <<-END_VERSIONS > versions.yml 
  "NFCORE_HIC:HIC:HICPRO:HICPRO_MAPPING:TRIM_READS": 
      gzip: $(echo $(gzip --version 2>&1) | head -1 | cut -d" " -f2) 
  END_VERSIONS 

Command exit status: 1 

Command output: 
  ##Fastq file: tmp.fastq 
  ##Restriction sites (n=1): 
  AAGCTAGCTT 
  ##Output File: HIC_MCF7_708_0_R2_trimmed.fastq 
  ##Trimmed reads: 1905 

Command error: 
  ##Fastq file: tmp.fastq 
  ##Restriction sites (n=1): 
  AAGCTAGCTT 
  ##Output File: HIC_MCF7_708_0_R2_trimmed.fastq 

  ##Trimmed reads: 1905 
  gzip: HIC_MCF7_708_0_R2_trimmed.fastq: Input/output error 

The restriction site in the error is correct, but there is no ^ to designate cut site. Is that the issue?

Thanks!

Command used and terminal output

Relevant files

No response

System information

No response

@aperreault aperreault added the bug Something isn't working label Nov 21, 2024
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working
Projects
None yet
Development

No branches or pull requests

1 participant