def _annotate(args, setclus):
logger.info("<STR_LIT>")<EOL>bedfile = generate_position_bed(setclus)<EOL>a = pybedtools.BedTool(bedfile, from_string=True)<EOL>beds = []<EOL>logger.info("<STR_LIT>")<EOL>if hasattr(args, '<STR_LIT>'):<EOL><INDENT>beds = args.list_files.split("<STR_LIT:U+002C>")<EOL>for filebed in beds:<EOL><INDENT>logger.info("<STR_LIT>" % filebed)<EOL>db = os.path.basename(filebed)<EOL>b = pybedtools.BedTool(filebed)<EOL>c = a.intersect(b, wo=True)<EOL>setclus = anncluster(c, setclus, db, args.type_ann, args.feature_id)<EOL><DEDENT><DEDENT>return setclus<EOL>
annotate transcriptional units with
    gtf/bed files provided by -b/g option