Description

Perform variant calling on a single tumor sample using mutect2 tumor only mode. Run the input bam file through getpileupsummarries and then calculatecontaminationto get the contamination and segmentation tables. Filter the mutect2 output vcf using filtermutectcalls and the contamination & segmentation tables for additional filtering.

Input

name
description
pattern

meta

Groovy Map containing sample information
e.g. [ id:‘test’ ]

input

list containing one BAM file, also able to take CRAM as an input

[ *.{bam/cram} ]

input_index

list containing one BAM file indexe, also able to take CRAM index as an input

[ *.{bam.bai/cram.crai} ]

fasta

The reference fasta file

*.fasta

fai

Index of reference fasta file

*.fasta.fai

dict

GATK sequence dictionary

*.dict

germline_resource

Population vcf of germline sequencing, containing allele fractions.

*.vcf.gz

germline_resource_tbi

Index file for the germline resource.

*.vcf.gz.tbi

panel_of_normals

vcf file to be used as a panel of normals.

*.vcf.gz

panel_of_normals_tbi

Index for the panel of normals.

*.vcf.gz.tbi

interval_file

File containing intervals.

*.interval_list

Output

name
description
pattern

versions

File containing software versions

versions.yml

mutect2_vcf

Compressed vcf file to be used for variant_calling.

[ *.vcf.gz ]

mutect2_tbi

Indexes of the mutect2_vcf file

[ *vcf.gz.tbi ]

mutect2_stats

Stats files for the mutect2 vcf

[ *vcf.gz.stats ]

pileup_table

File containing the pileup summary table.

*.pileups.table

contamination_table

File containing the contamination table.

*.contamination.table

segmentation_table

Output table containing segmentation of tumor minor allele fractions.

*.segmentation.table

filtered_vcf

file containing filtered mutect2 calls.

*.vcf.gz

filtered_tbi

tbi file that pairs with filtered vcf.

*.vcf.gz.tbi

filtered_stats

file containing statistics of the filtermutectcalls run.

*.filteringStats.tsv