comparison extract_aln_ends.py @ 50:0b9aab6aaebf draft

Uploaded 16cfcafe8b42055c5dd64e62c42b82b455027a40
author rnateam
date Tue, 26 Jan 2016 04:38:27 -0500
parents 570a7de9f151
children
comparison
equal deleted inserted replaced
49:303f6402a035 50:0b9aab6aaebf
1 #!/usr/bin/env python 1 #!/usr/bin/env python
2
3 import argparse
4 import logging
5 from sys import stdout
6 from shutil import rmtree
7 from tempfile import mkdtemp
8 from pybedtools import BedTool
9 import pysam
10 # avoid ugly python IOError when stdout output is piped into another program
11 # and then truncated (such as piping to head)
12 from signal import signal, SIGPIPE, SIG_DFL
13 signal(SIGPIPE, SIG_DFL)
2 14
3 tool_description = """ 15 tool_description = """
4 Extract alignment ends from sam file. 16 Extract alignment ends from sam file.
5 17
6 The resulting bed file contains the outer coordinates of the alignments. The 18 The resulting bed file contains the outer coordinates of the alignments. The
30 License: Apache 42 License: Apache
31 Email: maticzkd@informatik.uni-freiburg.de 43 Email: maticzkd@informatik.uni-freiburg.de
32 Status: Development 44 Status: Development
33 """ 45 """
34 46
35 import argparse
36 import logging
37 from sys import stdout
38 from shutil import rmtree
39 from tempfile import mkdtemp
40 from pybedtools import BedTool
41 import pysam
42
43 47
44 class DefaultsRawDescriptionHelpFormatter(argparse.ArgumentDefaultsHelpFormatter, 48 class DefaultsRawDescriptionHelpFormatter(argparse.ArgumentDefaultsHelpFormatter,
45 argparse.RawDescriptionHelpFormatter): 49 argparse.RawDescriptionHelpFormatter):
46 # To join the behaviour of RawDescriptionHelpFormatter with that of ArgumentDefaultsHelpFormatter 50 # To join the behaviour of RawDescriptionHelpFormatter with that of ArgumentDefaultsHelpFormatter
47 pass 51 pass
48
49 # avoid ugly python IOError when stdout output is piped into another program
50 # and then truncated (such as piping to head)
51 from signal import signal, SIGPIPE, SIG_DFL
52 signal(SIGPIPE, SIG_DFL)
53 52
54 # parse command line arguments 53 # parse command line arguments
55 parser = argparse.ArgumentParser(description=tool_description, 54 parser = argparse.ArgumentParser(description=tool_description,
56 epilog=epilog, 55 epilog=epilog,
57 formatter_class=DefaultsRawDescriptionHelpFormatter) 56 formatter_class=DefaultsRawDescriptionHelpFormatter)