add_new_project_from_vcf_and_ped.py

#!/usr/bin/env python2.7
import sys
import os
import argparse
import yaml
from datetime import datetime, date


p = argparse.ArgumentParser("")
p.add_argument("-i", "--project-id", help="Project id", required=True)
p.add_argument("-n", "--project-name", help="Project name", required=True)
p.add_argument("-v", "--vcf", help="The vcf file", required=True)
p.add_argument("-p", "--ped", help="The ped file. If not specified, the sample ids from the VCF will be used instead of the ped file.")
p.add_argument("--no-phenotips", action="store_true", help="don't add to phenotips")
p.add_argument("-f", "--force", help="Force annotation", action="store_true")
p.add_argument("-r", dest="run", action="store_true", help="Actually run the commands")
opts = p.parse_args()

vcf = opts.vcf
ped = opts.ped
project_id = opts.project_id
project_name = opts.project_name

if not os.path.isfile(vcf):
    p.error("Invalid vcf: " + vcf)
if ped and not os.path.isfile(ped):
    p.error("Invalid ped: " + ped)


commands = [
    "python2.7 -u manage.py add_project %(project_id)s '%(project_name)s' ",
    "python2.7 -u manage.py add_custom_population_to_project %(project_id)s gnomad-exomes",
    "python2.7 -u manage.py add_custom_population_to_project %(project_id)s gnomad-genomes",
    "python2.7 -u manage.py add_individuals_to_project %(project_id)s " +  ("--ped %(ped)s" if ped else "--vcf %(vcf)s"), 
    "python2.7 -u manage.py add_vcf_to_project %(project_id)s %(vcf)s", 
]

if not opts.no_phenotips:
    commands += [
        "python2.7 -u manage.py add_project_to_phenotips %(project_id)s '%(project_name)s' ",
        "python2.7 -u manage.py add_individuals_to_phenotips %(project_id)s --ped %(ped)s ",
    ]

commands += [
    "python2.7 -u manage.py generate_pedigree_images %(project_id)s",
    "python2.7 -u manage.py add_default_tags %(project_id)s",
    "python2.7 -u manage.py load_project %(project_id)s" + (" --force-annotations --force-clean " if opts.force else ""), 
    "python2.7 -u manage.py load_project_datastore %(project_id)s",
]

commands = map(lambda s: s % globals(), commands )

print(date.strftime(datetime.now(), "%m/%d/%Y %H:%M:%S  -- Will run: "))
for c in commands:
    print(c)
#i = input("Do you want to continue? [y/n] ")
if opts.run:
    for c in commands:
        print(date.strftime(datetime.now(), "%m/%d/%Y %H:%M:%S") + " -- Running: " + c)
        sys.stdout.flush()
        r = os.system(c)
        if r != 0:
            print(date.strftime(datetime.now(), "%m/%d/%Y %H:%M:%S") + " -- Command failed: " + c + "\nExiting.." )
            break