-
Notifications
You must be signed in to change notification settings - Fork 51
Subsitute samtools/bam2fq with samtools/fastq #274
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
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.
Great, thanks for making this change 🙂 Two small questions but they're not blockers.
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.
One last suggestion of a change to the emitted files but otherwise LGTM
@@ -54,7 +54,7 @@ workflow LONGREAD_HOSTREMOVAL { | |||
|
|||
emit: | |||
stats = SAMTOOLS_STATS.out.stats //channel: [val(meta), [reads ] ] | |||
reads = SAMTOOLS_BAM2FQ.out.reads // channel: [ val(meta), [ reads ] ] | |||
reads = SAMTOOLS_FASTQ.out.other // channel: [ val(meta), [ reads ] ] |
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.
Not that we use these files anywhere else, I think it would be safer to mix
all of teh outputs from FASTQ
, not just other? This would have the same behavior as the old module - otherwise here you might lose of the other reads I think (depending on how samtools fastq
interprets the input, no?)
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.
Do you suggest that we have
reads = SAMTOOLS_FASTQ.out.fastq?
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 think actually .out.fastq).mix(<module>.out.other)
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 see, that makes sense because samtools/bam2fq
had one emit with all reads.
I will test and update the PR.
I will merge this PR. |
This PR closes #257
PR checklist
nf-core lint
).nextflow run . -profile 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).