Mercurial > repos > petr-novak > dante
comparison dante.py @ 14:a6c55d1bdb6c draft
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) |