Hi Everyone,

I have troubles with the following code using the biomart perl API

 my $action='clean';
my $initializer = BioMart::Initializer->new('registryFile'=> $confFile, 'action'=>$action);
 my $registry = $initializer->getRegistry;
my $query = BioMart::Query->new('registry'=> $registry,'virtualSchemaName'=>'default');
 my $symbol = 'xxxx';
        
        $query->addAttribute("gene_flank");
        
        $query->addFilter("mgi_symbol", [$symbol]);
        $query->addFilter("upstream_flank",["1000"]);
        $query->addFilter("transcript_status", ["KNOWN"]);
        $query->addFilter("with_entrezgene", ["Only"]);

        my $query_runner = BioMart::QueryRunner->new();

        $query_runner->execute($query);
        $query_runner->printHeader();
        $query_runner->printResults();
        $query_runner->printFooter();

# I do not get any data back. If i change the attribute it to anything else (eg "cdna" or "peptide" or "5utr") I get actual sequence back. # "gene_flank" is listed as an attribute of dataset "mmusculus_gene_ensembl".

I have tried the code with different genes and on different datasets. No success.

Has anybody had a similar experience?

Thanks, Josef

Reply via email to