def find_coordinates(hmms, bit_thresh):
<EOL>seq2hmm = parse_hmm(hmms, bit_thresh)<EOL>seq2hmm = best_model(seq2hmm)<EOL>group2hmm = {} <EOL>for seq, info in list(seq2hmm.items()):<EOL><INDENT>group2hmm[seq] = {}<EOL>for group_num, group in enumerate(hit_groups(info[<NUM_LIT:1>])):<EOL><INDENT>best = sorted(group, reverse = True, key = itemgetter(-<NUM_LIT:1>))[<NUM_LIT:0>]<EOL>strand = best[<NUM_LIT:5>]<EOL>coordinates = [i[<NUM_LIT:0>] for i in group] + [i[<NUM_LIT:1>] for i in group]<EOL>coordinates = [min(coordinates), max(coordinates), strand]<EOL>matches = [i for i in group if i[<NUM_LIT:5>] == strand]<EOL>gaps = check_gaps(matches)<EOL>group2hmm[seq][group_num] = [info[<NUM_LIT:0>], strand, coordinates, matches, gaps]<EOL><DEDENT><DEDENT>return group2hmm<EOL>
find 16S rRNA gene sequence coordinates