Warning

This document is for an old release of Galaxy. You can alternatively view this page in the latest release if it exists or view the top of the latest release's documentation.

Source code for galaxy.datatypes.converters.wiggle_to_array_tree_converter

#!/usr/bin/env python

from __future__ import division

import sys

from bx.arrays.array_tree import array_tree_dict_from_reader, FileArrayTreeDict
from bx.arrays.wiggle import WiggleReader

BLOCK_SIZE = 100


def main():

    input_fname = sys.argv[1]
    out_fname = sys.argv[2]

    reader = WiggleReader(open(input_fname))

    # Fill array from reader
    d = array_tree_dict_from_reader(reader, {}, block_size=BLOCK_SIZE)

    for array_tree in d.values():
        array_tree.root.build_summary()

    FileArrayTreeDict.dict_to_file(d, open(out_fname, "w"))


if __name__ == "__main__":
    main()