[Py014] 剔除含N序列

剔除所有含有N碱基的序列

from Bio import SeqIO

no_N_records = [record for record in SeqIO.parse('test.fa','fasta') if 'N' not in record.seq]

SeqIO.write(no_N_records, 'out.fa','fasta')

剔除全是N碱基的序列

from Bio import SeqIO

not_only_N_records = [record for record in SeqIO.parse('test.fa','fasta') if ''.join(set(record.seq))!='N']

SeqIO.write(not_only_N_records, 'out.fa','fasta')
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

相关阅读更多精彩内容

友情链接更多精彩内容