Hi Alex

Thanks for the great program.

I've been getting a segmentation fault error which happens with both the newest patch (2.7.10a_alpha_220601) as well as previous patches for 2.7.10a_alpha and version 2.7.9a. I've narrowed the problem down to when I include --soloFeatures Velocyto as it works fine when I use "--soloFeatures Gene GeneFull SJ". Just noticed this was also mentioned as a side issue in another thread: #1366

As a work around am I right that if I subtract Gene counts from GeneFull counts I would effectively end up with un-spliced counts, and with Gene counts representing spliced+ambiguous counts? If so this may be the best approach anyway as the He et al. 2022 Alevin-fry paper (https://www.nature.com/articles/s41592-022-01408-3) shows that RNA velocity may be improved by combining ambiguous with spliced reads. Or are ambiguous reads not counted in Gene counts?

Thanks
Georgette