Skip to content

Commit

Permalink
copy pileup vcf to merge_output vcf if no full-alignment candidate
Browse files Browse the repository at this point in the history
  • Loading branch information
zhengzhenxian committed May 22, 2023
1 parent 42eaaee commit e8947ff
Showing 1 changed file with 5 additions and 5 deletions.
10 changes: 5 additions & 5 deletions preprocess/SelectCandidates.py
Original file line number Diff line number Diff line change
Expand Up @@ -273,11 +273,11 @@ def SelectCandidates(args):
['\t'.join([contig_name, str(x[0] - 1), str(x[1] - 1), ]) for x in
split_output]) + '\n') # bed format

if len(all_full_aln_regions) > 0:
all_full_aln_regions_path = os.path.join(split_folder, 'FULL_ALN_FILE_{}'.format(contig_name))
with open(all_full_aln_regions_path, 'w') as output_file:
output_file.write('\n'.join(all_full_aln_regions) + '\n')
else:

all_full_aln_regions_path = os.path.join(split_folder, 'FULL_ALN_FILE_{}'.format(contig_name))
with open(all_full_aln_regions_path, 'w') as output_file:
output_file.write('\n'.join(all_full_aln_regions) + '\n')
if len(all_full_aln_regions) == 0:
print_calling_step(output_fn=pileup_vcf_fn)
return

Expand Down

0 comments on commit e8947ff

Please sign in to comment.