Skip to content

Instantly share code, notes, and snippets.

@wenliangz

wenliangz/rosalind_hamm

Created February 15, 2015 23:56
Show Gist options

  • Select an option

    Learn more about clone URLs
  • Save wenliangz/6a7391d8d09ea79f75ae to your computer and use it in GitHub Desktop.

Select an option

Learn more about clone URLs
Save wenliangz/6a7391d8d09ea79f75ae to your computer and use it in GitHub Desktop.
Download ZIP

Revisions

  1. wenliangz created this gist Feb 15, 2015.
    31 changes: 31 additions & 0 deletions rosalind_hamm
    View file
    Original file line number Diff line number Diff line change
    @@ -0,0 +1,31 @@
    seq=[]
    with open('rosalind_hamm.txt') as f:
    for each in f:
    each=each.strip()
    seq.append(each)

    seq1=seq[0]
    seq2=seq[1]

    seq1_dict={}
    seq2_dict={}
    for position,base in enumerate(seq1):
    seq1_dict[position]=base

    for position,base in enumerate(seq2):
    seq2_dict[position]=base

    Hamm=set(seq1_dict.iteritems())-set(seq2_dict.iteritems())
    dH=len(Hamm)
    print dH

    seq1 ="ABCD"
    seq2 = "ABXD"
    dist = 0
    for i in range(0, min(len(seq1), len(seq2))):
    if seq1[i] != seq2[i]:
    dist += 1
    dist += abs(len(seq1) - len(seq2))

    #what the hell is going on