[ensembl-dev] Why do I get duplicated Variation Features

Benoît Ballester benoit at ebi.ac.uk
Wed Jul 4 14:46:06 BST 2012


Hi, 

I am trying to fetch some variants for some slices but get duplicated variation features. I don't understand why, as I would expect one variant per slice.
 
eg: 
7:Y:347772:347784:1  rs71900610      GA/-    Y       347778  347779  deletion        dbSNP   INTERGENIC
7:Y:347772:347784:1  rs71900610      GA/-    Y       347778  347779  deletion        dbSNP   INTERGENIC

or 

48:Y:386863:386872:1     rs10600708      ACAC/-  Y       386864  386867  deletion        dbSNP   INTERGENIC
48:Y:386863:386872:1     rs10600708      ACAC/-  Y       386864  386867  deletion        dbSNP   INTERGENIC

or 

22:Y:116772:116785:-1   rs36189917      G/A     Y       116775  116775  SNP     dbSNP   UPSTREAM
22:Y:116772:116785:-1   rs36189917      G/A     Y       116775  116775  SNP     dbSNP   INTERGENIC
(here UPSTREAM/INTERGENIC difference)


I am sure I am missing something obvious somewhere, but so far I couldn't put my finger on it.  


My code is pretty straightforward :

 my $vfs = $vfa->fetch_all_by_Slice($slice);
	foreach my $vf (@$vfs) {
	    my $v = $vf->variation();
	#print into on slice/variant/variantion-feature
        }
  }


Any feedback appreciated,

Ben

--
Benoit Ballester, PhD



More information about the Dev mailing list