生物信息脚本练习（1） 找出fasta文件中大于500的序列

最近做了一些生物信息的脚本练习。
这是第一个例子。
找出一个fasta文件中大于500的序列，并重定向到另一个新的文件中。
这个文件每条序列是如下的样子。

c100027.graph_c0|orf3 type=complete len=150nt loc=c100027.graph_c0:123-272:-
ATGAGGATCTTTACGCCAAATGAGGGCCTTGTTGTTGATCTTTTCAGTAAGAGACGTTGC
GGAAATATTGGCGAAAATTTAAGAAACCTAGTTAGTTTTAGCAGCCACATAGTTGGCAAG
AATTATACTATTCAAGCCATGTGGCACTGA

这是我一开始的解法：
思路是用正则匹配第一行中len的数值。找到他们之后，根据这数字算出这条序列有多少行，然后把这么些行的信息输出。

import re
output = open('data.txt', 'w')
seq = []
element = []
row = []
with open ("d:/Zanthoxylum_Bungeanum_Maxim.Unigene.cds.fa","r") as f:
    for l in f:
        seq.append(l)
for i in seq:
    if re.match(">.*len=[5-9][0-9][0-9]nt.*",i) or re.match(">.*len=[0-9][0-9][0-9][0-9]nt.*",i):
        element_number =