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