vcf to vranges
1
0
Entering edit mode
arsala521 ▴ 10
@arsala521-10976
Last seen 7.7 years ago

Hello

I made vranges object from my VCF file usingreadVcfAsVRanges. There are 37974 variants in the VCF file and in vranges there are 14682016 ranges. Can someone please explain this difference?

I guess that if a variant is in more than one sample then there will be more than one range for that variant in vranges object but even then there are too many number of ranges.

Regards

 

 

 

readvcfasvranges vcf to vranges • 1.1k views
ADD COMMENT
0
Entering edit mode
@valerie-obenchain-4275
Last seen 2.2 years ago
United States

The data are expanded (replicated into rows) by the number of ALT values as well as the samples which may explain the difference. Reading in a small subset with a 'param' would allow you to inspect manually to see what's going on, e.g.,

param <- ScanVcfParam(which=GRanges("chr6", IRanges(5000, 6000)))

readVcfAsVRanges(myfile, param=param)

Valerie

ADD COMMENT
0
Entering edit mode

thank you so much.

 

ADD REPLY

Login before adding your answer.

Traffic: 796 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