Hprc banner tamu.png


Revision as of 17:39, 9 January 2019 by Cmdickens (talk | contribs) (POTION)
Jump to: navigation, search


GCATemplates available: no

POTION homepage


A user has identified that the columns for the final_results.out file do not match the columns in the positive.out file.

The POTION developer has been notified and the developer has acknowledged that the columns in the final_results.out file are out of order.

For Example:
   In the file final_results.out there may be 3 P's (in column 2) for model m78
   and 76 P's for m12 (in column 3). However, in the positive.out file it is reversed, there
   are 76 sequences listed for model 8 (m78) and only 3 for model 2 (m12)

POTION (POsitive selecTION) is an open source, modular and end-to-end software for genomic scale detection of positive Darwinian selection in groups of homologous coding sequences through estimation of dN/dS ratios.

module spider POTION

POTION requires multiple steps to be run on the command line prior to running a job script:

module load POTION/1.1.2-intel-2015B-Perl-5.20.0

rsync -r $EBROOTPOTION/bin ./
rsync -r $EBROOTPOTION/config_files ./

ln -s $EBROOTPAML/bin/codeml bin/
ln -s $EBROOTPHYLIP/bin/consense bin/
ln -s $EBROOTMUSCLE/bin/muscle bin/
ln -s $EBROOTPHIPACK/bin/Phi bin/
ln -s $EBROOTPRANK/bin/prank bin/
ln -s $EBROOTPHYLIP/bin/proml bin/
ln -s $EBROOTPHYLIP/bin/dnaml bin/
ln -s $EBROOTPHYLIP/bin/seqboot bin/
ln -s $EBROOTTRIMAL/bin/trimal bin/
ln -s $EBROOTMAFFT/bin/mafft bin/
ln -s $EBROOTPAGAN/bin/pagan bin/
ln -s $EBROOTPHYML/bin/phyml bin/

cp ./config_files/dummy_potion_user_config_file_site_standard my_config_file

Edit the following variables in your my_config_file:

CDS_dir_path = /path/to/CDS                         # path to folder containing CDS data; you create the CDS directory and put your data there
homology_file_path = /path/to/all_orthomcl.out      # path to the ORTHOMCL 1.4 main output file
project_dir_path = /path/to/potion_results_out      # path to the main directory where results will be created.Parent directory must exist.

Then edit only the potion_dir variable in the ./config_files/potion_config file and set it to the full path of your current working directory

potion_dir = $EBROOTPOTION/           # change this to the full path of your working directory

OrthoMCL version 1.4 is required if generating the all_orthomcl.out input file for POTION.

Sample job script run from your working directory:

#BSUB -L /bin/bash
#BSUB -J potion_job
#BSUB -o stdout.%J
#BSUB -e stderr.%J
#BSUB -n 20
#BSUB -R "span[ptile=20]"
#BSUB -R "rusage[mem=2700]"
#BSUB -M 2700
#BSUB -W 96:00

module load POTION/1.1.2-intel-2015B-Perl-5.20.0

./bin/potion.pl --conf_file_path /absolute/pato/to/working/directory/my_config_file