Пример #1
0
#!/usr/bin/python
'''
!!!!NOT FINISHED, DO NOT USE!!!!

for similary functions, use fastaLen + fasta2tab + updateQuery + cut

1. fasta file name
2. asc/des
'''

from sys import argv,exit
from Bio import SeqIO

try:
    fname = argv[1]
    order = argv[2]
except:
    exit(__doc__)


f = open(fname,'r')
records = SeqIO.parst(f,'fasta')

for r in records:
    seqName = r.description
    sequence = r.seq
    

f.close()