marcelm / marcelm/cutadapt

Output cut offsets in separate file

Open
#281 7 comments 0 reactions 0 assignees View on GitHub

Nobody has claimed this yet.

feature-request
Dominant language
Python
Stars
587
Forks
144
Avg merge
2h 37m
Merged PRs (30d)
1

Description

Hi, I would like to use cutadapt in our pipeline, which uses a database to store data and requires traceability/reconstruction of the various steps, storing cutoff points for the adapter clipping along with the original sequences. To do that it would be practical to have an option to also produce a tab separated file holding the read identifier and the four cutoff offsets (two for the forward and two for the reverse reads in pairs where the adapter is found). Offsets for the right end should be given from the end of the sequence (i.e. number of bases cutted from the right end as opposed to the offset from the beginning of the read).
Would that be feasible? Where should I look in the source tree to implement such a feature?

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

No source files, tests, or entry points are named in the issue. Start by tracing paired-read adapter trimming and output generation, then define the tab-separated record format and verify that forward and reverse cut offsets can be reconstructed for reads where adapters are found.

Written by the indexing model from the issue text.

Assessment

Tech stack
python
Domain
bioinformatics
Issue type
Feature
Difficulty
5/5
Estimated time
Over a week
Activity status
Stale
Clarity
Mostly clear
Newbie friendliness
30/100

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.