Skip to content
Snippets Groups Projects
Commit 397e4149 authored by nina.marthe_ird.fr's avatar nina.marthe_ird.fr
Browse files

corrected the output of the gff for several genomes

parent 8fd17172
No related branches found
No related tags found
No related merge requests found
...@@ -27,7 +27,6 @@ def target_gff(first_seg,last_seg,feature_id,list_seg,max_diff): ...@@ -27,7 +27,6 @@ def target_gff(first_seg,last_seg,feature_id,list_seg,max_diff):
def transfer_on_target(pos_seg, gfa, out_once, out_var,out_clustal,target_genomes,max_diff): def transfer_on_target(pos_seg, gfa, out_once, out_var,out_clustal,target_genomes,max_diff):
[once,var,clustal,stats]=[True,True,True,False] [once,var,clustal,stats]=[True,True,True,False]
get_segments_positions_on_genome(pos_seg)
list_feature_to_transfer= Features.keys() list_feature_to_transfer= Features.keys()
if var: if var:
...@@ -35,11 +34,13 @@ def transfer_on_target(pos_seg, gfa, out_once, out_var,out_clustal,target_genome ...@@ -35,11 +34,13 @@ def transfer_on_target(pos_seg, gfa, out_once, out_var,out_clustal,target_genome
for genome in target_genomes: for genome in target_genomes:
print(f'generation of {genome} output') print(f'generation of {genome} output')
get_segments_positions_on_genome(pos_seg)
segments_list={} segments_list={}
if var: if var:
target_genome_path=paths[genome] target_genome_path=paths[genome]
if once: if once:
print(f'generation of {genome} gff')
file_out=open(out_once,'w') file_out=open(out_once,'w')
global output_target_gff global output_target_gff
output_target_gff=[0,"",file_out] output_target_gff=[0,"",file_out]
...@@ -73,6 +74,7 @@ def transfer_on_target(pos_seg, gfa, out_once, out_var,out_clustal,target_genome ...@@ -73,6 +74,7 @@ def transfer_on_target(pos_seg, gfa, out_once, out_var,out_clustal,target_genome
file_out.close() file_out.close()
if var: if var:
print(f'generation of {genome} gene variation details')
seg_seq=get_segments_sequence(gfa,segments_list) seg_seq=get_segments_sequence(gfa,segments_list)
file_out_var = open(out_var, 'w') file_out_var = open(out_var, 'w')
global output_variations global output_variations
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment