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

feat: add optional column for adapters #419

Merged
merged 16 commits into from
Jan 3, 2022
Merged

feat: add optional column for adapters #419

merged 16 commits into from
Jan 3, 2022

Conversation

thomasbtf
Copy link
Contributor

@thomasbtf thomasbtf commented Dec 19, 2021

Description

This PR adds the option to specify adapters sequence for trimming via the sample sheet. However, this is just an option. Per default the adapates specified in the config are used. Only for samples where the column adapters in the sample sheet these "custom adapters" are used.

It also provides some "predefined" updates sequences, which can be used by speceifing revelo-rna-seq or nimagen-easy-seq in the sample sheet.

Lastly, the docs are updated.

Related Issue

closes #356

Checklist

docs/user-guide/configuration.md Outdated Show resolved Hide resolved
docs/user-guide/configuration.md Outdated Show resolved Hide resolved
docs/user-guide/configuration.md Outdated Show resolved Hide resolved
docs/user-guide/configuration.md Outdated Show resolved Hide resolved
docs/user-guide/configuration.md Outdated Show resolved Hide resolved
docs/user-guide/configuration.md Outdated Show resolved Hide resolved
docs/user-guide/configuration.md Outdated Show resolved Hide resolved
docs/user-guide/configuration.md Outdated Show resolved Hide resolved
Comment on lines 848 to 852
if adapters == "revelo-rna-seq":
return "--adapter_sequence AGATCGGAAGAGCACACGTCTGAACTCCAGTCA --adapter_sequence_r2 AGATCGGAAGAGCGTCGTGTAGGGAAAGAGTGT"
# https://www.nimagen.com/shop/products/rc-cov096/easyseq-sars-cov-2-novel-coronavirus-whole-genome-sequencing-kit
elif adapters == "nimagen-easy-seq":
return "--adapter_sequence GCGAATTTCGACGATCGTTGCATTAACTCGCGAA --adapter_sequence_r2 AGATCGGAAGAGCGTCGTGTAGGGAAAGAGTGT"
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I'd still prefer to keep these in the config and not hard coded in the common.smk. Or am I missing something here?

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

These hardcoded sequences are just used if you enter "revelo-rna-seq" or "nimagen-easy-seq" in the new adapter column. Else the sequences in the config file are used.

However, for generating the ref data set report, we will have to use this new column.

@thomasbtf thomasbtf requested a review from alethomas December 23, 2021 18:12
@thomasbtf thomasbtf merged commit f883ce2 into master Jan 3, 2022
@thomasbtf thomasbtf deleted the feat-adapters branch January 3, 2022 13:17
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

Parallel Analysis of same Technology
2 participants