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

Allow choice between 3' and 5' read ends #42

Open
2 tasks
bmmalone opened this issue Feb 2, 2017 · 0 comments
Open
2 tasks

Allow choice between 3' and 5' read ends #42

bmmalone opened this issue Feb 2, 2017 · 0 comments
Labels

Comments

@bmmalone
Copy link
Contributor

bmmalone commented Feb 2, 2017

Currently, the code is hard-coded to use the 5' ends of riboseq reads. However, for some sequencing protocols (mostly for prokaryotes), it has been shown that using the 3' end is better. The logic for end selection is entirely in the riboseq-utils package. In particular, it appears to only be used in:

  • extract_metagene_profiles.py: around line 114, the call to bam_utils.get_five_prime_ends
  • ribo_utils.get_p_sites: around line 510, around the message "Discarding 3' end of reads"

It appears the choice could be added by passing a relevant flag to these locations in the code.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

No branches or pull requests

2 participants