Code:
#!/usr/bin/python from Bio import SeqIO import sys import os #usage: python long.seq.py in.fasta out.fasta 200 input_seq_iterator = SeqIO.parse(open(sys.argv[1], "rU"), "fasta") short_seq_iterator = (record for record in input_seq_iterator \ if len(record.seq) > int(sys.argv[3])) output_handle = open(sys.argv[2], "w") SeqIO.write(short_seq_iterator, output_handle, "fasta") output_handle.close()
Leave a comment: