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 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
|
||||
|
||||
usage() {
|
||||
echo "usage: generation.sh <number of sequences> <number_of_reads>"
|
||||
echo "usage: generation.sh <number of sequences> <number of reads>"
|
||||
exit 1
|
||||
}
|
||||
|
||||
@@ -17,6 +17,6 @@ filename="sequence"
|
||||
prefix="curesim_"
|
||||
|
||||
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
|
||||
rm "$data_directory/log.txt"
|
||||
|
||||
28
shell.nix
28
shell.nix
@@ -2,15 +2,35 @@
|
||||
|
||||
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 = [
|
||||
R
|
||||
rPackages.immuneSIM
|
||||
rPackages.Biostrings
|
||||
rPackages.stringr
|
||||
jdk
|
||||
# Development tools
|
||||
rPackages.languageserver
|
||||
rPackages.lintr
|
||||
CuReSim
|
||||
];
|
||||
}
|
||||
|
||||
Reference in New Issue
Block a user