[ensembl-dev] Variants seq_region_name excluding patches / seq_region_start for VCF
Genomeo Dev
genomeodev at gmail.com
Wed Dec 17 18:03:45 GMT 2014
Hi,
*Using Ensembl 75 VM*
I am trying to retrieve coordinates and sequence regions for some variants
using variant IDs:
my $var = $var_adaptor->fetch_by_name('rs7904594');
my $vfeatures = $var->get_all_VariationFeatures();
while (my $vf = shift @{$vfeatures}){
print $vf->seq_region_name(), "\t", $vf->seq_region_start(), "\n";
}
*Questions:*
1. Is there a way I can choose to print main chromosomes only excluding
patches?
2. I want to use the position information in a VCF. seq_region_start only
works for SNPs. How can I make use of that for indels?
Thanks!
--
G.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.ensembl.org/pipermail/dev_ensembl.org/attachments/20141217/38e86e63/attachment.html>
More information about the Dev
mailing list