secastel / phaser

phasing and Allele Specific Expression from RNA-seq
GNU General Public License v3.0
111 stars 36 forks source link

Issue in function split_read #84

Open chrbianchini opened 4 months ago

chrbianchini commented 4 months ago

Hi, I'm using phaser on the bams from the TCGA dataset and it crashes when the CIGAR string is *. In particular the instruction seq_len = int(number_build) at row 197 of file read_variant_map.py throws a ValueError because the string number_build is empty. I "fixed" this in my project by returning an empty list if cigar == "*" but maybe it can be handled in a more elegant way. Greetings,

Christian