-
Notifications
You must be signed in to change notification settings - Fork 419
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
First version of fermi2 / fermikit variant caller #1093
Conversation
|
<tool id="fermi2" name="fermi2" version="0.14.dev1"> | ||
<description>assembles Illumina reads into unitigs</description> | ||
<requirements> | ||
<requirement type="package" version="r193">fermi2</requirement> |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
indentaion
@@ -0,0 +1,65 @@ | |||
<tool id="fermi2" name="fermi2" version="0.14.dev1"> |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
is the version number correct? dev
seems wired to me.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Oh I see the versions are different between the tools:
+ <requirement type="package" version="0.14.dev1">fermikit</requirement>
So fermikit is a bunch of tools and fermi2 is a workflow/makefile thingy?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
There is a fermi2 binary, that focuses on the exploration of FMD-index as a graph. It relies on fermi or ropebwt2 for index construction.
fermi2.pl is a wrapper with which you can then do an assembly, and fermikit includes fermi2, bwa, trimadap etc so that you can go from raw reads to variant calls based on the assembly done by fermi2. I split the fermikit functionality into the assembly (using the fermi2
wrapper), mapping (using the standard bwa-mem wrapper with intractg
option) and variant calling (with the fermikit_variants
wrapper). I suppose we could rely on multiple exact conda packages instead of the fermikit package, if you prefer that.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
and I chose version 0.14dev1 because the dependencies in the fermikit conda package match the current master (see https://github.com/lh3/fermikit/blob/master/NEWS.md)
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Ok, makes somehow sense :)
Thanks! |
Thanks Bjoern! |
Fermikit can be used to call variants from high coverage short read sequencing
data.