Compare commits
5 Commits
36eb73b458
...
0.1.0
| Author | SHA1 | Date | |
|---|---|---|---|
|
ad8abcc4fc
|
|||
|
6440816a87
|
|||
|
0e005735bc
|
|||
|
4f0936718b
|
|||
|
1b6e2d13ea
|
46
README.org
46
README.org
@@ -1,3 +1,49 @@
|
|||||||
* locigenesis
|
* locigenesis
|
||||||
|
|
||||||
locigenesis is a tool that generates an immune repertoire and runs it through a sequence reader simulation tool, to generate sequencing errors.
|
locigenesis is a tool that generates an immune repertoire and runs it through a sequence reader simulation tool, to generate sequencing errors.
|
||||||
|
|
||||||
|
** Installation
|
||||||
|
|
||||||
|
This project uses [[https://nixos.org/][Nix]] to ensure reproducible builds.
|
||||||
|
|
||||||
|
1. Install Nix (compatible with MacOS, Linux and [[https://docs.microsoft.com/en-us/windows/wsl/about][WSL]]):
|
||||||
|
|
||||||
|
#+begin_src shell
|
||||||
|
curl -L https://nixos.org/nix/install | sh
|
||||||
|
#+end_src
|
||||||
|
|
||||||
|
1. Clone the repository:
|
||||||
|
|
||||||
|
#+begin_src shell
|
||||||
|
git clone https://git.coolneng.duckdns.org/coolneng/locigenesis
|
||||||
|
#+end_src
|
||||||
|
|
||||||
|
3. Change the working directory to the project:
|
||||||
|
|
||||||
|
#+begin_src shell
|
||||||
|
cd locigenesis
|
||||||
|
#+end_src
|
||||||
|
|
||||||
|
4. Enter the nix-shell:
|
||||||
|
|
||||||
|
#+begin_src shell
|
||||||
|
nix-shell
|
||||||
|
#+end_src
|
||||||
|
|
||||||
|
After running these commands, you will find yourself in a shell that contains all the needed dependencies.
|
||||||
|
|
||||||
|
** Usage
|
||||||
|
|
||||||
|
An execution script that accepts 2 parameters is provided, the following command invokes it:
|
||||||
|
|
||||||
|
#+begin_src shell
|
||||||
|
./generation.sh <number of sequences> <number of reads>
|
||||||
|
#+end_src
|
||||||
|
|
||||||
|
- <number of sequences>: an integer that specifies the number of different sequences to generate
|
||||||
|
- <number of reads>: an integer that specifies the number of reads to perform on each sequence
|
||||||
|
|
||||||
|
The script will generate 2 files under the data directory:
|
||||||
|
|
||||||
|
| HVR.fastq | Contains the original CDR3 sequence |
|
||||||
|
| CuReSim-HVR.fastq | Contains CDR3 after the read simulation, with sequencing errors |
|
||||||
|
|||||||
@@ -1,7 +1,7 @@
|
|||||||
#!/bin/sh
|
#!/bin/sh
|
||||||
|
|
||||||
usage() {
|
usage() {
|
||||||
echo "usage: generation.sh <number of sequences> <number_of_reads>"
|
echo "usage: generation.sh <number of sequences> <number of reads>"
|
||||||
exit 1
|
exit 1
|
||||||
}
|
}
|
||||||
|
|
||||||
@@ -17,6 +17,6 @@ filename="sequence"
|
|||||||
prefix="curesim_"
|
prefix="curesim_"
|
||||||
|
|
||||||
Rscript src/repertoire.r "$sequences" "$number_of_reads" &&
|
Rscript src/repertoire.r "$sequences" "$number_of_reads" &&
|
||||||
java -jar tools/CuReSim.jar -f "$data_directory$filename$fastq" -o "$data_directory$prefix$filename$fastq"
|
CuReSim -f "$data_directory$filename$fastq" -o "$data_directory$prefix$filename$fastq"
|
||||||
Rscript src/alignment.r
|
Rscript src/alignment.r
|
||||||
rm "$data_directory/log.txt"
|
rm "$data_directory/log.txt"
|
||||||
|
|||||||
28
shell.nix
28
shell.nix
@@ -2,15 +2,35 @@
|
|||||||
|
|
||||||
with pkgs;
|
with pkgs;
|
||||||
|
|
||||||
mkShell {
|
let
|
||||||
|
CuReSim = stdenv.mkDerivation rec {
|
||||||
|
name = "CuReSim";
|
||||||
|
version = "1.3";
|
||||||
|
src = fetchzip {
|
||||||
|
url =
|
||||||
|
"http://www.pegase-biosciences.com/wp-content/uploads/2015/08/${name}${version}.zip";
|
||||||
|
sha256 = "1hvlpgy4haqgqq52mkxhcl9i1fx67kgwi6f1mijvqzk0xff77hkp";
|
||||||
|
stripRoot = true;
|
||||||
|
extraPostFetch = ''
|
||||||
|
chmod go-w $out
|
||||||
|
'';
|
||||||
|
};
|
||||||
|
|
||||||
|
nativeBuildInputs = [ makeWrapper ];
|
||||||
|
|
||||||
|
installPhase = ''
|
||||||
|
mkdir -pv $out/share/java $out/bin
|
||||||
|
cp -r ${src} $out/share/java/${name}
|
||||||
|
makeWrapper ${pkgs.jdk}/bin/java $out/bin/CuReSim --add-flags "-jar $out/share/java/${name}/${name}.jar"
|
||||||
|
'';
|
||||||
|
};
|
||||||
|
in mkShell {
|
||||||
buildInputs = [
|
buildInputs = [
|
||||||
R
|
R
|
||||||
rPackages.immuneSIM
|
rPackages.immuneSIM
|
||||||
rPackages.Biostrings
|
rPackages.Biostrings
|
||||||
rPackages.stringr
|
rPackages.stringr
|
||||||
jdk
|
jdk
|
||||||
# Development tools
|
CuReSim
|
||||||
rPackages.languageserver
|
|
||||||
rPackages.lintr
|
|
||||||
];
|
];
|
||||||
}
|
}
|
||||||
|
|||||||
Reference in New Issue
Block a user