visualize.py 文件源码

def __check_input( opts, args, parser ):
    """
    Make sure the input is in the form of either a cmp.h5 file of aligned reads
    or a FOFN of unaligned bas.h5 files. Also make sure that a reference fasta 
    file is specified if 
    """
    if len(args)!=2:
        print "ERROR -- expecting two arguments: \
                 (1) <SEQ>_methyl_features.txt output from methylprofiles containing methylation features for mapping \
                 (2) <SEQ>_other_features.txt output from methylprofiles containing alternative sequence features for mapping"

    mfeats_fn  = args[0]
    ofeats_fn  = args[1]
    feature_type = None

    if not os.path.exists(mfeats_fn):
        parser.error("Can't find file of sequence features (methylprofiles output) for mapping: %s" % mfeats_fn)

    if not os.path.exists(ofeats_fn):
        parser.error("Can't find file of sequence features (methylprofiles output) for mapping: %s" % ofeats_fn)

    return mfeats_fn, ofeats_fn