Print column numbers that are 100% conserved

[CCGG26]

how can i write a Python program that prints the column numbers in a FASTA format multiple alignment that are 100% conserved. im really having trouble getting a grip of this concept...For example in the multiple alignment below

>human
ACC
>mouse
ACC
>cat
TCC
>dog
ACA

column 2 is 100% conserved but columns 1 and 3 are not 100% conserved.

[bvdet]

What is the output supposed to look like? I don't know what a Fasta format multiple alignment is. :(

[CCGG26]

the fasta format multiple alignment is just a name i believe, i do not think it is anything significant..

[bvdet]

What is the output supposed to look like?

[CCGG26]

output: column 2 = 100%

[bvdet]

What is column 2? 100% of what? Please be specific.

[CCGG26]

A 100% conserved column is
one that has the exact same nucleotide in every sequence. For example if the
user enters 1 and the the multiple alignment below is given as input

>human
ACC
>mouse
ACC
>cat
TCC
>dog
ACA

then the output should be "No". But if the user enters 2 then the output
should be "Yes".

[bvdet]

Following is a test script that show's how it can be done using set().

import random

data = '''>human
ACC
>mouse
ACC
>cat
TCC
>dog
ACA'''

def conserved(col, seq):
    colList = set([item[col] for item in seq])
    if len(colList) == 1:
        return True
    return False

dataList = data.split("\n")
sequences =[list(dataList[i]) for i in range(1, len(dataList), 2)]

column = random.choice([0,1,2])

result = conserved(column, sequences)
print "Column %s %s conserved" % (column, ["IS", "IS NOT"][not result or 0])

The three columns are 0, 1 and 2 which is consistent with a list index.

[dwblas]

Shouldn't it be true for 1 or 2, and false for 3 or 4?