-
Notifications
You must be signed in to change notification settings - Fork 0
sarahhp/rnaseq-basic
Folders and files
Name | Name | Last commit message | Last commit date | |
---|---|---|---|---|
Repository files navigation
First time use: Make sure you have snakemake, hisat2, samtools, deeptools and subread installed, or are using a conda virtual environment with these installed. Running $ conda env create -f envs/rnaseq-basic.yml should create an environment with the right packages installed. My conda environment is called rnaseq-basic, so before running the workflow I run: $ conda activate rnaseq-basic Then I do a practise run or "dry" run: $ snakemake -s scripts/rnaseq_with_replicates.py --configfile test_mar19.json -npr $ snakemake -s scripts/rnaseq2.py --configfile test_mar19.json -npr scripts/rnaseq2.py runs the "proper" way to make hisat indices which leads to better alignment but takes up a lot of memory. For ensembl I needed to use laminB and 330GB of RAM for this step. scripts/rnaseq_with_replicates.py runs the quick and dirty way to make hisat indices. Full usage: $ snakemake -s scripts/rnaseq_with_replicates.py --configfile test_mar19.json -pr $ snakemake -s scripts/rnaseq2.py --configfile test_mar19.json -pr Correct index creation in scripts/rnaseq2.py. For ensembl I needed to use laminB and 330GB of RAM for this step. $ snakemake -s scripts/rnaseq2.py --configfile rnaseq_mar19.json --cores 8 -pr build_hisat_index > build_logs/02.log 2>&1 &
About
No description, website, or topics provided.
Resources
Stars
Watchers
Forks
Releases
No releases published
Packages 0
No packages published