import Functions as fct # creates segments and features for the intersect between the graph for chr10 and the gff of IRGSP intersect_path='/home/nina/annotpangenome/align_genes/intersect_segments-genes_chr10.bed' fct.create_seg_feat(intersect_path) # outputs the gff of the graph for chr10 output_file='new_graph_chr10.gff' fct.graph_gff(output_file) pos_seg="seg_coord/AzucenaRS1_chromosome10.bed" gff="new_graph_chr10.gff" out="azucena_chr10.gff" fct.azu(pos_seg,gff,out)