samtools VS htslib

I don’t know about your first question, but for your second, use “bwa mem -p …” for smart pairing of an interleaved FASTQ/A with both paired end reads and singletons. BWA will recognize whether adjacent reads are paired if they have the same prefix (see here)

In this post, I'm going to add the command tanuki to Samtools to display a ASCII art(AA) tanuki.

Quite a number of bioinformatics related tool that were written in C/C++, eg. (https://github.com/samtools/samtools). But there are also a lot of modern packages for Python / R now

There is another nice multi-core gzip based library called BGZF[1]. It is commonly used in bioinformatics. BGZF has the added advantage that it is block compressed with built in indexing method to permit seeking in compressed files.

[1] https://github.com/samtools/htslib

In contrast, processing can start immediately and only transfer what is necessary if tooling can read bytes of data directly from Amazon S3. Tools based on htslib can do this, so you can run something like:

Another idea is add in reading cloud data natively e.g. htslib which samtools relies on can read s3 directly. https://github.com/samtools/htslib

Greetings! I'm a bioinformatics software dev in San Francisco and I'm looking for others interested in working on a high-performance genetic data analysis project. The project is in C, using https://github.com/samtools/htslib.

I'm working on a bioinformatics-related binding called ruby-htslib. htslib makes heavy use of bit fields throughout the library, so supporting bit fields is inevitable.