Mercurial > repos > petr-novak > dante
comparison dante.py @ 15:3151a72a6671 draft
Uploaded
| author | petr-novak |
|---|---|
| date | Tue, 03 Sep 2019 05:20:02 -0400 |
| parents | a6c55d1bdb6c |
| children | 039c45c01b47 |
comparison
equal
deleted
inserted
replaced
| 14:a6c55d1bdb6c | 15:3151a72a6671 |
|---|---|
| 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 == "+": | 588 if strand_gff == "+": |
| 589 feature_start = min(start_hit[regions_above_threshold])-1 | 589 feature_start = min(start_hit[regions_above_threshold]) + 1 |
| 590 feature_end = max(end_hit[regions_above_threshold]) | 590 feature_end = max(end_hit[regions_above_threshold]) |
| 591 else: | 591 else: |
| 592 feature_end = seq_len[region][0] - min(start_hit[regions_above_threshold]) - 1 | 592 feature_end = seq_len[region][0] - min(start_hit[regions_above_threshold]) |
| 593 feature_start = seq_len[region][0] - max(end_hit[regions_above_threshold]) + 1 | 593 feature_start = seq_len[region][0] - max(end_hit[regions_above_threshold]) + 1 |
| 594 create_gff3(domain_type, ann_substring, unique_annotations, | 594 create_gff3(domain_type, ann_substring, unique_annotations, |
| 595 ann_pos_counts, feature_start,feature_end, | 595 ann_pos_counts, feature_start,feature_end, |
| 596 step, best_idx, annotation_best, db_name_best, | 596 step, best_idx, annotation_best, db_name_best, |
| 597 db_starts_best, db_ends_best, strand_gff, score, | 597 db_starts_best, db_ends_best, strand_gff, score, |