2021-02-26 02:59:22 +01:00
|
|
|
#!/bin/sh
|
|
|
|
|
|
|
|
usage() {
|
2021-02-28 02:23:58 +01:00
|
|
|
echo "usage: generation.sh <number of sequences> <sequencing runs>"
|
2021-02-26 02:59:22 +01:00
|
|
|
exit 1
|
|
|
|
}
|
|
|
|
|
2021-02-28 02:23:58 +01:00
|
|
|
if [ $# != 2 ]; then
|
2021-02-26 02:59:22 +01:00
|
|
|
usage
|
|
|
|
fi
|
|
|
|
|
|
|
|
sequences=$1
|
2021-02-28 02:23:58 +01:00
|
|
|
sequencing_runs=$2
|
2021-02-26 02:59:22 +01:00
|
|
|
data_directory="data/"
|
2021-02-28 02:23:58 +01:00
|
|
|
file="sequence.fastq"
|
2021-02-26 02:59:22 +01:00
|
|
|
prefix="curesim_"
|
|
|
|
|
2021-02-28 02:23:58 +01:00
|
|
|
Rscript src/repertoire.r "$sequences" "$sequencing_runs"
|
|
|
|
java -jar tools/CuReSim.jar -f "$data_directory$file" -o "$data_directory$prefix$file"
|