Advertisement
Guest User

Untitled

a guest
Mar 26th, 2018
309
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Python 0.22 KB | None | 0 0
  1. from Bio import SeqIO
  2. import sys
  3.  
  4. ffile = SeqIO.parse(sys.argv[1], "fasta")
  5. header_set = {line.strip() for line in open(sys.argv[2])}
  6. SeqIO.write((r for r in ffile if r.name not in header_set), 'out_file.fa', 'fasta')
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement