def self_compare(fastas, id2desc, algorithm):
for fasta in fastas:<EOL><INDENT>blast = open(search(fasta, fasta, method = algorithm, alignment = '<STR_LIT>'))<EOL>for hit in best_blast(blast, <NUM_LIT:1>):<EOL><INDENT>id, bit = hit[<NUM_LIT:0>].split()[<NUM_LIT:0>], float(hit[-<NUM_LIT:1>])<EOL>id2desc[id].append(bit)<EOL><DEDENT><DEDENT>return id2desc<EOL>
compare genome to self to get the best possible bit score for each ORF