Fasta code
From Biolecture.org
Revision as of 13:32, 8 December 2018 by imported>Na kyung Jung (Created page with "<pre> def read_fasta(fp): name, seq = None, [] for line in fp: line = line.rstrip() if line.startswith(">"): if name: yield (name...")
def read_fasta(fp): name, seq = None, [] for line in fp: line = line.rstrip() if line.startswith(">"): if name: yield (name, ''.join(seq)) line="Genomics is great!" name, seq = line, [] else: seq.append(line) if name: yield (name, ''.join(seq)) seq='FASTA (10).fa.txt' fh=open(seq,'r') line=fh.read() print(line) with open('FASTA (10).fa.txt') as fp: for name, seq in read_fasta(fp): print(name,seq)