Skip to content

Commit

Permalink
Merge pull request #1091 from griffithlab/generate_protein_fasta_prox…
Browse files Browse the repository at this point in the history
…_vars

Generate protein fasta prox vars
  • Loading branch information
susannasiebert authored Jun 5, 2024
2 parents 7bc69d1 + f8ce28d commit ca00971
Showing 1 changed file with 3 additions and 0 deletions.
3 changes: 3 additions & 0 deletions pvactools/tools/pvacseq/generate_protein_fasta.py
Original file line number Diff line number Diff line change
Expand Up @@ -193,6 +193,9 @@ def main(args_input = sys.argv[1:]):

temp_dir = tempfile.mkdtemp()
proximal_variants_tsv = convert_vcf(args.input_vcf, temp_dir, args.sample_name, args.phased_proximal_variants_vcf, args.flanking_sequence_length, args.pass_only)
proximal_variants_file = "{}.proximal_variants.tsv".format(args.output_file)
if proximal_variants_tsv is not None:
shutil.copy(proximal_variants_tsv, proximal_variants_file)
generate_fasta(args.flanking_sequence_length, downstream_sequence_length, temp_dir, proximal_variants_tsv)
parse_files(args.output_file, temp_dir, args.mutant_only, args.input_tsv, args.aggregate_report_evaluation)
shutil.rmtree(temp_dir, ignore_errors=True)
Expand Down

0 comments on commit ca00971

Please sign in to comment.