Skip to content

Commit

Permalink
one particular sample cannot be analyzed by Nanoplot, skip
Browse files Browse the repository at this point in the history
  • Loading branch information
SHuang-Broad committed Feb 17, 2025
1 parent 863782f commit 9b9d449
Show file tree
Hide file tree
Showing 2 changed files with 15 additions and 14 deletions.
6 changes: 3 additions & 3 deletions wdl/pipelines/ONT/Preprocessing/ONTFlowcellWGSuBAM.wdl
Original file line number Diff line number Diff line change
Expand Up @@ -43,8 +43,8 @@ workflow ONTFlowcellWGSuBAM {
wgs_cov:
"whole genome mean coverage"

nanoplot_summ:
"Summary on alignment metrics provided by Nanoplot (todo: study the value of this output)"
# nanoplot_summ:
# "Summary on alignment metrics provided by Nanoplot (todo: study the value of this output)"
seqkit_stats:
"A few metrics output by seqkit stats"
Expand Down Expand Up @@ -102,7 +102,7 @@ workflow ONTFlowcellWGSuBAM {
File aligned_bai = FinalizeAlignedBai.gcs_path

Float wgs_cov = QCandMetrics.wgs_cov
Map[String, Float] nanoplot_summ = QCandMetrics.nanoplot_summ
# Map[String, Float] nanoplot_summ = QCandMetrics.nanoplot_summ
Map[String, Float] sam_flag_stats = QCandMetrics.sam_flag_stats

Map[String, Float]? seqkit_stats = FASTQstats.stats
Expand Down
23 changes: 12 additions & 11 deletions wdl/pipelines/TechAgnostic/Utility/AlignedBamQCandMetrics.wdl
Original file line number Diff line number Diff line change
Expand Up @@ -53,8 +53,8 @@ workflow Work {
wgs_cov:
"whole genome mean coverage"

nanoplot_summ:
"Summary on alignment metrics provided by Nanoplot (todo: study the value of this output)"
# nanoplot_summ:
# "Summary on alignment metrics provided by Nanoplot (todo: study the value of this output)"
sam_flag_stats:
"SAM flag stats"
Expand Down Expand Up @@ -100,7 +100,7 @@ workflow Work {

output {
Float wgs_cov = MosDepthWGS.wgs_cov
Map[String, Float] nanoplot_summ = NanoPlotFromBam.stats_map
# Map[String, Float] nanoplot_summ = NanoPlotFromBam.stats_map
Map[String, Float] sam_flag_stats = ParseFlagStatsJson.qc_pass_reads_SAM_flag_stats

# fingerprint
Expand Down Expand Up @@ -160,12 +160,12 @@ workflow Work {

################################
# nanoplot
call NP.NanoPlotFromBam { input: bam = bam, bai = bai, disk_type = disk_type }
FinalizationManifestLine b = object
{files_to_save: flatten([[NanoPlotFromBam.stats], NanoPlotFromBam.plots]),
is_singleton_file: false,
destination: metrics_output_dir + "/nanoplot",
output_attribute_name: "nanoplot"}
# call NP.NanoPlotFromBam { input: bam = bam, bai = bai, disk_type = disk_type }
# FinalizationManifestLine b = object
# {files_to_save: flatten([[NanoPlotFromBam.stats], NanoPlotFromBam.plots]),
# is_singleton_file: false,
# destination: metrics_output_dir + "/nanoplot",
# output_attribute_name: "nanoplot"}
###################################################################################
# OPTIONAL QC/METRICS
Expand Down Expand Up @@ -242,10 +242,11 @@ workflow Work {
# save results
###################################################################################
call SAVE.SaveFilestoDestination as FF { input:
instructions = select_all([a, b, c]),
# instructions = select_all([a, b, c]),
instructions = select_all([a, c]),
already_finalized = select_all([methyl_out]),
key_file = select_first(select_all([MosDepthWGS.summary_txt,
NanoPlotFromBam.stats,
# NanoPlotFromBam.stats,
fingerprint.fingerprint_summary,]))
}
}
Expand Down

0 comments on commit 9b9d449

Please sign in to comment.