Helper function for demultiplexingSource:
Helper function for demultiplexing sequencing reads, designed in a way to
allow for parallelization across barcodes (parallel extraction of reads by
barcode). This function takes a specific barcode (numeric index) from lists
of sample names/barcodes, a
Biostrings::DNAStringSet of barcodes by
sequence header, and a
Biostrings::QualityScaledXStringSet of reads
corresponding to the barcodes. Based on the barcode index given, it extracts
all reads for the indexed barcode and writes all the reads from that barcode
to a separate .fastq file.
extract_reads( barcodeIndex, barcodes, sampleNames, index, reads, location = "./demultiplex_fastq", rcBarcodes = TRUE, hDist = 0, quiet = TRUE )
Which barcode (integer number or index) in the barcodes or sample name to use for read extraction.
A list of all barcodes in the sequencing dataset. Correlates and in same order as
A list of sample names or identifiers associated with each barcode in the barcodes list.
Biostrings::DNAStringSetthat contains the read headers and barcode sequence for each header in the sequence slot.
Biostrings::QualityScaledXStringSetthat has the same headers and order as the index file, but contains the read sequences and their quality scores.
A directory location to store the demultiplexed read files. Defaults to generate a new subdirectory at './demultiplex_fastq'
Should the barcode indices in the barcodes list be reverse complemented to match the sequences in the index DNAStringSet? Defaults to
Uses a Hamming Distance or number of base differences to allow for inexact matches for the barcodes/indexes. Defaults to 0. Warning: if the Hamming Distance is >=1 and this leads to inexact index matches to more than one barcode, that read will be written to more than one demultiplexed read files.
Turns off most messages. Default is
Writes a single .fastq file that contains all reads whose index matches the barcode specified. This file will be written to the location directory, and will be named based on the specified sampleName and barcode, e.g. './demultiplex_fastq/SampleName1_GGAATTATCGGT.fastq.gz'
## Create temporary directory ref_temp <- tempfile() dir.create(ref_temp) ## Load example barcode, index, and read data into R session barcodePath <- system.file("extdata", "barcodes.txt", package = "MetaScope") bcFile <- read.table(barcodePath, sep = "\t", header = TRUE) indexPath <- system.file("extdata", "virus_example_index.fastq", package = "MetaScope") inds <- Biostrings::readDNAStringSet(indexPath, format = "fastq") readPath <- system.file("extdata", "virus_example.fastq", package = "MetaScope") reads <- Biostrings::readQualityScaledDNAStringSet(readPath) #> Warning: metadata columns on input DNAStringSet object were dropped ## Extract reads from the first barcode results <- extract_reads(1, bcFile[, 2], bcFile[, 1], inds, reads, rcBarcodes = FALSE, location = ref_temp) ## Extract reads from multiple barcodes more_results <- lapply(1:6, extract_reads, bcFile[, 2], bcFile[, 1], inds, reads, rcBarcodes = FALSE, location = ref_temp) ## Remove temporary directory unlink(ref_temp, recursive = TRUE)