diff --git a/Functions_output.py b/Functions_output.py index ca2f93c6ad868c31337eb61e6af8cb4cc1e1d9ec..d8af551852073318a827ad63fdbf826ada0485b3 100644 --- a/Functions_output.py +++ b/Functions_output.py @@ -51,7 +51,7 @@ def gff_one(first_seg,last_seg,feature_id,list_seg,max_diff): # writes the gff of azucena using the gff of the graph -def genome_gff(pos_seg, gff, gfa, out_once, out_detail, out_var,target_genome_name,max_diff): +def genome_gff(pos_seg, gfa, out_once, out_detail, out_var,target_genome_name,max_diff): print("generation of the genome's gff ") #Â create variables and open files @@ -83,7 +83,7 @@ def genome_gff(pos_seg, gff, gfa, out_once, out_detail, out_var,target_genome_na # no segment is missing, the feature is complete - feature_ok # total number of features, with missing segments or not - feature_total get_segments_positions_on_genome(pos_seg) - list_feature_to_transfer=get_all_features_in_gff(gff) # get the list of all the features to transfer from the gff + list_feature_to_transfer= Features.keys() for feat in list_feature_to_transfer: diff --git a/main.py b/main.py index dfa58ef6ef61d7924a38db8117ca23c5186ba74e..9f7eddad9ba19bc7792c1541adb50a3ec3544b3d 100755 --- a/main.py +++ b/main.py @@ -43,7 +43,7 @@ if run=="command_line": # outputs the gff of a genome for the chr10 max_diff=2 #Â maximum size difference (n times bigger or smaller) between the gene on the reference genome and the gene on the target genome for the gene to be transfered. - genome_gff(pos_seg,out_gff,gfa,out_once,out_detail,out_var,target_genome_name,max_diff) + genome_gff(pos_seg,gfa,out_once,out_detail,out_var,target_genome_name,max_diff)