[已解决]关于python选择性阅读文本的问题

其实我不帅

目的：给出一个fasta格式的文件，只要输入文件名就可以统计数该文件中碱基序列的A、T、C、G数量
我写的代码

1. def ATCG_counter(seq):
   
2.      a=t=c=g=0
   
3.      for each in seq:
   
4.           if each == 'A':
   
5.                a +=1
   
6.           elif each == 'T':
   
7.                t +=1
   
8.           elif each == 'G':
   
9.                g +=1
   
10.           elif each == 'C':
    
11.                c +=1
    
12. 
    
13.      print('A:%d'%(a))
    
14.      print('T:%d'%(t))
    
15.      print('C:%d'%(c))
    
16.      print('G:%d'%(g))
    
17.      print('All base:%d'%(a+t+c+g))
    
18.      
    
19. import os
    
20. file_name = input('请输入文件名：')
    
21. f = open(file_name)
    
22. sequence = f.read()
    
23. ATCG_counter(sequence)

复制代码

问题：fasta格式文件的第一行是这段序列的基本信息，之后才是序列碱基，但是使用这段代码会使第一行的代码信息中的字母被误以为是碱基，造成统计误差，该如何修改这段代码？
bty 求大神给看看这段代码有何改进之处

最佳答案

月排行榜 / 总排行榜

冬雪雪冬

2016-10-2 17:57:35

好在每行最后有换行符，可以把文件读取作为列表，再用列表切片，去掉第一行。

1. f = open('t.txt')
   
2. seq = list(f)
   
3. seq = seq[1:]
   
4. print(seq)
   
5. f.close()

复制代码

跳转到最佳答案楼层

冬雪雪冬

好在每行最后有换行符，可以把文件读取作为列表，再用列表切片，去掉第一行。

1. f = open('t.txt')
   
2. seq = list(f)
   
3. seq = seq[1:]
   
4. print(seq)
   
5. f.close()

复制代码

其实我不帅

冬雪雪冬 发表于 2016-10-2 17:57
好在每行最后有换行符，可以把文件读取作为列表，再用列表切片，去掉第一行。

好的 我去试试  谢谢