Confrontando sequenze di DNA in Python 3

Question

Il mio programma finora è:

seq_a = "TGGAGGCAATGGCGGCCAGCA" 
seq_b = "GACTCCTCCTCCTCCTGCTCA"  
len_a = len(seq_a)  
len_b = len(seq_b)  
print("Length of Sequence A: " + str(len_a))  
print()  
print("Length of Sequence B: " + str(len_b)) 
print() 

def sequence_compare(seq_a, seq_b): 
     len1 = len(seq_a) 
     len2 = len(seq_b) 
     mismatches = [] 
     for pos in range (0, min(len1, len2)) : 
       if seq_a[pos] != seq_b[pos]: 
        mismatches.append('|') 
       else: 
        mismatches.append(' ') 
     print (seq_a) 
     print (mismatches) 
     print (seq_b) 
sequence_compare(seq_a,seq_b) 

mi aspetto un'uscita lungo le linee di:

TGGAGGCAATGGCGGCCAGCA 
||||||||| ||||||||| 
GACTCCTCCTCCTCCTGCTCA 

Ma invece l'output è:

TGGAGGCAATGGCGGCCAGCA 
['|', '|', '|', '|', '|', '|', '|', '|', '|', ' ', '|', '|', '|', '|', '|', '|', '|', '|', '|', ' ', ' '] 
GACTCCTCCTCCTCCTGCTCA 

Answer 1

Sei abbastanza vicino ..

print ("".join(mismatches)) 

farà il trucco e stampare come richiesto che è

||||||||| |||||||||