-
Notifications
You must be signed in to change notification settings - Fork 3
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
Adapt vcf to scout #127
Adapt vcf to scout #127
Conversation
|
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.
Nice work lucia! It would be good to also make sure that the vcfs coming from the call_variants
sub workflow is normalised and decomposed.
ch_genome_chrsizes.flatten().map{chromsizes -> | ||
return [[id:variant_caller], chromsizes] | ||
} | ||
.set { ch_varcallerinfo } | ||
|
||
ADD_VARCALLER_TO_BED (ch_varcallerinfo).gz_tbi | ||
.map{meta,bed,tbi -> return [bed, tbi]} | ||
.set{ch_varcallerbed} | ||
|
||
REMOVE_DUPLICATES.out.vcf | ||
.join(TABIX_REMOVE_DUP.out.tbi) | ||
.combine(ch_varcallerbed) | ||
.combine(ch_foundin_header) | ||
.set { ch_annotate_in } | ||
|
||
BCFTOOLS_ANNOTATE(ch_annotate_in) | ||
|
||
TABIX_ANNOTATE(BCFTOOLS_ANNOTATE.out.vcf) |
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.
I would move this to the call_variants
subworkflow.
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.
Done!
Co-authored-by: Anders Jemt <jemten@users.noreply.github.com>
withName: '.*ALLELE_SPECIFIC_CALLING:ADD_VARCALLER_TO_BED' { | ||
ext.args2 = '-s 1 -b 2 -e 3' | ||
} | ||
|
||
withName: '.*ALLELE_SPECIFIC_CALLING:BCFTOOLS_ANNOTATE' { | ||
ext.args = "-c CHROM,FROM,TO,FOUND_IN --output-type z" | ||
ext.prefix = { "${meta.id}_split_rmdup_info" } |
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.
Let's remove this part since we now add the variant caller in the call_variants subworkflow
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.
Done!
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.
🎉
PR checklist
nf-core lint
).nextflow run . -profile test,docker --outdir <OUTDIR>
).nextflow run . -profile debug,test,docker --outdir <OUTDIR>
).docs/usage.md
is updated.docs/output.md
is updated.CHANGELOG.md
is updated.README.md
is updated (including new tool citations and authors/contributors).