coolneng
/
locigenesis
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

Rename output file to curesim-HVR.fastq

This commit is contained in:
coolneng 2021-05-14 20:01:56 +02:00
parent 91b3e37bd8
commit 2acec89f84
Signed by: coolneng
GPG Key ID: 9893DA236405AF57
2 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
README.md
View File

@ -63,6 +63,6 @@ command invokes it:
The script will generate 2 files under the data directory: The script will generate 2 files under the data directory:
|HVR.fastq | CuReSim-HVR.fastq | |HVR.fastq | curesim-HVR.fastq |
|:----:|:-----:| |:----:|:-----:|
|Contains the original CDR3 sequence|Contains CDR3 after the read simulation, with sequencing errors | |Contains the original CDR3 sequence|Contains CDR3 after the read simulation, with sequencing errors |

2
src/alignment.r
View File

@ -145,4 +145,4 @@ get_hvr_sequences <- function(sequences, vdj_segments, cores = detectCores()) {
data <- parse_data(file = "data/curesim_sequence.fastq") data <- parse_data(file = "data/curesim_sequence.fastq")
hvr <- get_hvr_sequences(sequences = data[[1]], vdj_segments = data[[2]]) hvr <- get_hvr_sequences(sequences = data[[1]], vdj_segments = data[[2]])
Biostrings::writeXStringSet(hvr, "data/CuReSim-HVR.fastq", format = "fastq") Biostrings::writeXStringSet(hvr, "data/curesim-HVR.fastq", format = "fastq")