check for genetic_map being totally outside the vcf region

odelaneau / shapeit4

Segmented HAPlotype Estimation and Imputation Tool

MIT License

89 stars 17 forks source link

check for genetic_map being totally outside the vcf region #13

Closed jaredo closed 4 years ago

jaredo commented 4 years ago

In some cases the genetic map can totally exclude the specified vcf region causing i_locus to exceed the length of vec_pos which results in a segfault.

For example chr21:5130578-10162681 on hg38 was problematic for me since the map starts at 10326676.

odelaneau commented 4 years ago

Just seen this! I do not get notified for pull requests...

ANy way, thank you!