getSNP: getSNP - finding SNPs in all samples for regions specified by...

Description Usage Arguments Value Author(s)

Description

Finding SNPs in all samples for regions specified by a BED file (sequential version) using Samtools mpileup functionality. Please make sure that samtools executables are installed and available from command line without specifing the full path(e.g. added to the PATH environment variable)

Usage

1
getSNP(covdesc, minQual, refSeqFile, bedFile = NULL,iParallel=TRUE)

Arguments

covdesc

Covdesc-like data frame - BAM files are read from row names. Similar to covdesc in simpleaffy and rnaSeqMap - tab delimited table of BAM files and groupings. This is the table of experimental design.

minQual

number specifing the minimum quality of the SNPs returned by samtools used to filter out these of poor quality. Please refer to Samtools documentation for further details on the definition of this measure.

refSeqFile

Reference sequence in FASTA format used for SNPs detection.

bedFile

BED file with same genomic coordinates as in BAMs (preferably standard hg19)

iParallel

boolean specifing whether calculations shoulb be run in parallel,true by default, The degree of parallelism is set up automagically and by default equals to the number of logical CPU cores.

Value

list of data frames, for one sample each containing coordinates of the SNP with its quality

Author(s)

Alicja Szabelska, Marek Wiewiorka, Michal Okoniewski


ampliQueso documentation built on May 2, 2019, 12:45 a.m.