GC content

计算核酸序列的GC含量（GC content）#!/usr/bin/env python# -*- coding: utf-8 -*-__author__ = 'shengwei ma'__author_email__ = 'shengweima@icloud.com'from Bio import SeqIOfrom Bio.SeqUtils import GCfor rec in S

msw521sg

1465人浏览 · 2017-03-03 15:36:27

msw521sg · 2017-03-03 15:36:27 发布

计算核酸序列的GC含量（GC content）

#!/usr/bin/env python
# -*- coding: utf-8 -*- 
__author__ = 'shengwei ma'
__author_email__ = 'shengweima@icloud.com'

from Bio import SeqIO
from Bio.SeqUtils import GC

for rec in SeqIO.parse("1.txt", "fasta"): #1.txt为fasta格式的文件
    print ' %s\t%2.2f%%' % (rec.name, GC(rec.seq))