getSNP {ampliQueso} | R Documentation

## getSNP - finding SNPs in all samples for regions specified by a BED file all regions if BED file not provided (sequential and parallel versions)

### 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

    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

[Package ampliQueso version 1.18.0 Index]