

find_format: Gets quality score encoding format from the FASTQ file. Return possibilities are Sanger(/Illumina1.8),

Solexa(/Illumina1.0), Illumina1.3, and Illumina1.5. This encoding is heuristic based and may not be 100 % accurate since there is overlap in the encodings used, so it is best if you already know the format.

Description

Gets quality score encoding format from the FASTQ file. Return possibilities are Sanger(/Illumina1.8), Solexa(/Illumina1.0), Illumina1.3, and Illumina1.5. This encoding is heuristic based and may not be 100 % accurate since there is overlap in the encodings used, so it is best if you already know the format.

Usage

find_format(infile, reads_used)

Arguments

Argument	Description
infile	A string giving the path for the fastq file
reads_used	int, the number of reads to use to determine the encoding format.

Value

A string denoting the read format. Possibilities are Sanger, Solexa, Illumina1.3, and Illumina1.5.

Examples

r infile <- system.file("extdata", "10^5_reads_test.fq.gz", package = "qckitfastq") find_format(infile,100)