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 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43
| import os
import sys
k=15
mon_fichier1=open(sys.argv[1],"r")
for line1 in mon_fichier1:
mon_fichier2=open(sys.argv[2],"r")
if line1[0]=='>':
seq1id=line1[:-2]
else:
seq1=line1
for line2 in mon_fichier2:
if line2[0]=='>':
seq2id=line2[:-2]
else:
seq2=line2
print ("seq1id :",seq1id)
print ("seq2id :",seq2id)
if seq1[0]!='>' and seq2[0]!='>':
L1=[]
L2=[]
for i in range(0, len(seq1)-k+1,k):
item=seq1[i:i+k]
if item not in L1:
L1.append(item)
for i in range(0, len(seq2)-k+1):
item=seq2[i:i+k]
if item not in L2:
L2.append(item)
C=0
for kmer1 in L1:
for kmer2 in L2:
if kmer1==kmer2:
C=C+1
print ("C= ",C)
if C>5:
print("il y a une ressemblence!!") #je vais ajouter d'autres traitement
mon_fichier1.close()
mon_fichier2.close() |
Partager