File: SeqAsm.py

package info (click to toggle)
pbsuite 15.8.24%2Bdfsg-3
  • links: PTS, VCS
  • area: main
  • in suites: buster
  • size: 14,512 kB
  • sloc: python: 10,962; sh: 147; xml: 21; makefile: 14
file content (27 lines) | stat: -rw-r--r-- 516 bytes parent folder | download | duplicates (4)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
import sys


from difflib import SequenceMatcher

USAGE= "%prog A small assembler that quickly guesses at consensus of small, completely overlapping sequences"""

#ss = ["CGCGTACGCG",
#     ["
a = "dsa jld lal"
b = "dsajld kll"
c = "dsc jle kal"
d = "dsd jlekal"

ss = [a,b,c,d]

s = SequenceMatcher()

for i in range(len(ss)):
    x = ss[i]
    s.set_seq1(x)
    for j in range(i+1,len(ss)):    
        y = ss[j]
        s.set_seq2(y)   
        print
        print s.ratio()
        print s.get_matching_blocks()