"dereplicate fastq file, and return fastq file"
"dereplicate fastq file, and return fastq tibble"
"dereplicate fastq tibble, and return fastq file"
"dereplicate fastq tibble, and return fastq tibble"
"dereplicate fasta file, and return fasta file"
"dereplicate fasta file, and return fasta tibble"
"dereplicate fasta tibble, and return fasta file"
"dereplicate fasta tibble, and return fasta tibble"
"dereplicate fastq file, and return fastq file with relabeling"
"dereplicate fastq file, and return fastq file with relbeling sha1"
"dereplicate fastq file, and return fastq file with fastq_qout_max"
"dereplicate fastq file, and return fastq tibble with vsearch_options"

dereplicate fastq file, and return fastq file

dereplicate fastq file, and return fastq tibble

dereplicate fastq tibble, and return fastq file

dereplicate fastq tibble, and return fastq tibble

dereplicate fasta file, and return fasta file

dereplicate fasta file, and return fasta tibble

dereplicate fasta tibble, and return fasta file

dereplicate fasta tibble, and return fasta tibble

dereplicate fastq file, and return fastq file with relabeling

dereplicate fastq file, and return fastq file with relbeling sha1

dereplicate fastq file, and return fastq file with fastq_qout_max

dereplicate fastq file, and return fastq tibble with vsearch_options

test_that("get merging lengths from merging two fastq files", {
merging_lengths_df <- vs_merging_lengths(fastq_input = fastq_input,
expected_df <- readRDS(test_path("testdata", "output", "merging_lengths_fq_files.rds"))
# Remove 'plot' attribute before comparison due to errors with ggplot
test_that("get merging lengths from merging two fastq tibbles", {
expected_df <- readRDS(test_path("testdata", "output", "merging_lengths_fq_tibbles.rds"))
test_that("vs_merging_lengths creates histograms when R1 and R2 lengths vary", {
# Artificially shorten half the sequences in both inputs to force multiple lengths
fastq_input$Sequence[1:n_half] <- substr(fastq_input$Sequence[1:n_half], 1, 20)
fastq_input$Quality[1:n_half] <- substr(fastq_input$Quality[1:n_half], 1, 20)
reverse$Sequence[1:n_half] <- substr(reverse$Sequence[1:n_half], 1, 20)
reverse$Quality[1:n_half] <- substr(reverse$Quality[1:n_half], 1, 20)
test_that("pe_df tibble can be merged", {
merging_lengths_df <- vs_merging_lengths(fastq_input = fastq_input)
expected_df <- readRDS(test_path("testdata", "output", "merging_lengths_pe_df.rds"))