计算碱基出现次数DNA

计算碱基出现次数

编写一个程序，从文本文件gene.dat中读取DNA序列，计算碱基出现频率

哈哈第一次编文案
这思路很简单：就是比较字母是否相等

import os
# 读取及打开文件
f=open(r路径)
h,a,t,c,g=0,0,0,0,0
str1=''
# `读取文件每行`
for line in f:
    s = line.split()
    #每行分割
    h+=1
    # 不要第一行
    if h !=1:
        for i in s:
            str1=str(i)
            #对每行里每个字母分割
            for x in list(str1):
            # 字母转换为大写
                X=x.upper()
            # 判别字母
                if X=='A':
                    a+=1
                elif X=='T':
                    t+=1
                elif X=='C':
                    c+=1
                elif X=='G':
                    g+=1
print("count：",a,t,c,g)
f.close()

这是一种简单的方式
还有一种引用库
使用biopython