Remove redundant sequencing runs argument

generation.sh

@@ -19,6 +19,5 @@ fastq=".fastq"
 filename="sequence"
 prefix="curesim_"
 
-Rscript src/repertoire.r "$sequences" "$sequencing_runs"
+Rscript src/repertoire.r "$sequences" && java -jar tools/CuReSim.jar -n "$sequencing_runs" -m "$read_mean_size" -sd "$read_variance_size" -f "$data_directory$filename$fasta" -o "$data_directory$prefix$filename$fastq"
-java -jar tools/CuReSim.jar -n "$sequencing_runs" -m "$read_mean_size" -sd "$read_variance_size" -f "$data_directory$filename$fasta" -o "$data_directory$prefix$filename$fastq"
 rm "$data_directory/log.txt"

src/repertoire.r

@@ -22,7 +22,8 @@ process_data <- function(repertoire) {
     "d_sequence_alignment", "j_sequence_alignment"
   )
   data <- repertoire[, columns]
-  data$sequence <- Biostrings::reverseComplement(data$sequence)
+  dna_sequence <- Biostrings::DNAStringSet(data$sequence)
+  data$sequence <- Biostrings::reverseComplement(dna_sequence)
   save_data(data)
 }
 
@@ -34,6 +35,6 @@ parse_cli_arguments <- function() {
   return(args[1])
 }
 
-arguments <- parse_cli_arguments(commandArgs(trailing))
+args <- parse_cli_arguments()
-repertoire <- generate_repertoire(number_of_sequences = arguments[1])
+repertoire <- generate_repertoire(number_of_sequences = as.integer(args[1]))
 process_data(repertoire)