Как увеличить счетчик при чтении строк скриптом?

Question

У меня есть эта таблица значений, и мне было интересно, как я могу позволить программе читать каждую строку. Для каждой строки с 'a', 'g', 'c' или 'u' я хочу увеличить количество на единицу. Для этого примера, когда я запускаю его, он должен иметь результат 12.

a  1    0.000 S
g  2    0.260 S
a  3    0.990 S
a  4    0.980 S
c  5    0.000 S
u  6    1.000 S
c  7    0.000 S
a  8    1.000 S
a  9    1.000 T
u 10    0.820 S
a 11    1.000 T
g 12    0.000 S
F 13    1.000 S
S 14    1.000 S
T 15    1.000 S

Код, который я попробовал ниже:

rna_residues = ['a','c','g','u']
count_dict = {}
        #Making the starting number 0
        rna_count = 0
        #if any lines of the file starts with one of the rna_residue
        if line.startswith(tuple(rna_residues)):
            for residue in line:
                if residue in rna_residues:
                    rna_count += 1
            count_dict[line] = [rna_count]  
            print count_dict    

Каким-то образом, когда я запускаю его, нет списка счетчиков:

{'a  1    0.000 S\n': [1]}
{'g  2    0.260 S\n': [1]}
{'a  3    0.990 S\n': [1]}
{'a  4    0.980 S\n': [1]}
{'c  5    0.000 S\n': [1]}
{'u  6    1.000 S\n': [1]}
{'c  7    0.000 S\n': [1]}
{'a  8    1.000 S\n': [1]}
{'a  9    1.000 T\n': [1]}
{'u 10    0.820 S\n': [1]}
{'a 11    1.000 T\n': [1]}
{'g 12    0.000 S\n': [1]}

Я знаю, что это много информации, но есть ли советы, которые могут мне помочь с этим? Большое спасибо!!

Answer 1

Вы используете целую строку в качестве ключа в словаре, поэтому, если у вас нет одинаковых строк, все значения будут равны 1. Зачем вам вообще нужен словарь? У меня сложилось впечатление, что вы хотите посчитать количество строк, начинающихся с любого из символов 'a','c','g','u',

Для этого достаточно следующего кода:

rna_residues = ['a','c','g','u']
rna_count = 0
with open('/path/to/file') as opened_file:    
    for line in opened_file:
        # or if line[0] in rna_residues
        if any(line.startswith(residue) for residue in rna_residues):
            rna_count += 1
print rna_count
# 12