There are various resource files and executables needed for this pipeline. If you are working on JUNO, you should be fine as default options will work fine for you. For other users, here are a list of resources needed in various steps in the pipeline, and their descriptions

Compile Reads

• Pooled bam directory

  • Directory containing list of donor bams (unfiltered) to be genotyped for systematic artifact filtering

  • Default:/work/access/production/resources/msk-access/current/novaseq_curated_duplex_bams_dmp/current/

• Fasta

  • Hg19 human reference fasta

  • Default:/work/access/production/resources/reference/current/Homo_sapiens_assembly19.fasta

• Genotyper

  • Path to the GBCMS genotyper executable

  • Default: /ifs/work/bergerm1/Innovation/software/maysun/GetBaseCountsMultiSample/GetBaseCountsMultiSample

• DMP IMPACT Github Repository

  • Repository of DMP IMPACT data updated daily through the cbio enterprise github

  • Default: /juno/work/access/production/resources/cbioportal/current/mskimpact

• DMP IMPACT raw data

  • Mirror bam directory

    • Directory containing list of DMP IMPACT bams

    • Default: /juno/res/dmpcollab/dmpshare/share/irb12_245/

  • Mirror bam key file -- ONLY 'IM' (SOLID TISSUE) SAMPLES ARE GENOTYPED

    • File containing DMP ID - BAM ID mapping

    • Default: /juno/res/dmpcollab/dmprequest/12-245/key.txt

  • Need to talk to Aijaz Syed about 12-245 access

>

Filter Calls

• CH list

  • list of signed out CH calls from DMP

  • Default: /juno/work/access/production/resources/dmp_signedout_CH/current/signedout_CH.txt

SV Incorporation

• DMP IMPACT Github Repository

  • Repository of DMP IMPACT data updated daily through the cbio enterprise github

  • Default: /juno/work/access/production/resources/cbioportal/current/mskimpact

1. Separating copy number output into individual files

Rscript cna_divide_by_sample.R -i /path/to/input/seg_clusp.txt -o /output/directory

The third step takes all the SV variants from all samples within each patient and present them in the same format as SNVs and incorporate SVs in the patient level table.

Usage

Default

What SV_incorporation.R does

1. 1. Only SVs implicating any ACCESS SV calling key genes are retained

4. 1. Not Called

   2. Not Covered -- none of the genes in key genes

   3. Called

5. Read in SNV table, row-bind with SV table, write out table

The pipeline aims to generate uniform and useful outputs for analyst in preliminary stage of analysis.

Example command to run through the pipeline:

Master reference file

An example of can be found in the data/ folder

Required Columns for maf file

1. Installing conda

Conda installation tutorial can be found here

2. Creating conda environment and installing R/python packages

conda create --name access_data_analysis python=3
conda activate access_data_analysis
conda install r-essentials r-base r-argparse r-ggpubr r-ggthemes r-plotly r-kableextra r-htmlwidgets r-dt
pip install genotype-variants

The final step takes the processed data from the previous steps and plots the genomic alterations over all samples of each patient. The report includes several sections with interactive plots:

1. Patient information

The first section displays the patient ID, DMP id (if provided), tumor type (if provided), and each sample. Any provided sample meta-information is also display for each sample.

2. Plot of SNV variant allele frequencies

The second section shows SNV/INDEL events are plotted out by VAFs over timepoints. Above the panel it also display sample timepoint annotation, such as treatment information (if provided). If you provide IMPACT sample information, it will segregate each mutation by whether it is known to be clonal in IMPACT, subclonal in IMPACT, or is present in ACCESS only. There are additional tabs that display a table of mutation data and methods description.

3. Plot of copy number alterations

The third section shows CNAs that are plotted by fold-change(fc) for each ACCESS sample and gene. If there are no CNAs, then this section is not displayed.

4. Plot of clonal SNV/INDEL VAFs adjusted for copy number

If you provided an IMPACT sample, this last section will show SNV/INDEL events that are plotted out by VAFs over timepoints. However, the VAFs are corrected for IMPACT copy number information. Details of the method are shown under the Description tab in this section. Similar to section 2, sample timepoint annotations are shown above the plot.

Usage

Rscript reports/create_report.R -h                                      
usage: reports/create_report.R [-h] -t TEMPLATE -p PATIENT_ID -r RESULTS -rc
                               CNA_RESULTS_DIR -tt TUMOR_TYPE -m METADATA
                               [-d DMP_ID] [-ds DMP_SAMPLE_ID] [-dm DMP_MAF]
                               [-o OUTPUT] [-ca] [-pi]

optional arguments:
  -h, --help            show this help message and exit
  -t TEMPLATE, --template TEMPLATE
                        Path to Rmarkdown template file.
  -p PATIENT_ID, --patient-id PATIENT_ID
                        Patient ID
  -r RESULTS, --results RESULTS
                        Path to CSV file containing mutation and genotype
                        results for the patient.
  -rc CNA_RESULTS_DIR, --cna-results-dir CNA_RESULTS_DIR
                        Path to directory containing CNA results for the
                        patient.
  -tt TUMOR_TYPE, --tumor-type TUMOR_TYPE
                        Tumor type
  -m METADATA, --metadata METADATA
                        Path to file containing meta data for each sample.
                        Should contain a 'cmo_sample_id_plasma', 'sex', and
                        'collection_date' columns. Can also optionally include
                        a 'timepoint' column (e.g. for treatment information).
  -d DMP_ID, --dmp-id DMP_ID
                        DMP patient ID (optional).
  -ds DMP_SAMPLE_ID, --dmp-sample-id DMP_SAMPLE_ID
                        DMP sample ID (optional).
  -dm DMP_MAF, --dmp-maf DMP_MAF
                        Path to DMP MAF file (optional).
  -o OUTPUT, --output OUTPUT
                        Output file
  -ca, --combine-access
                        Don't splite VAF plots by clonality.
  -pi, --plot-impact    Also plot VAFs from IMPACT samples.

This step generates a final CNA call set for plotting. This consists of:

• Calls passing de novo CNA calling threshold

  • Significant adjusted p value ( <= 0.05)

  • Significant fold change ( > 1.5 or < -1.5)

• Calls that can be rescued based on prior knowledge from IMPACT samples

  • Significant adjusted p value ( <= 0.05)

  • Lowered threshold for fold change ( < 1.2 or < -1.2)

Usage

Rscript R/CNA_processing.R -h                                       
usage: R/CNA_processing.R [-h] [-m MASTERREF] [-o RESULTSDIR] [-dmp DMPDIR]

optional arguments:
  -h, --help            show this help message and exit
  -m MASTERREF, --masterref MASTERREF
                        File path to master reference file
  -o RESULTSDIR, --resultsdir RESULTSDIR
                        Output directory
  -dmp DMPDIR, --dmpdir DMPDIR
                        Directory of clinical DMP IMPACT repository [default]

What CNA_processing.R does

1. Process DMP CNAs

2. Process ACCESS CNAs

3. CNA Calling (de novo and rescue)

Format of the final call set:

The second step takes all the genotypes generated from the first step and organized into a patient level variants table with VAFs and call status for each variant of each sample.

Each call is subjected to:

1. Read depth filter (hotspot vs non-hotspot)

2. Systematic artifact filter

3. Germline filters

   1. If any normal exist -- (buffy coat and DMP normal) 2:1 rule

   2. If not -- exac freq < 0.01% and VAF < 30%

4. CH tag

Usage

Rscript R/filter_calls.R -h                                         
usage: R/filter_calls.R [-h] [-m MASTERREF] [-o RESULTSDIR] [-dmpk DMPKEYPATH]
                        [-ch CHLIST] [-c CRITERIA]

optional arguments:
  -h, --help            show this help message and exit
  -m MASTERREF, --masterref MASTERREF
                        File path to master reference file
  -o RESULTSDIR, --resultsdir RESULTSDIR
                        Output directory
  -ch CHLIST, --chlist CHLIST
                        List of signed out CH calls [default]
  -c CRITERIA, --criteria CRITERIA
                        Calling criteria [default]

Default

Default options can be found here

What filter_calls.R does

Generate a reference of systematic artifacts -- any call with occurrence in more than or equal to 2 donor samples (occurrence defined as more than or equal to 2 duplex reads)

For each patient

1. Read in sample sheets -- reference for downstream analysis

2. Generate a preliminary patient level variants table

3. Read in and merging in hotspots, DMP signed out calls and occurrence in donor samples

4. Call status annotation

   1. All call passing read depth/genotype filter annotated as 'Called' or 'Genotyped'

   2. Any call not satisfying germline filters are overwritten with 'Not Called'

      1. Calls with zero coverage in plasma sample also annotated as 'Not Covered'

5. Final processing

   1. Combining duplex and simplex read counts

   2. CH tags

6. Write out table

Example of the patient level table:

There are two variantion:

• compile_reads.R : Works with Research ACCESS and Clinical IMPACT

• compile_reads_all.R: Works with Research ACCESS, Clinical ACCESS and Clinical IMPACT

The first step of the pipeline is to genotype all the variants of interest in the included samples (this means plasma, buffy coat, DMP tumor, DMP normal, and donor samples). Once we obtained the read counts at every loci of every sample, we then generate a table of VAFs and call status for each variant in all samples within a patient in the next step.

Usage compile_reads.R

Rscript R/compile_reads.R -h                                        
usage: R/compile_reads.R [-h] [-m MASTERREF] [-o RESULTSDIR]
                         [-pb POOLEDBAMDIR] [-fa FASTAPATH]
                         [-gt GENOTYPERPATH] [-dmp DMPDIR] [-mb MIRRORBAMDIR]
                         [-dmpk DMPKEYPATH]

optional arguments:
  -h, --help            show this help message and exit
  -m MASTERREF, --masterref MASTERREF
                        File path to master reference file
  -o RESULTSDIR, --resultsdir RESULTSDIR
                        Output directory
  -pb POOLEDBAMDIR, --pooledbamdir POOLEDBAMDIR
                        Directory for all pooled bams [default]
  -fa FASTAPATH, --fastapath FASTAPATH
                        Reference fasta path [default]
  -gt GENOTYPERPATH, --genotyperpath GENOTYPERPATH
                        Genotyper executable path [default]
  -dmp DMPDIR, --dmpdir DMPDIR
                        Directory of clinical DMP IMPACT repository [default]
  -mb MIRRORBAMDIR, --mirrorbamdir MIRRORBAMDIR
                        Mirror BAM file directory [default]
  -dmpk DMPKEYPATH, --dmpkeypath DMPKEYPATH
                        DMP mirror BAM key file [default]

Usage compile_reads_all.R

Rscript R/compile_reads_all.R -h
usage: R/compile_reads_all.R [-h] [-m MASTERREF] [-o RESULTSDIR]
                             [-pid PROJECTID] [-pb POOLEDBAMDIR]
                             [-fa FASTAPATH] [-gt GENOTYPERPATH] [-dmp DMPDIR]
                             [-mb MIRRORBAMDIR] [-mab MIRRORACCESSBAMDIR]
                             [-dmpk DMPKEYPATH] [-dmpak DMPACCESSKEYPATH]

optional arguments:
  -h, --help            show this help message and exit
  -m MASTERREF, --masterref MASTERREF
                        File path to master reference file
  -o RESULTSDIR, --resultsdir RESULTSDIR
                        Output directory
  -pid PROJECTID, --projectid PROJECTID
                        Project ID for submitted jobs involved in this run
  -pb POOLEDBAMDIR, --pooledbamdir POOLEDBAMDIR
                        Directory for all pooled bams [default]
  -fa FASTAPATH, --fastapath FASTAPATH
                        Reference fasta path [default]
  -gt GENOTYPERPATH, --genotyperpath GENOTYPERPATH
                        Genotyper executable path [default]
  -dmp DMPDIR, --dmpdir DMPDIR
                        Directory of clinical DMP repository [default]
  -mb MIRRORBAMDIR, --mirrorbamdir MIRRORBAMDIR
                        Mirror BAM file directory [default]
  -mab MIRRORACCESSBAMDIR, --mirroraccessbamdir MIRRORACCESSBAMDIR
                        Mirror BAM file directory for MSK-ACCESS [default]
  -dmpk DMPKEYPATH, --dmpkeypath DMPKEYPATH
                        DMP mirror BAM key file [default]
  -dmpak DMPACCESSKEYPATH, --dmpaccesskeypath DMPACCESSKEYPATH
                        DMP mirror BAM key file for MSK-ACCESS [default]

Default

Default options can be found here

What compile_reads does

For each patient

• Create a sample sheet -- similar to the one for genotype-variants

• Generate all variants of interests

  • DMP calls from cbio repo

  • ACCESS calls from SNV pipeline

• Generate unique variants list

• Tag hotspots on unique variants

• Genotype with genotype-variants

Afterwards, for donor bams

• Obtain all variants genotyped in any patient, generate a all unique list of variants

• Genotype with genotype-variants

Tool to do the following operations:

• Reads meta data file, and based on the timepoint information given, convert them to days for a samples belonging to a given patient_id

• Supports following date formats:

  • MM/DD/YY

  • M/D/YY

  • MM/D/YY

  • M/DD/YY

  • MM/DD/YYYY

  • YYYY/MM/DD

Requirements

• pandas

• typing

• arrow

Example command

Usage

There are intermediate files generated with each step in the /output/directory , here is a diagram for its organization

Table of Contents

get_cbioportal_variants

Requirement:

• pandas

• typing

• typer

• bed_lookup()

Example command

subset_cpt

subset_cst

subset_cna

subset_sv

subset_maf

Sub-modules

read_tsv

Read a tsv file

Arguments:

• maf File - Input MAF/tsv like format file

Returns:

• data_frame - Output a data frame containing the MAF/tsv

read_ids

make a list of ids

Arguments:

• sid tuple - Multiple ids as tuple

• ids File - File containing multiple ids

Returns:

• list - List containing all ids

filter_by_columns

Filter data by columns

Arguments:

• sid list - list of columns to subset over

• tsv_df data_frame - data_frame to subset from

Returns:

• data_frame - A copy of the subset of the data_frame

filter_by_rows

Filter the data by rows

Arguments:

• sid list - list of row names to subset over

• tsv_df data_frame - data_frame to subset from

• col_name string - name of the column to filter using names in the sid

Returns:

• data_frame - A copy of the subset of the data_frame

read_bed

Read BED file using bed_lookup

Arguments:

• bed file - File ins BED format to read

Returns:

object : bed file object to use for filtering

check_if_covered

Function to check if a variant is covered in a given bed file

Arguments:

• bedObj object - BED file object to check coverage

• mafObj data_frame - data frame to check coverage against coordinates using column 'Chromosome' and position column is 'Start_Position'

Returns:

• data_frame - description

get_row

Function to skip rows

Arguments:

• tsv_file file - file to be read

Returns:

• list - lines to be skipped

The Tool does the following operations:

• Read one or more files from the inputs

• Removes unwanted columns, modifying the column headers depending on the requirements

• Massaging the data frame to make it compatible with MAF format

• Write the data frame to a file in MAF format and Excel format

Requirements

• pandas

• openpyxl

• typing

• typer

Example command

Explicitly specifying files on command line

Specifying files in a text FileOfFiles

where FileOfFiles.txt

Keeping normal samples identified using "normal" string, by default they are filtered

Usage

Scripts for downstream analysis plotting of ACCESS variant calling pipeline output

This gitbook will walk you through:

• Requirements

• run_create_report

  • Main Script (run_create_report.py)

• Submodules

Requirements

access_data_analysis=>0.1.2 # works with this repo tag
typer==0.3.2
typing_extensions==3.10.0.0
pandas==1.2.5
rich==12.1.0

run_create_report

Main Script (run_create_report.py)

Usage: run_create_report.py [OPTIONS]

Options:
  -r, --repo PATH                 Base path to where the git repository is
                                  located for access_data_analysis

  -s, --script PATH               Path to the create_report.R script, fall
                                  back if `--repo` is not given

  -t, --template PATH             Path to the template.Rmd or
                                  template_days.Rmd to be used with
                                  create_report.R when `--repo` is not given

  -m, --manifest FILE             File containing meta information per sample.
                                  Require following columns in the header:
                                  cmo_patient_id, sample_id, dmp_patient_id,
                                  collection_date or collection_day,
                                  timepoint. If dmp_sample_id column is given
                                  and has information that will be used to run
                                  facets. If dmp_sample_id is not given and
                                  dmp_patient_id is given than it will be used
                                  to get the Tumor sample with lowest number.
                                  If dmp_sample_id or dmp_patient_id is not
                                  given then it will run without the facet maf
                                  file  [required]

  -v, --variant-results DIRECTORY
                                  Base path for all results of small variants
                                  as generated by filter_calls.R script in
                                  access_data_analysis (Make sure only High
                                  Confidence calls are included)  [required]

  -c, --cnv-results DIRECTORY     Base path for all results of CNV as
                                  generated by CNV_processing.R script in
                                  access_data_analysis  [required]

  -f, --facet-repo DIRECTORY      Base path for all results of facets on
                                  Clinical MSK-IMPACT samples  [default: /juno
                                  /work/ccs/shared/resources/impact/facets/all
                                  /]

  -bf, --best-fit                 If this is set to True then we will attempt
                                  to parse `facets_review.manifest` file to
                                  pick the best fit for a given dmp_sample_id
                                  [default: False]

  -l, --tumor-type TEXT           Tumor type label for the report  [required]
  -cfm, --copy-facet-maf          If this is set to True then we will copy the
                                  facet maf file in the directory specified in
                                  `copy_facet_dir`  [default: False]

  -cfd, --copy-facet-dir PATH     Directory path where the facet maf file
                                  should be copied.

  -d, --template-days             If the `--repo` option is specified and if
                                  this is set to True then we will use the
                                  template_days RMarkdown file as the template
                                  [default: False]

  -gm, --generate-markdown        If given, the create_report.R will be run
                                  with `-md` flag to generate markdown
                                  [default: False]

  -ff, --force                    If this is set to True then we will not stop
                                  if an error is encountered in a given sample
                                  while running create_report.R but keep on
                                  running for the next sample  [default:
                                  False]

  --install-completion            Install completion for the current shell.
  --show-completion               Show completion for the current shell, to
                                  copy it or customize the installation.

  --help                          Show this message and exit.

Wrapper script to run create_report.R

Arguments:

• repo_path Path, optional - "Base path to where the git repository is located for access_data_analysis".

• script_path Path, optional - "Path to the create_report.R script, fall back if --repo is not given".

• template_path Path, optional - "Path to the template.Rmd or template_days.Rmd to be used with create_report.R when --repo is not given".

• manifest Path, required - "File containing meta information per sample. Require following columns in the header: cmo_patient_id, sample_id, dmp_patient_id, collection_date or collection_day, timepoint. If dmp_sample_id column is given and has information that will be used to run facets. if dmp_sample_id is not given and dmp_patient_id is given than it will be used to get the Tumor sample with lowest number.If dmp_sample_id or dmp_patient_id is not given then it will run without the facet maf file".

• variant_path Path, required - "Base path for all results of small variants as generated by filter_calls.R script in access_data_analysis (Make sure only High Confidence calls are included)".

• cnv_path Path, required - "Base path for all results of CNV as generated by CNV_processing.R script in access_data_analysis".

• facet_repo Path, required - "Base path for all results of facets on Clinical MSK-IMPACT samples".

• best_fit bool, optional - "If this is set to True then we will attempt to parse facets_review.manifest file to pick the best fit for a given dmp_sample_id".

• tumor_type str, required - "Tumor type label for the report".

• copy_facet bool, optional - "If this is set to True then we will copy the facet maf file in the directory specified in copy_facet_dir".

• copy_facet_dir Path, optional - "Directory path where the facet maf file should be copied.".

• template_days bool, optional - "If the --repo option is specified and if this is set to True then we will use the template_days RMarkdown file as the template".

• markdown bool, optional - "If given, the create_report.R will be run with -md flag to generate markdown".

• force bool, optional - "If this is set to True then we will not stop if an error is encountered in a given sample but keep on running for the next sample".

Usage

• Using Generate Markdown, copy facet maf file, use template_days RMarkdown, force flag and best fit for facets

> python python/run_create_report/run_create_report.py \
-m /home/shahr2/bergerlab/Project_10619_D/small_variants/manifest_noDate_days.tsv \
-r /home/shahr2/github/access_data_analysis \
-v /home/shahr2/bergerlab/Project_10619_D/small_variants/results_20Jan2023/results_stringent/ \
-c /home/shahr2/bergerlab/Project_10619_D/small_variants/results_20Jan2023/CNA_final_call_set \
-l "Melanoma" -gm -d -cfm -ff -bf

• Using Generate Markdown, force flag and default fit for facets

> python python/run_create_report/run_create_report.py \
-m /home/shahr2/bergerlab/Project_10619_D/small_variants/manifest_noDate_days.tsv \
-r /home/shahr2/github/access_data_analysis \
-v /home/shahr2/bergerlab/Project_10619_D/small_variants/results_20Jan2023/results_stringent/ \
-c /home/shahr2/bergerlab/Project_10619_D/small_variants/results_20Jan2023/CNA_final_call_set \
-l "Melanoma" -gm -ff

Submodules

check_required_columns

check_required_columns

def check_required_columns(manifest, template_days=None)

Check if all required columns are present in the sample manifest file

Arguments:

• manifest data_frame - meta information file with information for each sample

• template_days bool - True|False if template days RMarkdown will be used

Raises:

• typer.Abort - if "cmo_patient_id" column not provided

• typer.Abort - if "cmo_sample_id/sample_id" column not provided

• typer.Abort - if "dmp_patient_id" column not provided

• typer.Abort - if "collection_date/collection_day" column not provided

• typer.Abort - if "timepoint" column not provided

Returns:

• list - column name for the manifest file

• data_frame - data_frame with unique ids to traverse over

generate_repo_paths

generate_repo_path

def generate_repo_path(repo_path=None, script_path=None, template_path=None, template_days=None)

Generate path to create_report.R and template RMarkdown file

Arguments:

• repo_path pathlib.Path, optional - Path to clone of git repo access_data_analysis. Defaults to None.

• script_path pathlib.Path, optional - Path to create_report.R. Defaults to None.

• template_path pathlib.Path, optional - Path to template RMarkdown file. Defaults to None.

• template_days bool, optional - True|False to use days template if using repo_path. Defaults to None.

Raises:

• typer.Abort - Abort if both repo_path and script_path are not given

• typer.Abort - Abort if both repo_path and template_path are not given

Returns:

• str - Path to create_report.R and path to template markdown file

read_manifest

read_manifest

def read_manifest(manifest)

Read manifest file

Arguments:

• manifest pathlib.PATH - description

Returns:

• data_frame - description

get_row

def get_row(tsv_file)

Function to skip rows

Arguments:

• tsv_file file - file to be read

Returns:

• list - lines to be skipped

get_small_variant_csv

get_small_variant_csv

def get_small_variant_csv(patient_id, csv_path)

Get the path to CSV file to be used for a given patient containing all variants

Arguments:

• patient_id str - patient id used to identify the csv file

• csv_path pathlib.path - base path where the csv file is expected to be present

Raises:

• typer.Abort - if no csv file is returned

• typer.Abort - if more then one csv file is returned

Returns:

• str - path to csv file containing the variants

run_cmd

run_cmd

def run_cmd(cmd)

Given a system command run it using subprocess

Arguments:

• cmd str - System command to be run as a string

run_multiple_cmd

def run_multiple_cmd(commands, parallel_process=None)

Given a system command run it using subprocess

Arguments:

• cmd list[str] - list of system commands to be run

generate_facet_maf_path

generate_facet_maf_path

def generate_facet_maf_path(facet_path, patient_id, sample_id=None)

Get path of maf associated with facet-suite output

Arguments:

• facet_path pathlib.PATH|str - path to search for the facet file

• patient_id str - patient id to be used to search, default is set to None

• sample_id str - sample id to be used to search, default is set to None

Returns:

• str - path of the facets maf

get_maf_path

def get_maf_path(maf_path, patient_id, sample_id)

Get the path to the maf file

Arguments:

• maf_path pathlib.Path - Base path of the maf file

• patient_id str: DMP Patient ID for facets

• sample_id str - DMP Sample ID if any for facets

Returns:

• str - Path to the maf file

get_best_fit_folder

def get_best_fit_folder(facet_manifest_path)

Get the best fit folder for the given facet manifest path

Arguments:

• facet_manifest_path str - manifest path to be used for determining best fit

Returns:

• pathlib.Path - path to the folder containing best fit maf files

generate_create_report_cmd

generate_create_report_cmd

def generate_create_report_cmd(script, markdown, template_file, cmo_patient_id, csv_file, manifest, cnv_path, dmp_patient_id, dmp_sample_id, dmp_facet_maf, tumor_type=None)

Create the system command that should be run for create_report.R

Arguments:

• script str - path for create_report.R

• markdown bool - True|False to generate markdown output

• template_file str - path for the template file

• cmo_patient_id str - patient id from CMO

• csv_file str - path to csv file containing variant information

• tumor_type str - tumor type label

• manifest pathlib.Path - path to the manifest containing meta data

• cnv_path pathlib.Path - path to directory having cnv files

• dmp_patient_id str - patient id of the clinical msk-impact sample

• dmp_sample_id str - sample id of the clinical msk-impact sample

• dmp_facet_maf str - path to the clinical msk-impact maf file annotated for facets results

Returns:

• cmd str - system command to run for create_report.R

• html_output pathlib.Path - where the output file should be written