Hello,
I'm trying to find a way to retrieve one read by position in a BAM file in a C program (using htslib).
Till now the best candidate seems to be *bam_itr_queryi*

e.g. if I want to fetch the read at position 20000 in reference 1 I'm trying to use *bam_itr_queryi* in this way

iterator = *bam_itr_queryi*(BAMindex, 1, 20000, 20000);

But iterator is empty...
Any idea?
thank you
Claudio

--
Claudio Alberti
----------------------------------------------
http://gramm.epfl.ch
EPFL SCI STI MM
ELG 140 (ELG Building)
Station 11
CH-1015 Lausanne - Switzerland
Tel. +41 21 6936869
----------------------------------------------

------------------------------------------------------------------------------
Transform Data into Opportunity.
Accelerate data analysis in your applications with
Intel Data Analytics Acceleration Library.
Click to learn more.
http://pubads.g.doubleclick.net/gampad/clk?id=278785231&iu=/4140
_______________________________________________
Samtools-help mailing list
Samtools-help@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/samtools-help

Reply via email to