Mercurial > repos > tiagoantao > clone_filter
annotate STACKS_clone_filter.py @ 7:2059c028c07a draft
planemo upload commit 4d00a0b7a4b859b642eceb3c7cd4adb1453903a0-dirty
author | tiagoantao |
---|---|
date | Thu, 04 Feb 2016 08:00:32 -0500 |
parents | 1a5ebb16aaea |
children | a8b87bc0c3ad |
rev | line source |
---|---|
0
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
1 #!/usr/bin/python |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
2 # -*- coding: utf-8 -*- |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
3 |
5
1a5ebb16aaea
planemo upload commit 4d00a0b7a4b859b642eceb3c7cd4adb1453903a0-dirty
tiagoantao
parents:
0
diff
changeset
|
4 import argparse |
1a5ebb16aaea
planemo upload commit 4d00a0b7a4b859b642eceb3c7cd4adb1453903a0-dirty
tiagoantao
parents:
0
diff
changeset
|
5 import glob |
1a5ebb16aaea
planemo upload commit 4d00a0b7a4b859b642eceb3c7cd4adb1453903a0-dirty
tiagoantao
parents:
0
diff
changeset
|
6 import os |
0
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
7 import subprocess |
5
1a5ebb16aaea
planemo upload commit 4d00a0b7a4b859b642eceb3c7cd4adb1453903a0-dirty
tiagoantao
parents:
0
diff
changeset
|
8 import shutil |
0
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
9 |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
10 |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
11 def __main__(): |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
12 parser = argparse.ArgumentParser() |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
13 parser.add_argument('-f', type=str) |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
14 parser.add_argument('-fwd', type=str) |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
15 parser.add_argument('-rev', type=str) |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
16 |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
17 #parser.add_argument('-k', action='store_true') |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
18 |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
19 options = parser.parse_args() |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
20 |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
21 # create the working dir |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
22 os.mkdir('job_outputs') |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
23 os.chdir('job_outputs') |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
24 |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
25 # create the populations command input line |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
26 cmd_line=['clone_filter'] |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
27 if options.f is not None: |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
28 cmd_line.extend(['-f', options.f]) |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
29 else: |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
30 cmd_line.extend(['-1', options.fwd]) |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
31 cmd_line.extend(['-2', options.rev]) |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
32 |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
33 print "[CMD]:"+' '.join(cmd_line) |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
34 subprocess.call(cmd_line) |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
35 |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
36 # postprocesses |
5
1a5ebb16aaea
planemo upload commit 4d00a0b7a4b859b642eceb3c7cd4adb1453903a0-dirty
tiagoantao
parents:
0
diff
changeset
|
37 clean_files = glob.glob('*.fq*') |
0
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
38 if options.f is not None: |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
39 shutil.copy(clean_files[0], 'clean.fq') |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
40 else: |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
41 shutil.copy(clean_files[0], 'fwd.fq') |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
42 shutil.copy(clean_files[1], 'rev.fq') |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
43 |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
44 |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
45 if __name__ == '__main__': |
42e64fbf4609
planemo upload commit 60a090d1dc89c2eb383feeb35cd6078b9130e0e5-dirty
tiagoantao
parents:
diff
changeset
|
46 __main__() |