star: Running Star two steps for variant calls

In kendomaniac/docker4seq: Running NGS computing demanding applications, e.g reads mapping and counting, wrapped in docker containers

Running Star two steps for variant calls

Description

This function executes the two steps STAR as sugested by best practice GATK for calling variants on RNAseq data only PE data are accepted

Usage

star(
  group = c("sudo", "docker"),
  fastq.folder = getwd(),
  scratch.folder = "/data/scratch",
  genome.folder,
  groupid,
  threads = 1,
  version = 2
)

Arguments

group	a character string. Two options: "sudo" or "docker", depending to which group the user belongs
fastq.folder	a character string indicating where gzip fastq files are located
scratch.folder	a character string indicating the scratch folder where docker container will be mounted
genome.folder	a character string indicating the folder where the indexed reference genome for STAR is located.
groupid	a character string to be inserted in the bam as identifier for the sample
threads	a number indicating the number of cores to be used from the application
version	a number indicating version in use

Value

three files: sorted_reads.bam, which is sorted and duplicates marked bam file, sort_reads.bai, and sort_reads.stats, which provides mapping statistics

Author(s)

Raffaele Calogero, raffaele.calogero [at] unito [dot] it, Bioinformatics and Genomics unit, University of Torino Italy

Examples

## Not run: 
    #downloading fastq files
    system("wget http://130.192.119.59/public/test_R1.fastq.gz")
    system("wget http://130.192.119.59/public/test_R2.fastq.gz")
    #running star nostrand pe
    star(group="docker",fastq.folder=getwd(), scratch.folder="/data/scratch",
    genome.folder="/data/scratch/hg38star", groupid="test", threads=12)


## End(Not run)

kendomaniac/docker4seq documentation built on Sept. 3, 2024, 6:42 p.m.