view mira.py @ 0:ddd2e3362c5e draft default tip

Uploaded v0.0.1, previously only on the Test Tool Shed
author peterjc
date Fri, 21 Nov 2014 06:38:46 -0500
parents
children
line wrap: on
line source

"""
MiraAssemblyFormat class for the 'mira' format within Galaxy
"""

from galaxy.datatypes.data import Text


class MiraAssemblyFormat(Text):
    """MIRA Assembly Format  data"""
    file_ext = "mira"

    def sniff( self, filename ):
        """Determines whether the file is a MIRA Assembly Format file.

        Note currently this only detects MIRA Assembly Format v2.0,
        as used in MIRA v3.9 and v4.0.

        It does not detect MIRA Assembly Format v1 as used in both
        MIRA v3.2 and v3.4.
        """
        h = open(filename)
        line = h.readline()
        if line.rstrip() != "@Version\t2\t0":
            h.close()
            return False
        line = h.readline()
        if line.rstrip() != "@Program\tMIRALIB":
            h.close()
            return False
        return True
    
    def merge(split_files, output_file):
        """Merging multiple MIRA files is non-trivial and may not be possible..."""
        if len(split_files) == 1:
            #For one file only, use base class method (move/copy)
            return Text.merge(split_files, output_file)
        if not split_files:
            raise ValueError("Given no MIRA, %r, to merge into %s" \
                             % (split_files, output_file))
        raise NotImplementedError("Merging MIRA Assembly Files has not been implemented")
    merge = staticmethod(merge)

    def split( cls, input_datasets, subdir_generator_function, split_params):
        """Split a MIRA Assembly File (not implemented for now)."""
        if split_params is None:
            return None
        raise NotImplementedError("Can't yet split a MIRA Assembly Format file")
    merge = staticmethod(merge)