-
Notifications
You must be signed in to change notification settings - Fork 186
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 wrapper for rbt vcf_fix_iupac_alleles #1209
base: master
Are you sure you want to change the base?
Conversation
|
||
with tempfile.TemporaryDirectory() as tmpdir: | ||
shell( | ||
"(rbt vcf-fix-iupac-alleles {extra} < {snakemake.input[0]} | bcftools sort --temp-dir {tmpdir} {bcftools_opts}) {log}" |
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.
Why the sort? rbt vcf-fix-iupac-alleles does not change the order of the records.
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.
The main idea was to allow for automatic output on any format (bcf
, vcf
, vcf.gz
).
Since I was using it to standardize resources (known variation from ensembl), I thought that sort
would be a nice-to-have (just to be sure the output was sorted) and the extra runtime would be ok since it would only be run once.
But I can also remove it (or switch to view
), if you prefer to keep the wrapper more general.
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.
Ping @johanneskoester
This PR was marked as stale because it has been open for 6 months with no activity. |
This PR was marked as stale because it has been open for 6 months with no activity. |
Description
QC
For all wrappers added by this PR,
input:
andoutput:
file paths in the resulting rule can be changed arbitrarily,threads: x
statement withx
being a reasonable default,map_reads
for a step that maps reads),environment.yaml
specifications follow the respective best practices,input:
oroutput:
),Snakefile
s and their entries are explained via comments (input:
/output:
/params:
etc.),stderr
and/orstdout
are logged correctly (log:
), depending on the wrapped tool,tempfile.gettempdir()
points to (see here; this also means that using any Pythontempfile
default behavior works),meta.yaml
contains a link to the documentation of the respective tool or command,Snakefile
s pass the linting (snakemake --lint
),Snakefile
s are formatted with snakefmt,