elif overlap_noncoding!="": return overlap_noncoding else: return "-" else: return "-" (coding_ref_input,noncoding_ref_input,out_dir,combined_gtf_input)=fetch_args() start_time=time.time() print "Classification start:" (cuffcompare_outdir,gtf_outdir)=prepare_dirs(out_dir) (coding_ref,noncoding_ref,combined_gtf)=fetch_gtfs( coding_ref_input, noncoding_ref_input, combined_gtf_input, cuffcompare_outdir) gtf=Gtf.simple_read(combined_gtf_input) tid_gid=gtf.get_tid_gid() gid_tid=tid_gid.key_by([2],[1]) all_gids=gid_tid.getCol(1) print "Input genes: ", len(all_gids) print "Running Cuffcompare ..." strand_coding_tmap=cuffcompare(coding_ref_input,combined_gtf_input,cuffcompare_outdir+'/strandc') strand_coding_tmap.update_col('cuff_gene_id',gtf.get_gid(strand_coding_tmap.getCol('cuff_id'))) strand_noncoding_tmap=cuffcompare(noncoding_ref_input,combined_gtf_input,cuffcompare_outdir+'/strandnc') strand_noncoding_tmap.update_col('cuff_gene_id',gtf.get_gid(strand_noncoding_tmap.getCol('cuff_id'))) coding_tmap=cuffcompare(coding_ref,combined_gtf,cuffcompare_outdir+'/coding') coding_tmap.update_col('cuff_gene_id',gtf.get_gid(coding_tmap.getCol('cuff_id')))
return overlap_noncoding else: return "-" else: return "-" (coding_ref_input, noncoding_ref_input, out_dir, combined_gtf_input) = fetch_args() start_time = time.time() print "Classification start:" (cuffcompare_outdir, gtf_outdir) = prepare_dirs(out_dir) (coding_ref, noncoding_ref, combined_gtf) = fetch_gtfs(coding_ref_input, noncoding_ref_input, combined_gtf_input, cuffcompare_outdir) gtf = Gtf.simple_read(combined_gtf_input) tid_gid = gtf.get_tid_gid() gid_tid = tid_gid.key_by([2], [1]) all_gids = gid_tid.getCol(1) print "Input genes: ", len(all_gids) print "Running Cuffcompare ..." strand_coding_tmap = cuffcompare(coding_ref_input, combined_gtf_input, cuffcompare_outdir + '/strandc') strand_coding_tmap.update_col( 'cuff_gene_id', gtf.get_gid(strand_coding_tmap.getCol('cuff_id'))) strand_noncoding_tmap = cuffcompare(noncoding_ref_input, combined_gtf_input, cuffcompare_outdir + '/strandnc') strand_noncoding_tmap.update_col( 'cuff_gene_id', gtf.get_gid(strand_noncoding_tmap.getCol('cuff_id')))