nf-core / nf-core/clipseq

Reads with / characters in the read names cause the pipeline to fail

Open
#57 1 comment 0 reactions 0 assignees View on GitHub

Nobody has claimed this yet.

bug
Dominant language
Nextflow
Stars
26
Forks
43
PR merge metrics
No merged PRs in 30d

Description

Check Documentation

I have checked the following places for your error:

Description of the bug

If / characters are in the read names, then the pipeline will fail at the dedup step. This is caused, however, by STAR truncating read names based on the / character. This removes the UMI information which was added earlier in the pipeline, causing the dedup step to fail.

Steps to reproduce

Steps to reproduce the behaviour:

Running the pipeline with reads with a / in the read name. Example

Expected behaviour

Either a verbose error message saying that the reads are incompatible, or automatic removal or replacement of the / characters.

Contributor guide

Open the contributing guide

First steps

  1. Read the whole issue, then the project's contributing guide.
  2. Comment on the issue to say you are picking it up — it saves two people doing the same work.
  3. Fork the repository and make your change on a branch.
  4. Open a pull request that references the issue number.

Research direction

Start by reproducing the failure with the ENA example linked in the issue, then inspect the STAR mapping step and the subsequent dedup step to confirm where read names are truncated. Done means the pipeline either handles read names containing / characters or reports a clear incompatibility error before deduplication.

Written by the indexing model from the issue text.

Assessment

Domain
bioinformatics
Issue type
Bug
Difficulty
4/5
Estimated time
3-5 days
Activity status
Stale
Clarity
Mostly clear
Newbie friendliness
35/100

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.