-
Notifications
You must be signed in to change notification settings - Fork 679
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
Fix #961 #963
Fix #961 #963
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.
I feel like the NfcoreTemplate.groovy
lib is getting more and more specific to rnaseq and less TEMPLATy, should we move things out? and keep this one file only for TEMPLATE stuff?
|
def samp_aln = '' | ||
for (samp in fail_percent_mapped) { | ||
samp_aln += " ${samp.value}%: ${samp.key}\n" | ||
if (fail_mapped_count > 0 || fail_trimmed_count > 0 || fail_strand_count > 0) { |
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.
Samples that fail strandedness check are not skipped, they are still processed but a warning is printed in the multiqc report. For example a run with all 6 samples with the wrong strandedness specified in the sampelsheet, the below output is printed to terminal:
-[nf-core/rnaseq] Pipeline completed successfully with skipped sampl(es)- -[nf-core/rnaseq] Please check MultiQC report: 6/6 samples failed strandedness check.-
Can this be changed to only report 'skipped samp(les)' for samples that fail map or trim count and not strandedness. The message that samples failed strandedness check is sufficient.
A more helpful long-term solution would be to run a QC and strand check BEFORE the whole rnaseq run with alignment and counting is performed. This is more consistent with what researchers want in real life.
Closes #961