Library Htsjdk to query a vcf file
0
0
Entering edit mode
8.8 years ago
onemoreuser ▴ 20

I'm using:

VCFFileReader file= new VCFFileReader(file.vcf, file.vcf.gz.tbi, true)

to open a file and then

CloseableIterator<VariantContext> vc = file.query("10", 77162995, 131048262)

to query specific regions of the file. The problem is that when I choose a region in the middle or the end of the chromosome the method doesn't return anything. For the region starting at position 77162995 and end at 131048262 it doesn't work.

java htsjdk vcf annotation • 1.8k views
ADD COMMENT
0
Entering edit mode

Are there any records in this interval? If there aren't any it might be normal to return nothing (I haven't tried...)

ADD REPLY

Login before adding your answer.

Traffic: 2110 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6