tests/testthat/test_arguments.r

In debar: A Post-Clustering Denoiser for COI-5P Barcode Data

test_that("Input arguments are handled correctly.", {
  
  #load in a set of example sequences  
  fastq_gz_real_file = system.file('extdata/coi_sequel_data_subset.fastq.gz', package = 'debar')
  data = read_fastq(fastq_gz_real_file)
  
  #
  ex1 = denoise(data$sequence[[1]], keep_phred = FALSE, to_file = FALSE)
  #x[['adjustment_count']] ==  1
  expect_equal(ex1[['adjustment_count']], 1)
  #is.null(ex1[['phred']]) == TRUE
  expect_equal(is.null(ex1[['phred']]), TRUE)
  
  super_short = 'GGTAGGAGCGTGTATACAGCGGCAGTCGAACATGTAGCTGACTCAGGTCACATTCAACAAATCATAAAGATAT'
  ex2 = denoise(super_short, keep_phred = FALSE, to_file = FALSE)
  
  #ex2$reject == TRUE
  expect_equal(ex2$reject, TRUE)
  
  ex3 = denoise(super_short, keep_phred = FALSE, terminate_rejects = FALSE, to_file = FALSE)
  #ex3$reject== FALSE
  expect_equal(ex3$reject, FALSE)
  #ex3$outseq == toupper(rev_comp(super_short))
  expect_equal(ex3$outseq, toupper(rev_comp(super_short)))
  
  #test that the warning is thrown if different lengths in the phred and DNA sequences.
  expect_error(denoise(data$sequence[[78]], phred = data$quality[[187]], to_file = FALSE), 
  "The length of the input DNA sequence and phred scores must match.")

})