这学期选了生信的选修课—perl/python在生物信息学中的应用
把结课作业的代码整理出来主要是python对FASTA文件的读取和数据处理
FASTA文件数据处理
FASTA文件读取:
只含一个基因序列
将FASTA文件的基因序列读取到一个列表中,列表中的每个元素为每一行基因序列构成的字符串
f=open('/home/miaoyr/perl_practice/test1_file/DTNBP1.fasta')
ls=[]
for line in f:
if not line.startswith('>'):
ls.append(line.replace('\n','&#