Skip to content

Commit 833c80b

Browse files
committed
Annotate END2 for SR_ONLY_LARGE_INVERSION
1 parent 0218a28 commit 833c80b

1 file changed

Lines changed: 8 additions & 0 deletions

File tree

src/svtk/svtk/cxsv/complex_sv.py

Lines changed: 8 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -67,11 +67,15 @@ def resolve(self, SR_only_cutoff):
6767
self.set_unresolved()
6868
self.vcf_record.info['SVTYPE'] = 'BND'
6969
self.vcf_record.info['UNRESOLVED_TYPE'] = 'SR_ONLY_LARGE_INVERSION'
70+
if 'END2' not in self.vcf_record.info:
71+
self.vcf_record.info['END2'] = self.vcf_record.stop
7072
self.cluster_type = 'SR_ONLY_LARGE_INVERSION'
7173
for r in self.records:
7274
r.info['UNRESOLVED'] = True
7375
r.info['UNRESOLVED_TYPE'] = 'SR_ONLY_LARGE_INVERSION'
7476
r.info['SVTYPE'] = 'BND'
77+
if 'END2' not in r.info:
78+
r.info['END2'] = r.stop
7579
if 'CPX_TYPE' in r.info.keys():
7680
r.info.pop('CPX_TYPE')
7781
if 'CPX_INTERVALS' in r.info.keys():
@@ -120,11 +124,15 @@ def _is_SR_only(record):
120124
self.set_unresolved()
121125
self.vcf_record.info['SVTYPE'] = 'BND'
122126
self.vcf_record.info['UNRESOLVED_TYPE'] = 'SR_ONLY_LARGE_INVERSION'
127+
if 'END2' not in self.vcf_record.info:
128+
self.vcf_record.info['END2'] = self.vcf_record.stop
123129
self.cluster_type = 'SR_ONLY_LARGE_INVERSION'
124130
for r in self.records:
125131
r.info['UNRESOLVED'] = True
126132
r.info['UNRESOLVED_TYPE'] = 'SR_ONLY_LARGE_INVERSION'
127133
r.info['SVTYPE'] = 'BND'
134+
if 'END2' not in r.info:
135+
r.info['END2'] = r.stop
128136
if 'CPX_TYPE' in r.info.keys():
129137
r.info.pop('CPX_TYPE')
130138
if 'CPX_INTERVALS' in r.info.keys():

0 commit comments

Comments
 (0)