Hello Peter,
thank you for the answer.
I actually tried several combinations of beginning and end, but it's not
clear to me how to fetch the BAM record if present.
I see in the iterator that */i/* is always = -1 which seems to indicate
that it's empty, but I'm not sure on how to use the iterator.
Is there any documentation about the HTSlib API apart from the code
comments?
How am I supposed to use the iterator to retrieve the fetched read(s)?
thanks
Claudio
On 17/03/2016 08:49, Peter Johansson wrote:
Hi Claudio,
On 03/17/2016 01:53 AM, Claudio Alberti wrote:
iterator = *bam_itr_queryi*(BAMindex, 1, 20000, 20000);
try
iterator = *bam_itr_queryi*(BAMindex, 1, 20000, 20001);
Cheers,
Peter
--
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
