@files( "pipeline.ini", PARAMS["orthology_pairwise"] )
def getPairwiseOrthologs( infile, outfile ):
    statement = '''psql -h db -U andreas -d postgres -F "," -A -c "select s.set_id, m.schema, m.gene_id, s.pattern 
                   from cgat_proj007_v2.ortholog_sets_members m, cgat_proj007_v2.ortholog_sets s 
                   where m.set_id=s.set_id and s.nspecies=2"
                   | sed s/,/\\\\t/g
                   > %(outfile)s '''
    P.run()
    
###################################################################
Export list of pairwise orthologous genes from all species comparisons from postgres database
dataset/ETHPy150Open CGATOxford/cgat/obsolete/pipeline_species_conservation.py/getPairwiseOrthologs