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

seqkit subseq ­- accept reverse coordinates in gtf #392

Closed
stas-malavin opened this issue Jun 9, 2023 · 3 comments
Closed

seqkit subseq ­- accept reverse coordinates in gtf #392

stas-malavin opened this issue Jun 9, 2023 · 3 comments

Comments

@stas-malavin
Copy link

stas-malavin commented Jun 9, 2023

Hello,
It makes sense to me if seqkit subseq would accept reverse coordinates (530 2) in gtf (and also probably bed, though I haven't tried that option).
In my scenario, I build gtf from hmmscan results, and some seqs in the initial DB are on the strand 2, hence reversed coordinates. I'm replacing them, but anyway, why not accepting as it is?
Thank you

@shenwei356
Copy link
Owner

I'll implement this later (>1-2 week(s)).

@shenwei356
Copy link
Owner

shenwei356 commented Jul 3, 2023

Oh, I forgot to tell you that it was implemented last week. Please use the binaries here: #393 (comment)

@stas-malavin
Copy link
Author

Perfect, thanks so much! :)

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

2 participants