diff --git a/subworkflows/nf-core/bam_qc_picard/main.nf b/subworkflows/nf-core/bam_qc_picard/main.nf
new file mode 100644
index 00000000..31edb2c3
--- /dev/null
+++ b/subworkflows/nf-core/bam_qc_picard/main.nf
@@ -0,0 +1,43 @@
+//
+// Run QC steps on BAM/CRAM files using Picard
+//
+
+params.options = [:]
+
+include { PICARD_COLLECTMULTIPLEMETRICS } from '../../../modules/picardcollectmultiplemetrics/main' addParams( options: params.options )
+include { PICARD_COLLECTWGSMETRICS      } from '../../../modules/picardcollectwgsmetrics/main'      addParams( options: params.options )
+include { PICARD_COLLECTHSMETRICS       } from '../../../modules/picardcollecthsmetrics/main'       addParams( options: params.options )
+
+workflow BAM_QC_PICARD {
+    take:
+    ch_bam_bai          // channel: [ val(meta), [ bam ], [bai/csi] ]
+    ch_fasta            // channel: [ fasta ]
+    ch_bait_interval    // channel: [ bait_interval ]
+    ch_target_interval  // channel: [ target_interval ]
+
+    main:
+    ch_versions = Channel.empty()
+
+    PICARD_COLLECTMULTIPLEMETRICS( ch_bam_bai, ch_fasta] )
+    ch_versions = ch_versions.mix(PICARD_COLLECTMULTIPLEMETRICS.out.versions.first())
+    if (!ch_bait_interval.isEmpty() || !ch_target_interval.isEmpty()) {
+        if (ch_bait_interval.isEmpty()) {
+            throw new Error("Bait interval channel is empty")
+        }
+        if (ch_target_interval.isEmpty()) {
+            throw new Error("Target interval channel is empty")
+        }
+        PICARD_COLLECTHSMETRICS( ch_bam_bai, ch_fasta, ch_bait_interval, ch_target_interval )
+        ch_versions = ch_versions.mix(PICARD_COLLECTHSMETRICS.out.versions.first())
+    } else {
+        PICARD_COLLECTWGSMETRICS( ch_bam_bai, ch_fasta )
+        ch_versions = ch_versions.mix(PICARD_COLLECTWGSMETRICS.out.versions.first())
+    }
+
+    emit:
+    hs_metrics          = PICARD_COLLECTHSMETRICS.out.hs_metrics    // channel: [ val(meta), [ hs_metrics ] ]
+    wgs_metrics         = PICARD_COLLECTWGSMETRICS.out.metrics      // channel: [ val(meta), [ wgs_metrics ] ]
+    multiple_metrics    = PICARD_COLLECTMULTIPLEMETRICS.out.metrics // channel: [ val(meta), [ multiple_metrics ] ]
+
+    versions = ch_versions                                          // channel: [ versions.yml ]
+}
diff --git a/subworkflows/nf-core/bam_qc_picard/meta.yml b/subworkflows/nf-core/bam_qc_picard/meta.yml
new file mode 100644
index 00000000..77104e82
--- /dev/null
+++ b/subworkflows/nf-core/bam_qc_picard/meta.yml
@@ -0,0 +1,64 @@
+name: bam_qc
+description: Produces comprehensive statistics from BAM file
+keywords:
+  - statistics
+  - counts
+  - hs_metrics
+  - wgs_metrics
+  - bam
+  - sam
+  - cram
+modules:
+  - picard/collectmultiplemetrics
+  - picard/collectwgsmetrics
+  - picard/collecthsmetrics
+input:
+  - meta:
+      type: map
+      description: |
+        Groovy Map containing sample information
+        e.g. [ id:'test', single_end:false ]
+  - bam:
+      type: file
+      description: BAM/CRAM/SAM file
+      pattern: "*.{bam,cram,sam}"
+  - bai:
+      type: file
+      description: Index for BAM/CRAM/SAM file
+      pattern: "*.{bai,crai,sai}"
+  - fasta:
+      type: optional file
+      description: Reference file the CRAM was created with
+      pattern: "*.{fasta,fa}"
+  - bait_intervals:
+      type: optional file
+      description: An interval list file that contains the locations of the baits used.
+      pattern: "baits.interval_list"
+  - target_intervals:
+      type: optional file
+      description: An interval list file that contains the locations of the targets.
+      pattern: "targets.interval_list"
+output:
+  - meta:
+      type: map
+      description: |
+        Groovy Map containing sample information
+        e.g. [ id:'test', single_end:false ]
+  - hs_metrics:
+      type: file
+      description: Alignment metrics files generated by picard CollectHsMetrics
+      pattern: "*_collecthsmetrics.txt"
+  - wgs_metrics:
+      type: file
+      description: Alignment metrics files generated by picard CollectWgsMetrics
+      pattern: "*_{metrics}"
+  - multiple_metrics:
+      type: file
+      description: Alignment metrics files generated by picard CollectMultipleMetrics
+      pattern: "*_{metrics}"
+  - versions:
+      type: file
+      description: File containing software versions
+      pattern: "versions.yml"
+authors:
+  - "@matthdsm"
diff --git a/tests/subworkflows/nf-core/bam_qc_picard/main.nf b/tests/subworkflows/nf-core/bam_qc_picard/main.nf
new file mode 100644
index 00000000..d88f2bf9
--- /dev/null
+++ b/tests/subworkflows/nf-core/bam_qc_picard/main.nf
@@ -0,0 +1,21 @@
+#!/usr/bin/env nextflow
+
+nextflow.enable.dsl = 2
+
+include { BAM_QC_PICARD } from '../../../../subworkflows/nf-core/bam_qc_picard/main' addParams([:])
+
+workflow test_bam_qc_picard_wgs {
+    input = [ [ id:'test', single_end:false ], // meta map
+                file(params.test_data['sarscov2']['illumina']['test_single_end_bam'], checkIfExists: true)
+            ]
+
+    BAM_QC_PICARD ( input, [], [], [] )
+}
+
+workflow test_bam_qc_picard_targetted {
+    input = [ [ id:'test', single_end:false ], // meta map
+                file(params.test_data['sarscov2']['illumina']['test_paired_end_bam'], checkIfExists: true)
+            ]
+
+    BAM_QC_PICARD ( input, [], file(params.test_data['sarscov2']['genome']['baits_interval_list'], checkIfExists: true), file(params.test_data['sarscov2']['genome']['targets_interval_list'], checkIfExists: true) )
+}
diff --git a/tests/subworkflows/nf-core/bam_qc_picard/test.yml b/tests/subworkflows/nf-core/bam_qc_picard/test.yml
new file mode 100644
index 00000000..af98230d
--- /dev/null
+++ b/tests/subworkflows/nf-core/bam_qc_picard/test.yml
@@ -0,0 +1,33 @@
+- name: bam qc picard wgs
+  command: nextflow run ./tests/subworkflows/nf-core/bam_qc_picard -entry test_bam_qc_picard_wgs -c tests/config/nextflow.config
+  tags:
+    - subworkflows
+    # - subworkflows/bam_qc_picard
+    # Modules
+    # - picard
+    # - picard/collectmultiplemetrics
+    # - picard/collectwgsmetrics
+  files:
+    - path: ./output/picard/test.CollectMultipleMetrics.alignment_summary_metrics
+    - path: ./output/picard/test.CollectMultipleMetrics.insert_size_metrics
+    - path: ./output/picard/test.CollectMultipleMetrics.base_distribution_by_cycle_metrics
+    - path: ./output/picard/test.CollectMultipleMetrics.quality_by_cycle_metrics
+    - path: ./output/picard/test.CollectMultipleMetrics.quality_distribution_metrics
+    - path: ./output/picard/test.CollectWgsMetrics.coverage_metrics
+
+- name: bam qc picard targetted
+  command: nextflow run ./tests/subworkflows/nf-core/bam_qc_picard -entry test_bam_qc_picard_targetted -c tests/config/nextflow.config
+  tags:
+    - subworkflows
+    # - subworkflows/bam_qc_picard
+    # Modules
+    # - picard
+    # - picard/collectmultiplemetrics
+    # - picard/collecthsmetrics
+  files:
+    - path: ./output/picard/test.CollectMultipleMetrics.alignment_summary_metrics
+    - path: ./output/picard/test.CollectMultipleMetrics.insert_size_metrics
+    - path: ./output/picard/test.CollectMultipleMetrics.base_distribution_by_cycle_metrics
+    - path: ./output/picard/test.CollectMultipleMetrics.quality_by_cycle_metrics
+    - path: ./output/picard/test.CollectMultipleMetrics.quality_distribution_metrics
+    - path: ./output/picard/test_collecthsmetrics.txt