-
Notifications
You must be signed in to change notification settings - Fork 35
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
Refactor mitochondrial analysis workflow #419
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.
Found one tiny thing, I will take another look before approving
docs/usage.md
Outdated
5. Variant calling - Structural variants (SV) (Tiddit & Manta) | ||
6. SNV annotation & ranking (rohcall, vcfanno, ensembl VEP, GENMOD) | ||
7. SV annotation & ranking (SVDB query, ensembl VEP, GENMOD) | ||
8. Mitochondrial analysis | ||
8. Mitochondrial annotation | ||
|
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.
Mitochondrial annotation is 8 here, but was 9 before (line 23) and is 9 later on (line 245), check if it's correct
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.
Massive work Ram. Some small comments and questions
.genotyped_intervals_vcf | ||
.collect{it[1]} | ||
.set { gcnvcaller_vcf } | ||
|
||
CALL_SV_MT (ch_mt_bam_bai, ch_genome_fasta) | ||
|
||
// CALL_SV_MT_SHIFT (ch_mtshift_bam_bai, ch_mtshift_fasta) |
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.
remove?
@@ -0,0 +1,111 @@ | |||
// | |||
// Merge and annotate MT |
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.
// Merge and annotate MT | |
// Merge MT variants |
ch_mt_vcf | ||
.combine(ANNOTATE_CADD.out.vcf.ifEmpty("null")) | ||
.branch { it -> | ||
merged: it[2].equals("null") | ||
return [it[0], it[1]] | ||
cadd: !(it[2].equals("null")) | ||
return [it[2], it[3]] | ||
} | ||
.set { ch_for_mix } | ||
ch_vep_in = ch_for_mix.merged.mix(ch_for_mix.cadd) |
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 there reason to use the same type of channel combinations as in #421 ?
Co-authored-by: Anders Jemt <jemten@users.noreply.github.com>
…ip ci] Co-authored-by: Anders Jemt <jemten@users.noreply.github.com>
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.
Looks good
PR checklist
This PR,
skip_mt_analysis
which was used to turn on/off the mitochondrial workflowskip_mt_annotation
which can be used to turn on/off annotation and ranking for mitochondrial SNVsWith changes in this PR, mitochondrial SNVs will no longer be called by deepvariant, and they will come from the GATKs best practices workflow for mitochondria.
nf-core lint
).nextflow run . -profile test,docker --outdir <OUTDIR>
).nextflow run . -profile test_one_sample,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).