Skip to content

Commit 3c48f89

Browse files
committed
Fix bug if reads jump over gene
If they a given read doesn't support any exons, than it shouldn't support excluding them either. Imagine a case like Gene A: >>>>>------------->>>> Gene B: <<<---<<< A read from gene A shouldn't tell you anything about the splicing of gene B.
1 parent d862119 commit 3c48f89

File tree

1 file changed

+5
-1
lines changed

1 file changed

+5
-1
lines changed

CalculatePSI.py

+5-1
Original file line numberDiff line numberDiff line change
@@ -55,7 +55,11 @@ def get_supported_exons(read, exons_in_gene):
5555
#print("StoppedIteration", cur_exon)
5656
pass
5757
finally:
58-
return supported_lens, supported_exons, excluded_exons
58+
return (
59+
supported_lens,
60+
supported_exons,
61+
(excluded_exons if supported_exons else [])
62+
)
5963

6064

6165
def get_exon_dictionary(gtf_file):

0 commit comments

Comments
 (0)