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
|
#!/usr/bin/python3
"""
Read a wiggle track and print out a series of lines containing
"chrom position score". Ignores track lines, handles bed, variableStep
and fixedStep wiggle lines.
"""
import sys
import bx.wiggle
if len(sys.argv) > 1:
in_file = open(sys.argv[1])
else:
in_file = sys.stdin
if len(sys.argv) > 2:
out_file = open(sys.argv[2], "w")
else:
out_file = sys.stdout
for fields in bx.wiggle.Reader(in_file):
print(" ".join(map(str, fields)))
in_file.close()
out_file.close()
|