SeqLib: a C plus plus API for rapid BAM manipulation, sequence alignment and sequence assembly
SeqLib: a C plus plus API for rapid BAM manipulation, sequence alignment and sequence assembly
复制标题
DOI:
10.1093/bioinformatics/btw741
复制
发表时间:
2017-03-01
期刊:
影响因子:
5.8
通讯作者:
Beroukhim, Rameen
中科院分区:
文献类型:
--
作者:
Wala, Jeremiah;Beroukhim, Rameen
We present SeqLib, a C++API and command line tool that provides a rapid and user-friendly interface to BAM/SAM/CRAM files, global sequence alignment operations and sequence assembly. Four C libraries perform core operations in SeqLib: HTSlib for BAM access, BWA-MEM and BLAT for sequence alignment and Fermi for error correction and sequence assembly. Benchmarking indicates that SeqLib has lower CPU and memory requirements than leading C++sequence analysis APIs. We demonstrate an example of how minimal SeqLib code can extract, error-correct and assemble reads from a CRAM file and then align with BWA-MEM. SeqLib also provides additional capabilities, including chromosome-aware interval queries and read plotting. Command line tools are available for performing integrated error correction, micro-assemblies and alignment.