We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
1 parent 7c2bc45 commit 0cd5e98Copy full SHA for 0cd5e98
1 file changed
strkit/call/call_locus.py
@@ -939,7 +939,10 @@ def get_locus_alignment_data_from_read(
939
)
940
# TODO: find a way of making use of these reads (e.g., C9orf72 expansion)
941
# if locus_read_coords.full_left_flank:
942
- # print(segment.query_sequence[locus_read_coords.left_flank_end:])
+ # print("left flank only", segment.query_sequence[locus_read_coords.left_flank_end:])
943
+ # print(repr(locus_read_coords))
944
+ # print(locus_with_ref_data.right_coord_adj, locus_with_ref_data.ref_right_flank_seq)
945
+ # print(segment.aligned_coords.pair_at_idx(len(segment.aligned_coords) - 1))
946
return None
947
948
if not (left := locus_read_coords.full_left_flank) or not locus_read_coords.full_right_flank:
0 commit comments