The post-TOGA processing pipeline.
- Adds
--outdir
to control where postoga output goes.- Adds
--isoforms
to allow the user specify external isoform tables.- Disables git features from
logger.py
[branch/commit] info.- Changes a small portion of
filter_query_annotation.py
&&make_query_table.py
to fix the--paralog
flag and allow a broader inclusion of isoforms.
To use postoga, just:
Clone the repository
# clone the repository
git clone --recursive https://github.com/alejandrogzi/postoga.git
cd postoga
Create a conda environment (recommended) OR set INSTALL_PYTHON in configure.sh to true
conda env create -f env.yml
conda activate postoga
Run configure.sh to install dependencies (Rust required)
./configure.sh
Run test to confirm functionality
./test.sh
If you see something like this at then end, postoga is ready!:
####################################
postoga: the post-TOGA processing pipeline
version: 0.5.0-devel
commit: 67ecb4e
branch: master
[XXXXXX] - INFO: postoga started!
[XXXXXX] - INFO: running in mode base with arguments: {'mode': 'base', 'path': './supply/test/', 'by_class': None, 'by_rel': None, 'threshold': '0.5', 'to': 'gff', 'assembly_qual': './supply/Ancestral_placental.txt'}
[XXXXXX] - INFO: found 10 projections, 10 unique transcripts, 10 unique genes
Note
postoga is dependent from TOGA. Any changes in TOGA will have a repercusion here. If you found any bug/errors, please report them here. This project is in constant development, any desired features are welcome!
Here is a descrption of postoga features:
usage: postoga.py [-h] {base,haplotype}
positional arguments:
{base,haplotype} Select mode
base Base mode
haplotype Haplotype mode
postoga.py base [-h] --outdir OUTDIR --togadir TOGADIR [-bc BY_CLASS] [-br BY_REL]
[-th THRESHOLD] -to {gtf,gff,bed} [-aq ASSEMBLY_QUAL]
[-sp {human,mouse,chicken}] [-src {ensembl,gene_name,entrez}]
[-phy {mammals,birds}] [-s] [-par PARALOG]
optional arguments:
-h, --help show this help message and exit
--outdir OUTDIR, -o OUTDIR
Path to posTOGA output directory
--togadir TOGADIR, --td TOGADIR
Path to TOGA results directory
-bc BY_CLASS, --by-class BY_CLASS
Filter parameter to only include certain orthology classes (I, PI, UL, M, PM,
L, UL)
-br BY_REL, --by-rel BY_REL
Filter parameter to only include certain orthology relationships (o2o, o2m,
m2m, m2m, o2z)
-th THRESHOLD, --threshold THRESHOLD
Filter parameter to preserve orthology scores greater or equal to a given
threshold (0.0 - 1.0)
-to {gtf,gff,bed}, --to {gtf,gff,bed}
Specify the conversion format for .bed (query_annotation/filtered) file (gtf,
gff3) or just keep it as .bed (bed)
-aq ASSEMBLY_QUAL, --assembly_qual ASSEMBLY_QUAL
Calculate assembly quality based on a list of genes provided by the user
(default: Ancestral_placental.txt)
-sp {human,mouse,chicken}, --species {human,mouse,chicken}
Species name to be used as a reference for the assembly quality calculation
(default: human)
-src {ensembl,gene_name,entrez}, --source {ensembl,gene_name,entrez}
Source of the ancestral gene names (default: ENSG)
-phy {mammals,birds}, --phylo {mammals,birds}
Phylogenetic group of your species (default: mammals)
-s, --skip Skip steps 2, 3, and 4 and only filter the .bed file
-par PARALOG, --paralog PARALOG
Filter parameter to preserve transcripts with paralog projection probabilities
less or equal to a given threshold (0.0 - 1.0)
-iso ISOFORMS, --isoforms ISOFORMS
Path to a custom isoform table (default: None)
postoga.py haplotype [-h] --outdir OUTDIR -hp HAPLOTYPE_DIR [-r RULE] [-s {query,loss}]
optional arguments:
-h, --help show this help message and exit
--outdir OUTDIR, -o OUTDIR
Path to posTOGA output directory
-hp HAPLOTYPE_DIR, --haplotype_dir HAPLOTYPE_DIR
Path to TOGA results directories separated by commas (path1,path2,path3)
-r RULE, --rule RULE Rule to merge haplotype assemblies (default: I>PI>UL>L>M>PM>PG>abs)
-s {query,loss}, --source {query,loss}
Source of the haplotype classes (query, loss)