Add v_call and j_call to sequence ID

src/repertoire.r

@@ -10,13 +10,19 @@ generate_repertoire <- function(number_of_sequences) {
   ))
 }
 
+save_data <- function(data, reads) {
+  vdj_sequences <- data[-1]
+  amplified_vdj <- vdj_sequences[rep(seq_len(nrow(vdj_sequences)), reads), ]
+  Biostrings::writeXStringSet(data$sequence, "data/sequence.fasta")
+}
+
 process_data <- function(repertoire, reads) {
   columns <- c("sequence", "v_call", "j_call")
   data <- repertoire[, columns]
   dna_sequence <- Biostrings::DNAStringSet(data$sequence)
   data$sequence <- Biostrings::reverseComplement(dna_sequence)
   names(data$sequence) <- paste(rownames(data), data$v_call, data$j_call)
-  Biostrings::writeXStringSet(data$sequence, "data/sequence.fasta")
+  save_data(data, reads)
 }
 
 parse_cli_arguments <- function() {