tests/testthat/test_get_tandem_id_from_genome_seq.R

In avidaR: A Computational Biologist’s Toolkit To Get Data From 'avidaDB'

test_that("get_tandem_id_from_genome_seq", {
  
  skip_if_offline(host = "graphdb.fortunalab.org")
  
  suppressMessages({
    triplestore <- triplestore_access$new()
    
    response <- triplestore$set_access_options(
      url = "https://graphdb.fortunalab.org",
      user = "public_avida",
      password = "public_avida",
      repository = "avidaDB_test"
    )
    
    sequence1 <- "acksdblwxnlhuxpjfmhcowivenwrcepkpwhcsdxzobuicloefcpsxpmovktqodgtpzursjelhdfxayfcvtdsqbloizzfgxypwoza"
    sequence2 <- "wakkoarhimcntsprxhcwpxvdvixzqkeezetmyvuvnwlxatoldiimqgtslvutvnquxefucqmkearisdaogqzbdywpbenlhdtwawtd"
  
    found1 <- get_tandem_id_from_genome_seq(
      genome_seq = sequence1,
      triplestore = triplestore
    )
  
    found2 <- get_tandem_id_from_genome_seq(
      genome_seq = c(sequence1, sequence2),
      seed_id = c(2,1),
      tandem_seq = TRUE,
      tandem_pos = FALSE,
      triplestore = triplestore
    )
  
    found3 <- get_tandem_id_from_genome_seq(
      genome_seq = c(sequence1, sequence2),
      seed_id = c(2,1),
      tandem_seq = FALSE,
      tandem_pos = TRUE,
      triplestore = triplestore
    )
    
    skip_if(is.null(response) || is.null(found1) || is.null(found2) || is.null(found3))
  
    expect_match(class(found1), "data.frame")
    expect_true(nrow(found1) > 0)
    expect_equal(colnames(found1), c("genome_id", "tandem_id"))
    expect_equal(found1$tandem_id[1], "tandem_4909824")
  
    expect_match(class(found2), "data.frame")
    expect_true(nrow(found2) > 0)
    expect_equal(colnames(found2), c("seed_id", "genome_id", "tandem_id", "tandem_seq"))
    expect_equal(found2$seed_id[1], "seed_2")
  
    expect_match(class(found3), "data.frame")
    expect_true(nrow(found3) > 0)
    expect_equal(colnames(found3), c("seed_id", "genome_id", "tandem_id", "tandem_pos"))
    expect_equal(grepl("|", found3$tandem_pos[1]), TRUE)    
  })

})