Skip to content

Commit

Permalink
added hound_isd_bed
Browse files Browse the repository at this point in the history
  • Loading branch information
Anya Korsakova committed Sep 25, 2024
1 parent 233043e commit 24c5212
Showing 1 changed file with 3 additions and 5 deletions.
8 changes: 3 additions & 5 deletions src/baskerville/scripts/hound_isd_bed.py
Original file line number Diff line number Diff line change
Expand Up @@ -68,7 +68,7 @@ def main():
parser.add_option(
"-o",
dest="out_dir",
default="sat_mut",
default="sat_del",
help="Output directory [Default: %default]",
)
parser.add_option(
Expand Down Expand Up @@ -269,7 +269,7 @@ def main():
# for each nucleotide
# copy and modify
alt_1hot = np.copy(ref_1hot)

# left-matched shift: delete 1 nucleotide at position mi
dna.hot1_delete(alt_1hot[0], mi, options.del_len)

Expand Down Expand Up @@ -300,9 +300,7 @@ def main():
ref_preds_stitch, alt_preds, options.snp_stats, None
)
for snp_stat in options.snp_stats:
scores_h5[snp_stat][si, mi - mut_start, 0] = ism_scores[
snp_stat
]
scores_h5[snp_stat][si, mi - mut_start, 0] = ism_scores[snp_stat]

# close output HDF5
scores_h5.close()
Expand Down

0 comments on commit 24c5212

Please sign in to comment.