msisensor-pro on Biowulf

MSIsensor-pro is an updated version of msisensor. MSIsensor-pro evaluates Microsatellite Instability (MSI) for cancer patients with next generation sequencing data. It accepts the whole genome sequencing, whole exome sequencing and target region (panel) sequencing data as input. MSIsensor-pro introduces a multinomial distribution model to quantify polymerase slippages for each tumor sample and a discriminative sites selection method to enable MSI detection without matched normal samples. For samples of various sequencing depths and tumor purities, MSIsensor-pro significantly outperformed the current leading methods in terms of both accuracy and computational cost. If you want to know more detail about MSIsensor-pro, please see the MSIsensor-pro Schematics and Internals MSIsensor-pro page.

References:

Documentation
Important Notes

Interactive job
Interactive jobs should be used for debugging, graphics, or applications that cannot be run as batch jobs.

Allocate an interactive session and run the program. Sample session:

[user@biowulf]$ sinteractive --mem=5g
salloc.exe: Pending job allocation 46116226
salloc.exe: job 46116226 queued and waiting for resources
salloc.exe: job 46116226 has been allocated resources
salloc.exe: Granted job allocation 46116226
salloc.exe: Waiting for resource configuration
salloc.exe: Nodes cn3144 are ready for job

module load msisensor-pro
[user@cn3144]$ cp -rp ${MSISENSORPRO_TEST_DATA} /data/$USER
[user@cn3144]$ cd /data/$USER/demo/scripts
[user@cn3144]$ ./1_test_scan.sh
[user@cn3144]$ ./2_test_msi.sh
[user@cn3144]$ ./3_test_baseline.sh
[user@cn3144]$ ./4_test_pro.sh
    

Batch job
Most jobs should be run as batch jobs.

Create a batch script file (e.g. msisensor-pro.sh). For example:

#!/bin/bash
module load msisensor-pro
cp -rp ${MSISENSORPRO_TEST_DATA} /data/$USER
cd /data/$USER/demo/scripts
./1_test_scan.sh
./2_test_msi.sh
./3_test_baseline.sh
./4_test_pro.sh

Submit this job using the Slurm sbatch command.

sbatch --mem=5g msisensor-pro.sh