Skip to content

Commit

Permalink
chore(call): rm outdated comment
Browse files Browse the repository at this point in the history
  • Loading branch information
davidlougheed committed Dec 4, 2024
1 parent ae541b3 commit ec3e5b7
Showing 1 changed file with 1 addition and 4 deletions.
5 changes: 1 addition & 4 deletions strkit/call/output/vcf.py
Original file line number Diff line number Diff line change
Expand Up @@ -146,7 +146,7 @@ def output_contig_vcf_lines(
peak_start_anchor_seqs_upper = tuple(iter_to_upper(peak_start_anchor_seqs))
common_anchor_prefix = commonprefix([ref_start_anchor, *peak_start_anchor_seqs_upper])
# anchor_offset = how many bases we can cut off from the front of the anchor
# since they're shared between all alleles.
# since they're shared between all alleles - yields a more compact representation.
# - we need to leave one base as an anchor for VCF compliance though, thus the min(...)
anchor_offset = min(len(common_anchor_prefix), params.vcf_anchor_size - 1)

Expand All @@ -167,9 +167,6 @@ def output_contig_vcf_lines(
key=lambda c: c[1] + c[0]
)

# If there's no anchor variation, we will use the last base for a compact representation.
# Otherwise, we include the full VCF anchor sized-sequence. TODO: common prefix strip

call = result["call"]
call_95_cis = result["call_95_cis"]

Expand Down

0 comments on commit ec3e5b7

Please sign in to comment.