radiaCalling: radia

In msubirana/ergWgsTools: What the Package Does (One Line, Title Case)

Description

Tumor-normal pair somatic variant calling using RADIA

Usage

radiaCalling(
  tumor_file,
  normal_file,
  sample_name,
  ref,
  out_path,
  radia_path,
  python_radia,
  samtools,
  threads,
  bcftools
)

Arguments

tumor_file	Tumor bam to file to perform the variant calling.
normal_file	Normal bam to file to perform the variant calling.
ref	Path for the reference genome to use for the alignment (fasta format) and the corresponding indexes generated with bwa index and a dictionary index file generated by CreateSequenceDictionary gatk tool.
out_path	Path where the output of the analysis will be saved.
python_radia	Path to the python binary with all the RADIA prerequisites.
samtools	Path of samtools binary.
threads	Number of threads to use in the analysis.

Examples

## Not run: 
radiaCalling(tumor_file = 'raw/sample_tumor.bam',
               normal_file = 'raw/sample_normal.bam',
               sample_name = 'sample'
               ref = 'ref/hg38.fa',
               out_path = 'rst',
               radia_path = '/bin/radia/scripts',
               python_radia = 'venv/radia/bin/python',
               samtools = '/bin/samtools')

## End(Not run)

msubirana/ergWgsTools documentation built on June 8, 2020, 8:07 a.m.