This is a test version of Biostars. For the public version, visit https://www.biostars.org.
htslib C API: getting reads from region

I want to use the htslib C API to perform this task

my $pos = "1:100-1000";
my $bam_file = "in.bam";
my @reads = get_reads_from_htslib($bam_file,$pos)

foreach(@reads) {
    # do stuff...
}

I wrote it in perl but you get the idea. I can't find an example anywhere online. Does anyone have experience with htslib and getting reads from a region?

samtools htslib c api
• 8,453 views • link
updated by Ram Ram • written by QVINTVS_FABIVS_MAXIMVS QVINTVS_FABIVS_MAXIMVS

1 answer

here is a C code:

Compilation:

gcc -Isamtools -Ihtslib -Lsamtools -Lhtslib biostar151053.c htslib/libhts.a samtools/libbam.a -lz -lpthread
• 1 views • link
updated by Ram Ram • written by Pierre Lindenbaum Pierre Lindenbaum

What function calls here come from samtools.h?

• 0 views • link
written by Danielk Danielk

Dear Pierre, I have another question. Now I have a sam file without the index file, and I want to iterate all the reads of the sam file, how can I make it? Thanks in advance.

• 0 views • link
written by lghust2011 lghust2011

It seems that I have to compress the sam file to bam file, and get the index file then I can make it..

• 0 views • link
written by lghust2011 lghust2011

Dear Pierre I used the exact code and bash line but I face an error. I've downloaded htslib and samtools and put them in the same folder as the .c file.

.... "_lzma_stream_decoder", referenced from: _cram_uncompress_block in libhts.a(cram_io.o) ld: symbol(s) not found for architecture x86_64 clang: error: linker command failed with exit code 1 (use -v to see invocation)

Should I build samtools and htslib as mentioned here ?

• 0 views • link
written by sina.majidian • sina.majidian

I think so, furthermore, the code above is ~5 years old, the htslib has since much changed.

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Log in to answer this question.

More posts like this