Implement HVR sequence alignment

2021-03-27 09:39:18 +01:00
1 changed files with 2 additions and 1 deletions
--- a/src/alignment.r
+++ b/src/alignment.r
@ -5,7 +5,8 @@ construct_dataframe <- function(data) {
  vdj_string_set <- lapply(data, FUN = Biostrings::DNAStringSet)
  vdj_dataframe <- as.data.frame(vdj_string_set)
  vdj_dataframe$hvr_region <- paste(vdj_dataframe$v_sequence,
-    vdj_dataframe$d_sequence, vdj_dataframe$j_sequence,
+    vdj_dataframe$d_sequence,
+    vdj_dataframe$j_sequence,
    sep = ""
  )
  return(vdj_dataframe)