pagoda2 VS alevin-fry

Try pagoda2 for this: https://github.com/kharchenkolab/pagoda2

I realize that this is a bit old, but I find pagoda2 faster: https://github.com/kharchenkolab/pagoda2

You can also get this information easily from the raw data using our alevin-fry tool. It produces a count matrix with spliced/unspliced/ambiguous counts that can easily be imported into python or R for velocity analysis. You can read more about alevin-fry here and check it out on github here. It's also easily installable via bioconda, and there's even a tool to simplify execution, called simpleaf. I'm happy to answer any questions you might have about using it!

I think it depends on where in the space one is working. Our most recent tool, for example, is alevin-fry, which is for efficient preprocessing of single cell RNA seq data. Rust is the natural replacement for tools in the space that otherwise would have been in C++ (aligners, assemblers, large scale indexing, and things like RNA seq quantification and preprocessing). I view julia as more appropriate to replace tools currently existing in Python or R. There, I believe it can offer some substantial benefits. For our applications however, garbage collection is usually a non-starter.

My lab has been trying to move most of our development over to rust (from C++). Our most recent paper presents a tool for single-cell and single-nucleus RNA-seq processing in rust https://rdcu.be/cIL35. The repo is https://github.com/COMBINE-lab/alevin-fry. It's one of our first "big" rust packages; so hopefully our design practices and rust skills will continue to improve.