comparison scripts/S01_find_orf_on_multiple_alignment.py @ 10:3d00be2d05f3 draft

planemo upload for repository https://github.com/abims-sbr/adaptsearch commit 7d45b8e07bf99420df9b195617515b8002b5eaf3

9:640ef4c06ed5

###### Reference: http://crazyhottommy.blogspot.fr/2013/10/python-code-for-getting-reverse.html
############################################################
def ReverseComplement2(seq):
    # too lazy to construct the dictionary manually, use a dict comprehension
    seq1 = 'ATCGN-TAGCN-atcgn-tagcn-'
    seq_dict = { seq1[i]:seq1[i+6] for i in range(24) if i < 6 or 12<=i<16 }
    return "".join([seq_dict[base] for base in reversed(seq)])
############################


#######################

10:3d00be2d05f3

###### Reference: http://crazyhottommy.blogspot.fr/2013/10/python-code-for-getting-reverse.html
############################################################
def ReverseComplement2(seq):
    # too lazy to construct the dictionary manually, use a dict comprehension
    seq1 = 'ATCGN-TAGCN-atcgn-tagcn-'
    seq_dict = { seq1[i]:seq1[i+6] for i in range(24) if i < 6 or 12<=i<=16 }
    return "".join([seq_dict[base] for base in reversed(seq)])
############################


#######################