Commit a5e8bfbb authored by nfontrod's avatar nfontrod
Browse files

src/bed_handler/get_other_exon_in_same_gene.py: modication of outpout file...

src/bed_handler/get_other_exon_in_same_gene.py: modication of outpout file name in get_other_exon_in_same_gene
parent e40ca2ed
Pipeline #143 passed with stages
in 51 seconds
......@@ -172,11 +172,12 @@ def get_other_exon_in_same_gene(bed: str,
final_df = get_final_bed(ctcf_dist, good_exon_df, distance_threshold)
final_df.to_csv(OutputBed.output / outfile, sep="\t", index=False)
gene_df = create_gene_bed4norm(BedConfig.gene_bed, final_df)
gene_df.to_csv(OutputBed.output / (outfile.replace(".bed", "") +
"-gene-dup.bed"), sep="\t", index=False)
gene_df.to_csv(OutputBed.output / outfile.replace("exon.bed",
"gene-dup.bed"),
sep="\t", index=False)
gene_df = gene_df.drop_duplicates()
gene_df.to_csv(OutputBed.output / (outfile.replace(".bed", "") +
"-gene.bed"), sep="\t", index=False)
gene_df.to_csv(OutputBed.output / outfile.replace("exon.bed", "gene.bed"),
sep="\t", index=False)
if __name__ == "__main__":
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment