EichlerLab / pav

Phased assembly variant caller
98 stars 8 forks source link

maybe a bug in call.snakefile #7

Closed PengJia6 closed 2 years ago

PengJia6 commented 2 years ago

https://github.com/EichlerLab/pav/blob/a9a15b9c61222c4f60f22eb1352d562466b7410c/rules/call.snakefile#L352

I need to changed this line to

lambda row: len(filter_tree[row['#CHROM']][row['POS']:row['END']]) == 0,

in my PAV run.

paudano commented 2 years ago

I about to push an update, and this code changed quite a bit. I verified this bug was fixed in those changes. Thank you for reporting this.