Commit c0e6014a authored by Jorge Navarro Muñoz's avatar Jorge Navarro Muñoz
Browse files

Don't stop if aligned sequences' lengths have a mismatch

parent 681d99a9
......@@ -371,12 +371,10 @@ def cluster_distance(A, B, A_domlist, B_domlist, anchor_domains):
# Sequences *should* be of the same length unless something went
# wrong elsewhere
if len(aligned_seqA) != len(aligned_seqB):
print("\tWARNING: mismatch in sequences' lengths while calculating sequence identity")
print("\t Domain: " + shared_domain)
print("\t Specific domain 1: " + aligned_seqA + " len: " + str(len(seq1)))
print("\t Specific domain 2: " + aligned_seqB + " len: " + str(len(seq2)))
sys.exit("debug stop")
print("\t trying to continue with shortest length...")
print("\tWARNING: mismatch in sequences' lengths while calculating sequence identity (" + shared_domain + ")")
if verbose:
print("\t Specific domain 1: " + aligned_seqA + " len: " + str(len(seq1)))
print("\t Specific domain 2: " + aligned_seqB + " len: " + str(len(seq2)))
seq_length = min(len(aligned_seqA), len(aligned_seqB))
else:
seq_length = len(aligned_seqA)
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment