rnasa

Gene Expression Level Calculator for RNA-seq

Installation

$ pip install -U rnasa

Dependent commands:

• pigz
• pbzip2
• bgzip
• samtools (and plot-bamstats)
• java
• fastqc
• trim_galore
• STAR
• rsem-prepare-reference
• rsem-refseq-extract-primary-assembly
• rsem-calculate-expression

Docker image

Pull the image from Docker Hub.

$ docker image pull dceoy/rnasa

Usage

1. Download and process resource data.

   $ rnasa download --genome=GRCh38 --dest-dir=/path/to/ref
   

2. Calculate TPM (transcripts per million) values.

   $ rnasa calculate \
       --workers=2 \
       --dest-dir=/path/to/output \
       /path/to/ref/GRCh38 \
       /path/to/sample1_fastq_prefix \
       /path/to/sample2_fastq_prefix \
       /path/to/sample3_fastq_prefix
   

   The command search for one (single-end) or two (paired-end) input FASTQ files by prefix.

3. Extract TPM values from RSEM results files.

   $ rnasa extract --dest-dir=. /path/to/output/rsem
   

Run rnasa --help for more information.