Warning

This document is for an in-development version 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()