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

Add 'existing paired alignment' to complex pipeline #247

Open
njrollins opened this issue Aug 7, 2020 · 2 comments
Open

Add 'existing paired alignment' to complex pipeline #247

njrollins opened this issue Aug 7, 2020 · 2 comments

Comments

@njrollins
Copy link

For some applications, such as protein-RNA complexes, I'd like to pair sequences by a custom method- and input already paired seqs (2 alignments where sequence X in ali 1 is paired with sequence X in ali 2) into pipeline.

This means skipping the 'genome distance' or 'best hit' step in the concatenation pipeline

@thomashopf
Copy link
Contributor

Yep definitely a good idea - pinging @aggreen re protein-protein pipeline

For protein-RNA, besides the obvious hack of putting the run through the protein complex pipeline, we will need to think this through in more detail, and maybe finally create the appropriate protocols and pipelines...

@aggreen
Copy link
Contributor

aggreen commented Aug 8, 2020

Yeah, the hack version of doing this now is if you use the best hit pipeline, use 'input existing alignment' for both of the monomer stages, and provide an input annotation file that has a unique identifier for each pair of sequences that you want to pair. Then best hit will just pair up all the sequences with the same annotation. But I agree it deserves a formal solution

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

No branches or pull requests

3 participants