Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- from Bio import SeqIO
- import sys
- ffile = SeqIO.parse(sys.argv[1], "fasta")
- header_set = {line.strip() for line in open(sys.argv[2])}
- 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