Update similiarity_to_reference.R

R/similiarity_to_reference.R

@@ -17,7 +17,7 @@ similiarity_to_reference <- function (seqs = named_DNAStringSet,
       as_tibble(rownames = "query") %>%
       pivot_longer(cols = -query, names_to = "seqnames",
                    values_to = "sim_to_ref") %>%
-      mutate(seq_var = if_else(query == names(ref)[1], "ref", "none")) %>%
+      mutate(seq_var = if_else(query %in% names(ref), "ref", "none")) %>%
       filter(seq_var == "ref") %>%
       select(-seq_var, -query)
     }