#maybe best for nucleic acid sequences not protein sequences
            #alignments = pairwise2.align.globalxx (target_sequence, query_sequence)
            #print "aligned:\n"
            #print format_alignment(*alignments[0])
            