seqDepth: Calculating sequencing depths from the output of 'seqkit...

View source: R/seqDepth.R

seqDepthR Documentation

Calculating sequencing depths from the output of 'seqkit stats'

Description

This function calculates the sequencing depth of each genome using the output of command 'seqkit stats' (https://github.com/shenwei356/seqkit)

Usage

seqDepth(reads, len, srt = TRUE)

Arguments

reads

A data frame imported from the output file of 'seqkit stats'. The first column of input file names must not contain paths. For example, sample_1.fastq.gz is valid whereas reads/filtered/sample_1.fastq.gz is invalid. Each FASTQ file should follow the pattern *_1, 2.fastq.gz.

len

Length (bp) of the reference genome

srt

Whether to sort rows by sequencing depths (default: TRUE)

Value

A data frame of isolate names and sequencing dapths


wanyuac/handyR documentation built on June 10, 2024, 1:24 a.m.