-
Notifications
You must be signed in to change notification settings - Fork 114
Expand file tree
/
Copy pathCRISPResso2_router.py
More file actions
35 lines (31 loc) · 1.46 KB
/
CRISPResso2_router.py
File metadata and controls
35 lines (31 loc) · 1.46 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
# /usr/bin/env python
# -*- coding: utf-8 -*-
"""CRISPResso2 - Kendell Clement and Luca Pinello 2018
Software pipeline for the analysis of genome editing outcomes from deep sequencing data
(c) 2018 The General Hospital Corporation. All Rights Reserved.
"""
import subprocess as sb
import sys
usage = '\n- CRISPResso2 Docker Container -\n\t- Possible commands: \n\t CRISPResso\n\t CRISPRessoBatch\n\t CRISPRessoPooled\n\t CRISPRessoWGS\n\t CRISPRessoCompare\n\t CRISPRessoPooledWGSCompare\n\t CRISPRessoAggregate\n\t License\n\t- this docker version should be run like this:\n\tdocker run -v ${PWD}:/DATA -w /DATA -i pinellolab/crispresso2 CRISPResso -r1 fastq1.fq -a AAAATTT \n'
if len(sys.argv) == 1:
print(usage)
sys.exit(1)
if sys.argv[1] == 'CRISPResso':
sb.call(["CRISPResso"] + sys.argv[2:])
elif sys.argv[1] == 'CRISPRessoBatch':
sb.call(["CRISPRessoBatch"] + sys.argv[2:])
elif sys.argv[1] == 'CRISPRessoCompare':
sb.call(["CRISPRessoCompare"] + sys.argv[2:])
elif sys.argv[1] == 'CRISPRessoPooled':
sb.call(["CRISPRessoPooled"] + sys.argv[2:])
elif sys.argv[1] == 'CRISPRessoWGS':
sb.call(["CRISPRessoWGS"] + sys.argv[2:])
elif sys.argv[1] == 'CRISPRessoPooledWGSCompare':
sb.call(["CRISPRessoPooledWGSCompare"] + sys.argv[2:])
elif sys.argv[1] == 'CRISPRessoAggregate':
sb.call(["CRISPRessoAggregate"] + sys.argv[2:])
elif sys.argv[1] == 'License':
with open("LICENSE.txt", 'r') as fin:
print(fin.read())
else:
print(usage)