From 91b3e37bd8117b686d9a235cd057edb31662a280 Mon Sep 17 00:00:00 2001 From: coolneng Date: Thu, 13 May 2021 19:06:58 +0200 Subject: [PATCH] Start j_alignment with the portion after the Cys --- src/alignment.r | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/alignment.r b/src/alignment.r index d28e028..c9cdacb 100644 --- a/src/alignment.r +++ b/src/alignment.r @@ -126,7 +126,7 @@ get_hvr_sequences <- function(sequences, vdj_segments, cores = detectCores()) { ) cys_coordinates <- parallel::mclapply(v_alignment, FUN = get_cys_coordinates) cys_df <- as.data.frame(do.call(rbind, cys_coordinates)) - remaining <- Biostrings::subseq(sequences, start = unlist(cys_df$end)) + remaining <- Biostrings::subseq(sequences, start = unlist(cys_df$end) + 1) j_alignment <- parallel::mcmapply(remaining, df$j_seq, FUN = align_sequence,