Compare commits
2 Commits
762791829a
...
bc666a37c7
Author | SHA1 | Date | |
---|---|---|---|
bc666a37c7 | |||
d5bf8f3b72 |
@ -21,3 +21,4 @@ prefix="curesim_"
|
|||||||
|
|
||||||
Rscript src/repertoire.r "$sequences" "$sequencing_runs"
|
Rscript src/repertoire.r "$sequences" "$sequencing_runs"
|
||||||
java -jar tools/CuReSim.jar -m "$read_mean_size" -sd "$read_variance_size" -f "$data_directory$filename$fasta" -o "$data_directory$prefix$filename$fastq"
|
java -jar tools/CuReSim.jar -m "$read_mean_size" -sd "$read_variance_size" -f "$data_directory$filename$fasta" -o "$data_directory$prefix$filename$fastq"
|
||||||
|
rm "$data_directory/log.txt"
|
||||||
|
@ -14,9 +14,9 @@ amplify_rows <- function(data, column, factor) {
|
|||||||
if (column == "sequence") {
|
if (column == "sequence") {
|
||||||
dna_string <- Biostrings::DNAStringSet(data)
|
dna_string <- Biostrings::DNAStringSet(data)
|
||||||
reverse_complement <- Biostrings::reverseComplement(dna_string)
|
reverse_complement <- Biostrings::reverseComplement(dna_string)
|
||||||
return(Biostrings::DNAStringSet(rep(reverse_complement, factor)))
|
return(rep(reverse_complement, factor))
|
||||||
}
|
}
|
||||||
return(as.character(rep(data, factor)))
|
return(rep(data, factor))
|
||||||
}
|
}
|
||||||
|
|
||||||
save_data <- function(data, name) {
|
save_data <- function(data, name) {
|
||||||
|
Loading…
Reference in New Issue
Block a user