Mercurial > repos > petr-novak > dante

comparison dante.py @ 14:a6c55d1bdb6c draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author petr-novak
date Wed, 28 Aug 2019 08:08:47 -0400
parents d0431a839606
children 3151a72a6671
comparison
equal deleted inserted replaced
13:ab56b34d67d8 14:a6c55d1bdb6c
583 db_name_best = db_names[best_idx_reg] 583 db_name_best = db_names[best_idx_reg]
584 db_starts_best = db_starts[best_idx_reg] 584 db_starts_best = db_starts[best_idx_reg]
585 db_ends_best = db_ends[best_idx_reg] 585 db_ends_best = db_ends[best_idx_reg]
586 if count_region == len(indices_plus): 586 if count_region == len(indices_plus):
587 strand_gff = "-" 587 strand_gff = "-"
588 if strand_gff == "+":
589 feature_start = min(start_hit[regions_above_threshold])-1
590 feature_end = max(end_hit[regions_above_threshold])
591 else:
592 feature_end = seq_len[region][0] - min(start_hit[regions_above_threshold]) - 1
593 feature_start = seq_len[region][0] - max(end_hit[regions_above_threshold]) + 1
588 create_gff3(domain_type, ann_substring, unique_annotations, 594 create_gff3(domain_type, ann_substring, unique_annotations,
589 ann_pos_counts, min(start_hit[regions_above_threshold])-1, 595 ann_pos_counts, feature_start,feature_end,
590 max(end_hit[regions_above_threshold]),
591 step, best_idx, annotation_best, db_name_best, 596 step, best_idx, annotation_best, db_name_best,
592 db_starts_best, db_ends_best, strand_gff, score, 597 db_starts_best, db_ends_best, strand_gff, score,
593 seq_id, db_seq, query_seq, domain_size, positions, gff, consensus) 598 seq_id, db_seq, query_seq, domain_size, positions, gff, consensus)
594 count_region += 1 599 count_region += 1
595 seq_ids.append(seq_id) 600 seq_ids.append(seq_id)