Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- import csv
- import sys
- import random
- import re
- # Ensure correct argv line
- if len(sys.argv)!= 3:
- print("Error !")
- sys.exit(1)
- with open(sys.argv[1]) as file:
- reader = csv.DictReader(file)
- sequences = next(reader)
- with open (sys.argv[2]) as dnafile:
- dna = csv.reader(dnafile)
- for row in dnafile:
- dnalist = row
- DNA = dnalist[0]
- dict_seq = {}
- search_re = re.compile('(?:' + STR + ')+')
- print(max((len(seq) for seq in search_re.findall(dna)), default=0) // len(STR))
- with open(sys.argv[1]) as file:
- reader = csv.DictReader(file)
- for row in reader:
- match = 0
- for DNA in dict_seq:
- if dict_seq[DNA] == int(row[DNA]):
- match += 1
- if match == len(dict_seq):
- print(row['name'])
- exit()
- #otherwise, no match
- print("No match")
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement